GC content elevates mutation and recombination rates in the yeast Saccharomyces cerevisiae

The recombination landscape of introgression in yeast

Meiotic recombination is an evolutionary force that acts by breaking up genomic linkage, increasing the efficacy of selection. Recombination is initiated with a double-strand break which is resolved via a crossover, which involves the reciprocal exchange of genetic material between homologous chromosomes, or a non-crossover, which results in small tracts of non-reciprocal exchange of genetic material. Crossover and non-crossover rates vary between species, populations, individuals, and across the genome. In recent years, recombination rate has been associated with the distribution of ancestry derived from past interspecific hybridization (introgression) in a variety of species. We explore this interaction of recombination and introgression by sequencing spores and detecting crossovers and non-crossovers from two crosses of the yeast Saccharomyces uvarum. One cross is between strains which each contain introgression from their sister species, S. eubayanus, while the other cross has no introgression present. We find that the recombination landscape is significantly different between S. uvarum crosses, and that some of these differences can be explained by the presence of introgression in one cross. Crossovers are reduced and non-crossovers are increased in heterozygous introgression compared to syntenic regions in the cross without introgression. This translates to reduced allele shuffling within introgressed regions, and an overall reduction of shuffling on most chromosomes with introgression compared to the syntenic regions and chromosomes without introgression. Our results suggest that hybridization can significantly influence the recombination landscape, and that the reduction in allele shuffling contributes to the initial purging of introgression in the generations following a hybridization event.

The Prediction of Recombination Hotspot Based on Automated Machine Learning

Meiotic recombination plays a pivotal role in genetic evolution. Genetic variation induced by recombination is a crucial factor in generating biodiversity and a driving force for evolution. At present, the development of recombination hotspot prediction methods has encountered challenges related to insufficient feature extraction and limited generalization capabilities. This paper focused on the research of recombination hotspot prediction methods. We explored deep learning-based recombination hotspot prediction and scrutinized the shortcomings of prevalent models in addressing the challenge of recombination hotspot prediction. To addressing these deficiencies, an automated machine learning approach was utilized to construct recombination hotspot prediction model. The model combined sequence information with physicochemical properties by employing TF-IDF-Kmer and DNA composition components to acquire more effective feature data. Experimental results validate the effectiveness of the feature extraction method and automated machine learning technology used in this study. The final model was validated on three distinct datasets and yielded accuracy rates of 97.14%, 79.71%, and 98.73%, surpassing the current leading models by 2%, 2.56%, and 4%, respectively. In addition, we incorporated tools such as SHAP and AutoGluon to analyze the interpretability of black-box models, delved into the impact of individual features on the results, and investigated the reasons behind misclassification of samples. Finally, an application of recombination hotspot prediction was established to facilitate easy access to necessary information and tools for researchers. The research outcomes of this paper underscore the enormous potential of automated machine learning methods in gene sequence prediction.

URA6 mutations provide an alternative mechanism for 5-FOA resistance in Saccharomyces cerevisiae

URA3 is frequently used in the yeast community as the mutation target for 5-fluoroorotic acid (5-FOA) resistance. We identified a novel class of ura6 mutants that can grow in the presence of 5-FOA. Unlike ura3 mutants, ura6 mutants remain prototrophic and grow in the absence of uracil. In addition to 5-FOA resistance, we found that mutations to URA6 also confer resistance to 5-fluorocytosine (5-FC) and 5-fluorouracil (5-FU). In total, we identified 50 unique missense mutations across 32 residues of URA6. We found that 28 out of the 32 affected residues are located in regions conserved between Saccharomyces cerevisiae and three clinically relevant pathogenic fungi. These findings suggest that mutations to URA6 present a second target for mutation screens utilizing 5-FOA as a selection marker as well as a potential mode of resistance to the antifungal therapeutic 5-FC.

Terminal regions of a protein are a hotspot for low complexity regions and selection

Volatile low complexity regions (LCRs) are a novel source of adaptive variation, functional diversification and evolutionary novelty. An interplay of selection and mutation governs the composition and length of low complexity regions. High %GC and mutations provide length variability because of mechanisms like replication slippage. Owing to the complex dynamics between selection and mutation, we need a better understanding of their coexistence. Our findings underscore that positively selected sites (PSS) and low complexity regions prefer the terminal regions of genes, co-occurring in most Tetrapoda clades. We observed that positively selected sites within a gene have position-specific roles. Central-positively selected site genes primarily participate in defence responses, whereas terminal-positively selected site genes exhibit non-specific functions. Low complexity region-containing genes in the Tetrapoda clade exhibit a significantly higher %GC and lower ω (dN/dS: non-synonymous substitution rate/synonymous substitution rate) compared with genes without low complexity regions. This lower ω implies that despite providing rapid functional diversity, low complexity region-containing genes are subjected to intense purifying selection. Furthermore, we observe that low complexity regions consistently display ubiquitous prevalence at lower purity levels, but exhibit a preference for specific positions within a gene as the purity of the low complexity region stretch increases, implying a composition-dependent evolutionary role. Our findings collectively contribute to the understanding of how genetic diversity and adaptation are shaped by the interplay of selection and low complexity regions in the Tetrapoda clade.

Global and local genomic features together modulate the spontaneous single nucleotide mutation rate

Spontaneous mutations are evolutionary engines as they generate variants for the evolutionary downstream processes that give rise to speciation and adaptation. Single nucleotide mutations (SNM) are the most abundant type of mutations among them. Here, we perform a meta-analysis to quantify the influence of selected global genomic parameters (genome size, genomic GC content, genomic repeat fraction, number of coding genes, gene count, and strand bias in prokaryotes) and local genomic features (local GC content, repeat content, CpG content and the number of SNM at CpG islands) on spontaneous SNM rates across the tree of life (prokaryotes, unicellular eukaryotes, multicellular eukaryotes) using wild-type sequence data in two different taxon classification systems. We find that the spontaneous SNM rates in our data are correlated with many genomic features in prokaryotes and unicellular eukaryotes irrespective of their sample sizes. On the other hand, only the number of coding genes was correlated with the spontaneous SNM rates in multicellular eukaryotes primarily contributed by vertebrates data. Considering local features, we notice that local GC content and CpG content significantly were correlated with the spontaneous SNM rates in the unicellular eukaryotes, while local repeat fraction is an important feature in prokaryotes and certain specific uni- and multi-cellular eukaryotes. Such predictive features of the spontaneous SNM rates often support non-linear models as the best fit compared to the linear model. We also observe that the strand asymmetry in prokaryotes plays an important role in determining the spontaneous SNM rates but the SNM spectrum does not.

Evolution of homologous recombination rates across bacteria

Bacteria are nonsexual organisms but are capable of exchanging DNA at diverse degrees through homologous recombination. Intriguingly, the rates of recombination vary immensely across lineages where some species have been described as purely clonal and others as "quasi-sexual." However, estimating recombination rates has proven a difficult endeavor and estimates often vary substantially across studies. It is unclear whether these variations reflect natural variations across populations or are due to differences in methodologies. Consequently, the impact of recombination on bacterial evolution has not been extensively evaluated and the evolution of recombination rate-as a trait-remains to be accurately described. Here, we developed an approach based on Approximate Bayesian Computation that integrates multiple signals of recombination to estimate recombination rates. We inferred the rate of recombination of 162 bacterial species and one archaeon and tested the robustness of our approach. Our results confirm that recombination rates vary drastically across bacteria; however, we found that recombination rate-as a trait-is conserved in several lineages but evolves rapidly in others. Although some traits are thought to be associated with recombination rate (e.g., GC-content), we found no clear association between genomic or phenotypic traits and recombination rate. Overall, our results provide an overview of recombination rate, its evolution, and its impact on bacterial evolution.

Absence of both MGME1 and POLG EXO abolishes mtDNA whereas absence of either creates unique mtDNA duplications

