Possible biased virulence attenuation in the Senegal strain of Ehrlichia ruminantium by ntrX gene conversion from an inverted segmental duplication.

Ehrlichia ruminantium is a tick-borne intracellular pathogen of ruminants that causes heartwater, a disease present in Sub-saharan Africa, islands in the Indian Ocean and the Caribbean, inducing significant economic losses. At present, three avirulent strains of E. ruminantium (Gardel, Welgevonden and Senegal isolates) have been produced by a process of serial passaging in mammalian cells in vitro, but unfortunately their use as vaccines do not offer a large range of protection against other strains, possibly due to the genetic diversity present within the species. So far no genetic basis for virulence attenuation has been identified in any E. ruminantium strain that could offer targets to facilitate vaccine production. Virulence attenuated Senegal strains have been produced twice independently, and require many fewer passages to attenuate than the other strains. We compared the genomes of a virulent and attenuated Senegal strain and identified a likely attenuator gene, ntrX, a global transcription regulator and member of a two-component system that is linked to environmental sensing. This gene has an inverted partial duplicate close to the parental gene that shows evidence of gene conversion in different E. ruminantium strains. The pseudogenisation of the gene in the avirulent Senegal strain occurred by gene conversion from the duplicate to the parent, transferring a 4 bp deletion which is unique to the Senegal strain partial duplicate amongst the wild isolates. We confirmed that the ntrX gene is not expressed in the avirulent Senegal strain by RT-PCR. The inverted duplicate structure combined with the 4 bp deletion in the Senegal strain can explain both the attenuation and the faster speed of attenuation in the Senegal strain relative to other strains of E. ruminantium. Our results identify nrtX as a promising target for the generation of attenuated strains of E. ruminantium by random or directed mutagenesis that could be used for vaccine production.

Compensatory Genetic and Transcriptional Cytonuclear Coordination in Allopolyploid Lager Yeast (Saccharomyces pastorianus).

Cytonuclear coordination between biparental-nuclear genomes and uniparental-cytoplasmic organellar genomes in plants is often resolved by genetic and transcriptional cytonuclear responses. Whether this mechanism also acts in allopolyploid members of other kingdoms is not clear. Additionally, cytonuclear coordination of interleaved allopolyploid cells/individuals within the same population is underexplored. The yeast Saccharomyces pastorianus provides the opportunity to explore cytonuclear coevolution during different growth stages and from novel dimensions. Using S. pastorianus cells from multiple growth stages in the same environment, we show that nuclear mitochondria-targeted genes have undergone both asymmetric gene conversion and growth stage-specific biased expression favoring genes from the mitochondrial genome donor (Saccharomyces eubayanus). Our results suggest that cytonuclear coordination in allopolyploid lager yeast species entails an orchestrated and compensatory genetic and transcriptional evolutionary regulatory shift. The common as well as unique properties of cytonuclear coordination underlying allopolyploidy between unicellular yeasts and higher plants offers novel insights into mechanisms of cytonuclear evolution associated with allopolyploid speciation.

Apricot kernels' extract and amygdalin alter bleomycin-induced Ty1 retrotransposition, mitotic gene conversion in the trp-5 locus and reverse point mutations in ilv1-92 allele in Saccharomyces cerevisiae.

The present study aims to analyze the effect of apricot kernels' extract (AKE) and amygdalin (AMY) on bleomycin-induced genetic alternations. Five endpoints were analyzed: cell survival, Ty1 retrotransposition, mitotic gene conversion in the trp-5 locus, reverse point mutations in ilv1-92 allele, and mitotic crossing-over in the ade2 locus. The present work provides the first experimental evidence that bleomycin induces Ty1 retrotransposition in Saccharomyces cerevisiae. New data is obtained that the degree of DNA protection of AMY and AKE depends on the studied genetic event. AKE has been found to provide significant protection against bleomycin-induced Ty1 retrotransposition due to better-expressed antioxidant potential. On the other side, AMY better-expressed protection against bleomycin-induced mitotic gene conversion and reverse mutations may be attributed to the activation of the repair enzymes.

DNA nicks induce mutational signatures associated with BRCA1 deficiency.

