Evolutionary Conservation Analysis Research Articles

A novel variant c.902C>A (p. A301D) in KCNQ4 associated with non-syndromic deafness 2A in a Chinese family.

Deafness autosomal dominant 2A (DFNA2A) is related to non-syndromic genetic hearing impairment. The KCNQ4 (Potassium Voltage-Gated Channel Subfamily Q Member 4) can lead to DFNA2A. In this study, we report a case of autosomal dominant non-syndromic hearing loss with six family members as caused by a novel variant in the KCNQ4 gene. The whole-exome sequencing (WES) and pure tone audiometry were performed on the proband of the family. Sanger sequencing was conducted on family members to determine if the novel variant in the KCNQ4 gene was present. Evolutionary conservation analysis and computational tertiary structure protein prediction of the wild-type KCNQ4 protein and its variant were then performed. In addition, voltage-gated channel activity of the wild-type KCNQ4 protein and its variant were tested using whole-cell patch clamp. It was observed that the proband had inherited autosomal dominant, non-syndromic sensorineural hearing loss as a trait. A novel co-segregating heterozygous missense variant (c.902C>A, p.Ala301Asp) of the KCNQ4 gene was identified in the proband and other five affected family members. This variant was predicted to cause an alanine-to-aspartic acid substitution at position 301 in the KCNQ4 protein. The alanine at position 301 is well conserved across different species. Whole-cell patch clamp showed that there was a significant difference between the WT protein currents and the mutant protein currents in the voltage-gated channel activity. In the present study, performing WES in conjunction with Sanger sequencing enhanced the detection of a novel, potentially causative variant (c301 A>G; p.Ala301Asp) in exon 6 of the KCNQ4 gene. Therefore, our findings contributed to the mutation spectrum of the KCNQ4 gene and may be useful in the diagnosis and gene therapy of deafness autosomal dominant 2A.

Identifying Oncogenic Missense Single Nucleotide Polymorphisms in Human SAT1 Gene Using Computational Algorithms and Molecular Dynamics Tools.

The human SAT1 gene encodes spermidine/spermine N1-acetyltransferase 1 (SSAT1), a regulatory biological catalyst of polyamine catabolism. Numerous essential biological processes, such as cellular proliferation, differentiation, and survival, depend on polyamines like spermidine and spermine. Thus, SSAT1 is involved in key cellular activities such as proliferation and survival of cells and mediates various diseases including cancer. A plethora of studies established the involvement of missense single nucleotide polymorphisms (SNPs) in numerous pathological conditions due to their ability to adversely affect the structure and subsequent function of the protein. To date, an in silico study to identify the pathogenic missense SNPs of the human SAT1 gene has not been accomplished yet. This study aimed to filter the missense SNPs that were functionally detrimental and pathogenic. The rs757435207 (I21N) was ascertained to be the most deleterious and pathogenic by all algorithmic tools. Stability and evolutionary conservation analysis tools also stated that I21N variant decreased the stability and was located in the highly conserved residue. Molecular dynamics simulation revealed that I21N caused substantial alterations in the conformational stability and dynamics of the SSAT1 protein. Consequently, the I21N variant could disrupt the native functional roles of the SSAT1 enzyme. Therefore, the I21N variant was identified and concluded to be an oncogenic missense variant of the human SAT1 gene. Overall, the findings of this study would be a great directory of future experimental research to develop personalized medicine.

Analysis of Evolutionary Conservation, Expression Level, and Genetic Association at a Genome-wide Scale Reveals Heterogeneity Across Polygenic Phenotypes.

Understanding the expression level and evolutionary rate of associated genes with human polygenic diseases provides crucial insights into their disease-contributing roles. In this work, we leveraged genome-wide association studies (GWASs) to investigate the relationship between the genetic association and both the evolutionary rate (dN/dS) and expression level of human genes associated with the two polygenic diseases of schizophrenia and coronary artery disease. Our findings highlight a distinct variation in these relationships between the two diseases. Genes associated with both diseases exhibit a significantly greater variance in evolutionary rate compared to those implicated in monogenic diseases. Expanding our analyses to 4,756 complex traits in the GWAS atlas database, we unraveled distinct trait categories with a unique interplay among the evolutionary rate, expression level, and genetic association of human genes. In most polygenic traits, highly expressed genes were more associated with the polygenic phenotypes compared to lowly expressed genes. About 69% of polygenic traits displayed a negative correlation between genetic association and evolutionary rate, while approximately 30% of these traits showed a positive correlation between genetic association and evolutionary rate. Our results demonstrate the presence of a spectrum among complex traits, shaped by natural selection. Notably, at opposite ends of this spectrum, we find metabolic traits being more likely influenced by purifying selection, and immunological traits that are more likely shaped by positive selection. We further established the polygenic evolution portal (evopolygen.de) as a resource for investigating relationships and generating hypotheses in the field of human polygenic trait evolution.