Both POLG and MGME1 are needed for mitochondrial DNA (mtDNA) maintenance in animal cells. POLG, the primary replicative polymerase of the mitochondria, has an exonuclease activity (3'→5') that corrects for the misincorporation of bases. MGME1 serves as an exonuclease (5'→3'), producing ligatable DNA ends. Although both have a critical role in mtDNA replication and elimination of linear fragments, these mechanisms are still not fully understood. Using digital PCR to evaluate and compare mtDNA integrity, we show that Mgme1 knock out (Mgme1 KK) tissue mtDNA is more fragmented than POLG exonuclease-deficient "Mutator" (Polg MM) or WT tissue. In addition, next generation sequencing of mutant hearts showed abundant duplications in/nearby the D-loop region and unique 100 bp duplications evenly spaced throughout the genome only in Mgme1 KK hearts. However, despite these unique mtDNA features at steady-state, we observed a similar delay in the degradation of mtDNA after an induced double strand DNA break in both Mgme1 KK and Polg MM models. Lastly, we characterized double mutant (Polg MM/Mgme1 KK) cells and show that mtDNA cannot be maintained without at least one of these enzymatic activities. We propose a model for the generation of these genomic abnormalities which suggests a role for MGME1 outside of nascent mtDNA end ligation. Our results highlight the role of MGME1 in and outside of the D-loop region during replication, support the involvement of MGME1 in dsDNA degradation, and demonstrate that POLG EXO and MGME1 can partially compensate for each other in maintaining mtDNA.

Epigenetic and Genetic Differentiation Between Coregonus Species Pairs

Phenotypic diversification is classically associated with genetic differentiation and gene expression variation. However, increasing evidence suggests that DNA methylation is involved in evolutionary processes due to its phenotypic and transcriptional effects. Methylation can increase mutagenesis and could lead to increased genetic divergence between populations experiencing different environmental conditions for many generations, though there has been minimal empirical research on epigenetically induced mutagenesis in diversification and speciation. Whitefish, freshwater members of the salmonid family, are excellent systems to study phenotypic diversification and speciation due to the repeated divergence of benthic-limnetic species pairs serving as natural replicates. Here we investigate whole genome genetic and epigenetic differentiation between sympatric benthic-limnetic species pairs in lake and European whitefish (Coregonus clupeaformis and Coregonus lavaretus) from four lakes (N = 64). We found considerable, albeit variable, genetic and epigenetic differences between species pairs. All SNP types were enriched at CpG sites supporting the mutagenic nature of DNA methylation, though C>T SNPs were most common. We also found an enrichment of overlaps between outlier SNPs with the 5% highest FST between species and differentially methylated loci. This could possibly represent differentially methylated sites that have caused divergent genetic mutations between species, or divergent selection leading to both genetic and epigenetic variation at these sites. Our results support the hypothesis that DNA methylation contributes to phenotypic divergence and mutagenesis during whitefish speciation.

Novel insights into the effects of 5-hydroxymethfurural on genomic instability and phenotypic evolution using a yeast model

5-Hydroxymethfurural (5-HMF) is naturally found in a variety of foods and beverages and represents a main inhibitor in the lignocellulosic hydrolysates used for fermentation. This study investigated the impact of 5-HMF on the genomic stability and phenotypic plasticity of the yeast Saccharomyces cerevisiae. Using next-generation sequencing technology, we examined the genomic alterations of diploid S. cerevisiae isolates that were subcultured on a medium containing 1.2 g/L 5-HMF. We found that in 5-HMF-treated cells, the rates of chromosome aneuploidy, large deletions/duplications, and loss of heterozygosity were elevated compared with that in untreated cells. 5-HMF exposure had a mild impact on the rate of point mutations but altered the mutation spectrum. Contrary to what was observed in untreated cells, more monosomy than trisomy occurred in 5-HMF-treated cells. The aneuploidy mutant with monosomic chromosome IX was more resistant to 5-HMF than the diploid parent strain because of the enhanced activity of alcohol dehydrogenase. Finally, we found that overexpression of ADH6 and ZWF1 effectively stabilized the yeast genome under 5-HMF stress. Our findings not only elucidated the global effect of 5-HMF on the genomic integrity of yeast but also provided novel insights into how chromosomal instability drives the environmental adaptability of eukaryotic cells.IMPORTANCESingle-cell microorganisms are exposed to a range of stressors in both natural and industrial settings. This study investigated the effects of 5-hydroxymethfurural (5-HMF), a major inhibitor found in baked foods and lignocellulosic hydrolysates, on the chromosomal instability of yeast. We examined the mechanisms leading to the distinct patterns of 5-HMF-induced genomic alterations and discovered that chromosomal loss, typically viewed as detrimental to cell growth under most conditions, can contribute to yeast tolerance to 5-HMF. Our results increased the understanding of how specific stressors stimulate genomic plasticity and environmental adaptation in yeast.

The Causes for Genomic Instability and How to Try and Reduce Them Through Rational Design of Synthetic DNA

Genetic engineering has revolutionized our ability to manipulate DNA and engineer organisms for various applications. However, this approach can lead to genomic instability, which can result in unwanted effects such as toxicity, mutagenesis, and reduced productivity. To overcome these challenges, smart design of synthetic DNA has emerged as a promising solution. By taking into consideration the intricate relationships between gene expression and cellular metabolism, researchers can design synthetic constructs that minimize metabolic stress on the host cell, reduce mutagenesis, and increase protein yield. In this chapter, we summarize the main challenges of genomic instability in genetic engineering and address the dangers of unknowingly incorporating genomically unstable sequences in synthetic DNA. We also demonstrate the instability of those sequences by the fact that they are selected against conserved sequences in nature. We highlight the benefits of using ESO, a tool for the rational design of DNA for avoiding genetically unstable sequences, and also summarize the main principles and working parameters of the software that allow maximizing its benefits and impact.

Narrow genetic diversity in germplasm from the Guinean and Sudano-Guinean zones in Benin indicates the need to broaden the genetic base of sweet fig banana (Musa acuminata cv Sotoumon)

Sweet fig (M. acuminata cv. Sotoumon) is an economically important dessert banana in Benin, with high nutritional, medicinal, and cultural values. Nevertheless, its productivity and yield are threatened by biotic and abiotic stresses. Relevant knowledge of the genetic diversity of this economically important crop is essential for germplasm conservation and the development of breeding programs. However, very little is known about the genetic makeup of this cultivar in Benin. To advance the understanding of genetic diversity in sweet fig banana germplasm, a Genotype-By-Sequencing (GBS) was performed on a panel of 273 accessions collected in different phytogeographical zones of Benin. GBS generated 8,457 quality SNPs, of which 1992 were used for analysis after filtering. The results revealed a low diversity in the studied germplasm (He = 0.0162). Genetic differentiation was overall very low in the collection as suggested by the negative differentiation index (Fstg = -0.003). The Analysis of Molecular Variance (AMOVA) indicated that the variation between accessions within populations accounted for 83.8% of the total variation observed (P < 0.001). The analysis of population structure and neighbor-joining tree partitioned the germplasm into three clusters out of which a predominant major one contained 98.1% of all accessions. These findings demonstrate that current sweet fig banana genotypes shared a common genetic background, which made them vulnerable to biotic and abiotic stress. Therefore, broadening the genetic base of the crop while maintaining its quality attributes and improving yield performance is of paramount importance. Moreover, the large genetic group constitutes an asset for future genomic selection studies in the crop and can guide the profiling of its conservation strategies.

Homologous Recombination and Repair Functions Required for Mutagenicity during Yeast Meiosis

Several meiotic events reshape the genome prior to its transfer (via gametes) to the next generation. The occurrence of new meiotic mutations is tightly linked to homologous recombination (HR) and firmly depends on Spo11-induced DNA breaks. To gain insight into the molecular mechanisms governing mutagenicity during meiosis, we examined the timing of mutation and recombination events in cells deficient in various DNA HR-repair genes, which represent distinct functions along the meiotic recombination process. Despite sequence similarities and overlapping activities of the two DNA translocases, Rad54 and Tid1, we observed essential differences in their roles in meiotic mutation occurrence: in the absence of Rad54, meiotic mutagenicity was elevated 8-fold compared to the wild type (WT), while in the tid1Δ mutant, there were few meiotic mutations, nine percent compared to the WT. We propose that the presence of Rad54 channels recombinational repair to a less mutagenic pathway, whereas repair assisted by Tid1 is more mutagenic. A 3.5-fold increase in mutation level was observed in dmc1∆ cells, suggesting that single-stranded DNA (ssDNA) may be a potential source for mutagenicity during meiosis. Taken together, we suggest that the introduction of de novo mutations also contributes to the diversification role of meiotic recombination. These rare meiotic mutations revise genomic sequences and may contribute to long-term evolutionary changes.

