A catalog of neutral and deleterious polymorphism in yeast

Initiating-clone analysis in patients with acute myeloid leukemia secondary to essential thrombocythemia

Most of essential thrombocythemia (ET) patients have the clone harboring a mutation in one of the JAK2, CALR, or MPL gene, and these clones generally acquire additional mutations at transformation to acute myeloid leukemia (AML). However, the proliferation of triple-negative clones has sometimes been observed at AML transformation. To clarify the clonal evolution of ET to AML, we analyzed paired samples at ET and AML transformation in eight patients. We identified that JAK2-unmutated AML clones proliferated at AML transformation in three patients in whom the JAK2-mutated clone was dominant at ET. In two patients, TET2-mutated, but not JAK2-mutated, clones might be common initiating clones for ET and transformed AML. In a patient with JAK2-mutated ET, SMARCC2, UBR4, and ZNF143, but not JAK2, -mutated clones proliferated at AML transformation. Precise analysis using single-cell sorted CD34+/CD38- fractions suggested that ET clone with JAK2-mutated and AML clone with TP53 mutation was derived from the common clone with these mutations. Although further study is required to clarify the biological significance of SMARCC2, UBR4, and ZNF143 mutations during disease progression of ET and AML transformation, the present results demonstrate the possibility of a common initial clone involved in both ET and transformed AML.

Application of whole exome sequencing in the diagnosis of muscular disorders: a study of Taiwanese pediatric patients

Background

Muscular dystrophies and congenital myopathies encompass various inherited muscular disorders that present diagnostic challenges due to clinical complexity and genetic heterogeneity.

Methods

This study aimed to investigate the use of whole exome sequencing (WES) in diagnosing muscular disorders in pediatric patients in Taiwan. Out of 161 pediatric patients suspected to have genetic/inherited myopathies, 115 received a molecular diagnosis through conventional tests, single gene testing, and gene panels. The remaining 46 patients were divided into three groups: Group 1 (multiplex ligation-dependent probe amplification-negative Duchenne muscular dystrophy) with three patients (6.5%), Group 2 (various forms of muscular dystrophies) with 21 patients (45.7%), and Group 3 (congenital myopathies) with 22 patients (47.8%).

Results

WES analysis of these groups found pathogenic variants in 100.0% (3/3), 57.1% (12/21), and 68.2% (15/22) of patients in Groups 1 to 3, respectively. WES had a diagnostic yield of 65.2% (30 patients out of 46), detecting 30 pathogenic or potentially pathogenic variants across 28 genes.

Conclusion

WES enables the diagnosis of rare diseases with symptoms and characteristics similar to congenital myopathies and muscular dystrophies, such as muscle weakness. Consequently, this approach facilitates targeted therapy implementation and appropriate genetic counseling.

From beer to breadboards: yeast as a force for biological innovation

The history of yeast Saccharomyces cerevisiae, aka brewer's or baker's yeast, is intertwined with our own. Initially domesticated 8,000 years ago to provide sustenance to our ancestors, for the past 150 years, yeast has served as a model research subject and a platform for technology. In this review, we highlight many ways in which yeast has served to catalyze the fields of functional genomics, genome editing, gene-environment interaction investigation, proteomics, and bioinformatics-emphasizing how yeast has served as a catalyst for innovation. Several possible futures for this model organism in synthetic biology, drug personalization, and multi-omics research are also presented.

In silico analysis of missense SNPs in TNFR1a and their possible therapeutic or pathogenic role in immune diseases

Tumor necrosis factor alpha (TNFa) is an inflammatory cytokine that is involved in the pathogenesis of various inflammatory disorders including rheumatoid arthritis. TNF-alpha receptor I (TNFR1a) is one of the receptors TNFa binds with for its activation. Any variation in this receptor might affect the role of TNFa in successive events. Amino acid residue substitutions might happen in TNFR1a through non-synonymous single nucleotide polymorphisms (nsSNPs) which may alter the functioning of TNFa, hence, identifying any such substitutions is of paramount significance. In this study, six nsSNPs at five different evolutionary conserved regions are predicted to be detrimental to the structure and/or function of TNFR1a by using numerous computational tools. Their 3D models are also proposed in this study. Besides, they were found to reduce the stability and affect the molecular mechanisms of this protein. Two contrasting possibilities might happen because of these substitutions. One, they might reduce the production of TNFa which is overexpressed in inflammatory diseases, hence can play therapeutic role in such diseases. Second, they might possibly hinder the apoptosis to occur which can effectuate the uncontrolled division of cells, hence can be pathogenic in diseases like cancer. Further investigations on these nsSNPs using animal models and at cellular level will open doors to understand the underlying mechanisms behind various diseases.

Gene-by-environment interactions influence the fitness cost of gene copy-number variation in yeast

Variation in gene copy number can alter gene expression and influence downstream phenotypes; thus copy-number variation provides a route for rapid evolution if the benefits outweigh the cost. We recently showed that genetic background significantly influences how yeast cells respond to gene overexpression, revealing that the fitness costs of copy-number variation can vary substantially with genetic background in a common-garden environment. But the interplay between copy-number variation tolerance and environment remains unexplored on a genomic scale. Here, we measured the tolerance to gene overexpression in four genetically distinct Saccharomyces cerevisiae strains grown under sodium chloride stress. Overexpressed genes that are commonly deleterious during sodium chloride stress recapitulated those commonly deleterious under standard conditions. However, sodium chloride stress uncovered novel differences in strain responses to gene overexpression. West African strain NCYC3290 and North American oak isolate YPS128 are more sensitive to sodium chloride stress than vineyard BC187 and laboratory strain BY4743. Consistently, NCYC3290 and YPS128 showed the greatest sensitivities to overexpression of specific genes. Although most genes were deleterious, hundreds were beneficial when overexpressed-remarkably, most of these effects were strain specific. Few beneficial genes were shared between the sodium chloride-sensitive isolates, implicating mechanistic differences behind their sodium chloride sensitivity. Transcriptomic analysis suggested underlying vulnerabilities and tolerances across strains, and pointed to natural copy-number variation of a sodium export pump that likely contributes to strain-specific responses to overexpression of other genes. Our results reveal extensive strain-by-environment interactions in the response to gene copy-number variation, raising important implications for the accessibility of copy-number variation-dependent evolutionary routes under times of stress.

Genome-wide probing of eukaryotic nascent RNA structure elucidates cotranscriptional folding and its antimutagenic effect

The transcriptional intermediates of RNAs fold into secondary structures with multiple regulatory roles, yet the details of such cotranscriptional RNA folding are largely unresolved in eukaryotes. Here, we present eSPET-seq (Structural Probing of Elongating Transcripts in eukaryotes), a method to assess the cotranscriptional RNA folding in Saccharomyces cerevisiae. Our study reveals pervasive structural transitions during cotranscriptional folding and overall structural similarities between nascent and mature RNAs. Furthermore, a combined analysis with genome-wide R-loop and mutation rate approximations provides quantitative evidence for the antimutator effect of nascent RNA folding through competitive inhibition of the R-loops, known to facilitate transcription-associated mutagenesis. Taken together, we present an experimental evaluation of cotranscriptional folding in eukaryotes and demonstrate the antimutator effect of nascent RNA folding. These results suggest genome-wide coupling between the processing and transmission of genetic information through RNA folding.

Role of genetic architecture in phenotypic plasticity

Phenotypic plasticity, the ability of an organism to display different phenotypes across environments, is widespread in nature. Plasticity aids survival in novel environments. Herein, we review studies from yeast that allow us to start uncovering the genetic architecture of phenotypic plasticity. Genetic variants and their interactions impact the phenotype in different environments, and distinct environments modulate the impact of genetic variants and their interactions on the phenotype. Because of this, certain hidden genetic variation is expressed in specific genetic and environmental backgrounds. A better understanding of the genetic mechanisms of phenotypic plasticity will help to determine short- and long-term responses to selection and how wide variation in disease manifestation occurs in human populations.