Analysis of human cancer genome sequences has revealed specific mutational signatures associated with BRCA1-deficient tumors, but the underlying mechanisms remain poorly understood. Here, we show that one-ended DNA double strand breaks (DSBs) converted from CRISPR/Cas9-induced nicks by DNA replication, not two-ended DSBs, cause more characteristic chromosomal aberrations and micronuclei in Brca1-deficient cells than in wild-type cells. BRCA1 is required for efficient homologous recombination of these nick-converted DSBs and suppresses bias towards long tract gene conversion and tandem duplication (TD) mediated by two-round strand invasion in a replication strand asymmetry. However, aberrant repair of these nick-converted one-ended DSBs, not that of two-ended DSBs in Brca1-deficient cells, generates mutational signatures such as small indels with microhomology (MH) at the junctions, translocations and small MH-mediated TDs, resembling those in BRCA1-deficient tumors. These results suggest a major contribution of DNA nicks to mutational signatures associated with BRCA1 deficiency in cancer and the underlying mechanisms.

Stop Codon Usage as a Window into Genome Evolution: Mutation, Selection, Biased Gene Conversion and the TAG Paradox.

Protein coding genes terminate with one of three stop codons (TAA, TGA, or TAG) that, like synonymous codons, are not employed equally. With TGA and TAG having identical nucleotide content, analysis of their differential usage provides an unusual window into the forces operating on what are ostensibly functionally identical residues. Across genomes and between isochores within the human genome, TGA usage increases with G + C content but, with a common G + C → A + T mutation bias, this cannot be explained by mutation bias-drift equilibrium. Increased usage of TGA in G + C-rich genomes or genomic regions is also unlikely to reflect selection for the optimal stop codon, as TAA appears to be universally optimal, probably because it has the lowest read-through rate. Despite TAA being favored by selection and mutation bias, as with codon usage bias G + C pressure is the prime determinant of between-species TGA usage trends. In species with strong G + C-biased gene conversion (gBGC), such as mammals and birds, the high usage and conservation of TGA is best explained by an A + T → G + C repair bias. How to explain TGA enrichment in other G + C-rich genomes is less clear. Enigmatically, across bacterial and archaeal species and between human isochores TAG usage is mostly unresponsive to G + C pressure. This unresponsiveness we dub the TAG paradox as currently no mutational, selective, or gBGC model provides a well-supported explanation. That TAG does increase with G + C usage across eukaryotes makes the usage elsewhere yet more enigmatic. We suggest resolution of the TAG paradox may provide insights into either an unknown but common selective preference (probably at the DNA/RNA level) or an unrecognized complexity to the action of gBGC.

Estimating the rates of crossover and gene conversion from individual genomes.

Recombination can occur either as a result of crossover or gene conversion events. Population genetic methods for inferring the rate of recombination from patterns of linkage disequilibrium generally assume a simple model of recombination that only involves crossover events and ignore gene conversion. However, distinguishing the 2 processes is not only necessary for a complete description of recombination, but also essential for understanding the evolutionary consequences of inversions and other genomic partitions in which crossover (but not gene conversion) is reduced. We present heRho, a simple composite likelihood scheme for coestimating the rate of crossover and gene conversion from individual diploid genomes. The method is based on analytic results for the distance-dependent probability of heterozygous and homozygous states at 2 loci. We apply heRho to simulations and data from the house mouse Mus musculus castaneus, a well-studied model. Our analyses show (1) that the rates of crossover and gene conversion can be accurately coestimated at the level of individual chromosomes and (2) that previous estimates of the population scaled rate of recombination ρ=4Ner under a pure crossover model are likely biased.

Gene Conversion Explains Elevated Diversity in the Immunity Modulating APL1 Gene of the Malaria Vector Anopheles funestus.

Leucine-rich repeat proteins and antimicrobial peptides are the key components of the innate immune response to Plasmodium and other microbial pathogens in Anopheles mosquitoes. The APL1 gene of the malaria vector Anopheles funestus has exceptional levels of non-synonymous polymorphism across the range of An. funestus, with an average πn of 0.027 versus a genome-wide average of 0.002, and πn is consistently high in populations across Africa. Elevated APL1 diversity was consistent between the independent pooled-template and target-enrichment datasets, however no link between APL1 diversity and insecticide resistance was observed. Although lacking the diversity of APL1, two further mosquito innate-immunity genes of the gambicin anti-microbial peptide family had πn/πs ratios greater than one, possibly driven by either positive or balancing selection. The cecropin antimicrobial peptides were expressed much more highly than other anti-microbial peptide genes, a result discordant with current models of anti-microbial peptide activity. The observed APL1 diversity likely results from gene conversion between paralogues, as evidenced by shared polymorphisms, overlapping read mappings, and recombination events among paralogues. In conclusion, we hypothesize that higher gene expression of APL1 than its paralogues is correlated with a more open chromatin formation, which enhances gene conversion and elevated diversity at this locus.