Identification of the Keratin-Associated Protein 22-2 Gene in the Capra hircus and Association of Its Variation with Cashmere Traits

The Cashmere goat is an excellent local goat breed in Gansu Province of China, and it is expected to improve cashmere production and cashmere quality through selection and breeding to enhance its commercial value. Keratin-associated proteins (KAPs) play an important role in maintaining wool structure. The gene encoding the keratin-associated protein 22-2 (KAP22-2) gene has been identified in selected species other than goats, such as humans, mice, and sheep. In this study, the sequence of the sheep KAP22-2 gene (KRTAP22-2) was aligned into the goat genome, and the sequence with the highest homology was assumed to be the goat KRTAP22-2 sequence and used to design primers to amplify the goat gene sequence. A total of 356 Longdong Cashmere goats (Gansu Province, China) were used for screening of genetic variants. Four specific bands were detected by polymerase chain reaction-single-stranded conformational polymorphism (PCR-SSCP) analysis, and they formed a total of six band types individually or in combination. Four alleles were identified by DNA sequencing of PCR amplification products. A total of four single nucleotide polymorphic sites (SNPs) were detected in the four sequenced KRTAP22-2 alleles. Two of them are in the 5'UTR region and the other two are in the coding region, and the variants in the coding region are all non-synonymous mutations. In addition, there was a 6 bp length variation in allele C. The gene was expressed in the cortical layer of primary and secondary hair follicles, the inner root sheath, as well as hair papillae and hair maternal cells in goats. The results of the correlation analysis between genotypes and cashmere traits showed that after excluding genotypes with a gene frequency of less than 5%, the mean fiber diameter (MFD) of cashmere was significantly higher in the AB genotype than in the AA and AC genotypes. That is, the KRTAP22-2 gene variants are associated with mean fiber diameter in cashmere. The above results suggest that the goat KRTAP22-2 variant can be utilized as a molecular marker candidate gene for cashmere traits.

Short Interrupted Repeat Cassette (SIRC)-Novel Type of Repetitive DNA Element Found in Arabidopsis thaliana

Short interrupted repeat cassette (SIRC)-a novel DNA element found throughout the A. thaliana nuclear genome. SIRCs are represented by short direct repeats interrupted by diverse DNA sequences. The maxima of SIRC's distribution are located within pericentromeric regions. We suggest that originally SIRC was a special case of the complex internal structure of the miniature inverted repeat transposable element (MITE), and further MITE amplification, transposition, and loss of terminal inverted repeats gave rise to SIRC as an independent DNA element. SIRC sites were significantly enriched with several histone modifications associated with constitutive heterochromatin and mobile genetic elements. The majority of DNA-binding proteins, strongly associated with SIRC, are related to histone modifications for transcription repression. A part of SIRC was found to overlap highly inducible protein-coding genes, suggesting a possible regulatory role for these elements, yet their definitive functions need further investigation.

Increased mutation and gene conversion within human segmental duplications

Single-nucleotide variants (SNVs) in segmental duplications (SDs) have not been systematically assessed because of the limitations of mapping short-read sequencing data1,2. Here we constructed 1:1 unambiguous alignments spanning high-identity SDs across 102 human haplotypes and compared the pattern of SNVs between unique and duplicated regions3,4. We find that human SNVs are elevated 60% in SDs compared to unique regions and estimate that at least 23% of this increase is due to interlocus gene conversion (IGC) with up to 4.3 megabase pairs of SD sequence converted on average per human haplotype. We develop a genome-wide map of IGC donors and acceptors, including 498 acceptor and 454 donor hotspots affecting the exons of about 800 protein-coding genes. These include 171 genes that have 'relocated' on average 1.61 megabase pairs in a subset of human haplotypes. Using a coalescent framework, we show that SD regions are slightly evolutionarily older when compared to unique sequences, probably owing to IGC. SNVs in SDs, however, show a distinct mutational spectrum: a 27.1% increase in transversions that convert cytosine to guanine or the reverse across all triplet contexts and a 7.6% reduction in the frequency of CpG-associated mutations when compared to unique DNA. We reason that these distinct mutational properties help to maintain an overall higher GC content of SD DNA compared to that of unique DNA, probably driven by GC-biased conversion between paralogous sequences5,6.

The fine-scale recombination rate variation and associations with genomic features in a butterfly

Recombination is a key molecular mechanism that has profound implications on both micro- and macroevolutionary processes. However, the determinants of recombination rate variation in holocentric organisms are poorly understood, in particular in Lepidoptera (moths and butterflies). The wood white butterfly (Leptidea sinapis) shows considerable intraspecific variation in chromosome numbers and is a suitable system for studying regional recombination rate variation and its potential molecular underpinnings. Here, we developed a large whole-genome resequencing data set from a population of wood whites to obtain high-resolution recombination maps using linkage disequilibrium information. The analyses revealed that larger chromosomes had a bimodal recombination landscape, potentially caused by interference between simultaneous chiasmata. The recombination rate was significantly lower in subtelomeric regions, with exceptions associated with segregating chromosome rearrangements, showing that fissions and fusions can have considerable effects on the recombination landscape. There was no association between the inferred recombination rate and base composition, supporting a limited influence of GC-biased gene conversion in butterflies. We found significant but variable associations between the recombination rate and the density of different classes of transposable elements, most notably a significant enrichment of short interspersed nucleotide elements in genomic regions with higher recombination rate. Finally, the analyses unveiled significant enrichment of genes involved in farnesyltranstransferase activity in recombination coldspots, potentially indicating that expression of transferases can inhibit formation of chiasmata during meiotic division. Our results provide novel information about recombination rate variation in holocentric organisms and have particular implications for forthcoming research in population genetics, molecular/genome evolution, and speciation.

The enormous repetitive Antarctic krill genome reveals environmental adaptations and population insights

Antarctic krill (Euphausia superba) is Earth's most abundant wild animal, and its enormous biomass is vital to the Southern Ocean ecosystem. Here, we report a 48.01-Gb chromosome-level Antarctic krill genome, whose large genome size appears to have resulted from inter-genic transposable element expansions. Our assembly reveals the molecular architecture of the Antarctic krill circadian clock and uncovers expanded gene families associated with molting and energy metabolism, providing insights into adaptations to the cold and highly seasonal Antarctic environment. Population-level genome re-sequencing from four geographical sites around the Antarctic continent reveals no clear population structure but highlights natural selection associated with environmental variables. An apparent drastic reduction in krill population size 10 mya and a subsequent rebound 100 thousand years ago coincides with climate change events. Our findings uncover the genomic basis of Antarctic krill adaptations to the Southern Ocean and provide valuable resources for future Antarctic research.

Mitogenome-based phylogenomics provides insights into the positions of the enigmatic sinensis group and the sanguinolenta group in Selaginellaceae (Lycophyte)

Spikemoss (Selaginellaceae) is one of the basal lineages of vascular plants. This family has a single genus Selaginella which consists of about 750 extant species. The phylogeny of Selaginellaceae has been extensively studied mainly based on plastid DNA and a few nuclear sequences. However, the placement of the enigmatic sinensis group is a long-term controversy because of the long branch in the plastid DNA phylogeny. The sanguinolenta group is also a phylogenetically problematic clade owing to two alternative positions resulted from different datasets. Here, we newly sequenced 34 mitochondrial genomes (mitogenomes) of individuals representing all seven subgenera and major clades in Selaginellaceae. We assembled the draft mitogenomes and annotated the genes and performed phylogenetic analyses based on the shared 17 mitochondrial genes. Our major results include: (1) all the assembled mitogenomes have complicated structures, unparalleled high GC content and a small gene content set, and the positive correlations among GC content, substitution rates and the number of RNA editing sites hold; (2) the sinensis group was well supported as a member of subg. Stachygynandrum; (3) the sanguinolenta group was strongly resolved as sister to all other Selaginella species except for subg. Selaginella. This study demonstrates the potential of mitogenome data in providing novel insights into phylogenetically recalcitrant problems.