Identification genetic variations in some heat shock protein genes of Tali goat breed and study their structural and functional effects on relevant proteins

BACKGROUND

Animals of different regions have adapted to adverse environmental conditions by modifying their phenotypic and genotypic characteristics in the long run.

OBJECTIVES

In this study, the effect of genetic variations of 10 heat shock protein (HSP) genes (HSP70A4, HSP70A9, HSP40C17, HSP40C27, HSP90AA1, HSP90AB1, HSPB7, HSPB11, HSPD1 and HSPE1) on the three-dimensional protein structure and function of proteins in Tali goat (a tropical breed) were studied and were compared with Saanen goat (as a sensitive breed).

METHODS

A pooled DNA of 15 samples from blood was sequenced and mapped to the goat reference sequence. The bioinformatics analysis was used to identify nsSNPs in the Tali breed and was compared with the Saanen goat. Four online bioinformatics tools (Sorting Intolerant from Tolerant, Protein Variation Effect Analyzer, Polymorphism Phenotyping version2 and Single Nucleotide Polymorphism Database and Gene Ontology) showed three deleterious missense nsSNPs and seven natural missense SNPs in these HSPs genes of Tali goat.

RESULTS

Out of 10 reported nsSNPs, 5 nsSNPs in HSP70A4, 1 nsSNP inHSP70A9, 2 nsSNPs in HSP40C17, 1 nsSNP in HSP40C27 and 1 nsSNP in HSPD1 were detected. ConSurf tools showed that the majority of the predicted nsSNPs occur in conserved sites. Moreover, several post-translational modification (PTM) predictors computed the probability of post-translation change of nsSNPs. The putative phosphorylation and glycosylation sites in HSPs proteins were substitutions rs669769139 and rs666336692 of the Tali goat breed.

CONCLUSION

These results on the effect of type of genetic variants on the function of HSP proteins will assist to predict the resistance to hard conditions in goat breeds. Considering that the identified SNPid rs669769139 (S248) which is located on the N-terminal ATPase domain of HSP70A4 is a PTM site with a highly conserved score and a natural substitution on changing the stability and benign protein that can affect the functional and structural characterization of HSPs protein for adaptation to the local climate.

Gene-by-environment interactions influence the fitness cost of gene copy-number variation in yeast

Variation in gene copy number can alter gene expression and influence downstream phenotypes; thus copy-number variation (CNV) provides a route for rapid evolution if the benefits outweigh the cost. We recently showed that genetic background significantly influences how yeast cells respond to gene over-expression (OE), revealing that the fitness costs of CNV can vary substantially with genetic background in a common-garden environment. But the interplay between CNV tolerance and environment remains unexplored on a genomic scale. Here we measured the tolerance to gene OE in four genetically distinct Saccharomyces cerevisiae strains grown under sodium chloride (NaCl) stress. OE genes that are commonly deleterious during NaCl stress recapitulated those commonly deleterious under standard conditions. However, NaCl stress uncovered novel differences in strain responses to gene OE. West African strain NCYC3290 and North American oak isolate YPS128 are more sensitive to NaCl stress than vineyard BC187 and laboratory strain BY4743. Consistently, NCYC3290 and YPS128 showed the greatest sensitivities to gene OE. Although most genes were deleterious, hundreds were beneficial when overexpressed - remarkably, most of these effects were strain specific. Few beneficial genes were shared between the NaCl-sensitive isolates, implicating mechanistic differences behind their NaCl sensitivity. Transcriptomic analysis suggested underlying vulnerabilities and tolerances across strains, and pointed to natural CNV of a sodium export pump that likely contributes to strain-specific responses to OE of other genes. Our results reveal extensive strain-by-environment interaction in the response to gene CNV, raising important implications for the accessibility of CNV-dependent evolutionary routes under times of stress.

In-silico assessment of high-risk non-synonymous SNPs in ADAMTS3 gene associated with Hennekam syndrome and their impact on protein stability and function

Hennekam Lymphangiectasia-Lymphedema Syndrome 3 (HKLLS3) is a rare genetical disorder caused by mutations in a few genes including ADAMTS3. It is characterized by lymphatic dysplasia, intestinal lymphangiectasia, severe lymphedema and distinctive facial appearance. Up till now, no extensive studies have been conducted to elucidate the mechanism of the disease caused by various mutations. As a preliminary investigation of HKLLS3, we sorted out the most deleterious nonsynonymous single nucleotide polymorphisms (nsSNPs) that might affect the structure and function of ADAMTS3 protein by using a variety of in silico tools. A total of 919 nsSNPs in the ADAMTS3 gene were identified. 50 nsSNPs were predicted to be deleterious by multiple computational tools. 5 nsSNPs (G298R, C567Y, A370T, C567R and G374S) were found to be the most dangerous and can be associated with the disease as predicted by different bioinformatics tools. Modelling of the protein shows it can be divided into segments 1, 2 and 3, which are connected by short loops. Segment 3 mainly consists of loops without substantial secondary structures. With prediction tools and molecular dynamics simulation, some SNPs were found to significantly destabilize the protein structure and disrupt the secondary structures, especially in segment 2. The deleterious effects of mutations in segment 1 are possibly not from destabilization but from other factors such as the change in phosphorylation as suggested by post-translational modification (PTM) studies. This is the first-ever study of ADAMTS3 gene polymorphism, and the predicted nsSNPs in ADAMST3, some of which have not been reported yet in patients, will serve for diagnostic purposes and further therapeutic implications in Hennekam syndrome, contributing to better diagnosis and treatment.

In Silico Analysis of nsSNPs of Human KRAS Gene and Protein Modeling Using Bioinformatic Tools

The KRAS gene belongs to the RAS family and codes for 188 amino acid residues of KRAS protein, with a molecular mass of 21.6 kD. Non-synonymous single-nucleotide polymorphisms (nsSNPs) have been identified within the coding region in which some are associated with different diseases. However, structural changes are not well defined yet. In this study, we first categorized SNPs in the KRAS coding area and then used computational methods to determine their impact on the protein structure and stability. In addition, the three-dimensional model of KRAS was taken from the Protein Data Bank for structural modeling. Furthermore, genomic data were extracted from a variety of sources, including the 1000 Genome Project, dbSNPs, and ENSEMBLE, and assessed through in silico methods. Based on various tools used in this study, 10 out of 48 missense SNPs with rsIDs were found deleterious. The substitution of alanine for proline at position 146 pushed several residues toward the center of the protein. Arginine instead of leucine has a minor effect on protein structure and stability. In addition, the substitution of proline for leucine at the 34th position disrupted the structure and led to a bigger size than the wild-type protein, hence interrupting the protein interaction. Using the well-intended computational approach and applying several bioinformatic tools, we characterized and identified most damaging nsSNPs and further explored the structural dynamics and stability of KRAS protein.

UGDR: a generic pipeline to detect recombined regions in polyploid and complex hybrid yeast genomes

MOTIVATION

In eukaryotes, homologous recombination between the parental genomes frequently occurs during the evolutionary conserved process of meiosis, generating the genetic diversity transmitted by the gametes. The genome-wide determination of the frequency and location of the recombination events can now be efficiently performed by genotyping the offspring's polymorphic markers. However, genotyping recombination in complex hybrid genomes with existing methods remains challenging because of their strain and ploidy specificity and the degree of diversity and complexity of the parental genomes, especially in [Formula: see text] polyploids.

RESULTS

We present UGDR, a pipeline to genotype the polymorphisms of complex hybrid yeast genomes. It is based on optimal mapping strategies of NGS reads, comparative analyses of the allelic ratio variation and read depth coverage. We tested the UGDR pipeline with sequencing reads from recombined hybrid diploid yeast strains and various clinical strains exhibiting different degrees of ploidy. UGDR allows to plot the markers distribution and recombination profile per chromosome.

CONCLUSION

UGDR detects and plots recombination events in haploids and polyploid yeasts, which facilitates the discovery and understanding of the yeast genetic recombination map and identify new out-performing recombinants.