Key amino acid residues govern the substrate selectivity of the transporter Xltr1p from Trichoderma reesei for glucose, mannose, and galactose

This research identified four amino acid residues (Leu174, Asn297, Tyr301, and Gln291) that contribute to substrate recognition by the high-affinity glucose transporter Xltr1p from Trichoderma reesei. Potential hotspots affecting substrate specificity were selected through homology modeling, evolutionary conservation analyses, and substrate-docking modeling of Xltr1p. Variants carrying mutations at these hotspots were subsequently obtained via in silico screening. Replacement of Leu174 or Asn297 in Xltr1p with alanine resulted in loss of hexose transport activity, indicating that Leu174 and Asn297 play essential roles in hexose transport. The Y301W variant exhibited accelerated mannose transport, but lost galactose transport capacity, and mutation of Gln291 to alanine greatly accelerated mannose transport. These results suggest that amino acids located in transmembrane α-helix 7 (Asn297, Tyr301, and Gln291) play critical roles in substrate recognition by the hexose transporter Xltr1p. Our results will help expand the potential applications of this transporter and provide insights into the mechanisms underlying its function and specificity.

An evolutionary view of the Fusarium core genome

Fusarium, a member of the Ascomycota fungi, encompasses several pathogenic species significant to plants and animals. Some phytopathogenic species have received special attention due to their negative economic impact on the agricultural industry around the world. Traditionally, identification and taxonomic analysis of Fusarium have relied on morphological and phenotypic features, including the fungal host, leading to taxonomic conflicts that have been solved using molecular systematic technologies. In this work, we applied a phylogenomic approach that allowed us to resolve the evolutionary history of the species complexes of the genus and present evidence that supports the F. ventricosum species complex as the most basal lineage of the genus. Additionally, we present evidence that proposes modifications to the previous hypothesis of the evolutionary history of the F. staphyleae, F. newnesense, F. nisikadoi, F. oxysporum, and F. fujikuroi species complexes. Evolutionary analysis showed that the genome GC content tends to be lower in more modern lineages, in both, the whole-genome and core-genome coding DNA sequences. In contrast, genome size gain and losses are present during the evolution of the genus. Interestingly, core genome duplication events positively correlate with genome size. Evolutionary and genome conservation analysis supports the F3 hypothesis of Fusarium as a more compact and conserved group in terms of genome conservation. By contrast, outside of the F3 hypothesis, the most basal clades only share 8.8% of its genomic sequences with the F3 clade.

Growth-Coupled Evolutionary Pressure Improving Epimerases for D-Allulose Biosynthesis Using a Biosensor-Assisted In Vivo Selection Platform.

Fast screening strategies that enable high-throughput evaluation and identification of desired variants from diversified enzyme libraries are crucial to tailoring biocatalysts for the synthesis of D-allulose, which is currently limited by the poor catalytic performance of ketose 3-epimerases (KEases). Here, the study designs a minimally equipment-dependent, high-throughput, and growth-coupled in vivo screening platform founded on a redesigned D-allulose-dependent biosensor system. The genetic elements modulating regulator PsiR expression levels undergo systematic optimization to improve the growth-responsive dynamic range of the biosensor, which presents ≈30-fold facilitated growth optical density with a high signal-to-noise ratio (1.52 to 0.05) toward D-allulose concentrations from 0 to 100mm. Structural analysisand evolutionary conservation analysis of Agrobacterium sp. SUL3 D-allulose 3-epimerase (ADAE) reveal a highly conserved catalytic active site and variable hydrophobic pocket, which together regulate substrate recognition. Structure-guided rational design and directed evolution are implemented using the growth-coupled in vivo screening platform to reprogram ADAE, in which a mutant M42 (P38N/V102A/Y201L/S207N/I251R) is identified with a 6.28-fold enhancement of catalytic activity and significantly improved thermostability with a 2.5-fold increase of the half-life at 60 °C. The research demonstrates that biosensor-assisted growth-coupled evolutionary pressure combined with structure-guided rational design provides a universal route for engineering KEases.