Very Low Rates of Spontaneous Gene Deletions and Gene Duplications in Dictyostelium discoideum

The study of spontaneous mutation rates has revealed a wide range of heritable point mutation rates across species, but there are comparatively few estimates for large-scale deletion and duplication rates. The handful of studies that have directly calculated spontaneous rates of deletion and duplication using mutation accumulation lines have estimated that genes are duplicated and deleted at orders of magnitude greater rates than the spontaneous point mutation rate. In our study, we tested whether spontaneous gene deletion and gene duplication rates are also high in Dictyostelium discoideum, a eukaryote with among the lowest point mutation rates (2.5 × 10-11 per site per generation) and an AT-rich genome (GC content of 22%). We calculated mutation rates of gene deletions and duplications using whole-genome sequencing data originating from a mutation accumulation experiment and determined the association between the copy number mutations and GC content. Overall, we estimated an average of 3.93 × 10-8 gene deletions and 1.18 × 10-8 gene duplications per gene per generation. While orders of magnitude greater than their point mutation rate, these rates are much lower compared to gene deletion and duplication rates estimated from mutation accumulation lines in other organisms (that are on the order of ~ 10-6 per gene/generation). The deletions and duplications were enriched in regions that were AT-rich even compared to the genomic background, in contrast to our expectations if low GC content was contributing to low mutation rates. The low deletion and duplication mutation rates in D. discoideum compared to other eukaryotes mirror their low point mutation rates, supporting previous work suggesting that this organism has high replication fidelity and effective molecular machinery to avoid the accumulation of mutations in their genome.

Shuffling the yeast genome using CRISPR/Cas9-generated DSBs that target the transposable Ty1 elements

Although homologous recombination between transposable elements can drive genomic evolution in yeast by facilitating chromosomal rearrangements, the details of the underlying mechanisms are not fully clarified. In the genome of the yeast Saccharomyces cerevisiae, the most common class of transposon is the retrotransposon Ty1. Here, we explored how Cas9-induced double-strand breaks (DSBs) directed to Ty1 elements produce genomic alterations in this yeast species. Following Cas9 induction, we observed a significant elevation of chromosome rearrangements such as deletions, duplications and translocations. In addition, we found elevated rates of mitotic recombination, resulting in loss of heterozygosity. Using Southern analysis coupled with short- and long-read DNA sequencing, we revealed important features of recombination induced in retrotransposons. Almost all of the chromosomal rearrangements reflect the repair of DSBs at Ty1 elements by non-allelic homologous recombination; clustered Ty elements were hotspots for chromosome rearrangements. In contrast, a large proportion (about three-fourths) of the allelic mitotic recombination events have breakpoints in unique sequences. Our analysis suggests that some of the latter events reflect extensive processing of the broken ends produced in the Ty element that extend into unique sequences resulting in break-induced replication. Finally, we found that haploid and diploid strain have different preferences for the pathways used to repair double-stranded DNA breaks. Our findings demonstrate the importance of DNA lesions in retrotransposons in driving genome evolution.

The mutational dynamics of short tandem repeats in large, multigenerational families

BACKGROUND

Short tandem repeats (STRs) compose approximately 3% of the genome, and mutations at STR loci have been linked to dozens of human diseases including amyotrophic lateral sclerosis, Friedreich ataxia, Huntington disease, and fragile X syndrome. Improving our understanding of these mutations would increase our knowledge of the mutational dynamics of the genome and may uncover additional loci that contribute to disease. To estimate the genome-wide pattern of mutations at STR loci, we analyze blood-derived whole-genome sequencing data for 544 individuals from 29 three-generation CEPH pedigrees. These pedigrees contain both sets of grandparents, the parents, and an average of 9 grandchildren per family.

RESULTS

We use HipSTR to identify de novo STR mutations in the 2nd generation of these pedigrees and require transmission to the third generation for validation. Analyzing approximately 1.6 million STR loci, we estimate the empirical de novo STR mutation rate to be 5.24 × 10^-5 mutations per locus per generation. Perfect repeats mutate about 2 × more often than imperfect repeats. De novo STRs are significantly enriched in Alu elements.

CONCLUSIONS

Approximately 30% of new STR mutations occur within Alu elements, which compose only 11% of the genome, but only 10% are found in LINE-1 insertions, which compose 17% of the genome. Phasing these mutations to the parent of origin shows that parental transmission biases vary among families. We estimate the average number of de novo genome-wide STR mutations per individual to be approximately 85, which is similar to the average number of observed de novo single nucleotide variants.

Recombinant Adeno-Associated Virus-Mediated Editing of the G551D Cystic Fibrosis Transmembrane Conductance Regulator Mutation in Ferret Airway Basal Cells

Mutations in the cystic fibrosis transmembrane conductance regulator (CFTR) gene cause cystic fibrosis (CF), a chronic disease that affects multiple organs, including the lung. We developed a CF ferret model of a scarless G551→D substitution in CFTR (CFTRG551D-KI), enabling approaches to correct this gating mutation in CF airways via gene editing. Homology-directed repair (HDR) was tested in Cas9-expressing CF airway basal cells (Cas9-GKI) from this model, as well as reporter basal cells (Y66S-Cas9-GKI) that express an integrated nonfluorescent Y66S-EGFP (enhanced green fluorescent protein) mutant gene to facilitate rapid assessment of HDR by the restoration of fluorescence. Recombinant adeno-associated virus (rAAV) vectors were used to deliver two DNA templates and sgRNAs for dual-gene editing at the EGFP and CFTR genes, followed by fluorescence-activated cell sorting of EGFPY66S-corrected cells. When gene-edited airway basal cells were polarized at an air-liquid interface, unsorted and EGFPY66S-corrected sorted populations gave rise to 26.0% and 70.4% CFTR-mediated Cl- transport of that observed in non-CF cultures, respectively. The consequences of gene editing at the CFTRG551D locus by HDR and nonhomologous end joining (NHEJ) were assessed by targeted gene next-generation sequencing (NGS) against a specific amplicon. NGS revealed HDR corrections of 3.1% of G551 sequences in the unsorted population of rAAV-infected cells, and 18.4% in the EGFPY66S-corrected cells. However, the largest proportion of sequences had indels surrounding the CRISPR (clustered regularly interspaced short palindromic repeats) cut site, demonstrating that NHEJ was the dominant repair pathway. This approach to simultaneously coedit at two genomic loci using rAAV may have utility as a model system for optimizing gene-editing efficiencies in proliferating airway basal cells through the modulation of DNA repair pathways in favor of HDR.

Recent insights into the evolution of mutation rates in yeast

Mutation is the origin of all genetic variation, good and bad. The mutation process can evolve in response to mutations, positive or negative selection, and genetic drift, but how these forces contribute to mutation-rate variation is an unsolved problem at the heart of genetics research. Mutations can be challenging to measure, but genome sequencing and other tools have allowed for the collection of larger and more detailed datasets, particularly in the yeast-model system. We review key hypotheses for the evolution of mutation rates and describe recent advances in understanding variation in mutational properties within and among yeast species. The multidimensional spectrum of mutations is increasingly recognized as holding valuable clues about how this important process evolves.

Host translesion polymerases are required for viral genome integrity

Human cells encode up to 15 DNA polymerases with specialized functions in chromosomal DNA synthesis and damage repair. In contrast, complex DNA viruses, such as those of the herpesviridae family, encode a single B-family DNA polymerase. This disparity raises the possibility that DNA viruses may rely on host polymerases for synthesis through complex DNA geometries. We tested the importance of error-prone Y-family polymerases involved in translesion synthesis (TLS) to human cytomegalovirus (HCMV) infection. We find most Y-family polymerases involved in the nucleotide insertion and bypass of lesions restrict HCMV genome synthesis and replication. In contrast, other TLS polymerases, such as the polymerase ζ complex, which extends past lesions, was required for optimal genome synthesis and replication. Depletion of either the polζ complex or the suite of insertion polymerases demonstrate that TLS polymerases suppress the frequency of viral genome rearrangements, particularly at GC-rich sites and repeat sequences. Moreover, while distinct from HCMV, replication of the related herpes simplex virus type 1 is impacted by host TLS polymerases, suggesting a broader requirement for host polymerases for DNA virus replication. These findings reveal an unexpected role for host DNA polymerases in ensuring viral genome stability.