In silico investigation of nonsynonymous single nucleotide polymorphisms in BCL2 apoptosis regulator gene to design novel protein-based drugs against cancer

BCL2 apoptosis regulator gene encodes Bcl-2 pro-survival protein, which plays an important role to evade apoptosis in various cancers. Moreover, single nucleotide polymorphisms (SNPs) in the BCL2 gene can be nonsynonymous (nsSNPs), which might affect the protein stability and probably its function. Therefore, we implement cutting-edge computational techniques based on the Spherical Polar Fourier and Monte-Carlo algorithms to investigate the impact of these SNPs on the B cell lymphoma-2 (Bcl-2) stability and therapeutic potential of protein-based molecules to inhibit this protein. As a result, we identified two nsSNPs (Q118R and R129C) to be deleterious and highly conserved, having a negative effect on protein stability. Additionally, molecular docking and molecular dynamics simulations confirmed the decreased binding affinity of mutated Bcl-2 variants to bind three-helix bundle protein inhibitor as these mutations occurred in the protein-protein binding site. Overall, this computational approach investigating nsSNPs provides a useful basis for designing novel molecules to inhibit Bcl-2 pro-survival pathway in malignant cells.

Diagnosis of a Single-Nucleotide Variant in Whole-Exome Sequencing Data for Patients With Inherited Diseases: Machine Learning Study Using Artificial Intelligence Variant Prioritization

BACKGROUND

In recent years, thanks to the rapid development of next-generation sequencing (NGS) technology, an entire human genome can be sequenced in a short period. As a result, NGS technology is now being widely introduced into clinical diagnosis practice, especially for diagnosis of hereditary disorders. Although the exome data of single-nucleotide variant (SNV) can be generated using these approaches, processing the DNA sequence data of a patient requires multiple tools and complex bioinformatics pipelines.

OBJECTIVE

This study aims to assist physicians to automatically interpret the genetic variation information generated by NGS in a short period. To determine the true causal variants of a patient with genetic disease, currently, physicians often need to view numerous features on every variant manually and search for literature in different databases to understand the effect of genetic variation.

METHODS

We constructed a machine learning model for predicting disease-causing variants in exome data. We collected sequencing data from whole-exome sequencing (WES) and gene panel as training set, and then integrated variant annotations from multiple genetic databases for model training. The model built ranked SNVs and output the most possible disease-causing candidates. For model testing, we collected WES data from 108 patients with rare genetic disorders in National Taiwan University Hospital. We applied sequencing data and phenotypic information automatically extracted by a keyword extraction tool from patient's electronic medical records into our machine learning model.

RESULTS

We succeeded in locating 92.5% (124/134) of the causative variant in the top 10 ranking list among an average of 741 candidate variants per person after filtering. AI Variant Prioritizer was able to assign the target gene to the top rank for around 61.1% (66/108) of the patients, followed by Variant Prioritizer, which assigned it for 44.4% (48/108) of the patients. The cumulative rank result revealed that our AI Variant Prioritizer has the highest accuracy at ranks 1, 5, 10, and 20. It also shows that AI Variant Prioritizer presents better performance than other tools. After adopting the Human Phenotype Ontology (HPO) terms by looking up the databases, the top 10 ranking list can be increased to 93.5% (101/108).

CONCLUSIONS

We successfully applied sequencing data from WES and free-text phenotypic information of patient's disease automatically extracted by the keyword extraction tool for model training and testing. By interpreting our model, we identified which features of variants are important. Besides, we achieved a satisfactory result on finding the target variant in our testing data set. After adopting the HPO terms by looking up the databases, the top 10 ranking list can be increased to 93.5% (101/108). The performance of the model is similar to that of manual analysis, and it has been used to help National Taiwan University Hospital with a genetic diagnosis.

Gene loss and compensatory evolution promotes the emergence of morphological novelties in budding yeast

Deleterious mutations are generally considered to be irrelevant for morphological evolution. However, they could be compensated by conditionally beneficial mutations, thereby providing access to new adaptive paths. Here we use high-dimensional phenotyping of laboratory-evolved budding yeast lineages to demonstrate that new cellular morphologies emerge exceptionally rapidly as a by-product of gene loss and subsequent compensatory evolution. Unexpectedly, the capacities for invasive growth, multicellular aggregation and biofilm formation also spontaneously evolve in response to gene loss. These multicellular phenotypes can be achieved by diverse mutational routes and without reactivating the canonical regulatory pathways. These ecologically and clinically relevant traits originate as pleiotropic side effects of compensatory evolution and have no obvious utility in the laboratory environment. The extent of morphological diversity in the evolved lineages is comparable to that of natural yeast isolates with diverse genetic backgrounds and lifestyles. Finally, we show that both the initial gene loss and subsequent compensatory mutations contribute to new morphologies, with their synergistic effects underlying specific morphological changes. We conclude that compensatory evolution is a previously unrecognized source of morphological diversity and phenotypic novelties.

Idiosyncratic epistasis leads to global fitness-correlated trends

Epistasis can markedly affect evolutionary trajectories. In recent decades, protein-level fitness landscapes have revealed extensive idiosyncratic epistasis among specific mutations. By contrast, other work has found ubiquitous and apparently nonspecific patterns of global diminishing-returns and increasing-costs epistasis among mutations across the genome. Here, we used a hierarchical CRISPR gene drive system to construct all combinations of 10 missense mutations from across the genome in budding yeast and measured their fitness in six environments. We show that the resulting fitness landscapes exhibit global fitness-correlated trends but that these trends emerge from specific idiosyncratic interactions. We thus provide experimental validation of recent theoretical work arguing that fitness-correlated trends can emerge as the generic consequence of idiosyncratic epistasis.

Novel Disease-Associated Missense Single-Nucleotide Polymorphisms Variants Predication by Algorithms Tools and Molecular Dynamics Simulation of Human TCIRG1 Gene Causing Congenital Neutropenia and Osteopetrosis