Frequent somatic gene conversion as a mechanism for loss of heterozygosity in tumor suppressor genes.

The major processes in carcinogenesis include the inactivation of tumor-suppressor genes (TSGs). Although Knudson's two-hit model requires two independent inactivating mutations, perhaps more frequently, a TSG inactivation can occur through a loss of heterozygosity (LOH) of an inactivating mutation. Deletion and uniparental disomy (UPD) have been well documented as LOH mechanisms, but the role of gene conversion is poorly understood. Here, we developed a simple algorithm to detect somatic gene conversion from short-read sequencing data. We applied it to 6285 cancer patient samples, from which 4978 somatic mutations that underwent gene conversion to achieve LOH were found. This number accounted for 14.8% of the total LOH mutations. We further showed that LOH by gene conversion was enriched in TSGs compared with non-TSG genes, showing a significant contribution of gene conversion to carcinogenesis.

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.

Somatic Diversification of Rearranged Antibody Gene Segments by Intra- and Interchromosomal Templated Mutagenesis.

The ability of the humoral immune system to generate Abs capable of specifically binding a myriad of Ags is critically dependent on the somatic hypermutation program. This program induces both templated mutations (i.e., gene conversion) and untemplated mutations. In humans, somatic hypermutation is widely believed to result in untemplated point mutations. In this study, we demonstrate detection of large-scale templated events that occur in human memory B cells and circulating plasmablasts. We find that such mutations are templated intrachromosomally from IGHV genes and interchromosomally from IGHV pseudogenes as well as other homologous regions unrelated to IGHV genes. These same donor regions are used in multiple individuals, and they predominantly originate from chromosomes 14, 15, and 16. In addition, we find that exogenous sequences placed at the IgH locus, such as LAIR1, undergo templated mutagenesis and that homology appears to be the major determinant for donor choice. Furthermore, we find that donor tracts originate from areas in proximity with open chromatin, which are transcriptionally active, and are found in spatial proximity with the IgH locus during the germinal center reaction. These donor sequences are inserted into the Ig gene segment in association with overlapping activation-induced cytidine deaminase hotspots. Taken together, these studies suggest that diversity generated during the germinal center response is driven by untemplated point mutations as well as templated mutagenesis using local and distant regions of the genome.

Gene conversion: a non-Mendelian process integral to meiotic recombination.

Meiosis is undoubtedly the mechanism that underpins Mendelian genetics. Meiosis is a specialised, reductional cell division which generates haploid gametes (reproductive cells) carrying a single chromosome complement from diploid progenitor cells harbouring two chromosome sets. Through this process, the hereditary material is shuffled and distributed into haploid gametes such that upon fertilisation, when two haploid gametes fuse, diploidy is restored in the zygote. During meiosis the transient physical connection of two homologous chromosomes (one originally inherited from each parent) each consisting of two sister chromatids and their subsequent segregation into four meiotic products (gametes), is what enables genetic marker assortment forming the core of Mendelian laws. The initiating events of meiotic recombination are DNA double-strand breaks (DSBs) which need to be repaired in a certain way to enable the homologous chromosomes to find each other. This is achieved by DSB ends searching for homologous repair templates and invading them. Ultimately, the repair of meiotic DSBs by homologous recombination physically connects homologous chromosomes through crossovers. These physical connections provided by crossovers enable faithful chromosome segregation. That being said, the DSB repair mechanism integral to meiotic recombination also produces genetic transmission distortions which manifest as postmeiotic segregation events and gene conversions. These processes are non-reciprocal genetic exchanges and thus non-Mendelian.

Globally Relaxed Selection and Local Adaptation in Boechera stricta.