Seq-InSite: sequence supersedes structure for protein interaction site prediction.

Proteins accomplish cellular functions by interacting with each other, which makes the prediction of interaction sites a fundamental problem. As experimental methods are expensive and time consuming, computational prediction of the interaction sites has been studied extensively. Structure-based programs are the most accurate, while the sequence-based ones are much more widely applicable, as the sequences available outnumber the structures by two orders of magnitude. Ideally, we would like a tool that has the quality of the former and the applicability of the latter. We provide here the first solution that achieves these two goals. Our new sequence-based program, Seq-InSite, greatly surpasses the performance of sequence-based models, matching the quality of state-of-the-art structure-based predictors, thus effectively superseding the need for models requiring structure. The predictive power of Seq-InSite is illustrated using an analysis of evolutionary conservation for four protein sequences. Seq-InSite is freely available as a web server at http://seq-insite.csd.uwo.ca/ and as free source code, including trained models and all datasets used for training and testing, at https://github.com/lucian-ilie/Seq-InSite.

Whole exome sequencing identifies a novel mutation in Annexin A4 that is associated with recurrent spontaneous abortion.

Recurrent spontaneous abortion (RSA) is a multifactorial disease, the exact causes of which are still unknown. Environmental, maternal, and genetic factors have been shown to contribute to this condition. The aim of this study was to investigate the presence of mutations in the ANXA4 gene in patients with RSA. Genomic DNA was extracted from 325 patients with RSA and 941 control women with a normal reproductive history for whole-exome sequencing (WES). The detected variants were annotated and filtered, and the pathogenicity of the variants was predicted through the SIFT online tool, functional enrichment analyses, Sanger sequencing validation, prediction of changes in protein structure, and evolutionary conservation analysis. Furthermore, plasmid construction, Western blotting, RT-qPCR, and cell migration, invasion and adhesion assays were used to detect the effects of ANXA4 mutations on protein function. An ANXA4 mutation (p.G8D) in 1 of the 325 samples from patients with RSA (RSA-219) was identified through WES. This mutation was not detected in 941 controls or included in public databases. Evolutionary conservation analysis revealed that the amino acid residue affected by the mutation (p.G8D) was highly conserved among 13 vertebrate species, and the SIFT program and structural modeling analysis predicted that this mutation was harmful. Furthermore, functional assays revealed that this mutation could inhibit cell migration, invasion and adhesion. Our study suggests that an unreported novel ANXA4 mutation (p.G8D) plays an important role in the pathogenesis of RSA and may contribute to the genetic diagnosis of RSA.

Unmasking crucial residues in adipose triglyceride lipase for coactivation with comparative gene identification-58

Lipolysis is an essential metabolic process that releases unesterified fatty acids from neutral lipid stores to maintain energy homeostasis in living organisms. Adipose triglyceride lipase (ATGL) plays a key role in intracellular lipolysis and can be coactivated upon interaction with the protein comparative gene identification-58 (CGI-58). The underlying molecular mechanism of ATGL stimulation by CGI-58 is incompletely understood. Based on analysis of evolutionary conservation, we used site directed mutagenesis to study a C-terminally truncated variant and full-length mouse ATGL providing insights in the protein coactivation on a per-residue level. We identified the region from residues N209-N215 in ATGL as essential for coactivation by CGI-58. ATGL variants with amino acids exchanges in this region were still able to hydrolyze triacylglycerol at the basal level and to interact with CGI-58, yet could not be activated by CGI-58. Our studies also demonstrate that full-length mouse ATGL showed higher tolerance to specific single amino acid exchanges in the N209-N215 region upon CGI-58 coactivation compared to C-terminally truncated ATGL variants. The region is either directly involved in protein-protein interaction or essential for conformational changes required in the coactivation process. Three-dimensional models of the ATGL/CGI-58 complex with the artificial intelligence software AlphaFold demonstrated that a large surface area is involved in the protein-protein interaction. Mapping important amino acids for coactivation of both proteins, ATGL and CGI-58, onto the 3D model of the complex locates these essential amino acids at the predicted ATGL/CGI-58 interface thus strongly corroborating the significance of these residues in CGI-58–mediated coactivation of ATGL.