Is the mutation rate lower in genomic regions of stronger selective constraints?

A study of the plant Arabidopsis thaliana detected lower mutation rates in genomic regions where mutations are more likely to be deleterious, challenging the principle that mutagenesis is blind to its consequence. To examine the generality of this finding, we analyze large mutational data from baker's yeast and humans. The yeast data do not exhibit this trend, whereas the human data show an opposite trend that disappears upon the control of potential confounders. We find that the Arabidopsis study identified substantially more mutations than reported in the original data-generating studies and expected from Arabidopsis' mutation rate. These extra mutations are enriched in polynucleotide tracts and have relatively low sequencing qualities so are likely sequencing errors. Furthermore, the polynucleotide “mutations” can produce the purported mutational trend in Arabidopsis. Together, our results do not support lower mutagenesis of genomic regions of stronger selective constraints in the plant, fungal, and animal models examined.

Plastid genome evolution of a monophyletic group in the subtribe Lauriineae (Laureae, Lauraceae)

Litsea, a non-monophyletic group of the tribe Laureae (Lauraceae), plays important roles in the tropical and subtropical forests of Asia, Australia, Central and North America, and the islands of the Pacific. However, intergeneric relationships between Litsea and Laurus, Lindera, Parasassafras and Sinosassafras of the tribe Laureae remain unresolved. In this study, we present phylogenetic analyses of seven newly sequenced Litsea plastomes, together with 47 Laureae plastomes obtained from public databases, representing six genera of the Laureae. Our results highlight two highly supported monophyletic groups of Litsea taxa. One is composed of 16 Litsea taxa and two Lindera taxa. The 18 plastomes of these taxa were further compared for their gene structure, codon usage, contraction and expansion of inverted repeats, sequence repeats, divergence hotspots, and gene evolution. The complete plastome size of newly sequenced taxa varied between 152,377 bp (Litsea auriculata) and 154,117 bp (Litsea pierrei). Seven of the 16 Litsea plastomes have a pair of insertions in the IRa (trnL-trnH) and IRb (ycf2) regions. The 18 plastomes of Litsea and Lindera taxa exhibit similar gene features, codon usage, oligonucleotide repeats, and inverted repeat dynamics. The codons with the highest frequency among these taxa favored A/T endings and each of these plastomes had nine divergence hotspots, which are located in the same regions. We also identified six protein coding genes (accD, ndhJ, rbcL, rpoC2, ycf1 and ycf2) under positive selection in Litsea; these genes may play important roles in adaptation of Litsea species to various environments.

Patterns and distribution of de novo mutations in multiplex Middle Eastern families

While de novo mutations (DNMs) are key to genetic diversity, they are also responsible for a high number of rare disorders. To date, no study has systematically examined the rate and distribution of DNMs in multiplex families in highly consanguineous populations. Leveraging WGS profiles of 645 individuals in 146 families, we implemented a combinatorial approach using 3 complementary tools for DNM discovery in 353 unique trio combinations. We found a total of 27,168 DNMs (median: 70 single-nucleotide and 6 insertion-deletions per individual). Phasing revealed around 80% of DNMs were paternal in origin. Notably, using whole-genome methylation data of spermatogonial stem cells, these DNMs were significantly more likely to occur at highly methylated CpGs (OR: 2.03; p value = 6.62 × 10⁻¹¹). We then examined the effects of consanguinity and ethnicity on DNMs, and found that consanguinity does not seem to correlate with DNM rate, and special attention has to be considered while measuring such a correlation. Additionally, we found that Middle-Eastern families with Arab ancestry had fewer DNMs than African families, although not significant (p value = 0.16). Finally, for families with diseased probands, we examined the difference in DNM counts and putative impact across affected and unaffected siblings, but did not find significant differences between disease groups, likely owing to the enrichment for recessive disorders in this part of the world, or the small sample size per clinical condition. This study serves as a reference for DNM discovery in multiplex families from the globally under-represented populations of the Middle-East.

Unusual mammalian usage of TGA stop codons reveals that sequence conservation need not imply purifying selection

The assumption that conservation of sequence implies the action of purifying selection is central to diverse methodologies to infer functional importance. GC-biased gene conversion (gBGC), a meiotic mismatch repair bias strongly favouring GC over AT, can in principle mimic the action of selection, this being thought to be especially important in mammals. As mutation is GC→AT biased, to demonstrate that gBGC does indeed cause false signals requires evidence that an AT-rich residue is selectively optimal compared to its more GC-rich allele, while showing also that the GC-rich alternative is conserved. We propose that mammalian stop codon evolution provides a robust test case. Although in most taxa TAA is the optimal stop codon, TGA is both abundant and conserved in mammalian genomes. We show that this mammalian exceptionalism is well explained by gBGC mimicking purifying selection and that TAA is the selectively optimal codon. Supportive of gBGC, we observe (i) TGA usage trends are consistent at the focal stop codon and elsewhere (in UTR sequences); (ii) that higher TGA usage and higher TAA→TGA substitution rates are predicted by a high recombination rate; and (iii) across species the difference in TAA <-> TGA substitution rates between GC-rich and GC-poor genes is largest in genomes that possess higher between-gene GC variation. TAA optimality is supported both by enrichment in highly expressed genes and trends associated with effective population size. High TGA usage and high TAA→TGA rates in mammals are thus consistent with gBGC’s predicted ability to “drive” deleterious mutations and supports the hypothesis that sequence conservation need not be indicative of purifying selection. A general trend for GC-rich trinucleotides to reside at frequencies far above their mutational equilibrium in high recombining domains supports the generality of these results.

Recessive cerebellar and afferent ataxias - clinical challenges and future directions

Cerebellar and afferent ataxias present with a characteristic gait disorder that reflects cerebellar motor dysfunction and sensory loss. These disorders are a diagnostic challenge for clinicians because of the large number of acquired and inherited diseases that cause cerebellar and sensory neuron damage. Among such conditions that are recessively inherited, Friedreich ataxia and RFC1-associated cerebellar ataxia, neuropathy, vestibular areflexia syndrome (CANVAS) include the characteristic clinical, neuropathological and imaging features of ganglionopathies, a distinctive non-length-dependent type of sensory involvement. In this Review, we discuss the typical and atypical phenotypes of Friedreich ataxia and CANVAS, along with the features of other recessive ataxias that present with a ganglionopathy or polyneuropathy, with an emphasis on recently described clinical features, natural history and genotype-phenotype correlations. We review the main developments in understanding the complex pathology that affects the sensory neurons and cerebellum, which seem to be most vulnerable to disorders that affect mitochondrial function and DNA repair mechanisms. Finally, we discuss disease-modifying therapeutic advances in Friedreich ataxia, highlighting the most promising candidate molecules and lessons learned from previous clinical trials.

Evolutionary Stability Optimizer (ESO): A Novel Approach to Identify and Avoid Mutational Hotspots in DNA Sequences While Maintaining High Expression Levels

Modern synthetic biology procedures rely on the ability to generate stable genetic constructs that keep their functionality over long periods of time. However, maintenance of these constructs requires energy from the cell and thus reduces the host’s fitness. Natural selection results in loss-of-functionality mutations that negate the expression of the construct in the population. Current approaches for the prevention of this phenomenon focus on either small-scale, manual design of evolutionary stable constructs or the detection of mutational sites with unstable tendencies. We designed the Evolutionary Stability Optimizer (ESO), a software tool that enables the large-scale automatic design of evolutionarily stable constructs with respect to both mutational and epigenetic hotspots and allows users to define custom hotspots to avoid. Furthermore, our tool takes the expression of the input constructs into account by considering the guanine-cytosine (GC) content and codon usage of the host organism, balancing the trade-off between stability and gene expression, allowing to increase evolutionary stability while maintaining the high expression. In this study, we present the many features of the ESO and show that it accurately predicts the evolutionary stability of endogenous genes. The ESO was created as an easy-to-use, flexible platform based on the notion that directed genetic stability research will continue to evolve and revolutionize current applications of synthetic biology. The ESO is available at the following link: https://www.cs.tau.ac.il/~tamirtul/ESO/.