The strength of selection varies among populations and across the genome, but the determinants of efficacy of selection remain unclear. In this study, we used whole-genome sequencing data from 467 Boechera stricta accessions to quantify the strength of selection and characterize the pattern of local adaptation. We found low genetic diversity on 0-fold degenerate sites and conserved non-coding sites, indicating functional constraints on these regions. The estimated distribution of fitness effects and the proportion of fixed substitutions suggest relaxed negative and positive selection in B. stricta. Among the four population groups, the NOR and WES groups have smaller effective population size (Ne), higher proportions of effectively neutral sites, and lower rates of adaptive evolution compared with UTA and COL groups, reflecting the effect of Ne on the efficacy of natural selection. We also found weaker selection on GC-biased sites compared with GC-conservative (unbiased) sites, suggested that GC-biased gene conversion has affected the strength of selection in B. stricta. We found mixed evidence for the role of the recombination rate on the efficacy of selection. The positive and negative selection was stronger in high-recombination regions compared with low-recombination regions in COL but not in other groups. By scanning the genome, we found different subsets of selected genes suggesting differential adaptation among B. stricta groups. These results show that differences in effective population size, nucleotide composition, and recombination rate are important determinants of the efficacy of selection. This study enriches our understanding of the roles of natural selection and local adaptation in shaping genomic variation.

Transcription Dependent Loss of an Ectopically Expressed Variant Surface Glycoprotein during Antigenic Variation in Trypanosoma brucei.

In the mammalian host, Trypanosoma brucei is coated in a single-variant surface glycoprotein (VSG) species. Stochastic switching of the expressed VSG allows the parasite to escape detection by the host immune system. DNA double-strand breaks (DSB) trigger VSG switching, and repair via gene conversion results in an antigenically distinct VSG being expressed from the single active bloodstream-form expression site (BES). The single active BES is marked by VSG exclusion 2 (VEX2) protein. Here, we have disrupted monoallelic VSG expression by stably expressing a second telomeric VSG from a ribosomal locus. We found that cells expressing two VSGs contained one VEX2 focus that was significantly larger in size than the wild-type cells; this therefore suggests the ectopic VSG is expressed from the same nuclear position as the active BES. Unexpectedly, we report that in the double VSG-expressing cells, the DNA sequence of the ectopic copy is lost following a DSB in the active BES, despite it being spatially separated in the genome. The loss of the ectopic VSG is dependent on active transcription and does not disrupt the number or variety of templates used to repair a BES DSB and elicit a VSG switch. We propose that there are stringent mechanisms within the cell to reinforce monoallelic expression during antigenic variation. IMPORTANCE The single-cell parasite Trypanosoma brucei causes the fatal disease human African trypanosomiasis and is able to colonize the blood, fat, skin, and central nervous system. Trypanosomes survive in the mammalian host owing to a dense protective protein coat that consists of a single-variant surface glycoprotein species. Stochastic switching of one VSG for an immunologically distinct one enables the parasite to escape recognition by the host immune system. We have disrupted monoallelic antigen expression by expressing a second VSG and report that following DSB-triggered VSG switching, the DNA sequence of the ectopic VSG is lost in a transcription-dependent manner. We propose that there are strict requirements to ensure that only one variant antigen is expressed following a VSG switch, which has important implications for understanding how the parasite survives in the mammalian host.

Evolution of CCR5 and CCR2 Genes in Bats Showed Multiple Independent Gene Conversion Events.

Chemokine receptors are an important determinant for the infectiousness of different pathogens, which are able to target the host cells by binding to the extracellular domains of these proteins. This is the mechanism of infection of HIV-1, among other concerning human diseases. Over the past years, it has been shown that two chemokine receptors, CCR2 and CCR5, have been shaped by events of gene conversion in different mammalian lineages, which has been linked to a possible selective advantage against pathogens. Here, by taking advantage of available bat genomes, we present the first insight of CCR2 and CCR5 evolution within the Chiroptera order. In total, four independent events of recombination between CCR2 and CCR5 were detected: two in a single species, Miniopterus natalensis; one in two species from the Rhinolophoidea superfamily; and one in four species from the Pteropodidae family. The regions affected by the gene conversions were generally extensive and always encompassed extracellular domains. Overall, we demonstrate that CCR2 and CCR5 have been subject to extensive gene conversion in multiple species of bats. Considering that bats are known to be large reservoirs of virus in nature, these results might indicate that chimeric CCR2-CCR5 genes might grant some bat species a selective advantage against viruses that rely in the extracellular portions of either CCR2 or CCR5 as gateways into the cell.

Factors That Affect the Rates of Adaptive and Nonadaptive Evolution at the Gene Level in Humans and Chimpanzees.