Comparative analysis of PDZ-binding motifs in the diacylglycerol kinase family.

Diacylglycerol kinases (DGKs) control local and temporal amounts of diacylglycerol (DAG) and phosphatidic acid (PA) by converting DAG to PA through phosphorylation in cells. Certain DGK enzymes possess C-terminal sequences that encode potential PDZ-binding motifs (PBMs), which could be involved in their recruitment into supramolecular signaling complexes. In this study, we used two different interactomic approaches, quantitative native holdup (nHU) and qualitative affinity purification (AP), both coupled to mass spectrometry (MS) to investigate the PDZ partners associated with the potential PBMs of DGKs. Complementing these results with site-specific affinity interactomic data measured on isolated PDZ domain fragments and PBM motifs, as well as evolutionary conservation analysis of the PBMs of DGKs, we explored functional differences within different DGK groups. All our results indicate that putative PBM sequences of type II enzymes, namely DGKδ, DGKη, and DGKκ, are likely to be nonfunctional. In contrast, type IV enzymes, namely DGKζ and DGKι, possess highly promiscuous PBMs that interact with a set of PDZ proteins with very similar affinity interactomes. The combination of various interactomic assays and evolutionary analyses provides a useful strategy for identifying functional domains and motifs within diverse enzyme families.

Occurrence and evolutionary conservation analysis of α-helical cationic amphiphilic segments in the human proteome.

The existence of encrypted fragments with antimicrobial activity in human proteins has been thoroughly demonstrated in the literature. Recently, algorithms for the large-scale identification of these segments in whole proteomes were developed, and the pervasiveness of this phenomenon was stated. These algorithms typically mine encrypted cationic and amphiphilic segments of proteins, which, when synthesized as individual polypeptide sequences, exert antimicrobial activity by membrane disruption. In the present report, the human reference proteome was submitted to the software kamal for the uncovering of protein segments that correspond to putative intragenic antimicrobial peptides (IAPs). The assessment of the identity of these segments, frequency, functional classes of parent proteins, structural relevance, and evolutionary conservation of amino acid residues within their corresponding proteins was conducted in silico. Additionally, the antimicrobial and anticancer activity of six selected synthetic peptides was evaluated. Our results indicate that cationic and amphiphilic segments can be found in 2% of all human proteins, but are more common in transmembrane and peripheral membrane proteins. These segments are surface-exposed basic patches whose amino acid residues present similar conservation scores to other residues with similar solvent accessibility. Moreover, the antimicrobial and anticancer activity of the synthetic putative IAP sequences was irrespective to whether these are associated to membranes in the cellular setting. Our study discusses these findings in light of the current understanding of encrypted peptide sequences, offering some insights into the relevance of these segments to the organism in the context of their harboring proteins or as separate polypeptide sequences.

Molecular dynamics simulations corroborate recombinant expression studies carried out on three αIIb β-propeller mutations reported in Indian Glanzmann thrombasthenia patients.

Mutations in the αIIb β-propeller domain have long been known to disrupt heterodimerization and intracellular trafficking of αIIbβ3 complexes leading to diminished surface expression and/or function, resulting in Glanzmann thrombasthenia. Our previous study on three β-propeller mutations, namely G128S, S287L, and G357S, showed variable defects in protein transport correlated with the patient's clinical phenotypes. Pulse-chase experiments revealed differences in αIIbβ3 complex maturation among the three mutations. Hence, the current study aims to correlate conformational changes caused by each one of them. Evolutionary conservation analysis, stability analysis, and molecular dynamics simulations of the three mutant structures were carried out. Stability analysis revealed that, while G128S and G357S mutations destabilized the β-propeller structure, S287L retained the stability. Wild-type and mutant β-propeller structures, when subjected to molecular dynamics simulations, confirmed that G128S and G357S were both destabilizing in nature when compared with the wild-type and S287L based on several parameters studied, like RMSD, RMSF, Rg, FEL, PCA, secondary structure, and hydrogen bonds. In our previous study, we demonstrated that mutant S287L αIIbβ3 complexes were more stable than the wild-type αIIbβ3 complexes, as evidenced in pulse-chase experiments. These findings corroborate variable intracellular fates of mutant αIIbβ3 complexes as a result of these β-propeller mutations.

