Search

Organelle genomes are known to have large sizes and substantial non-coding content, despite conserved coding regions and low substitution rates. Notably, volvocine green algae exhibit significant variation in plastid genome size, with some species harboring ptDNA ten times larger than the average. To explain this variability, my thesis explores two hypotheses. The first proposes that genetic divergence accumulates due to weak negative selection and genetic drift, resulting in similar evolution rates for coding and non-coding regions. The second suggests high evolution rates in non-coding sequences are due to error-prone repair mechanisms. Analyzing new plastid genomes from volvocine green algae, I found a potential for high silent-site substitution rates in intergenic regions. My analysis shows that these hypotheses can be applied to plastid genomes of close relatives to advance our understanding of the mechanisms of sequence evolution specific to non-coding DNA accumulation within the volvocine green algae.

In recent years, the association between synonymous mutations (SMs) and human diseases has been uncovered in many studies. It is a challenge for identifying deleterious SMs in the field of medical genomics. Although there are several computational methods proposed in the past years, the precise prediction of deleterious SMs is still challenging. In this work, we proposed a predictor named as EnDSM, which is an accurate method based on the ensemble framework. We explored multimodal features across four groups including functional score, conservation, splicing, and sequence features, and we then trained eight conceptually different machine learning classifiers for each of them, resulting in 32 base classification models. We further selected four base models referring to their prediction performance and the predictive probabilities of these base classification models were subsequently used as the input feature vectors of logistic regression classifier to construct the ensemble learning model. The results suggested that EnDSM achieved better performance comparing with other state-of-the-art predictors on the training and independent test datasets. We anticipate that our ensemble predictor EnDSM will become a valuable tool for deleterious SM prediction. The EnDSM server interface along with the benchmarking data sets are freely available at http://bioinfo.ahu.edu.cn/EnDSM.

