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.

CRISPR/Cas9-mediated gene knockout and interallelic gene conversion in human induced pluripotent stem cells using non-integrative bacteriophage-chimeric retrovirus-like particles

Background The application of CRISPR/Cas9 technology in human induced pluripotent stem cells (hiPSC) holds tremendous potential for basic research and cell-based gene therapy. However, the fulfillment of these promises relies on the capacity to efficiently deliver exogenous nucleic acids and harness the repair mechanisms induced by the nuclease activity in order to knock-out or repair targeted genes. Moreover, transient delivery should be preferred to avoid persistent nuclease activity and to decrease the risk of off-target events. We recently developed bacteriophage-chimeric retrovirus-like particles that exploit the properties of bacteriophage coat proteins to package exogenous RNA, and the benefits of lentiviral transduction to achieve highly efficient, non-integrative RNA delivery in human cells. Here, we investigated the potential of bacteriophage-chimeric retrovirus-like particles for the non-integrative delivery of RNA molecules in hiPSC for CRISPR/Cas9 applications. Results We found that these particles efficiently convey RNA molecules for transient expression in hiPSC, with minimal toxicity and without affecting the cell pluripotency and subsequent differentiation. We then used this system to transiently deliver in a single step the CRISPR-Cas9 components (Cas9 mRNA and sgRNA) to generate gene knockout with high indel rate (up to 85%) at multiple loci. Strikingly, when using an allele-specific sgRNA at a locus harboring compound heterozygous mutations, the targeted allele was not altered by NHEJ/MMEJ, but was repaired at high frequency using the homologous wild type allele, i.e., by interallelic gene conversion. Conclusions Our results highlight the potential of bacteriophage-chimeric retrovirus-like particles to efficiently and safely deliver RNA molecules in hiPSC, and describe for the first time genome engineering by gene conversion in hiPSC. Harnessing this DNA repair mechanism could facilitate the therapeutic correction of human genetic disorders in hiPSC.

Linked-Read Sequencing of Eight Falcons Reveals a Unique Genomic Architecture in Flux

Falcons are diverse birds of cultural and economic importance. They have undergone major lineage specific chromosomal rearrangements, resulting in greatly-reduced chromosome counts relative to other birds. Here, we use 10X Genomics linked reads to provide new high-contiguity genomes for two gyrfalcons, a saker falcon, a lanner falcon, three subspecies of peregrine falcons, and the common kestrel. Assisted by a transcriptome sequenced from 22 gyrfalcon tissues, we annotate these genomes for a variety of genomic features, estimate historical demography, and then investigate genomic equilibrium in the context of falcon-specific chromosomal rearrangements. We find that falcon genomes are not in AT-GC equilibrium with a bias in mutations towards higher AT content; this bias is predominantly driven by, but not dependent on, hypermutability of CpG sites. Small indels and large structural variants were also biased towards insertions rather than deletions. Patterns of disequilibrium were linked to chromosomal rearrangements: falcons have lost GC content in regions that have fused to larger chromosomes from microchromosomes and gained GC content in regions of macrochromosomes that have translocated to microchromosomes. Inserted bases have accumulated on regions ancestrally belonging to microchromosomes, consistent with insertion-biased gene conversion. We also find an excess of interspersed repeats on regions of microchromosomes that have fused to macrochromosomes. Our results reveal that falcon genomes are in a state of flux. They further suggest that many of the key differences between microchromosomes and macrochromosomes are driven by differences in chromosome size, and indicate a clear role for recombination and biased gene conversion in determining genomic equilibrium.

Meiotic Cas9 expression mediates gene conversion in the male and female mouse germline

Highly efficient gene conversion systems have the potential to facilitate the study of complex genetic traits using laboratory mice and, if implemented as a “gene drive,” to limit loss of biodiversity and disease transmission caused by wild rodent populations. We previously showed that such a system of gene conversion from heterozygous to homozygous after a sequence targeted CRISPR/Cas9 double-strand DNA break (DSB) is feasible in the female mouse germline. In the male germline, however, all DSBs were instead repaired by end joining (EJ) mechanisms to form an “insertion/deletion” (indel) mutation. These observations suggested that timing Cas9 expression to coincide with meiosis I is critical to favor conditions when homologous chromosomes are aligned and interchromosomal homology-directed repair (HDR) mechanisms predominate. Here, using a Cas9 knock-in allele at the Spo11 locus, we show that meiotic expression of Cas9 does indeed mediate gene conversion in the male as well as in the female germline. However, the low frequency of both HDR and indel mutation in both male and female germlines suggests that Cas9 may be expressed from the Spo11 locus at levels too low for efficient DSB formation. We suggest that more robust Cas9 expression initiated during early meiosis I may improve the efficiency of gene conversion and further increase the rate of “super-mendelian” inheritance from both male and female mice.

A Role for Synaptonemal Complex in Meiotic Mismatch Repair