Computational identification of a multi-peptide vaccine candidate in E2 glycoprotein against diverse Hepatitis C virus genotypes

Hepatitis C Virus (HCV) is estimated to affect nearly 180 million people worldwide, culminating in ∼0.7 million yearly casualties. However, a safe vaccine against HCV is not yet available. This study endeavored to identify a multi-genotypic, multi-epitopic, safe, and globally competent HCV vaccine candidate. We employed a consensus epitope prediction strategy to identify multi-epitopic peptides in all known envelope glycoprotein (E2) sequences, belonging to diverse HCV genotypes. The obtained peptides were screened for toxicity, allergenicity, autoimmunity and antigenicity, resulting in two favorable peptides viz., P2 (VYCFTPSPVVVG) and P3 (YRLWHYPCTV). Evolutionary conservation analysis indicated that P2 and P3 are highly conserved, supporting their use as part of a designed multi-genotypic vaccine. Population coverage analysis revealed that P2 and P3 are likely to be presented by >89% Human Leukocyte Antigen (HLA) molecules from six geographical regions. Indeed, molecular docking predicted the physical binding of P2 and P3 to various representative HLAs. We designed a vaccine construct using these peptides and assessed its binding to toll-like receptor 4 (TLR-4) by molecular docking and simulation. Subsequent analysis by energy-based and machine learning tools predicted high binding affinity and pinpointed the key binding residues (i.e. hotspots) in P2 and P3. Also, a favorable immunogenic profile of the construct was predicted by immune simulations. We encourage the scientific community to validate our vaccine construct in vitro and in vivo. Communicated by Ramaswamy H. Sarma

Mutational analysis of minichromosome maintenance complex component (MCM) family genes in Chinese Han women with polycystic ovarian syndrome

Purpose To investigate whether mutations in the minichromosome maintenance complex component (MCM) family genes were present in patients with polycystic ovary syndrome (PCOS) of Chinese descent. Methods A total of 365 Chinese patients with PCOS and 860 women without PCOS as control who underwent with assisted reproductive technology were enrolled. Genomic DNA was extracted from the peripheral blood of these patients for PCR and Sanger sequencing. The potential damage of these mutations/rare variants was analyzed through evolutionary conservation analysis and bioinformatic programs. Results Twenty-nine missense or nonsense mutations/rare variants in the MCM genes were identified in 365 patients with PCOS (7.9%, 29/365), all these mutations/rare variants were predicted to be ‘disease causing’ by SIFT and PolyPhen2 programs. Among those, four mutations were reported here for the first time, p.S7C (c.20C > G) in MCM2 (NM_004526.3), p.K350R (c.1049A > G) in MCM5 (NM_006739.3), p.K283N (c.849G > T) in MCM10 (NM_182751.2), and p.S1708F (c.5123C > T) in MCM3AP (NM_003906.4). All of these novel mutations were not found in our 860 control women, or also absent in public databases. In addition, the evolutionary conservation analysis results suggested that these novel mutations caused highly conserved amino acid substitutions among 10 vertebrate species. Conclusion This study identified a high frequency of potential pathogenic rare variants/mutations in MCM family genes in Chinese women with PCOS, which further expands the genotype spectrum in PCOS.

Novel and recurrent FBN1 mutations causing Marfan syndrome in two Chinese families.

To explore the genetic defects of two families with autosomal dominant Marfan syndrome (MFS). Two families with MFS were enrolled in this study. The detailed ocular presentations of the patients were recorded. Whole exome sequencing was performed to explore the pathogenic variants and Sanger sequencing was performed to confirm the gene mutations. Segregation analysis among the family members was made and bioinformatics analysis was performed to predict the functional impact of the mutations. The main ocular presentations of the probands were increased axial length and ectopia lentis. Using whole exome sequencing and Sanger sequencing, a novel heterozygous missense mutation (c.5060G > C, p.Cys1687Ser) and a recurrent missense mutation (c.2168A > T, p.Asp723Val) were identified within FBN1, which were co-segregated with the MFS phenotype in the families. Evolutionary conservation analysis showed that codons 723 and 1,687 were highly conserved among several species. Functional impact predictions made using several online programs suggested that the mutations were pathogenic. We identified a novel and a recurrent missense mutation in FBN1 in two Chinese families with MFS using whole exome sequencing, and our bioinformatics analysis indicated that the mutations were disease-causing. Our results expand the mutation spectrum of FBN1 and could help us better understand the genetic defects of the patients with MFS.