Unusual clinical manifestations and predominant stopgain ATM gene variants in a single centre cohort of ataxia telangiectasia from North India

Germline ATM gene variations result in phenotypic heterogeneity characterized by a variable degree of disease severity. We retrospectively collected clinical, genetic, and immunological data of 26 cases with A-T. Clinical manifestations included oculocutaneous telangiectasia (100%), ataxia (100%), fever, loose stools or infection (67%), cerebellar atrophy (50%), nystagmus (8%), dysarthria (15.38%), and visual impairment (8%). Genetic analysis confirmed ATM gene variations in 16 unrelated cases. The most common type of variation was stopgain variants (56%). Immunoglobulin profile indicated reduced IgA, IgG, and IgM in 94%, 50%, and 20% cases, respectively. T cell lymphopenia was observed in 80% of cases among those investigated. Unusual presentations included an EBV-associated smooth muscle tumour located in the liver in one case and Hyper IgM syndrome-like presentation in two cases. Increased immunosenescence was observed in T-cell subsets (CD4+CD57+ and CD8+CD57+). T-cell receptor excision circles (TRECs) were reduced in 3/8 (37.50%) cases.

Sex-specific splicing of Z- and W-borne nr5a1 alleles suggests sex determination is controlled by chromosome conformation

Pogona vitticeps has female heterogamety (ZZ/ZW), but the master sex-determining gene is unknown, as it is for all reptiles. We show that nr5a1 (Nuclear Receptor Subfamily 5 Group A Member 1), a gene that is essential in mammalian sex determination, has alleles on the Z and W chromosomes (Z-nr5a1 and W-nr5a1), which are both expressed and can recombine. Three transcript isoforms of Z-nr5a1 were detected in gonads of adult ZZ males, two of which encode a functional protein. However, ZW females produced 16 isoforms, most of which contained premature stop codons. The array of transcripts produced by the W-borne allele (W-nr5a1) is likely to produce truncated polypeptides that contain a structurally normal DNA-binding domain and could act as a competitive inhibitor to the full-length intact protein. We hypothesize that an altered configuration of the W chromosome affects the conformation of the primary transcript generating inhibitory W-borne isoforms that suppress testis determination. Under this hypothesis, the genetic sex determination (GSD) system of P. vitticeps is a W-borne dominant female-determining gene that may be controlled epigenetically.

Cloning and codon optimization of a novel feline interferon omega gene for production by Pichia pastoris and its antiviral efficacy in polyethylene glycol-modified form

Feline viral diseases, such as feline panleukopenia, feline infectious peritonitis, and feline coronaviral enteritis, seriously endanger the health of cats, and restrict the development of pet industry. Meanwhile, there is a current lack of effective vaccines to protect against feline viral diseases. Thus, effective therapeutic agents are highly desirable. Interferons (IFNs) are important mediators of the antiviral host defense in animals, particularly type I IFNs. In this study, a novel feline IFN omega (feIFN-ω) gene was extracted from the cat stimulated with feline parvovirus (FPV) combined with poly(I:C), and following codon optimization encoding the feIFN-ω, the desired gene (feIFN-ω’) fragment was inserted into plasmid pPICZαA, and transformed into Pichia pastoris GS115, generating a recombinant P. pastoris GS115 strain expressing the feIFN-ω’. After induction, we found that the expression level of the feIFN-ω’ was two times more than that of feIFN-ω (p < 0.01). Subsequently, the feIFN-ω’ was purified and modified with polyethylene glycol, and its antiviral efficacy was evaluated in vitro and in vivo, using vesicular stomatitis virus (VSV) and FPV as model virus. Our results clearly demonstrated that the feIFN-ω’ had significant antiviral activities on both homologous and heterologous animal cells in vitro. Importantly, the feIFN-ω’ can effectively promote the expression of antiviral proteins IFIT3, ISG15, Mx1, and ISG56, and further enhance host defense to eliminate FPV infection in vivo, suggesting a potential candidate for the development of therapeutic agent against feline viral diseases.

From Genome Variation to Molecular Mechanisms: What we Have Learned From Yeast Mitochondrial Genomes?

Analysis of genome variation provides insights into mechanisms in genome evolution. This is increasingly appreciated with the rapid growth of genomic data. Mitochondrial genomes (mitogenomes) are well known to vary substantially in many genomic aspects, such as genome size, sequence context, nucleotide base composition and substitution rate. Such substantial variation makes mitogenomes an excellent model system to study the mechanisms dictating mitogenome variation. Recent sequencing efforts have not only covered a rich number of yeast species but also generated genomes from abundant strains within the same species. The rich yeast genomic data have enabled detailed investigation from genome variation into molecular mechanisms in genome evolution. This mini-review highlights some recent progresses in yeast mitogenome studies.

Transposon-induced epigenetic silencing in the X chromosome as a novel form of dmrt1 expression regulation during sex determination in the fighting fish

Background

Fishes are the one of the most diverse groups of animals with respect to their modes of sex determination, providing unique models for uncovering the evolutionary and molecular mechanisms underlying sex determination and reversal. Here, we have investigated how sex is determined in a species of both commercial and ecological importance, the Siamese fighting fish Betta splendens.

Results

We conducted association mapping on four commercial and two wild populations of B. splendens. In three of the four commercial populations, the master sex determining (MSD) locus was found to be located in a region of ~ 80 kb on LG2 which harbours five protein coding genes, including dmrt1, a gene involved in male sex determination in different animal taxa. In these fish, dmrt1 shows a male-biased gonadal expression from undifferentiated stages to adult organs and the knockout of this gene resulted in ovarian development in XY genotypes. Genome sequencing of XX and YY genotypes identified a transposon, drbx1, inserted into the fourth intron of the X-linked dmrt1 allele. Methylation assays revealed that epigenetic changes induced by drbx1 spread out to the promoter region of dmrt1. In addition, drbx1 being inserted between two closely linked cis-regulatory elements reduced their enhancer activities. Thus, epigenetic changes, induced by drbx1, contribute to the reduced expression of the X-linked dmrt1 allele, leading to female development. This represents a previously undescribed solution in animals relying on dmrt1 function for sex determination. Differentiation between the X and Y chromosomes is limited to a small region of ~ 200 kb surrounding the MSD gene. Recombination suppression spread slightly out of the SD locus. However, this mechanism was not found in the fourth commercial stock we studied, or in the two wild populations analysed, suggesting that it originated recently during domestication.

Conclusions

Taken together, our data provide novel insights into the role of epigenetic regulation of dmrt1 in sex determination and turnover of SD systems and suggest that fighting fish are a suitable model to study the initial stages of sex chromosome evolution.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01205-y.

Deviation from baseline mutation burden provides powerful and robust rare-variants association test for complex diseases

Identifying rare variants that contribute to complex diseases is challenging because of the low statistical power in current tests comparing cases with controls. Here, we propose a novel and powerful rare variants association test based on the deviation of the observed mutation burden of a gene in cases from a baseline predicted by a weighted recursive truncated negative-binomial regression (RUNNER) on genomic features available from public data. Simulation studies show that RUNNER is substantially more powerful than state-of-the-art rare variant association tests and has reasonable type 1 error rates even for stratified populations or in small samples. Applied to real case-control data, RUNNER recapitulates known genes of Hirschsprung disease and Alzheimer's disease missed by current methods and detects promising new candidate genes for both disorders. In a case-only study, RUNNER successfully detected a known causal gene of amyotrophic lateral sclerosis. The present study provides a powerful and robust method to identify susceptibility genes with rare risk variants for complex diseases.

Distinct sequence features underlie microdeletions and gross deletions in the human genome