Herein we detected single nucleotide polymorphisms in MEF2B and UCP3 by DNA sequencing and the KASPar technology and analyzed their association with sheep growth traits. Two synonymous mutations, g.1826 C > T and g.10266 G > C, were detected, respectively, and they were found to be significantly associated with sheep growth traits (p T, the average body weight and chest and cannon circumference of sheep with the CC genotype were significantly higher than those of sheep with the CT and TT genotypes (p C, the average body weight and chest and cannon circumference of sheep with the GG genotype were significantly higher than those of sheep with the GC and CC genotypes (p < 0.05). Moreover, the average body weight of sheep with the CC/GG genotype was higher compared with those of other genotype combinations. We also assessed MEF2B and UCP3 expression in different sheep tissues, confirming their expression in all examined tissues. To summarize, we believe that the polymorphisms identified in MEF2B and UCP3 can serve as molecular markers for sheep growth traits.

Genetic alterations under positive selection in healthy tissues have implications for cancer risk. However, total levels of positive selection across the genome remain unknown. Passenger mutations are influenced by all driver mutations, regardless of type or location in the genome. Therefore, the total number of passengers can be used to estimate the total number of drivers—including unidentified drivers outside of cancer genes that are traditionally missed. Here we analyze the variant allele frequency spectrum of synonymous mutations from healthy blood and esophagus to quantify levels of missing positive selection. In blood, we find that only 30% of passengers can be explained by single-nucleotide variants in driver genes, suggesting high levels of positive selection for mutations elsewhere in the genome. In contrast, more than half of all passengers in the esophagus can be explained by just the two driver genes NOTCH1 and TP53, suggesting little positive selection elsewhere. Synonymous passenger mutations are used to measure levels of positive selection in healthy blood and esophagus. This approach can quantify missing selection due to unidentified drivers.

Cancer progression has been attributed to somatic changes in single-nucleotide variants, copy-number aberrations, loss of heterozygosity, chromosomal instability, epistatic interactions, and the tumor microenvironment. It is not entirely clear which of these changes are essential and which are ancillary to cancer. The dynamic nature of cancer evolution in a patient can be illuminated using several concepts and tools from classical evolutionary biology. Neutral mutation rates in cancer cells are calculable from genomic data such as synonymous mutations, and selective pressures are calculable from rates of fixation occurring beyond the expectation by neutral mutation and drift. However, these cancer effect sizes of mutations are complicated by epistatic interactions that can determine the likely sequence of gene mutations. In turn, longitudinal phylogenetic analyses of somatic cancer progression offer an opportunity to identify key moments in cancer evolution, relating the timing of driver mutations to corresponding landmarks in the clinical timeline. These analyses reveal temporal aspects of genetic and phenotypic change during tumorigenesis and across clinical timescales. Using a related framework, clonal deconvolution, physical locations of clones, and their phylogenetic relations can be used to infer tumor migration histories. Additionally, genetic interactions with the tumor microenvironment can be analyzed with longstanding approaches applied to organismal genotype-by-environment interactions. Fitness landscapes for cancer evolution relating to genotype, phenotype, and environment could enable more accurate, personalized therapeutic strategies. An understanding of the trajectories underlying the evolution of neoplasms, primary, and metastatic tumors promises fundamental advances toward accurate and personalized predictions of therapeutic response.

Autosomal dominant polycystic kidney disease (ADPKD) is mainly caused by the pathogenic mutation of PKD1 or PKD2 gene and usually affects bilateral kidneys. Synonymous mutations are generally assumed to be neutral as they do not alter amino acids. Herein, we described an extremely rare ADPKD child caused by a heterozygous synonymous mutation of PKD2 gene accompanied by massive proteinuria and congenital solitary kidney. Clinical characteristics of the patients were summarized. Whole-exome sequencing was performed to screen the disease-causing gene mutation, and reverse transcription polymerase chain reaction (RT-PCR) and Sanger sequencing were applied to analyze the impact of the identified mutation on gene transcription and splicing. Polycystic changes were found in the solitary kidney of a girl initially presented with nephrotic-range proteinuria. Thereafter her mother and 2 other family members were diagnosed to be ADPKD. Whole-exome sequencing of the proband identified a heterozygous synonymous mutation (c.1716G>A, p.Lys572=) located in the splicing site of exon 7 in PKD2 gene, which was co-segregated with the PKD phenotype in the family. RT-PCR and direct sequencing of amplified products revealed that this heterozygous synonymous mutation led to exon7 skipping in PKD2 gene. We reported an extremely rare child case of ADPKD2 in combination with solitary kidney and nephrotic-range proteinuria, and firstly confirmed the pathogenicity of a heterozygous synonymous mutation (c.1716G>A) in PKD2 gene. The results indicate that synonymous mutations should not be excluded from disease-causing if they are located in splicing site of an exon.

Objectives Hypomagnesemia 1, intestinal (HOMG1) is characterized by neurological symptoms that occur due to hypocalcemia and hypomagnesemia and caused by mutations in the TRPM6. Most of the identified variants in TRPM6 lead to premature termination: nonsense, frameshift, deletion, and splice site mutations. Case presentation Herein, we report a 1.5 month-old case who presented with convulsion due to hypocalcemia and hypomagnesemia in the early infancy. Sequencing of TRPM6 revealed a novel homozygous synonymous variant [c.2538G > A (p.Thr846Thr)] in the last codon of exon 19, which is most likely to affect the splicing. We report a novel homozygous synonymous variant in the TRPM6 leading to HOMG1, expanding the mutational spectrum. Conclusions Synonymous mutations that were previously considered as harmless should be evaluated at the nucleotide level, keeping in mind that they may affect splicing and cause to the disease.

The number of ribs is an important economic trait in the sheep industry when the sheep are raised for mutton. However, in sheep, the genetic mechanisms regulating rib number are poorly understood. In the present study, we aimed to identify important candidate genes that affect the increase in rib number in sheep. Whole-genome resequencing of 36 Hu sheep with an increased number of ribs (R14) and 36 sheep with normal (R13) rib numbers was carried out. Analysis using three methods (fixation index (FST), Fisher's exact test, and Chi-squared test) showed that 219 single nucleotide polymorphism sites overlapped among the results of the three methods, which represented 206 genes. Gene ontology and Kyoto Encyclopedia of Genes and Genomes analyses indicated that the genes were mainly associated with regulation of developmental process, inorganic anion transport, cellular biosynthetic process, tight junction, the oxytocin signaling pathway, and arrhythmogenic right ventricular cardiomyopathy. Four mutations were selected according to the significantly selected genomic regions and important pathways for genotyping and association analysis. The result demonstrated that three synonymous mutations correlated significantly with the rib number. Importantly, we revealed that the CPOX (encoding coproporphyrinogen oxidase), KCNH1 (encoding potassium voltage-gated channel subfamily H member 1), and CPQ (encoding carboxypeptidase Q) genes have a combined effect on rib number in Hu sheep. Our results identified candidate molecular markers for rib number in sheep breeding.

β-ketoacyl-CoA synthase is a key enzyme in the biosynthesis of over-long-chain fattty acids; thus, it plays a crucial role in plant resistance to stress. Herein, 33 OsKCSs were classified into five subfamilies based on phylogenetic analysis and domain structure in rice. Nonsynonymous to synonymous substitution ratio (Ka/Ks) analysis showed that duplicated genes were possibly under strong purifying selection pressure. Cis-element analysis revealed that OsKCSs may be regulated by different transcription factors. SNP distribution of OsKCSs in 295 germplasm accessions seemed associated with functional disparity under environmental stress. Additionally, 10 of the 33 OsKCSs were not expressed in different rice tissues. RNA-Seq and qRT-PCR analysis demonstrated that seven OsKCSs (OsKCS5, OsKCS6, OsKCS14, OsKCS15, OsKCS4, OsKCS30, and OsKCS21) had significantly different expression levels under 0.5% NaCl, 0.5% Na2CO3, and a 10-day incubation at 17 °C. Our results provide a solid foundation for understanding the roles of OsKCS genes and the mechanism by which they regulate stress resistance.

Background and Aim: Myostatin (MSTN), a member of the transforming growth factor-β family, is a negative regulator of muscle mass. This study aimed to detect the genetic variation of the 1160 bp fragment of exon 1 and part of intron 1 of the MSTN gene in several cattle populations raised in Indonesia. Materials and Methods: Polymerase chain reaction products of the MSTN gene amplified from 92 animals representing 10 cattle populations (Peranakan Ongole [PO], Belgian Blue x PO cross, Rambon, PO x Bali cross, Jabres, Galekan, Sragen, Donggala, Madura, and Bali) were sequenced, compared, and aligned with bovine MSTN of Bos taurus (GenBank Acc. No. AF320998.1) and Bos indicus (GenBank Acc. No. AY794986.1). Results: Four nucleotide substitutions (nt 1045 and 1066 in intron 1; nt 262 and 418 in exon 1) and two indels (nt 807 and 869 in intron 1) were synonymous mutations. Among these substitutions, only the nt 262G>C and nt 418A>G loci were polymorphic in all populations, except Bali cattle. The frequencies of the nt 262C (0.82) and nt 418A (0.65) alleles were highest. For the nt 262G>C locus, the CC genotype had the highest frequency (0.66) followed by GC (0.30) and CC (0.03). For the nt 418A>G locus, the AG genotype had the highest frequency (0.52) followed by AA (0.39) and GG (0.09). Conclusion: The results, showing genetic variations in exon 1 and intron 1 of the MSTN gene, might be helpful for future association studies.

Synopsis Extant anurans (frogs and toads) exhibit reduced dentition, ranging from a lack of mandibular teeth to complete edentulation, as observed in the true toads of the family Bufonidae. The evolutionary time line of these reductions remains vague due to a poor fossil record. Previous studies have demonstrated an association between the lack of teeth in edentulous vertebrates and the pseudogenization of the major tooth enamel gene amelogenin (AMEL) through accumulation of deleterious mutations and the disruption of its coding sequence. In this study, we have harnessed the pseudogenization of AMEL as a molecular dating tool to correlate loss of dentition with genomic mutation patterns during the rise of the family Bufonidae. Specifically, we have utilized AMEL pseudogenes in three members of the family as a tool to estimate the putative date of edentulation in true toads. Comparison of AMEL sequences from Rhinella marina, Bufo gargarizans and Bufo bufo, with nine extant, dentulous frogs, revealed mutations confirming AMEL inactivation in Bufonidae. AMEL pseudogenes in modern bufonids also exhibited remarkably high 86–93% sequence identity among each other, with only a slight increase in substitution rate and relaxation of selective pressure, in comparison with functional copies in other anurans. Moreover, using selection intensity estimates and synonymous substitution rates, analysis of functional and pseudogenized AMEL resulted in an estimated inactivation window of 46–60 million years ago in the lineage leading to modern true toads, a time line that coincides with the rise of the family Bufonidae.

Veronica nakaiana Ohwi (Plantaginaceae) is an endemic taxon on Ulleungdo Island, Korea. We report the second complete chloroplast genome sequence of V. nakaiana. Its genome size is 152,319 bp in length, comprising a large single-copy of 83,195 bp, a small single-copy of 17,702 bp, and a pair of inverted repeat regions of 25,711 bp. The complete genome contains 115 genes, including 51 protein-coding genes, four rRNA genes, and 31 tRNA genes. When comparing the two chloroplast genomes of V. nakaiana, 11 variable sites are recognized: seven SNPs and four indels. Two substitutions in the coding regions are recognized: rpoC2 (synonymous substitution) and rpl22 (nonsynonymous substitution). In nine noncoding regions, one is in the tRNA gene (trnK-UUU), one is in the intron of atpF, and seven are in the intergenic spacers (trnH-GUG~psbA, trnK-UUU, rps16~trnQ-UUG, trnC-GCA~petN, psbZ~trnG-GCC, ycf3~trnS-GGA, ycf4~cemA, and psbB~psbT). The data provide the level of genetic variation in V. nakaiana. This result will be a useful resource to formulate conservation strategies for V. nakaiana, which is a rare endemic species in Korea.

PURPOSE: To identify the contribution of mutations in the Desert Hedgehog (DHH) gene to the disorders of sexual differentiation (DSD) and male infertility. METHODS: The study included a total 430 subjects, including 47 gonadal dysgenesis cases, 6 patients with undescended testis and infertility characterized by azoospermia, 125 infertile male patients characterized by oligoasthenozoospermia, 24 patients with oligoasthenoteratozoospermia, and 200 ethnically matched normozoospermic fertile men who had fathered a child in the last two years. Sequencing of the complete coding region of the DHH gene was undertaken to find its contribution to the DSD and male infertility. RESULTS: We observed four novel mutations in the DHH gene in the cases with different reproductive anomalies. A synonymous substitution, c. 543C>T (p.His181His) was observed in 6.6% oligoasthenozoospermic infertile males and 1.5% normozoospermic fertile control samples (RR = 4.4077, 95%CI 1.19–16.29). Another synonymous substitution, c.990G>A (p.Ala330Ala) was observed in an infertile patient with unilateral undescended testis (case #12). Insertion of G at c.1156insG (p.Arg385fs) was observed in a case with bilateral undescended testis and azoospermia (case #23). In gonadal dysgenesis category, two mutations, insertion of G at c.1156insG (p.Arg385fs) and c.997A>G (p.Thr333Ala) substitution were observed in one case (case #34). These mutations were completely absent in control samples. CONCLUSION: Mutations in the DHH gene impact reproduction with mild mutations affecting fertility, and severe or multiple mutations resulting in gonadal dysgenesis.

Genetic predisposition has been suggested to play a role in the pathogenesis of endometriosis. The most common pathogenetic hypotheses is that there is possible dysregulation of the ligand-receptorial signaling involving the main sex hormones, eostrogen and progesterone. This hypothesis indicates the need for studies to identify the genes and genomic variants involved in the pathogenesis of endometriosis. The present study investigated the genomic variants in eostrogen receptor 1 and progesterone receptor genes in baboons with induced endometriosis. Female adult olive baboons (n = 10) were induced with endometriosis by intraperitoneal inoculation of autologous menstrual endometrium. DNA was extracted from thirty plasma samples and analyzed by direct sequencing using gene specific primers to determine polymorphisms in eostrogen receptor 1 and progesterone receptor genes. To determine the single nucleotide polymorphisms, the DNA sequences were analyzed using Bioedit software. Twenty-six point mutations in eostrogen receptor 1 gene resulted to a change in amino acid (non-synonymous mutations) while thirteen had no effect on the amino acid sequence (synonymous mutations). There were two point mutations in progesterone gene although one had no effect on the amino acid sequence. Transition substitutions were more predominant than transversions in the eostrogen receptor 1 gene. These findings suggest that genetic variants of eostrogen receptor 1 gene are related to susceptibility to endometriosis and therefore this warrants further investigation on how they influence the disease pathogenesis.

Identifying patient-specific flow of signal transduction perturbed by multiple single-nucleotide alterations is critical for improving patient outcomes in cancer cases. However, accurate estimation of mutational effects at the pathway level for such patients remains an open problem. While probabilistic pathway topology methods are gaining interest among the scientific community, the overwhelming majority do not account for network perturbation effects from multiple single-nucleotide alterations. Here we present an improvement of the mutational forks formalism to infer the patient-specific flow of signal transduction based on multiple single-nucleotide alterations, including non-synonymous and synonymous mutations. The lung adenocarcinoma and skin cutaneous melanoma datasets from TCGA Pan-Cancer Atlas have been employed to show the utility of the proposed method. We have comprehensively characterized six mutational forks. The number of mutated nodes ranged from one to four depending on the topological characteristics of a fork. Transitional confidences (TCs) have been computed for every possible combination of single-nucleotide alterations in the fork. The performed analysis demonstrated the capacity of the mutational forks formalism to follow a biologically explainable logic in the identification of high-likelihood signaling routes in lung adenocarcinoma and skin cutaneous melanoma patients. The findings have been largely supported by the evidence from the biomedical literature. We conclude that the formalism has a great chance to enable an assessment of patient-specific flow by leveraging information from multiple single-nucleotide alterations to adjust the transitional likelihoods that are solely based on the canonical view of a disease.

Background The high incidence of Plasmodium vivax infections associated with clinical severity and the emergence of chloroquine (CQ) resistance has posed a challenge to control efforts aimed at eliminating this disease. Despite conflicting evidence regarding the role of mutations of P. vivax multidrug resistance 1 gene (pvmdr1) in drug resistance, this gene can be a tool for molecular surveillance due to its variability and spatial patterns. Methods Blood samples were collected from studies conducted between 2006 and 2015 in the Northern and Southern Amazon Basin and the North Coast of Peru. Thick and thin blood smears were prepared for malaria diagnosis by microscopy and PCR was performed for detection of P. vivax monoinfections. The pvmdr1 gene was subsequently sequenced and the genetic data was used for haplotype and diversity analysis. Results A total of 550 positive P. vivax samples were sequenced; 445 from the Northern Amazon Basin, 48 from the Southern Amazon Basin and 57 from the North Coast. Eight non-synonymous mutations and three synonymous mutations were analysed in 4,395 bp of pvmdr1. Amino acid changes at positions 976F and 1076L were detected in the Northern Amazon Basin (12.8%) and the Southern Amazon Basin (4.2%) with fluctuations in the prevalence of both mutations in the Northern Amazon Basin during the course of the study that seemed to correspond with a malaria control programme implemented in the region. A total of 13 pvmdr1 haplotypes with non-synonymous mutations were estimated in Peru and an overall nucleotide diversity of π = 0.00054. The Northern Amazon Basin was the most diverse region (π = 0.00055) followed by the Southern Amazon and the North Coast (π = 0.00035 and π = 0.00014, respectively). Conclusion This study showed a high variability in the frequencies of the 976F and 1076L polymorphisms in the Northern Amazon Basin between 2006 and 2015. The low and heterogeneous diversity of pvmdr1 found in this study underscores the need for additional research that can elucidate the role of this gene on P. vivax drug resistance as well as in vitro and clinical data that can clarify the extend of CQ resistance in Peru.

Mitochondrial genome provides useful information about species concerning its evolution and phylogenetics. We have taken the advantage of high throughput next-generation sequencing technique to sequence the complete mitogenome of Yellow-billed babbler (Turdoides affinis), a species endemic to Peninsular India and Sri Lanka. Both, reference-based and de-novo assemblies of mitogenome were performed and observed that de-novo assembled mitogenome was most appropriate. The complete mitogenome of yellow-billed babbler (assembled de-novo) was 17,672 bp in length with 53.2% AT composition. Thirteen protein-coding genes along with two rRNAs and 22 tRNAs were detected. The arrangement pattern of these genes was found conserved among Leiothrichidae family mitogenomes. Duplicated control regions were found in the newly sequenced mitogenome. Downstream bioinformatics analysis revealed the effect of translational efficiency and purifying selection pressure over thirteen protein-coding genes in yellow-billed babbler mitogenome. Ka/Ks analysis indicated the highest synonymous substitution rate in the nad6 gene. Evolutionary analysis revealed the conserved nature of all the protein-coding genes across Leiothrichidae family mitogenomes. Our limited phylogeny results placed T. affinis in a separate group, a sister group of Garrulax. Overall, our results provide a useful information for future studies on the evolutionary and adaptive mechanisms of birds belong to the Leiothrichidae family.

Intracellular pathogens are responsible for an enormous amount of worldwide morbidity and mortality, and each has evolved specialized strategies to establish and maintain their replicative niche. Listeria monocytogenes is a facultative intracellular pathogen that secretes a pore-forming cytolysin called listeriolysin O (LLO), which disrupts the phagosomal membrane and, thereby, allows the bacteria access to their replicative niche in the cytosol. Nonsynonymous and synonymous mutations in a PEST-like domain near the LLO N terminus cause enhanced LLO translation during intracellular growth, leading to host cell death and loss of virulence. Here, we explore the mechanism of translational control and show that there is extensive codon restriction within the PEST-encoding region of the LLO messenger RNA (mRNA) (hly). This region has considerable complementarity with the 5' UTR and is predicted to form an extensive secondary structure that overlaps the ribosome binding site. Analysis of both 5' UTR and synonymous mutations in the PEST-like domain that are predicted to disrupt the secondary structure resulted in up to a 10,000-fold drop in virulence during mouse infection, while compensatory double mutants restored virulence to WT levels. We showed by dynamic protein radiolabeling that LLO synthesis was growth phase-dependent. These data provide a mechanism to explain how the bacteria regulate translation of LLO to promote translation during starvation in a phagosome while repressing it during growth in the cytosol. These studies also provide a molecular explanation for codon bias at the 5' end of this essential determinant of pathogenesis.

Background Plasmodium malariae is a widely spread but neglected human malaria parasite, which causes chronic infections. Studies on genetic polymorphisms of anti-malarial drug target genes in P. malariae are limited. Previous reports have shown polymorphisms in the P. malariae dihydrofolate reductase gene associated with pyrimethamine resistance and linked to pyrimethamine drug pressure. This study investigated polymorphisms of the P. malariae homologous genes, chloroquine resistant transporter and multidrug resistant 1, associated with chloroquine and mefloquine resistance in Plasmodium falciparum. Methods The orthologous P. malariae crt and mdr1 genes were studied in 95 patients with P. malariae infection between 2002 and 2016 from Thailand (N = 51) and Myanmar (N = 44). Gene sequences were analysed using BioEdit, MEGA7, and DnaSP programs. Mutations and gene amplifications were compared with P. falciparum and Plasmodium vivax orthologous genes. Protein topology models derived from the observed pmcrt and pmmdr1 haplotypes were constructed and analysed using Phyre2, SWISS MODEL and Discovery Studio Visualization V 17.2. Results Two non-synonymous mutations were observed in exon 2 (H53P, 40%) and exon 8 (E278D, 44%) of pmcrt. The topology model indicated that H53P and E278D were located outside of the transmembrane domain and were unlikely to affect protein function. Pmmdr1 was more diverse than pmcrt, with 10 non-synonymous and 3 synonymous mutations observed. Non-synonymous mutations were located in the parasite cytoplasmic site, transmembrane 11 and nucleotide binding domains 1 and 2. Polymorphisms conferring amino acid changes in the transmembrane and nucleotide binding domains were predicted to have some effect on PmMDR1 conformation, but were unlikely to affect protein function. All P. malariae parasites in this study contained a single copy of the mdr1 gene. Conclusions The observed polymorphisms in pmcrt and pmmdr1 genes are unlikely to affect protein function and unlikely related to chloroquine drug pressure. Similarly, the absence of pmmdr1 copy number variation suggests limited mefloquine drug pressure on the P. malariae parasite population, despite its long time use in Thailand for the treatment of falciparum malaria.

Synonymous substitution rates in plant mitochondrial genomes vary by orders of magnitude among species, whereas synonymous rates among genes within a genome are generally consistent. Exceptionally, genes within the Ajuga reptans (Lamiaceae) mitochondrial genome exhibit unprecedented intragenomic heterogeneity in synonymous sequence divergence, but the biological mechanisms underlying this rate variation remain unclear. We tracked the origin and evolutionary trajectory of mitochondrial rate variations by dense sampling in Ajugoideae and found differences in the timing and magnitude of rate acceleration for particular genes. The most divergent genes accelerated earlier, retained a high rate across Ajugoideae, and are generally devoid of RNA editing, whereas moderately diverged genes accelerated later and retained relatively higher RNA editing frequency. The acceleration of mutation rates correlates with increased guanine-cytosine (GC) content, suggesting a key role for GC-biased gene conversion and/or repair after the breakage of ancestral gene clusters.

Genetic variation of the multi-zinc finger BTB domain transcription factor ZBTB18 can cause a spectrum of human neurodevelopmental disorders, but the underlying mechanisms are not well understood. Recently, we reported that pathogenic, de novo ZBTB18 missense mutations alter its DNA-binding specificity and gene regulatory functions, leading to human neurodevelopmental disease. However, the functional impact of the general population ZBTB18 missense variants is unknown. Here, we investigated such variants documented in the Genome Aggregation Database (gnomAD) to discover that ZBTB gene family members are intolerant to loss-of-function and missense mutations, but not synonymous mutations. We studied ZBTB18 as a protein-DNA complex to find that general population missense variants are rare, and disproportionately map to non-DNA-contact residues, in contrast to the majority of disease-associated variants that map to DNA-contact residues, essential to motif binding. We studied a selection of variants (n = 12), which spans the multi-zinc finger region to find 58.3% (7/12) of variants displayed altered DNA binding, 41.6% (5/12) exhibited altered transcriptional activity in a luciferase reporter assay, 33.3% (4/12) exhibited altered DNA binding and transcriptional activity, whereas 33.3% (4/12) displayed a negligible functional impact. Our results demonstrate that general population variants, while rare, can influence ZBTB18 function, with potential consequences for neurodevelopment, homeostasis, and disease.

The analysis of the nucleotide sequences of exons and adjacent non-coding regions of the GC gene in 108 representatives of various ethnic groups of aboriginal population of Siberia was carried out. Polymorphism was found in four nucleotide positions: non-synonymous substitutions at the rs4588 and rs7041 loci, a synonymous substitution at the rs4752 locus, and a replacement in the non-coding region at the rs3733359 locus. Seven haplotypes of the GC gene were identified. Of these, 4 haplotypes encode the Gc1F isoform, 2 haplotypes encode the Gc1S isoform, and 1 haplotype encodes the Gc2 isoform. Between-regional differences were found in the distribution of variants of the GC gene: in the northeast and in the central part of Siberia, the highest prevalence of the Gc1F and Gc1F/Gc1F variants is observed, and in the south and west of Siberia, the Gc2, Gc1S/Gc2 and Gc2/Gc2 variants are most common. In the case of the GC gene, gene-environment interactions are apparently aimed at creating a balance between the activity of vitamin D-binding protein and the level of 25-hydroxyvitamin D in the blood serum.

There are known limitations in methods of detecting positive selection. Common methods do not enable differentiation between positive selection and compensatory covariation, a major limitation. Further, the traditional method of calculating the ratio of nonsynonymous to synonymous substitutions (dN/dS) does not take into account the 3D structure of biomacromolecules nor differences between amino acids. It also does not account for saturation of synonymous mutations (dS) over long evolutionary time that renders codon-based methods ineffective for older divergences. This work aims to address these shortcomings for detecting positive selection through the development of a statistical model that examines clusters of substitutions in clusters of variable radii. Additionally, it uses a parametric bootstrapping approach to differentiate positive selection from compensatory processes. A previously reported case of positive selection in the leptin protein of primates was reexamined using this methodology.

Deficiency of the transcription factor GATA2 is a highly penetrant genetic disorder predisposing to myelodysplastic syndromes (MDS) and immunodeficiency. It has been recognized as the most common cause underlying primary MDS in children. Triggered by the discovery of a recurrent synonymous GATA2 variant, we systematically investigated 911 patients with phenotype of pediatric MDS or cellular deficiencies for the presence of synonymous alterations in GATA2. In total, we identified nine individuals with five heterozygous synonymous mutations: c.351C>G, p.T117T (N = 4); c.649C>T, p.L217L; c.981G>A, p.G327G; c.1023C>T, p.A341A; and c.1416G>A, p.P472P (N = 2). They accounted for 8.2% (9/110) of cases with GATA2 deficiency in our cohort and resulted in selective loss of mutant RNA. While for the hotspot mutation (c.351C>G) a splicing error leading to RNA and protein reduction was identified, severe, likely late stage RNA loss without splicing disruption was found for other mutations. Finally, the synonymous mutations did not alter protein function or stability. In summary, synonymous GATA2 substitutions are a new common cause of GATA2 deficiency. These findings have broad implications for genetic counseling and pathogenic variant discovery in Mendelian disorders.

Recombination increases the local GC-content in genomic regions through GC-biased gene conversion (gBGC). The recent discovery of a large genomic region with extreme GC-content in the fat sand rat Psammomys obesus provides a model to study the effects of gBGC on chromosome evolution. Here, we compare the GC-content and GC-to-AT substitution patterns across protein-coding genes of four gerbil species and two murine rodents (mouse and rat). We find that the known high-GC region is present in all the gerbils, and is characterized by high substitution rates for all mutational categories (AT-to-GC, GC-to-AT, and GC-conservative) both at synonymous and nonsynonymous sites. A higher AT-to-GC than GC-to-AT rate is consistent with the high GC-content. Additionally, we find more than 300 genes outside the known region with outlying values of AT-to-GC synonymous substitution rates in gerbils. Of these, over 30% are organized into at least 17 large clusters observable at the megabase-scale. The unusual GC-skewed substitution pattern suggests the evolution of genomic regions with very high recombination rates in the gerbil lineage, which can lead to a runaway increase in GC-content. Our results imply that rapid evolution of GC-content is possible in mammals, with gerbil species providing a powerful model to study the mechanisms of gBGC.

Genetic variations have often been detected in plants regenerated from tissue culture, which is an essential step of molecular breeding and tissue culture techniques. The incidence of unintended genetic variations should be kept low to ensure the stability of biotechnological crops. In this study, regenerated Chinese cabbage plants (Brassica rapa ssp. pekinensis) were produced and sequenced to identify the genetic variation generated during tissue culture. The sequencing data of each regenerated plants were compared with reference genome sequences. A total of 6644 genetic variations, including single-nucleotide polymorphisms as well as insertions or deletions in the regenerated plants, were detected using bioinformatic analysis, and 109 mutation candidates that are expected to consistently occur in the regenerated R0 group were identified. To investigate the inheritance of the variations, the R1 and R2 lines of the two regenerated plants were also sequenced and analyzed. The results indicated that 39 of the regenerant-specific mutation candidates were expected to be maintained despite the progress of the generation. Among those, five of the SNP candidates occurred in exons, and sequence analysis was performed on these mutations. For detailed analysis, we validated the nucleotide sequences and deduced the amino acid sequences of the genes with five base substitution candidates in the exonic regions. Among the five SNP candidates, three were found to be synonymous mutations that resulted in no changes in the amino acid; however, two mutations were nonsynonymous mutations that altered the amino acid compositions. Our findings suggested a genetic diversity of the regenerated plants and revealed the mutation candidates that consistently occurred in the regenerated plants and were maintained in the progeny lines.

Purpose Several recent studies have documented CTNNB1 and BRAF mutations which are mutually exclusive for adamantinomatous craniopharyngioma (ACP) and papillary craniopharyngioma (PCP) tumors. This discovery is helpful in the development of novel targeted therapies in successful clinical trials with BRAF mutations in PCP cases. However, no such targeted therapy is available yet for ACP. Here, we report novel mutations, which are not previously reported, in a case of an adult ACP using NGS analysis. Results Patient DNA was sequenced using Ion PI v3 chip on Ion Proton. A total of 16 variants were identified in this tumor by NGS analysis, out of which four were missense mutations, seven were synonymous mutations, and five were intronic variants. In CTNNB1 gene a known missense mutation in c.101G>T; in TP53 a known missense mutation in c.215C>G; and two known missense variants in PIK3CA, viz., in c.1173A>G; in exon 7, and in c.3128T>C; in exon 21, were found, respectively. Seven synonymous mutations were detected in this tumor, viz., in IDH1 (rs11554137), in FGFR3 (rs7688609), in PDGFRA (rs1873778), in APC (COSM3760869), in EGFR (rs1050171), in MET (rs35775721), and in RET (rs1800861), respectively. Three known, intronic variants were found in genes, such as PIK3CA, KDR, and JAK3, respectively. Also, a 3'-UTR and a splice site acceptor site variant in CSF1R and FLT3 genes were found in this tumor. We have shown allele coverage, allele ratio, and p-value, for all these mutations. The p-values and Phred quality score were significantly high for these variants. Conclusion As reported in previous studies, in ACP tumors we found a CTNNB1 mutation by NGS analysis. The PIK3CA variants we detected were not known previously in ACP tumors. Finding the PIK3CA mutations in the ACP tumors may help develop targeted therapy for a subset of craniopharyngiomas with PIK3CA activating mutations. Clinical trials are in progress with specific PIK3CA inhibitors in advanced stages of many cancers.

Background von Hippel-Lindau (VHL) disease is a familial neoplasia syndrome that results from the germline mutation of VHL. Pathogenic VHL mutations include deletion, frameshift, nonsense and missense mutations. Synonymous mutations are expected to be phenotypically silent and their role in VHL disease remains poorly understood. Case presentation We report a Caucasian male with a family history of pheochromocytoma and the synonymous VHL mutation c.414A > G (p.Pro138Pro). At 47-years, MRI revealed pheochromocytoma in the left adrenal gland and hemangioblastomas in the spine and brain. Pheochromocytoma was treated by adrenalectomy. Radiotherapy, followed by craniotomy and resection were needed to reduce hemangioblastomas to residual lesions. Two of three of the proband’s children inherited the mutation and both presented with retinal hemangioblastomas without pheochromocytoma at age 7: one twin needed four laser treatments. Primary skin fibroblasts carrying the heterozygous mutation or wild type VHL were established from the family. Mutant fibroblasts downregulated full-length VHL mRNA and protein, and upregulated the short VHL mRNA isoform (a result of exon 2 skipping in splicing) at the mRNA level but not at the protein level. Conclusions Our study shows that the synonymous VHL mutation c.414A > G can within 7 years induce pediatric retinal hemangioblastoma in absence of pheochromocytoma. This highlights the need to include splicing-altering synonymous mutations into the screening for VHL disease. This is also the first report on detecting and validating a synonymous VHL mutation using patient-derived fibroblasts. The mutation c.414A > G translates to p.Pro138Pro, yet it is not functionally silent, because it causes aberrant splicing by skipping exon 2. The reduced but not completely abolished pVHL protein in a loss-of-heterozygosity genetic backdrop may underlie the etiology of VHL disease.

Background Biological monitoring of antiretroviral treatment (ART) in human immunodeficiency virus (HIV)-infected pediatric population remains challenging. The aim of the present study was to assess the long-term HIV-1 genetic diversity in pol gene in HIV-1-infected children in virological failure under antiretroviral regimen adapted according to the successive World Health Organization (WHO) guidelines for resource-constrained settings. Methods HIV-1 diversity in pol gene was assessed in HIV-1-infected children and adolescents born from HIV-infected mothers (median age at follow-up: 13.8 years) in virological failure (VF+) despite long-term regimen recommended by the WHO. The numbers of nonsynonymous substitutions per potential nonsynonymous site (dN) and of synonymous substitutions at potential synonymous sites (dS) in HIV-1 pol gene and the dN/dS ratios were used to estimate the selective pressure on circulating HIV-1. Results The immunological responses to ART basically corresponded to: 1) Full therapeutic failure with immunological (I-) and virological nonresponses in one-quarter (24.6%) of study children ((I-, VF+) subgroup); 2) Discordant immunovirological responses with paradoxical high CD4 T cell counts (I+) and high HIV-1 RNA load in the remaining cohort patients (75.4%) ((I+, VF+) subgroup). The mean dS was 1.8-fold higher in (I+, VF+) than (I-, VF+) subgroup (25.9 ± 18.4 vs. 14.3 ± 10.8). In the (I+, VF+) subgroup, the mean dS was 1.6-fold higher than the mean dN. Finally, the mean dN/dS ratio was 2.1-fold lower in (I+, VF+) than (I-, VF+) subgroup (0.6 ± 0.3 vs. 1.3 ± 0.7), indicating purifying selection in the immunovirological discordant (I+, VF+) subgroup and positive selection in the immunovirological failure (I-, VF+) subgroup. Conclusions Children and adolescents in immunovirological therapeutic failure harbor positive selection of HIV-1 strains favoring diversifying in pol-encoded amino acids. In contrast, children with persistent discordant immunovirological responses show accumulation of mutations and purifying selection in pol gene sequences, indicating limited genetic evolution and likely suggesting genetic adaptation of viruses to host functional constraints.