A novel EDAR variant identified in non-syndromic tooth agenesis: Insights from molecular dynamics

ObjectiveThis study aims to investigate a novel pathogenic variant in a Chinese family of non-syndromic tooth agenesis (NSTA) and study the impact of the variant on related protein and pathway. DesignOne NSTA family was collected. Whole exome sequencing and Sanger sequencing were performed on the proband with NSTA and his 5 family members. The pathogenic influence of the mutant is evaluated by bioinformatics analyses including evolutionary conservation analysis and secondary structure prediction. Molecular dynamics (MD) simulations and binding free energy calculations were then performed to explore changes in the tertiary structure and binding ability of the protein. ResultsWe found a novel missense ectodysplasin A receptor (EDAR) variant (c .1292 T > G; p.Ile431Arg) in all affected family members. The results of bioinformatics analyses revealed that the EDAR had harmful changes after mutation. MD simulations and the binding free energy calculations results showed that the mutant EDAR protein and EDAR/ectodysplasin-A receptor-associated adapter (EDARADD) complex displayed tertiary structural change, and EDAR possessed a lower affinity to EDARADD after mutation. ConclusionsWe found a novel EDAR variant (c.1292 T > G; p.Ile431Arg) in one NSTA family, which affects the binding of EDAR and EDARADD.

Evolutionary Conservation and Transcriptome Analyses Attribute Perenniality and Flowering to Day-Length Responsive Genes in Bulbous Barley (Hordeum bulbosum).

Rapid population growth and dramatic climatic turnovers are challenging global crop production. These challenges are spurring plant breeders to enhance adaptation and sustainability of major crops. One intriguing approach is to turn annual systems into perennial ones, yet long-term classical breeding efforts to induce perenniality have achieved limited success. Here, we report the results of our investigation of the genetic basis of bulb formation in the nonmodel organism Hordeum bulbosum, a perennial species closely related to barley. To identify candidate genes that regulate bulb formation in H. bulbosum, we applied two complementary approaches. First, we explored the evolutionary conservation of expressed genes among annual Poaceae species. Next, we assembled a reference transcriptome for H. bulbosum and conducted a differential expression (DE) analysis before and after stimulating bulb initiation. Low conservation was identified in genes related to perenniality in H. bulbosum compared with other species, including bulb development and sugar accumulation genes. We also inspected these genes using a DE analysis, which enabled identification of additional genes responsible for bulb initiation and flowering regulation. We propose a molecular model for the regulation of bulb formation involving storage organ development and starch biosynthesis genes. The high conservation observed along a major part of the pathway between H. bulbosum and barley suggests a potential for the application of biotechnological techniques to accelerate breeding toward perenniality in barley.

Identification of a novel ESR1 mutation in a Chinese PCOS woman with estrogen insensitivity in IVF treatment

BackgroundPolycystic ovary syndrome (PCOS) is a complex reproductive disorder, that affects approximately 5–10% of women of reproductive age. The disease is complex because its evolution may be impacted by genetic, lifestyle and environmental factors. Previous studies have emphasized the important roles of estrogen receptors in the pathogenesis of PCOS.ObjectiveTo use whole exome sequencing (WES) to assess possible pathogenic factors in a PCOS patient who exhibited estrogen insensitivity during hormone replacement therapy (HRT) treatment.MethodsGenome sequencing and variant filtering via WES were performed in a patient with PCOS. DNA extraction from 364 unrelated female controls without PCOS was followed by PCR amplification, Sanger sequencing and sequence alignment. Evolutionary conservation analysis, protein structural modelling and in silico prediction were applied to analyse the potential pathogenicity of the novel ESR1 mutation.Result(s)During the controlled ovarian hyperstimulation (COH) period of an IVF cycle, the patient experienced markedly prolonged ovarian stimulation due to a poor response to gonadotropins (Gn) and elevated serum FSH. A novel heterozygous ESR1 mutation, c.619G > A/p.A207T, leading to the replacement of a highly conserved alanine with a threonine, was identified in this patient, via WES analysis. This novel variant was not identified in 364 unrelated female controls without PCOS, or in the Exome Aggregation Consortium (ExAC) or 1000 Genome Project.Conclusion(s)We identified a novel heterozygous ESR1 mutation in a Han Chinese PCOS woman exhibiting clinical signs of estrogen insensitivity. This study may provide new strategies for IVF therapy, especially for patients who exhibit estrogen insensitivity during IVF cycle.