Microdeletions and gross deletions are important causes (~20%) of human inherited disease and their genomic locations are strongly influenced by the local DNA sequence environment. This notwithstanding, no study has systematically examined their underlying generative mechanisms. Here, we obtained 42,098 pathogenic microdeletions and gross deletions from the Human Gene Mutation Database (HGMD) that together form a continuum of germline deletions ranging in size from 1 to 28,394,429 bp. We analyzed the DNA sequence within 1 kb of the breakpoint junctions and found that the frequencies of non‐B DNA‐forming repeats, GC‐content, and the presence of seven of 78 specific sequence motifs in the vicinity of pathogenic deletions correlated with deletion length for deletions of length ≤30 bp. Further, we found that the presence of DR, GQ, and STR repeats is important for the formation of longer deletions (>30 bp) but not for the formation of shorter deletions (≤30 bp) while significantly (χ², p < 2E−16) more microhomologies were identified flanking short deletions than long deletions (length >30 bp). We provide evidence to support a functional distinction between microdeletions and gross deletions. Finally, we propose that a deletion length cut‐off of 25–30 bp may serve as an objective means to functionally distinguish microdeletions from gross deletions.

The Complexity of the Ovine and Caprine Keratin-Associated Protein Genes

Sheep (Ovis aries) and goats (Capra hircus) have, for more than a millennia, been a source of fibres for human use, be it for use in clothing and furnishings, for insulation, for decorative and ceremonial purposes, or for combinations thereof. While use of these natural fibres has in some respects been superseded by the use of synthetic and plant-based fibres, increased accounting for the carbon and water footprint of these fibres is creating a re-emergence of interest in fibres derived from sheep and goats. The keratin-associated proteins (KAPs) are structural components of wool and hair fibres, where they form a matrix that cross-links with the keratin intermediate filaments (KIFs), the other main structural component of the fibres. Since the first report of a complete KAP protein sequence in the late 1960s, considerable effort has been made to identify the KAP proteins and their genes in mammals, and to ascertain how these genes and proteins control fibre growth and characteristics. This effort is ongoing, with more and more being understood about the structure and function of the genes. This review consolidates that knowledge and suggests future directions for research to further our understanding.

GC-biased gene conversion in X-chromosome palindromes conserved in human, chimpanzee, and rhesus macaque

Gene conversion is GC-biased across a wide range of taxa. Large palindromes on mammalian sex chromosomes undergo frequent gene conversion that maintains arm-to-arm sequence identity greater than 99%, which may increase their susceptibility to the effects of GC-biased gene conversion. Here, we demonstrate a striking history of GC-biased gene conversion in 12 palindromes conserved on the X chromosomes of human, chimpanzee, and rhesus macaque. Primate X-chromosome palindrome arms have significantly higher GC content than flanking single-copy sequences. Nucleotide replacements that occurred in human and chimpanzee palindrome arms over the past 7 million years are one-and-a-half times as GC-rich as the ancestral bases they replaced. Using simulations, we show that our observed pattern of nucleotide replacements is consistent with GC-biased gene conversion with a magnitude of 70%, similar to previously reported values based on analyses of human meioses. However, GC-biased gene conversion since the divergence of human and rhesus macaque explains only a fraction of the observed difference in GC content between palindrome arms and flanking sequence, suggesting that palindromes are older than 29 million years and/or had elevated GC content at the time of their formation. This work supports a greater than 2:1 preference for GC bases over AT bases during gene conversion and demonstrates that the evolution and composition of mammalian sex chromosome palindromes is strongly influenced by GC-biased gene conversion.

Molecular landscape of DYSF mutations in dysferlinopathy: From a Chinese multicenter analysis to a worldwide perspective

Dysferlinopathy is one of the most common subgroup of autosomal recessive limb-girdle muscular dystrophies that is caused by mutations in DYSF gene. However, there is currently no worldwide comprehensive genetic analysis of DYSF variants. Through a national multicenter collaborative effort in China, we identified 222 DYSF variants with 40 novel variants from 245 patients. We then integrated DYSF variants from disease-related genetic databases including LOVD (n = 1020) and Clinvar (n = 1179), to depict the global landscape of disease-related DYSF variants. Normal-population-derived DSYF variants from gnomAD (n = 4318) and ChinaMAP (n = 13,330) were also analyzed in comparison. In Chinese patients, gender instead of genotype showed influence on the onset age of dysferlinopathy, with males showing an earlier age of onset. After integrative analysis, we identified two hotspot DYSF mutations, c.2997G>T in world patients and c.1375dup in Chinese patients, respectively. Both the pathogenic and likely pathogenic variants scattered on the whole gene length of DYSF. However, three specific domains (C2F-C2G-TM, DysF, and C2B-Ferl-C2C) contained variants at higher frequencies than reported in both the databases and Chinese patients. This study comprehensively collected available DYSF variant data, which may pave way for genetic counselling and future clinical trial design for gene therapies in dysferlinopathy.

Gene expression levels modulate germline mutation rates through the compound effects of transcription-coupled repair and damage

Of all mammalian organs, the testis has long been observed to have the most diverse gene expression profile. To account for this widespread gene expression, we have proposed a mechanism termed 'transcriptional scanning', which reduces germline mutation rates through transcription-coupled repair (TCR). Our hypothesis contrasts with an earlier observation that mutation rates are overall positively correlated with gene expression levels in yeast, implying that transcription is mutagenic due to effects dominated by transcription-coupled damage (TCD). Here we report evidence that the compound effects of both TCR and TCD during spermatogenesis modulate human germline mutation rates, with TCR dominating in most genes, thus supporting the transcriptional scanning hypothesis. Our analyses address potentially confounding factors, distinguish the differential mutagenic effects acting on the highly expressed genes and the low-to-moderately expressed genes, and resolve concerns relating to the validation of the results using a de novo mutation dataset. We also discuss the theoretical possibility of transcriptional scanning hypothesis from an evolutionary perspective. Together, these analyses support a model by which the coupling of transcription-coupled repair and damage establishes the pattern of germline mutation rates and provide an evolutionary explanation for widespread gene expression during spermatogenesis.

The fidelity of DNA replication, particularly on GC-rich templates, is reduced by defects of the Fe-S cluster in DNA polymerase δ

Iron-sulfur clusters (4Fe–4S) exist in many enzymes concerned with DNA replication and repair. The contribution of these clusters to enzymatic activity is not fully understood. We identified the MET18 (MMS19) gene of Saccharomyces cerevisiae as a strong mutator on GC-rich genes. Met18p is required for the efficient insertion of iron-sulfur clusters into various proteins. met18 mutants have an elevated rate of deletions between short flanking repeats, consistent with increased DNA polymerase slippage. This phenotype is very similar to that observed in mutants of POL3 (encoding the catalytic subunit of Pol δ) that weaken binding of the iron-sulfur cluster. Comparable mutants of POL2 (Pol ϵ) do not elevate deletions. Further support for the conclusion that met18 strains result in impaired DNA synthesis by Pol δ are the observations that Pol δ isolated from met18 strains has less bound iron and is less processive in vitro than the wild-type holoenzyme.

Rtt105 promotes high-fidelity DNA replication and repair by regulating the single-stranded DNA-binding factor RPA

Single-stranded DNA (ssDNA) covered with the heterotrimeric Replication Protein A (RPA) complex is a central intermediate of DNA replication and repair. How RPA is regulated to ensure the fidelity of DNA replication and repair remains poorly understood. Yeast Rtt105 is an RPA-interacting protein required for RPA nuclear import and efficient ssDNA binding. Here, we describe an important role of Rtt105 in high-fidelity DNA replication and recombination and demonstrate that these functions of Rtt105 primarily depend on its regulation of RPA. The deletion of RTT105 causes elevated spontaneous DNA mutations with large duplications or deletions mediated by microhomologies. Rtt105 is recruited to DNA double-stranded break (DSB) ends where it promotes RPA assembly and homologous recombination repair by gene conversion or break-induced replication. In contrast, Rtt105 attenuates DSB repair by the mutagenic single-strand annealing or alternative end joining pathway. Thus, Rtt105-mediated regulation of RPA promotes high-fidelity replication and recombination while suppressing repair by deleterious pathways. Finally, we show that the human RPA-interacting protein hRIP-α, a putative functional homolog of Rtt105, also stimulates RPA assembly on ssDNA, suggesting the conservation of an Rtt105-mediated mechanism.

Insight into Codon Utilization Pattern of Tumor Suppressor Gene EPB41L3 from Different Mammalian Species Indicates Dominant Role of Selection Force

Simple Summary