T Cell Immune Regulator 1, ATPase H + Transporting V0 Subunit A3 (TCIRG1 gene provides instructions for making one part, the a3 subunit, of a large protein complex known as a vacuolar H + -ATPase (V-ATPase). V-ATPases are a group of similar complexes that act as pumps to move positively charged hydrogen atoms (protons) across membranes. Single amino acid changes in highly conserved areas of the TCIRG1 protein have been linked to autosomal recessive osteopetrosis and severe congenital neutropenia. We used multiple computational approaches to classify disease-prone single nucleotide polymorphisms (SNPs) in TCIRG1. We used molecular dynamics analysis to identify the deleterious nsSNPs, build mutant protein structures, and assess the impact of mutation. Our results show that fifteen nsSNPs (rs199902030, rs200149541, rs372499913, rs267605221, rs374941368, rs375717418, rs80008675, rs149792489, rs116675104, rs121908250, rs121908251, rs121908251, rs149792489 and rs116675104) variants are likely to be highly deleterious mutations as by incorporating them into wild protein they destabilize the wild protein structure and function. They are also located in the V-ATPase I domain, which may destabilize the structure and impair TCIRG1 protein activation, as well as reduce its ATPase effectiveness. These mutants have not yet been identified in patients suffering from CN and osteopetrosis while (G405R, R444L, and D517N) reported in our study are already associated with osteopetrosis. Mutation V52L reported in our study was identified in a patient suspected for CN. Finally, these mutants can help to further understand the broad pool of illness susceptibilities associated with TCIRG1 catalytic kinase domain activation and aid in the development of an effective treatment for associated diseases.

In Silico Analysis Identified Putative Pathogenic Missense nsSNPs in Human SLITRK1 Gene

Human DNA contains several variations, which can affect the structure and normal functioning of a protein. These variations could be single nucleotide polymorphisms (SNPs) or insertion-deletions (InDels). SNPs, as opposed to InDels, are more commonly present in DNA and may cause genetic disorders. In the current study, several bioinformatic tools were used to prioritize the pathogenic variants in the SLITRK1 gene. Out of all of the variants, 16 were commonly predicted to be pathogenic by these tools. All the variants had very low frequency, i.e., <0.0001 in the global population. The secondary structure of all filtered variants was predicted, but no structural change was observed at the site of variation in any variant. Protein stability analysis of these variants was then performed, which determined a decrease in protein stability of 10 of the variants. Amino acid conservation analysis revealed that all the amino acids were highly conserved, indicating their structural and functional importance. Protein 3D structure of wildtype SLITRK1 and all of its variants was predicted using I-TASSER, and the effect of variation on 3D structure of the protein was observed using the Missense3D tool, which presented the probable structural loss in three variants, i.e., Asn529Lys, Leu496Pro and Leu94Phe. The wildtype SLITRK1 protein and these three variants were independently docked with their close interactor protein PTPRD, and remarkable differences were observed in the docking sites of normal and variants, which will ultimately affect the functional activity of the SLITRK1 protein. Previous studies have shown that mutations in SLITRK1 are involved in Tourette syndrome. The present study may assist a molecular geneticist in interpreting the variant pathogenicity in research as well as diagnostic setup.

Non-genetic and genetic rewiring underlie adaptation to hypomorphic alleles of an essential gene

Adaptive evolution to cellular stress is a process implicated in a wide range of biological and clinical phenomena. Two major routes of adaptation have been identified: non-genetic changes, which allow expression of different phenotypes in novel environments, and genetic variation achieved by selection of fitter phenotypes. While these processes are broadly accepted, their temporal and epistatic features in the context of cellular evolution and emerging drug resistance are contentious. In this manuscript, we generated hypomorphic alleles of the essential nuclear pore complex (NPC) gene NUP58. By dissecting early and long-term mechanisms of adaptation in independent clones, we observed that early physiological adaptation correlated with transcriptome rewiring and upregulation of genes known to interact with the NPC; long-term adaptation and fitness recovery instead occurred via focal amplification of NUP58 and restoration of mutant protein expression. These data support the concept that early phenotypic plasticity allows later acquisition of genetic adaptations to a specific impairment. We propose this approach as a genetic model to mimic targeted drug therapy in human cells and to dissect mechanisms of adaptation.

Computational algorithmic and molecular dynamics study of functional and structural impacts of non-synonymous single nucleotide polymorphisms in human DHFR gene

Human dihydrofolate reductase (DHFR) is a conserved enzyme that is central to folate metabolism and is widely targeted in pathogenic diseases as well as cancers. Although studies have reported the fact that genetic mutations in DHFR leads to a rare autosomal recessive inborn error of folate metabolism and drug resistance, there is a lack of an extensive study on how the deleterious non-synonymous SNPs (nsSNPs) disrupt its phenotypic effects. In this study, we aim at discovering the structural and functional consequences of nsSNPs in DHFR by employing a combined computational approach consisting of ten recently developed in silico tools for identification of damaging nsSNPs and molecular dynamics (MD) simulation for getting deeper insights into the magnitudes of damaging effects. Our study revealed the presence of 12 most deleterious nsSNPs affecting the native phenotypic effects, with three (R71T, G118D, Y122D) identified in the co-factor and ligand binding active sites. MD simulations also suggested that these three SNPs particularly Y122D, alter the overall structural flexibility and dynamics of the native DHFR protein which can provide more understandings into the crucial roles of these mutants in influencing the loss of DHFR function.

Pathogenic genetic variants from highly connected cancer susceptibility genes confer the loss of structural stability

Genetic polymorphisms in DNA damage repair and tumor suppressor genes have been associated with increasing the risk of several types of cancer. Analyses of putative functional single nucleotide polymorphisms (SNP) in such genes can greatly improve human health by guiding choice of therapeutics. In this study, we selected nine genes responsible for various cancer types for gene enrichment analysis and found that BRCA1, ATM, and TP53 were more enriched in connectivity. Therefore, we used different computational algorithms to classify the nonsynonymous SNPs which are deleterious to the structure and/or function of these three proteins. The present study showed that the major pathogenic variants (V1687G and V1736G of BRCA1, I2865T and V2906A of ATM, V216G and L194H of TP53) might have a greater impact on the destabilization of the proteins. To stabilize the high-risk SNPs, we performed mutation site-specific molecular docking analysis and validated using molecular dynamics (MD) simulation and molecular mechanics/Poisson Boltzmann surface area (MM/PBSA) studies. Additionally, SNPs of untranslated regions of these genes affecting miRNA binding were characterized. Hence, this study will assist in developing precision medicines for cancer types related to these polymorphisms.

Cis-regulatory variants affect gene expression dynamics in yeast

Evolution of cis-regulatory sequences depends on how they affect gene expression and motivates both the identification and prediction of cis-regulatory variants responsible for expression differences within and between species. While much progress has been made in relating cis-regulatory variants to expression levels, the timing of gene activation and repression may also be important to the evolution of cis-regulatory sequences. We investigated allele-specific expression (ASE) dynamics within and between Saccharomyces species during the diauxic shift and found appreciable cis-acting variation in gene expression dynamics. Within-species ASE is associated with intergenic variants, and ASE dynamics are more strongly associated with insertions and deletions than ASE levels. To refine these associations, we used a high-throughput reporter assay to test promoter regions and individual variants. Within the subset of regions that recapitulated endogenous expression, we identified and characterized cis-regulatory variants that affect expression dynamics. Between species, chimeric promoter regions generate novel patterns and indicate constraints on the evolution of gene expression dynamics. We conclude that changes in cis-regulatory sequences can tune gene expression dynamics and that the interplay between expression dynamics and other aspects of expression is relevant to the evolution of cis-regulatory sequences.

Genetic interaction network has a very limited impact on the evolutionary trajectories in continuous culture-grown populations of yeast

Background

The impact of genetic interaction networks on evolution is a fundamental issue. Previous studies have demonstrated that the topology of the network is determined by the properties of the cellular machinery. Functionally related genes frequently interact with one another, and they establish modules, e.g., modules of protein complexes and biochemical pathways. In this study, we experimentally tested the hypothesis that compensatory evolutionary modifications, such as mutations and transcriptional changes, occur frequently in genes from perturbed modules of interacting genes.

Results

Using Saccharomyces cerevisiae haploid deletion mutants as a model, we investigated two modules lacking COG7 or NUP133, which are evolutionarily conserved genes with many interactions. We performed laboratory evolution experiments with these strains in two genetic backgrounds (with or without additional deletion of MSH2), subjecting them to continuous culture in a non-limiting minimal medium. Next, the evolved yeast populations were characterized through whole-genome sequencing and transcriptome analyses. No obvious compensatory changes resulting from inactivation of genes already included in modules were identified. The supposedly compensatory inactivation of genes in the evolved strains was only rarely observed to be in accordance with the established fitness effect of the genetic interaction network. In fact, a substantial majority of the gene inactivations were predicted to be neutral in the experimental conditions used to determine the interaction network. Similarly, transcriptome changes during continuous culture mostly signified adaptation to growth conditions rather than compensation of the absence of the COG7, NUP133 or MSH2 genes. However, we noticed that for genes whose inactivation was deleterious an upregulation of transcription was more common than downregulation.

Conclusions

Our findings demonstrate that the genetic interactions and the modular structure of the network described by others have very limited effects on the evolutionary trajectory following gene deletion of module elements in our experimental conditions and has no significant impact on short-term compensatory evolution. However, we observed likely compensatory evolution in functionally related (albeit non-interacting) genes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01830-9.

Pervasive Phenotypic Impact of a Large Nonrecombining Introgressed Region in Yeast

To explore the origin of the diversity observed in natural populations, many studies have investigated the relationship between genotype and phenotype. In yeast species, especially in Saccharomyces cerevisiae, these studies are mainly conducted using recombinant offspring derived from two genetically diverse isolates, allowing to define the phenotypic effect of genetic variants. However, large genomic variants such as interspecies introgressions are usually overlooked even if they are known to modify the genotype–phenotype relationship. To have a better insight into the overall phenotypic impact of introgressions, we took advantage of the presence of a 1-Mb introgressed region, which lacks recombination and contains the mating-type determinant in the Lachancea kluyveri budding yeast. By performing linkage mapping analyses in this species, we identified a total of 89 loci affecting growth fitness in a large number of conditions and 2,187 loci affecting gene expression mostly grouped into two major hotspots, one being the introgressed region carrying the mating-type locus. Because of the absence of recombination, our results highlight the presence of a sexual dimorphism in a budding yeast for the first time. Overall, by describing the phenotype–genotype relationship in the Lachancea kluyveri species, we expanded our knowledge on how genetic characteristics of large introgression events can affect the phenotypic landscape.

Autosomal dominant diabetes associated with a novel ZYG11A mutation resulting in cell cycle arrest in beta-cells

Diabetes is a genetically heterogeneous disease, for which we are aiming to identify causative genes. Here, we report a missense mutation (c.T1424C:p.L475P) in ZYG11A identified by exome sequencing as segregating with hyperglycemia in a Thai family with autosomal dominant diabetes. ZYG11A functions as a target recruitment subunit of an E3 ubiquitin ligase complex that plays an important role in the regulation of cell cycle. We demonstrate an increase in cells arrested at G₂/mitotic phase among beta-cells deficient for ZYG11A or overexpressing L475P-ZYG11A, which is associated with a decreased growth rate. This is the first evidence linking a ZYG11A mutation to hyperglycemia, and suggesting ZYG11A as a cell cycle regulator required for beta-cell growth. Since most family members were either overweight or obese, but only mutation carriers developed hyperglycemia, our data also suggests the ZYG11A mutation as a genetic factor predisposing obese individuals to beta-cell failure in maintenance of glucose homeostasis.

Molecular and evolutionary processes generating variation in gene expression

Heritable variation in gene expression is common within and between species. This variation arises from mutations that alter the form or function of molecular gene regulatory networks that are then filtered by natural selection. High-throughput methods for introducing mutations and characterizing their cis- and trans-regulatory effects on gene expression (particularly, transcription) are revealing how different molecular mechanisms generate regulatory variation, and studies comparing these mutational effects with variation seen in the wild are teasing apart the role of neutral and non-neutral evolutionary processes. This integration of molecular and evolutionary biology allows us to understand how the variation in gene expression we see today came to be and to predict how it is most likely to evolve in the future.

Comprehensive in-silico analysis of damage associated SNPs in hOCT1 affecting Imatinib response in chronic myeloid leukemia

Non-synonymous single nucleotide polymorphisms (nsSNPs) in hOCT1 (encoded by SLC22A1 gene) are expected to affect Imatinib uptake in chronic myeloid leukemia (CML). In this study, sequence homology-based genetic analysis of a set of 270 coding SNPs identified 18 nsSNPs to be putatively damaging/deleterious using eight different algorithms. Subsequently, based on conservation of amino acid residues, stability analysis, posttranscriptional modifications, and solvent accessibility analysis, the possible structural-functional relationship was established for high-confidence nsSNPs. Furthermore, based on the modeling results, some dissimilarities of mutant type amino acids from wild-type amino acids such as size, charge, interaction and hydrophobicity were revealed. Three highly deleterious mutations consisting of P283L, G401S and R402G in SLC22A1 gene that may alter the protein structure, function and stability were identified. These results provide a filtered data to explore the effect of uncharacterized nsSNP and find their association with Imatinib resistance in CML.

In silico analysis predicting effects of deleterious SNPs of human RASSF5 gene on its structure and functions

Ras association domain-containing protein 5 (RASSF5), one of the prospective biomarkers for tumors, generally plays a crucial role as a tumor suppressor. As deleterious effects can result from functional differences through SNPs, we sought to analyze the most deleterious SNPs of RASSF5 as well as predict the structural changes associated with the mutants that hamper the normal protein-protein interactions. We adopted both sequence and structure based approaches to analyze the SNPs of RASSF5 protein. We also analyzed the putative post translational modification sites as well as the altered protein-protein interactions that encompass various cascades of signals. Out of all the SNPs obtained from the NCBI database, only 25 were considered as highly deleterious by six in silico SNP prediction tools. Among them, upon analyzing the effect of these nsSNPs on the stability of the protein, we found 17 SNPs that decrease the stability. Significant deviation in the energy minimization score was observed in P350R, F321L, and R277W. Besides this, docking analysis confirmed that P350R, A319V, F321L, and R277W reduce the binding affinity of the protein with H-Ras, where P350R shows the most remarkable deviation. Protein-protein interaction analysis revealed that RASSF5 acts as a hub connecting two clusters consisting of 18 proteins and alteration in the RASSF5 may lead to disassociation of several signal cascades. Thus, based on these analyses, our study suggests that the reported functional SNPs may serve as potential targets for different proteomic studies, diagnosis and therapeutic interventions.

Chromosomal-level genome assembly of the scimitar-horned oryx: Insights into diversity and demography of a species extinct in the wild

Captive populations provide a valuable insurance against extinctions in the wild. However, they are also vulnerable to the negative impacts of inbreeding, selection and drift. Genetic information is therefore considered a critical aspect of conservation management. Recent developments in sequencing technologies have the potential to improve the outcomes of management programmes; however, the transfer of these approaches to applied conservation has been slow. The scimitar-horned oryx (Oryx dammah) is a North African antelope that has been extinct in the wild since the early 1980s and is the focus of a large-scale and long-term reintroduction project. To enable the selection of suitable founder individuals, facilitate post-release monitoring and improve captive breeding management, comprehensive genomic resources are required. Here, we used 10X Chromium sequencing together with Hi-C contact mapping to develop a chromosomal-level genome assembly for the species. The resulting assembly contained 29 chromosomes with a scaffold N50 of 100.4 Mb, and displayed strong chromosomal synteny with the cattle genome. Using resequencing data from six additional individuals, we demonstrated relatively high genetic diversity in the scimitar-horned oryx compared to other mammals, despite it having experienced a strong founding event in captivity. Additionally, the level of diversity across populations varied according to management strategy. Finally, we uncovered a dynamic demographic history that coincided with periods of climate variation during the Pleistocene. Overall, our study provides a clear example of how genomic data can uncover valuable insights into captive populations and contributes important resources to guide future management decisions of an endangered species.

Resistance Mechanisms of Saccharomyces cerevisiae to Commercial Formulations of Glyphosate Involve DNA Damage Repair, the Cell Cycle, and the Cell Wall Structure

The use of glyphosate-based herbicides is widespread and despite their extensive use, their effects are yet to be deciphered completely. The additives in commercial formulations of glyphosate, though labeled inert when used individually, have adverse effects when used in combination with other additives along with the active ingredient. As a species, Saccharomyces cerevisiae has a wide range of resistance to glyphosate-based herbicides. To investigate the underlying genetic differences between sensitive and resistant strains, global changes in gene expression were measured, when yeast were exposed to a glyphosate-based herbicide (GBH). Expression of genes involved in numerous pathways crucial to the cell's functioning, such as DNA replication, MAPK signaling, meiosis, and cell wall synthesis changed. Because so many diverse pathways were affected, these strains were then subjected to in-lab-evolutions (ILE) to select mutations that confer increased resistance. Common fragile sites were found to play a role in adaptation to resistance to long-term exposure of GBHs. Copy number increased in approximately 100 genes associated with cell wall proteins, mitochondria, and sterol transport. Taking ILE and transcriptomic data into account it is evident that GBHs affect multiple biological processes in the cell. One such component is the cell wall structure which acts as a protective barrier in alleviating the stress caused by exposure to inert additives in GBHs. Sed1, a GPI-cell wall protein, plays an important role in tolerance of a GBH. Hence, a detailed study of the changes occurring at the genome and transcriptome levels is essential to better understand the effects of an environmental stressor such as a GBH, on the cell as a whole.

Evolution of Ty1 copy number control in yeast by horizontal transfer and recombination

Transposable elements constitute a large fraction of most eukaryotic genomes. Insertion of mobile DNA sequences typically has deleterious effects on host fitness, and thus diverse mechanisms have evolved to control mobile element proliferation. Mobility of the Ty1 retrotransposon in Saccharomyces yeasts is regulated by copy number control (CNC) mediated by a self-encoded restriction factor derived from the Ty1 gag capsid gene that inhibits virus-like particle function. Here, we survey a panel of wild and human-associated strains of S. cerevisiae and S. paradoxus to investigate how genomic Ty1 content influences variation in Ty1 mobility. We observe high levels of mobility for a tester element with a gag sequence from the canonical Ty1 subfamily in permissive strains that either lack full-length Ty1 elements or only contain full-length copies of the Ty1' subfamily that have a divergent gag sequence. In contrast, low levels of canonical Ty1 mobility are observed in restrictive strains carrying full-length Ty1 elements containing a canonical gag sequence. Phylogenomic analysis of full-length Ty1 elements revealed that Ty1' is the ancestral subfamily present in wild strains of S. cerevisiae, and that canonical Ty1 in S. cerevisiae is a derived subfamily that acquired gag from S. paradoxus by horizontal transfer and recombination. Our results provide evidence that variation in the ability of S. cerevisiae and S. paradoxus strains to repress canonical Ty1 transposition via CNC is regulated by the genomic content of different Ty1 subfamilies, and that self-encoded forms of transposon control can spread across species boundaries by horizontal transfer.

pCADD: SNV prioritisation in Sus scrofa

BACKGROUND

In animal breeding, identification of causative genetic variants is of major importance and high economical value. Usually, the number of candidate variants exceeds the number of variants that can be validated. One way of prioritizing probable candidates is by evaluating their potential to have a deleterious effect, e.g. by predicting their consequence. Due to experimental difficulties to evaluate variants that do not cause an amino-acid substitution, other prioritization methods are needed. For human genomes, the prediction of deleterious genomic variants has taken a step forward with the introduction of the combined annotation dependent depletion (CADD) method. In theory, this approach can be applied to any species. Here, we present pCADD (p for pig), a model to score single nucleotide variants (SNVs) in pig genomes.

RESULTS

To evaluate whether pCADD captures sites with biological meaning, we used transcripts from miRNAs and introns, sequences from genes that are specific for a particular tissue, and the different sites of codons, to test how well pCADD scores differentiate between functional and non-functional elements. Furthermore, we conducted an assessment of examples of non-coding and coding SNVs, which are causal for changes in phenotypes. Our results show that pCADD scores discriminate between functional and non-functional sequences and prioritize functional SNVs, and that pCADD is able to score the different positions in a codon relative to their redundancy. Taken together, these results indicate that based on pCADD scores, regions with biological relevance can be identified and distinguished according to their rate of adaptation.

CONCLUSIONS

We present the ability of pCADD to prioritize SNVs in the pig genome with respect to their putative deleteriousness, in accordance to the biological significance of the region in which they are located. We created scores for all possible SNVs, coding and non-coding, for all autosomes and the X chromosome of the pig reference sequence Sscrofa11.1, proposing a toolbox to prioritize variants and evaluate sequences to highlight new sites of interest to explain biological functions that are relevant to animal breeding.

Recombination-induced revertant mosaicism in ichthyosis with confetti and loricrin keratoderma

Revertant mosaicism refers to a condition in which a pathogenic germline mutation is spontaneously corrected in somatic cells, resulting in the presence of two or more cell populations with different genotypes in an organism arising from a single fertilized egg. If the revertant cells are clonally expanded due to a survival advantage over the surrounding mutant cells, patients benefit from this self-healing phenomenon which leads to the development of milder-than-expected clinical phenotypes; in genetic skin diseases, patients with revertant mosaicism present with small islands of healthy skin. To date, revertant mosaicism has been reported in ∼50 genetic diseases involving the skin, blood, liver, muscle, and brain. In this review, I briefly summarize current knowledge on revertant mosaicism in two particular skin diseases, ichthyosis with confetti (IWC) and loricrin keratoderma (LK), both of which develop numerous revertant skin patches. Notably, homologous recombination (HR) is the only mechanism underlying the reversion of pathogenic mutations in IWC and LK, and this was identified following the analysis of ∼50 revertant epidermis samples. All the samples showed long-tract loss of heterozygosity (LOH) that originated at regions centromeric to pathogenic mutations and extended to the telomere of the mutation-harboring chromosomes. Elucidating the molecular mechanisms underlying revertant mosaicism in IWC and LK-especially how mutant proteins induce long-tract LOH-would potentially expand the possibility of manipulating HR to induce the reversion of disease-causing mutations and help devising novel therapies not only for IWC and LK but also for other intractable genetic diseases.

Bioinformatics classification of mutations in patients with Mucopolysaccharidosis IIIA

Mucopolysaccharidosis (MPS) IIIA, also known as Sanfilippo syndrome type A, is a severe, progressive disease that affects the central nervous system (CNS). MPS IIIA is inherited in an autosomal recessive manner and is caused by a deficiency in the lysosomal enzyme sulfamidase, which is required for the degradation of heparan sulfate. The sulfamidase is produced by the N-sulphoglucosamine sulphohydrolase (SGSH) gene. In MPS IIIA patients, the excess of lysosomal storage of heparan sulfate often leads to mental retardation, hyperactive behavior, and connective tissue impairments, which occur due to various known missense mutations in the SGSH, leading to protein dysfunction. In this study, we focused on three mutations (R74C, S66W, and R245H) based on in silico pathogenic, conservation, and stability prediction tool studies. The three mutations were further subjected to molecular dynamic simulation (MDS) analysis using GROMACS simulation software to observe the structural changes they induced, and all the mutants exhibited maximum deviation patterns compared with the native protein. Conformational changes were observed in the mutants based on various geometrical parameters, such as conformational stability, fluctuation, and compactness, followed by hydrogen bonding, physicochemical properties, principal component analysis (PCA), and salt bridge analyses, which further validated the underlying cause of the protein instability. Additionally, secondary structure and surrounding amino acid analyses further confirmed the above results indicating the loss of protein function in the mutants compared with the native protein. The present results reveal the effects of three mutations on the enzymatic activity of sulfamidase, providing a molecular explanation for the cause of the disease. Thus, this study allows for a better understanding of the effect of SGSH mutations through the use of various computational approaches in terms of both structure and functions and provides a platform for the development of therapeutic drugs and potential disease treatments.

Prediction of Deleterious Non-synonymous SNPs of Human STK11 Gene by Combining Algorithms, Molecular Docking, and Molecular Dynamics Simulation

Serine-threonine kinase11 (STK11) is a tumor suppressor gene which plays a key role in regulating cell growth and apoptosis. It is widely known as a multitasking kinase and engaged in cell polarity, cell cycle arrest, chromatin remodeling, energy metabolism, and Wnt signaling. The substitutions of single amino acids in highly conserved regions of the STK11 protein are associated with Peutz-Jeghers syndrome (PJS), which is an autosomal dominant inherited disorder. The abnormal function of the STK11 protein is still not well understood. In this study, we classified disease susceptible single nucleotide polymorphisms (SNPs) in STK11 by using different computational algorithms. We identified the deleterious nsSNPs, constructed mutant protein structures, and evaluated the impact of mutation by employing molecular docking and molecular dynamics analysis. Our results show that W239R and W308C variants are likely to be highly deleterious mutations found in the catalytic kinase domain, which may destabilize structure and disrupt the activation of the STK11 protein as well as reduce its catalytic efficiency. The W239R mutant is likely to have a greater impact on destabilizing the protein structure compared to the W308C mutant. In conclusion, these mutants can help to further realize the large pool of disease susceptibilities linked with catalytic kinase domain activation of STK11 and assist to develop an effective drug for associated diseases.

Multiple Changes Underlie Allelic Divergence of CUP2 Between Saccharomyces Species

Under the model of micromutationism, phenotypic divergence between species is caused by accumulation of many small-effect changes. While mapping the causal changes to single nucleotide resolution could be difficult for diverged species, genetic dissection via chimeric constructs allows us to evaluate whether a large-effect gene is composed of many small-effect nucleotide changes. In a previously described non-complementation screen, we found an allele difference of CUP2, a copper-binding transcription factor, underlies divergence in copper resistance between Saccharomyces cerevisiae and S. uvarum. Here, we tested whether the allele effect of CUP2 was caused by multiple nucleotide changes. By analyzing chimeric constructs containing four separate regions in the CUP2 gene, including its distal promoter, proximal promoter, DNA binding domain and transcriptional activation domain, we found that all four regions of the S. cerevisiae allele conferred copper resistance, with the proximal promoter showing the largest effect, and that both additive and epistatic effects are likely involved. These findings support a model of multiple changes underlying evolution and suggest an important role of both protein coding and cis-regulatory changes in evolution.

In silico analysis of nonsynonymous single-nucleotide polymorphisms (nsSNPs) of the SMPX gene

Mutations in the SMPX gene can disrupt the regular activity of the SMPX protein, which is involved in the hearing process. Recent reports showing a link between nonsynonymous single-nucleotide polymorphisms (nsSNPs) in SMPX and hearing loss, thus classifying deleterious SNPs in SMPX will be an uphill task before designing a more extensive population study. In this study, damaging nsSNPs of SMPX from the dbSNP database were identified by using 13 bioinformatics tools. Initially, the impact of nsSNPs in the SMPX gene were evaluated through different in silico predictors; and the deleterious convergent changes were analyzed by energy-minimization-guided residual network analysis. In addition, the pathogenic effects of mutations in SMPX-mediated protein-protein interactions were also characterized by structural modeling and binding energy calculations. A total of four mutations (N19D, A29T, K54N, and S71L) were found to be highly deleterious by all the tools, which are located at highly conserved regions. Furthermore, all four mutants showed structural alterations, and the communities of amino acids for mutant proteins were readily changed, compared to the wild-type. Among them, A29T (rs772775896) was revealed as the most damaging nsSNP, which caused significant structural deviation of the SMPX protein, as a result reducing the binding affinity to other functional partners. These findings reflect the computational insights into the deleterious role of nsSNPs in SMPX, which might be helpful for subjecting wet-lab confirmatory analysis.

Adaptation by Loss of Heterozygosity in Saccharomyces cerevisiae Clones Under Divergent Selection

Loss of heterozygosity (LOH) is observed during vegetative growth and reproduction of diploid genotypes through mitotic crossovers, aneuploidy caused by nondisjunction, and gene conversion. We aimed to test the role that LOH plays during adaptation of two highly heterozygous Saccharomyces cerevisiae genotypes to multiple environments over a short time span in the laboratory. We hypothesized that adaptation would be observed through parallel LOH events across replicate populations. Using genome resequencing of 70 clones, we found that LOH was widespread with 5.2 LOH events per clone after ∼500 generations. The most common mode of LOH was gene conversion (51%) followed by crossing over consistent with either break-induced replication or double Holliday junction resolution. There was no evidence that LOH involved nondisjunction of whole chromosomes. We observed parallel LOH in both an environment-specific and environment-independent manner. LOH largely involved recombining existing variation between the parental genotypes, but also was observed after de novo, presumably beneficial, mutations occurred in the presence of canavanine, a toxic analog of arginine. One highly parallel LOH event involved the ENA salt efflux pump locus on chromosome IV, which showed repeated LOH to the allele from the European parent, an allele originally derived by introgression from S. paradoxus Using CRISPR-engineered LOH we showed that the fitness advantage provided by this single LOH event was 27%. Overall, we found extensive evidence that LOH could be adaptive and is likely to be a greater source of initial variation than de novo mutation for rapid evolution of diploid genotypes.

Computational analysis of high-risk SNPs in human CHK2 gene responsible for hereditary breast cancer: A functional and structural impact

Nowadays CHK2 mutation is studied frequently in hereditary breast and ovarian cancer patients in addition to BRCA1/BRCA2. CHK2 is a tumor suppressor gene that encodes a serine/threonine kinase, also involved in pathways such as DNA repair, cell cycle regulation and apoptosis in response to DNA damage. CHK2 is a well-studied moderate penetrance gene that correlates with third high risk susceptibility gene with an increased risk for breast cancer. Hence before planning large population study, it is better to scrutinize putative functional SNPs of CHK2 using different computational tools. In this study, we have used various computational approaches to identify nsSNPs which are deleterious to the structure and/or function of CHK2 protein that might be causing this disease. Computational analysis was performed by different in silico tools including SIFT, Align GVGD, SNAP-2, PROVEAN, Poly-Phen-2, PANTHER, PhD-SNP, MUpro, iPTREE-STAB, Consurf, InterPro, NCBI Conserved Domain Search tool, ModPred, SPARKS-X, RAMPAGE, Verify-3D, FT Site, COACH and PyMol. Out of 78 nsSNP of human CHK2 gene, seven nsSNPs were predicted functionally most significant SNPs. Among these seven nsSNP, p.Arg160Gly, p.Gly210Arg and p.Ser415Phe are highly conserved residues with conservation score of 9 and three nsSNP were predicted to be involved in post translational modification. The p.Arg160Gly and p.Gly210Arg may interfere in phosphopeptide binding site on FHA conserved domain. The p.Ser415Phe may interfere in formation of activation loop of protein-kinase domain and might interfere in interactions of CHK2 with ligand. The study concludes that mutation of serine to phenylalanine at position 415 is a major mutation in native CHK2 protein which might contribute to its malfunction, ultimately causing disease. This is the first comprehensive study, where CHK2 gene variants are analyzed using in silico tools hence it will be of great help while considering large scale studies and also in developing precision medicines related to these polymorphisms in the era of personalized medicine.

Identification of most damaging nsSNPs in human CCR6 gene: In silico analyses

Single nucleotide polymorphisms in CCR6 (C-C chemokine receptor type 6) gene have been found to be the possible cause of many diseases like rheumatoid arthritis, psoriasis, lupus nephritis and systemic sclerosis and other autoimmune diseases. Therefore, identification of structurally and functionally important polymorphisms in CCR6 is important in order to study its potential malfunctioning and discovering therapeutic targets. Several bioinformatics tools were used to identify most damaging nsSNPs that might be vital for CCR6 structure and function. The in silico tools included PROVEAN, SIFT, SNP&GO and PolyPhen2 followed by I-Mutant MutPred and ConSurf. Phyre2 and I-TASSER were used for protein 3-D Modelling while gene-gene interaction was predicted by STRING and GeneMANIA. Our study suggested that three nsSNPs rs1376162684, rs751102128 and rs1185426631 are the most damaging in CCR6 gene while 7 missense SNPs rs1438637216, rs139697820, rs768420505, rs1282264186, rs1394647982, rs769360638 and rs1263402382 are found to revert into stop codons. Prediction of post-transcriptional modifications highlighted the significance of rs1376162684 because it effected potential phosphorylation site. Gene-gene interactions showed relation of CCR6 with other genes depicting its importance in several pathways and co-expressions. In future, studying diseases related to CCR6 should include investigation of these 10 nsSNPs. Being the first of its type, this study also proposes future perspectives that will help in precision medicines. For such purposes, CCR6 proteins from patients of autoimmune diseases should be explored. Animal models can also be of significance find out the effects of CCR6 in diseases.

Quantitative Trait Nucleotides Impacting the Technological Performances of Industrial Saccharomyces cerevisiae Strains

The budding yeast Saccharomyces cerevisiae is certainly the prime industrial microorganism and is related to many biotechnological applications including food fermentations, biofuel production, green chemistry, and drug production. A noteworthy characteristic of this species is the existence of subgroups well adapted to specific processes with some individuals showing optimal technological traits. In the last 20 years, many studies have established a link between quantitative traits and single-nucleotide polymorphisms found in hundreds of genes. These natural variations constitute a pool of QTNs (quantitative trait nucleotides) that modulate yeast traits of economic interest for industry. By selecting a subset of genes functionally validated, a total of 284 QTNs were inventoried. Their distribution across pan and core genome and their frequency within the 1,011 Saccharomyces cerevisiae genomes were analyzed. We found that 150 of the 284 QTNs have a frequency lower than 5%, meaning that these variants would be undetectable by genome-wide association studies (GWAS). This analysis also suggests that most of the functional variants are private to a subpopulation, possibly due to their adaptive role to specific industrial environment. In this review, we provide a literature survey of their phenotypic impact and discuss the opportunities and the limits of their use for industrial strain selection.

A Quasi-Domesticate Relic Hybrid Population of Saccharomyces cerevisiae × S. paradoxus Adapted to Olive Brine

The adaptation of the yeast Saccharomyces cerevisiae to man-made environments for the fermentation of foodstuffs and beverages illustrates the scientific, social, and economic relevance of microbe domestication. Here we address a yet unexplored aspect of S. cerevisiae domestication, that of the emergence of lineages harboring some domestication signatures but that do not fit completely in the archetype of a domesticated yeast, by studying S. cerevisiae strains associated with processed olives, namely table olives, olive brine, olive oil, and alpechin. We confirmed earlier observations that reported that the Olives population results from a hybridization between S. cerevisiae and S. paradoxus. We concluded that the olive hybrids form a monophyletic lineage and that the S. cerevisiae progenitor belonged to the wine population of this species. We propose that homoploid hybridization gave rise to a diploid hybrid genome, which subsequently underwent the loss of most of the S. paradoxus sub-genome. Such a massive loss of heterozygosity was probably driven by adaptation to the new niche. We observed that olive strains are more fit to grow and survive in olive brine than control S. cerevisiae wine strains and that they appear to be adapted to cope with the presence of NaCl in olive brine through expansion of copy number of ENA genes. We also investigated whether the S. paradoxus HXT alleles retained by the Olives population were likely to contribute to the observed superior ability of these strains to consume sugars in brine. Our experiments indicate that sugar consumption profiles in the presence of NaCl are different between members of the Olives and Wine populations and only when cells are cultivated in nutritional conditions that support adaptation of their proteome to the high salt environment, which suggests that the observed differences are due to a better overall fitness of olives strains in the presence of high NaCl concentrations. Although relic olive hybrids exhibit several characteristics of a domesticated lineage, tangible benefits to humans cannot be associated with their phenotypes. These strains can be seen as a case of adaptation without positive or negative consequences to humans, that we define as a quasi-domestication.

A new chromosomal rearrangement improves the adaptation of wine yeasts to sulfite

Sulfite-generating compounds are widely used during winemaking as preservatives because of its antimicrobial and antioxidant properties. Thus, wine yeast strains have developed different genetic strategies to increase its sulfite resistance. The most efficient sulfite detoxification mechanism in Saccharomyces cerevisiae uses a plasma membrane protein called Ssu1 to efflux sulfite. In wine yeast strains, two chromosomal translocations (VIIItXVI and XVtXVI) involving the SSU1 promoter region have been shown to upregulate SSU1 expression and, as a result, increase sulfite tolerance. In this study, we have identified a novel chromosomal rearrangement that triggers wine yeast sulfite adaptation. An inversion in chromosome XVI (inv-XVI) probably due to sequence microhomology, which involves SSU1 and GCR1 regulatory regions, increases the expression of SSU1 and the sulfite resistance of a commercial wine yeast strain. A detailed dissection of this chimeric SSU1 promoter indicates that both the removed SSU1 promoter sequence and the relocated GCR1 sequence contribute to SSU1 upregulation and sulfite tolerance. However, no relevant function has been attributed to the SSU1-promoter-binding transcription factor Fzf1. These results unveil a new genomic event that confers an evolutive advantage to wine yeast strains.

A polyploid admixed origin of beer yeasts derived from European and Asian wine populations

Strains of Saccharomyces cerevisiae used to make beer, bread, and wine are genetically and phenotypically distinct from wild populations associated with trees. The origins of these domesticated populations are not always clear; human-associated migration and admixture with wild populations have had a strong impact on S. cerevisiae population structure. We examined the population genetic history of beer strains and found that ale strains and the S. cerevisiae portion of allotetraploid lager strains were derived from admixture between populations closely related to European grape wine strains and Asian rice wine strains. Similar to both lager and baking strains, ale strains are polyploid, providing them with a passive means of remaining isolated from other populations and providing us with a living relic of their ancestral hybridization. To reconstruct their polyploid origin, we phased the genomes of two ale strains and found ale haplotypes to both be recombinants between European and Asian alleles and to also contain novel alleles derived from extinct or as yet uncharacterized populations. We conclude that modern beer strains are the product of a historical melting pot of fermentation technology.

Identification and structural characterization of deleterious non-synonymous single nucleotide polymorphisms in the human SKP2 gene

In SCF (Skp, Cullin, F-box) ubiquitin-protein ligase complexes, S-phase kinase 2 (SKP2) is one of the major players of F-box family, that is responsible for the degradation of several important cell regulators and tumor suppressor proteins. Despite of having significant evidence for the role of SKP2 on tumorgenesis, there is a lack of available data regarding the effect of non-synonymous polymorphisms. In this communication, the structural and functional consequences of non-synonymous single nucleotide polymorphisms (nsSNPs) of SKP2 have been reported by employing various computational approaches and molecular dynamics simulation. Initially, several computational tools like SIFT, PolyPhen-2, PredictSNP, I-Mutant 2.0 and ConSurf have been implicated in this study to explore the damaging SNPs. In total of 172 nsSNPs, 5 nsSNPs were identified as deleterious and 3 of them were predicted to be decreased the stability of protein. Guided from ConSurf analysis, P101L (rs761253702) and Y346C (rs755010517) were categorized as the highly conserved and functional disrupting mutations. Therefore, these mutations were subjected to three dimensional model building and molecular dynamics simulation study for the detailed structural consequences upon the mutations. The study revealed that P101L and Y346C mutations increased the flexibility and changed the structural dynamics. As both these mutations are located in the most functional regions of SKP2 protein, these computational insights might be helpful to consider these nsSNPs for wet-lab confirmatory analysis as well as in rationalizing future population based studies and structure based drug design against SKP2.

Deleterious Mutation Burden and Its Association with Complex Traits in Sorghum (Sorghum bicolor)

Sorghum (Sorghum bicolor (L.) Moench) is a major staple food cereal for millions of people worldwide. Valluru et al. identify putative deleterious mutations among ∼5.5M segregating variants of 229 diverse sorghum...

Sorghum (Sorghum bicolor L.) is a major food cereal for millions of people worldwide. The sorghum genome, like other species, accumulates deleterious mutations, likely impacting its fitness. The lack of recombination, drift, and the coupling with favorable loci impede the removal of deleterious mutations from the genome by selection. To study how deleterious variants impact phenotypes, we identified putative deleterious mutations among ∼5.5 M segregating variants of 229 diverse biomass sorghum lines. We provide the whole-genome estimate of the deleterious burden in sorghum, showing that ∼33% of nonsynonymous substitutions are putatively deleterious. The pattern of mutation burden varies appreciably among racial groups. Across racial groups, the mutation burden correlated negatively with biomass, plant height, specific leaf area (SLA), and tissue starch content (TSC), suggesting that deleterious burden decreases trait fitness. Putatively deleterious variants explain roughly one-half of the genetic variance. However, there is only moderate improvement in total heritable variance explained for biomass (7.6%) and plant height (average of 3.1% across all stages). There is no advantage in total heritable variance for SLA and TSC. The contribution of putatively deleterious variants to phenotypic diversity therefore appears to be dependent on the genetic architecture of traits. Overall, these results suggest that incorporating putatively deleterious variants into genomic models slightly improves prediction accuracy because of extensive linkage. Knowledge of deleterious variants could be leveraged for sorghum breeding through either genome editing and/or conventional breeding that focuses on the selection of progeny with fewer deleterious alleles.