Novel rare mutation in a conserved site of PTPRB causes human hypoplastic left heart syndrome.

Hypoplastic left heart syndrome (HLHS) is a rare but fatal birth defect in which the left side of the heart is underdeveloped. HLHS accounts for 2% to 4% of congenital heart anomalies. Whole genome sequencing (WGS) was conducted for a family trio consisting of a proband and his parents. A homozygous rare variant was detected in the PTPRB (Protein Tyrosine Phosphatase Receptor Type B) gene of the proband by functional annotation and co-segregation analysis. Sanger sequencing was used to confirm genotypes of the variant. The in silico prediction tools, including Mutation Taster, SpliceAI, and CADD, were used to predict the impact of the mutation. The allele frequencies across populations were compared based on multiple databases, including "1000 genomes" and "gnomAD". We used two vectors (pcMINI and pcDNA3.1) to generate a minigene construct to validate the mutational effect at the transcriptional level. Family-based WGS analyses showed that only a homozygous splice acceptor variant (NC_000012.12: g.70636068T>G, NM_001109754.4: c.56-2A>C, NG_029940.2: g.6373A>C) at the exon-intron border of PTPRB gene associates with HLHS. This variant is also within the region with the enhancer activity based on UCSC genome annotation. Genotyping and Sanger sequencing revealed that the proband's parents are heterozygous for this variant. Evolutionary conservation analysis revealed that the site (NC_000012.12: g.70636068) is extremely conserved across species, supporting the evolutionary functional constraints of the ancestral wild type (T). In silico tools universally predicted a deleterious or disease-causing impact of the mutation from T to G. The mutation was not found in the 1000 genomes and gnomAD databases, which indicates that this mutation is very rare in most human populations. A splicing assay indicated that the mutated minigene caused aberrant splicing of mRNA, in which a 3bp missing in the second exon resulted in the deletion of one amino acid (NP_001103224.1:p.Glu19del) compared to the normal protein of PRPTB (also the VE-PTP). Structure prediction revealed that the deletion occurred within the C-region of the signal peptide of VE-PTP, suggesting signal peptide-related defects as a potential mechanism for the HLHS cellular pathogeny. We report a rare homozygous variant with splicing error in PTPRB associated with HLHS. Previous model species studies revealed conserved functions of PTPRB in cardiovascular and heart development in mice and zebrafish. Our study is the first report to show the association between PTPRB and HLHS in humans.

Clinical characteristics and in silico analysis of congenital pseudarthrosis of the tibia combined with neurofibromatosis type 1 caused by a novel NF1 mutation.

Congenital pseudarthrosis of the tibia (CPT) is a rare congenital bone malformation, which has a strong relationship with Neurofibromatosis type 1 (NF1). NF1 is an autosomal dominant disease leading to multisystem disorders. Here, we presented the genotypic and phenotypic characteristics of one unique case of a five-generation Chinese family. The proband was CPT accompanied with NF1 due to NF1 mutation. The proband developed severe early-onset CPT combined with NF1 after birth. Appearance photos and X-ray images of the left limb of the proband showed significant bone malformation. Slit-lamp examination showed Lisch nodules in both eyes of the proband. Whole-exome sequencing (WES) and Sanger sequencing confirmed the truncation variant of NF1 (c.871G>T, p. E291*). Sequence conservative and evolutionary conservation analysis indicated that the novel mutation (p.E291*) was highly conserved. The truncated mutation led to the loss of functional domains, including CSRD, GRD, TBD, SEC14-PH, CTD, and NLS. It may explain why the mutation led to a severe clinical feature. Our report expands the genotypic spectrum of NF1 mutations and the phenotypic spectrum of CPT combined with NF1.