A large subset of meiotic recombination intermediates form within the physical context of synaptonemal complex (SC), but the functional relationship between SC structure and homologous recombination remains obscure. Our prior analysis of strains deficient for SC central element proteins demonstrated that tripartite SC is dispensable for interhomolog recombination in S. cerevisiae. Here we report that while dispensable for recombination per se, SC proteins promote efficient mismatch repair at interhomolog recombination sites. Failure to repair mismatches within heteroduplex-containing meiotic recombination intermediates leads to genotypically sectored colonies (post meiotic segregation events). We discovered increased post-meiotic segregation at THR1 in cells lacking Ecm11 or Gmc2, or in the SC-deficient but recombination-proficient zip1[Δ21-163] mutant. High-throughput sequencing of octad meiotic products furthermore revealed a genome-wide increase in recombination events with unrepaired mismatches in ecm11 mutants relative to wild type. Meiotic cells missing Ecm11 display longer gene conversion tracts, but tract length alone does not account for the higher frequency of unrepaired mismatches. Interestingly, the per-nucleotide mismatch frequency is elevated in ecm11 when analyzing all gene conversion tracts, but is similar between wild type and ecm11 if considering only those events with unrepaired mismatches. Thus, in both wild type and ecm11 strains a subset of recombination events is susceptible to a similar degree of inefficient mismatch repair, but in ecm11 mutants a larger fraction of events fall into this inefficient repair category. Finally, we observe elevated post-meiotic segregation at THR1 in mutants with a dual deficiency in MutSγ crossover recombination and SC assembly, but not in the mlh3 mutant, which lacks MutSγ crossovers but has abundant SC. We propose that SC structure promotes efficient mismatch repair of joint molecule recombination intermediates, and that absence of SC is the molecular basis for elevated post-meiotic segregation in both MutSγ crossover-proficient (ecm11, gmc2) and MutSγ crossover-deficient (msh4, zip3) strains.

Illegitimate Recombination between Duplicated Genes Generated from Recursive Polyploidizations Accelerated the Divergence of the Genus Arachis

The peanut (Arachis hypogaea L.) is the leading oil and food crop among the legume family. Extensive duplicate gene pairs generated from recursive polyploidizations with high sequence similarity could result from gene conversion, caused by illegitimate DNA recombination. Here, through synteny-based comparisons of two diploid and three tetraploid peanut genomes, we identified the duplicated genes generated from legume common tetraploidy (LCT) and peanut recent allo-tetraploidy (PRT) within genomes. In each peanut genome (or subgenomes), we inferred that 6.8–13.1% of LCT-related and 11.3–16.5% of PRT-related duplicates were affected by gene conversion, in which the LCT-related duplicates were the most affected by partial gene conversion, whereas the PRT-related duplicates were the most affected by whole gene conversion. Notably, we observed the conversion between duplicates as the long-lasting contribution of polyploidizations accelerated the divergence of different Arachis genomes. Moreover, we found that the converted duplicates are unevenly distributed across the chromosomes and are more often near the ends of the chromosomes in each genome. We also confirmed that well-preserved homoeologous chromosome regions may facilitate duplicates’ conversion. In addition, we found that these biological functions contain a higher number of preferentially converted genes, such as catalytic activity-related genes. We identified specific domains that are involved in converted genes, implying that conversions are associated with important traits of peanut growth and development.

Gene conversion explains elevated diversity in the immunity modulating APL1 gene of the malaria vector Anopheles funestus

The leucine rich repeat gene APL1 is a key component of immunity to Plasmodium and other microbial pathogens in Anopheles mosquitoes. In the malaria vector Anopheles funestus the APL1 gene has four paralogues which occur along the same chromosome arm. We show that APL1 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 (and πs) is consistently high in populations across Africa. The pattern of APL1 diversity was consistent between independent pooled-template and target-enrichment datasets, however no link between APL1 diversity and insecticide-resistance was observed with the phenotyped target-enrichment dataset. Two further 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. Cecropin antimicrobial peptides were expressed much more highly than other anti-microbial peptide genes, an observation discordant with current models of anti-microbial peptide activity. The observed APL1 diversity likely results from gene conversion between paralogs, as evidenced by shared polymorphisms, overlapping read mappings, and recombination events among paralogues. Gene conversion at APL1 versus alternative explanations is concordant with similarly elevated diversity in APL1 and TEP1 loci in An. gambiae. In contrast, the more closely related An. stephensi which also encodes a single-copy of APL1 does not show this elevated diversity. We hypothesise that a more open chromatin formation at the APL1 locus due to higher gene expression than its paralogues enhances gene conversion, and therefore increased polymorphism, at APL1.

Muller’s ratchet of the Y chromosome with gene conversion

Muller’s ratchet is a process in which deleterious mutations are fixed irreversibly in the absence of recombination. The degeneration of the Y chromosome, and the gradual loss of its genes, can be explained by Muller’s ratchet. However, most theories consider single-copy genes, and may not be applicable to Y chromosomes, which have a number of duplicated genes in many species, which are probably undergoing concerted evolution by gene conversion. We developed a model of Muller’s ratchet to explore the evolution of the Y chromosome. The model assumes a non-recombining chromosome with both single-copy and duplicated genes. We used analytical and simulation approaches to obtain the rate of gene loss in this model, with special attention to the role of gene conversion. Homogenization by gene conversion makes both duplicated copies either mutated or intact. The former promotes the ratchet, and the latter retards, and we ask which of these counteracting forces dominates under which conditions. We found that the effect of gene conversion is complex, and depends upon the fitness effect of gene duplication. When duplication has no effect on fitness, gene conversion accelerates the ratchet of both single-copy and duplicated genes. If duplication has an additive fitness effect, the ratchet of single-copy genes is accelerated by gene duplication, regardless of the gene conversion rate, whereas gene conversion slows the degeneration of duplicated genes. Our results suggest that the evolution of the Y chromosome involves several parameters, including the fitness effect of gene duplication by increasing dosage and gene conversion rate.