The present study envisaged the codon usage pattern analysis of tumor suppressor gene EPB41L3 for the human, brown rat, domesticated cattle, and Sumatran orangutan. Most amino acids are coded by more than one synonymous codon, but they are used in a biased manner. The codon usage bias results from multiple factors like compositional properties, dinucleotide abundance, neutrality, parity, tRNA pool, etc. Understanding codon bias is central to fields as diverse as molecular evolution, gene expressivity, protein translation, and protein folding. This kind of studies is important to see the effects of various evolutionary forces on codon usage. The present study indicated that the selection force is dominant over other forces shaping codon usage in the envisaged organisms.

Abstract

Uneven codon usage within genes as well as among genomes is a usual phenomenon across organisms. It plays a significant role in the translational efficiency and evolution of a particular gene. EPB41L3 is a tumor suppressor protein-coding gene, and in the present study, the pattern of codon usage was envisaged. The full-length sequences of the EPB41L3 gene for the human, brown rat, domesticated cattle, and Sumatran orangutan available at the NCBI were retrieved and utilized to analyze CUB patterns across the selected mammalian species. Compositional properties, dinucleotide abundance, and parity analysis showed the dominance of A and G whilst RSCU analysis indicated the dominance of G/C-ending codons. The neutrality plot plotted between GC12 and GC3 to determine the variation between the mutation pressure and natural selection indicated the dominance of selection pressure (R = 0.926; p < 0.00001) over the three codon positions across the gene. The result is in concordance with the codon adaptation index analysis and the ENc-GC3 plot analysis, as well as the translational selection index (P2). Overall selection pressure is the dominant pressure acting during the evolution of the EPB41L3 gene.

Cerebellar ataxia, neuropathy, vestibular areflexia syndrome: genetic and clinical insights

PURPOSE OF REVIEW

This review aims to summarise the present cerebellar ataxia, neuropathy, vestibular ataxia syndrome (CANVAS) literature, providing both clinical and genetic insights that might facilitate the timely clinical and genetic diagnosis of this disease.

RECENT FINDINGS

Recent advancements in the range of the clinical features of CANVAS have aided the development of a broader, more well-defined clinical diagnostic criteria. Additionally, the identification of a biallelic repeat expansion in RFC1 as the cause of CANVAS and a common cause of late-onset ataxia has opened the door to the potential discovery of a pathogenic mechanism, which in turn, may lead to therapeutic advancements and improved patient care.

SUMMARY

The developments in the clinical and genetic understanding of CANVAS will aid the correct and timely diagnosis of CANVAS, which continues to prove challenging within the clinic. The insights detailed within this review will raise the awareness of the phenotypic spectrum and currently known genetics. We also speculate on the future directions of research into CANVAS.

DNA Topoisomerase 3α Is Involved in Homologous Recombination Repair and Replication Stress Response in Trypanosoma cruzi

DNA topoisomerases are enzymes that modulate DNA topology. Among them, topoisomerase 3α is engaged in genomic maintenance acting in DNA replication termination, sister chromatid separation, and dissolution of recombination intermediates. To evaluate the role of this enzyme in Trypanosoma cruzi, the etiologic agent of Chagas disease, a topoisomerase 3α knockout parasite (TcTopo3α KO) was generated, and the parasite growth, as well as its response to several DNA damage agents, were evaluated. There was no growth alteration caused by the TcTopo3α knockout in epimastigote forms, but a higher dormancy rate was observed. TcTopo3α KO trypomastigote forms displayed reduced invasion rates in LLC-MK2 cells when compared with the wild-type lineage. Amastigote proliferation was also compromised in the TcTopo3α KO, and a higher number of dormant cells was observed. Additionally, TcTopo3α KO epimastigotes were not able to recover cell growth after gamma radiation exposure, suggesting the involvement of topoisomerase 3α in homologous recombination. These parasites were also sensitive to drugs that generate replication stress, such as cisplatin (Cis), hydroxyurea (HU), and methyl methanesulfonate (MMS). In response to HU and Cis treatments, TcTopo3α KO parasites showed a slower cell growth and was not able to efficiently repair the DNA damage induced by these genotoxic agents. The cell growth phenotype observed after MMS treatment was similar to that observed after gamma radiation, although there were fewer dormant cells after MMS exposure. TcTopo3α KO parasites showed a population with sub-G1 DNA content and strong γH2A signal 48 h after MMS treatment. So, it is possible that DNA-damaged cell proliferation due to the absence of TcTopo3α leads to cell death. Whole genome sequencing of MMS-treated parasites showed a significant reduction in the content of the multigene families DFG-1 and RHS, and also a possible erosion of the sub-telomeric region from chromosome 22, relative to non-treated knockout parasites. Southern blot experiments suggest telomere shortening, which could indicate genomic instability in TcTopo3α KO cells owing to MMS treatment. Thus, topoisomerase 3α is important for homologous recombination repair and replication stress in T. cruzi, even though all the pathways in which this enzyme participates during the replication stress response remains elusive.

Using GC Content to Compare Recombination Patterns on the Sex Chromosomes and Autosomes of the Guppy, Poecilia reticulata, and Its Close Outgroup Species

Genetic and physical mapping of the guppy (Poecilia reticulata) have shown that recombination patterns differ greatly between males and females. Crossover events occur evenly across the chromosomes in females, but in male meiosis they are restricted to the tip furthest from the centromere of each chromosome, creating very high recombination rates per megabase, as in pseudoautosomal regions of mammalian sex chromosomes. We used GC content to indirectly infer recombination patterns on guppy chromosomes, based on evidence that recombination is associated with GC-biased gene conversion, so that genome regions with high recombination rates should be detectable by high GC content. We used intron sequences and third positions of codons to make comparisons between sequences that are matched, as far as possible, and are all probably under weak selection. Almost all guppy chromosomes, including the sex chromosome (LG12), have very high GC values near their assembly ends, suggesting high recombination rates due to strong crossover localization in male meiosis. Our test does not suggest that the guppy XY pair has stronger crossover localization than the autosomes, or than the homologous chromosome in the close relative, the platyfish (Xiphophorus maculatus). We therefore conclude that the guppy XY pair has not recently undergone an evolutionary change to a different recombination pattern, or reduced its crossover rate, but that the guppy evolved Y-linkage due to acquiring a male-determining factor that also conferred the male crossover pattern. We also identify the centromere ends of guppy chromosomes, which were not determined in the genome assembly.

Evidence for a synergistic effect of post-translational modifications and genomic composition of eEF-1α on the adaptation of Phytophthora infestans

Genetic variation plays a fundamental role in pathogen's adaptation to environmental stresses. Pathogens with low genetic variation tend to survive and proliferate more poorly due to their lack of genotypic/phenotypic polymorphisms in responding to fluctuating environments. Evolutionary theory hypothesizes that the adaptive disadvantage of genes with low genomic variation can be compensated for structural diversity of proteins through post-translation modification (PTM) but this theory is rarely tested experimentally and its implication to sustainable disease management is hardly discussed. In this study, we analyzed nucleotide characteristics of eukaryotic translation elongation factor-1α (eEF-lα) gene from 165 Phytophthora infestans isolates and the physical and chemical properties of its derived proteins. We found a low sequence variation of eEF-lα protein, possibly attributable to purifying selection and a lack of intra-genic recombination rather than reduced mutation. In the only two isoforms detected by the study, the major one accounted for >95% of the pathogen collection and displayed a significantly higher fitness than the minor one. High lysine representation enhances the opportunity of the eEF-1α protein to be methylated and the absence of disulfide bonds is consistent with the structural prediction showing that many disordered regions are existed in the protein. Methylation, structural disordering, and possibly other PTMs ensure the ability of the protein to modify its functions during biological, cellular and biochemical processes, and compensate for its adaptive disadvantage caused by sequence conservation. Our results indicate that PTMs may function synergistically with nucleotide codes to regulate the adaptive landscape of eEF-1α, possibly as well as other housekeeping genes, in P. infestans. Compensatory evolution between pre- and post-translational phase in eEF-1α could enable pathogens quickly adapting to disease management strategies while efficiently maintaining critical roles of the protein playing in biological, cellular, and biochemical activities. Implications of these results to sustainable plant disease management are discussed.