The rate of amino acid substitution has been shown to be correlated to a number of factors including the rate of recombination, the age of the gene, the length of the protein, mean expression level, and gene function. However, the extent to which these correlations are due to adaptive and nonadaptive evolution has not been studied in detail, at least not in hominids. We find that the rate of adaptive evolution is significantly positively correlated to the rate of recombination, protein length and gene expression level, and negatively correlated to gene age. These correlations remain significant when each factor is controlled for in turn, except when controlling for expression in an analysis of protein length; and they also generally remain significant when biased gene conversion is taken into account. However, the positive correlations could be an artifact of population size contraction. We also find that the rate of nonadaptive evolution is negatively correlated to each factor, and all these correlations survive controlling for each other and biased gene conversion. Finally, we examine the effect of gene function on rates of adaptive and nonadaptive evolution; we confirm that virus-interacting proteins (VIPs) have higher rates of adaptive and lower rates of nonadaptive evolution, but we also demonstrate that there is significant variation in the rate of adaptive and nonadaptive evolution between GO categories when removing VIPs. We estimate that the VIP/non-VIP axis explains about 5-8 fold more of the variance in evolutionary rate than GO categories.

Charge Carrier Transport Mechanism in Ta2 O5 , TaON, and Ta3 N5 Studied by Applying Polaron Hopping and Bandlike Models.

TaON and Ta3 N5 are considered promising materials for photocatalytic and photoelectrochemical water splitting. In contrast, their counterpart Ta2 O5 does not exhibit good photocatalytic performance. This may be explained with the different charge carrier transport mechanisms in these materials, which are not well understood yet. Herein, we investigate the charge transport properties in Ta2 O5 , TaON, and Ta3 N5 by polaron hopping and bandlike models. First, the polaron binding energies were calculated to evaluate whether the small polaron occurs in these materials. Then we performed calculations to localize the excess carriers as small polarons using a hybrid density functional. We find that the small polaron hopping is the charge transfer mechanism in Ta2 O5, whereas our calculations indicate that this mechanism may not occur in TaON and Ta3 N5 . We also investigated the bandlike model mechanism by calculating the charge carrier mobility of these materials using the effective mass approximation, but the calculated mobility is not consistent with experimental results. This study is a first step towards understanding charge transport in oxynitrides and nitrides and furthermore establishes a simple rule to determine whether a small polaron occurs in a material.

Skim-Based Genotyping by Sequencing Using a Double Haploid Population to Call SNPs, Infer Gene Conversions, and Improve Genome Assemblies.

Genotyping by sequencing (GBS) is an emerging technology to rapidly call an abundance of single nucleotide polymorphisms (SNPs) using genome sequencing technology. Several different methodologies and approaches have recently been established, most of these relying on a specific preparation of data. Here we describe our GBS pipeline, which uses high coverage reads from two parents and low coverage reads from their double haploid offspring to call SNPs on a large scale. The upside of this approach is the high resolution and scalability of the method.

GC content of plant genes is linked to past gene duplications.

The frequency of G and C nucleotides in genomes varies from species to species, and sometimes even between different genes in the same genome. The monocot grasses have a bimodal distribution of genic GC content absent in dicots. We categorized plant genes from 5 dicots and 4 monocot grasses by synteny to related species and determined that syntenic genes have significantly higher GC content than non-syntenic genes at their 5`-end in the third position within codons for all 9 species. Lower GC content is correlated with gene duplication, as lack of synteny to distantly related genomes is associated with past interspersed gene duplications. Two mutation types can account for biased GC content, mutation of methylated C to T and gene conversion from A to G. Gene conversion involves non-reciprocal exchanges between homologous alleles and is not detectable when the alleles are identical or heterozygous for presence-absence variation, both likely situations for genes duplicated to new loci. Gene duplication can cause production of siRNA which can induce targeted methylation, elevating mC→T mutations. Recently duplicated plant genes are more frequently methylated and less likely to undergo gene conversion, each of these factors synergistically creating a mutational environment favoring AT nucleotides. The syntenic genes with high GC content in the grasses compose a subset that have undergone few duplications, or for which duplicate copies were purged by selection. We propose a "biased gene duplication / biased mutation" (BDBM) model that may explain the origin and trajectory of the observed link between duplication and genic GC bias. The BDBM model is supported by empirical data based on joint analyses of 9 angiosperm species with their genes categorized by duplication status, GC content, methylation levels and functional classes.