Evolutionary and biomedical insights from a marmoset diploid genome assembly

The accurate and complete assembly of both haplotype sequences of a diploid organism is essential to understanding the role of variation in genome functions, phenotypes, and diseases1. Here, using a trio-binning approach, we present a high-quality, diploid reference genome, with both haplotypes assembled independently at the chromosome level, for the common marmoset (Callithrix jacchus), an important primate model system widely used in biomedical research2,3. The full heterozygosity spectrum between the two haplotypes involves 1.36% of the genome, much higher than the 0.13% indicated by the standard single nucleotide heterozygosity estimation alone. The de novo mutation rate is 0.43 × 10-8 per site per generation, where the paternal inherited genome acquired twice as many mutations as the maternal. Our diploid assembly enabled us to discover a recent expansion of the sex differentiated region and unique evolutionary changes in the marmoset Y chromosome. Additionally, we identified many genes with signatures of positive selection that might have contributed to the evolution of Callithrix biological features. Brain related genes were highly conserved between marmosets and humans, though several genes experienced lineage-specific copy number variations or diversifying selection, providing important implications for the application of marmosets as a model system.

Fast effect size shrinkage software for beta-binomial models of allelic imbalance

Allelic imbalance occurs when the two alleles of a gene are differentially expressed within a diploid organism and can indicate important differences in cis-regulation and epigenetic state across the two chromosomes. Because of this, the ability to accurately quantify the proportion at which each allele of a gene is expressed is of great interest to researchers. This becomes challenging in the presence of small read counts and/or sample sizes, which can cause estimators for allelic expression proportions to have high variance. Investigators have traditionally dealt with this problem by filtering out genes with small counts and samples. However, this may inadvertently remove important genes that have truly large allelic imbalances. Another option is to use pseudocounts or Bayesian estimators to reduce the variance. To this end, we evaluated the accuracy of four different estimators, the latter two of which are Bayesian shrinkage estimators: maximum likelihood, adding a pseudocount to each allele, approximate posterior estimation of GLM coefficients (apeglm) and adaptive shrinkage (ash). We also wrote C++ code to quickly calculate ML and apeglm estimates and integrated it into the apeglm package. The four methods were evaluated on two simulations and one real data set. Apeglm consistently performed better than ML according to a variety of criteria, and generally outperformed use of pseudocounts as well. Ash also performed better than ML in one of the simulations, but in the other performance was more mixed. Finally, when compared to five other packages that also fit beta-binomial models, the apeglm package was substantially faster and more numerically reliable, making our package useful for quick and reliable analyses of allelic imbalance. Apeglm is available as an R/Bioconductor package at http://bioconductor.org/packages/apeglm.

PLK-1 promotes the merger of the parental genome into a single nucleus by triggering lamina disassembly

Life of sexually reproducing organisms starts with the fusion of the haploid egg and sperm gametes to form the genome of a new diploid organism. Using the newly fertilized Caenorhabditis elegans zygote, we show that the mitotic Polo-like kinase PLK-1 phosphorylates the lamin LMN-1 to promote timely lamina disassembly and subsequent merging of the parental genomes into a single nucleus after mitosis. Expression of non-phosphorylatable versions of LMN-1, which affect lamina depolymerization during mitosis, is sufficient to prevent the mixing of the parental chromosomes into a single nucleus in daughter cells. Finally, we recapitulate lamina depolymerization by PLK-1 in vitro demonstrating that LMN-1 is a direct PLK-1 target. Our findings indicate that the timely removal of lamin is essential for the merging of parental chromosomes at the beginning of life in C. elegans and possibly also in humans, where a defect in this process might be fatal for embryo development.

Fast effect size shrinkage software for beta-binomial models of allelic imbalance

Allelic imbalance occurs when the two alleles of a gene are differentially expressed within a diploid organism, and can indicate important differences in cis-regulation and epigenetic state across the two chromosomes. Because of this, the ability to accurately quantify the proportion at which each allele of a gene is expressed is of great interest to researchers. This becomes challenging in the presence of small read counts and/or sample sizes, which can cause estimates for allelic expression proportions to have high variance. Investigators have traditionally dealt with this problem by filtering out genes with small counts and samples. However, this may inadvertently remove important genes that have truly large allelic imbalances. Another option is to use Bayesian estimators to reduce the variance. To this end, we evaluated the accuracy of three different estimators, the latter two of which are Bayesian shrinkage estimators: maximum likelihood, approximate posterior estimation of GLM coefficients (apeglm) and adaptive shrinkage (ash). We also wrote C++ code to quickly calculate ML and apeglm estimates, and integrated it into the apeglm package. The three methods were evaluated on both simulated and real data. Apeglm consistently performed better than ML according to a variety of criteria, including mean absolute error and concordance at the top. While ash had lower error and greater concordance than ML on the simulations, it also had a tendency to over-shrink large effects, and performed worse on the real data according to error and concordance. Furthermore, when compared to five other packages that also fit beta-binomial models, the apeglm package was substantially faster, making our package useful for quick and reliable analyses of allelic imbalance. Apeglm is available as an R/Bioconductor package at http://bioconductor.org/packages/apeglm.

Killer meiotic drive and dynamic evolution of the wtf gene family

Natural selection works best when the two alleles in a diploid organism are transmitted to offspring at equal frequencies. Despite this, selfish loci known as meiotic drivers that bias their own transmission into gametes are found throughout eukaryotes. Drive is thought to be a powerful evolutionary force, but empirical evolutionary analyses of drive systems are limited by low numbers of identified meiotic drive genes. Here, we analyze the evolution of the wtf gene family of Schizosaccharomyces pombe that contains both killer meiotic drive genes and suppressors of drive. We completed assemblies of all wtf genes for two S. pombe strains, as well as a subset of wtf genes from over 50 strains. We find that wtf copy number can vary greatly between strains, and that amino acid substitutions, expansions and contractions of DNA sequence repeats, and nonallelic gene conversion between family members all contribute to dynamic wtf gene evolution. This work demonstrates the power of meiotic drive to foster rapid evolution and identifies a recombination mechanism through which transposons can indirectly mobilize meiotic drivers.

Cross-species fertilization: the hamster egg receptor, Juno, binds the human sperm ligand, Izumo1

Fertilization is the culminating event in sexual reproduction and requires the recognition and fusion of the haploid sperm and egg to form a new diploid organism. Specificity in these recognition events is one reason why sperm and eggs from different species are not normally compatible. One notable exception is the unusual ability of zona-free eggs from the Syrian golden hamster ( Mesocricetus auratus ) to recognize and fuse with human sperm, a phenomenon that has been exploited to assess sperm quality in assisted fertility treatments. Following our recent finding that the interaction between the sperm and egg recognition receptors Izumo1 and Juno is essential for fertilization, we now demonstrate concordance between the ability of Izumo1 and Juno from different species to interact, and the ability of their isolated gametes to cross-fertilize each other in vitro . In particular, we show that Juno from the golden hamster can directly interact with human Izumo1. These data suggest that the interaction between Izumo1 and Juno plays an important role in cross-species gamete recognition, and may inform the development of improved prognostic tests that do not require the use of animals to guide the most appropriate fertility treatment for infertile couples.

Sperm competition and sperm cooperation: the potential role of diploid and haploid expression

Sperm competition is a powerful selective force driving the evolution of sperm shape and function. Recent findings suggest that sperm cooperation is a potential evolutionary response to sperm competition. Sperm cooperation may enhance the performance of the ejaculate increasing a male's chance to outcompete rival males in competition for fertilisation. Whether and how sperm cooperation may evolve is the focal point of this review. The relative importance of haploid and diploid gene expression for the evolution of sperm cooperation and the potential conflict of interest between (i) haploid sperm and diploid male and (ii) among sibling sperm, since sibling sperm only share an average of 50% of their genes in a diploid organism, are discussed. Furthermore, sperm cooperation is defined and the literature for empirical evidence of sperm cooperation is reviewed in light of the author's definitions.

A Telomere-mediated Chromosome Fragmentation Approach to Assess Mitotic Stability and Ploidy Alterations ofLeishmaniaChromosomes

We have used a telomere-associated chromosome fragmentation strategy to induce internal chromosome-specific breakage ofLeishmaniachromosomes. The integration of telomeric repeats from the kinetoplastidTrypanosoma bruceiinto defined positions of theLeishmaniagenome by homologous recombination can induce chromosome breakage accompanied by the deletion of the chromosomal part that is distal to the site of the break. The cloned telomeric DNA at the end of the truncated chromosomes is functional and it can seed the formation of new telomeric repeats. We found that genome ploidy is often altered upon telomere-mediated chromosome fragmentation events resulting in large chromosomal deletions. In most cases diploidy is either preserved, or partial trisomic cells are observed, but interestingly we report here the generation of partial haploid mutants in this diploid organism. Partial haploidLeishmaniamutants should facilitate studies on the function of chromosome-assigned genes. We also present several lines of evidence for the presence of sequences involved in chromosome mitotic stability and segregation during cell cycle in this parasitic protozoan. Telomere-directed chromosome fragmentation studies inLeishmaniamay constitute a useful tool to assay for centromere function.

Recombination and the evolution of diploidy.

With two copies of every gene, a diploid organism is able to mask recessive deleterious mutations. In this paper we present the analysis of a two-locus model designed to determine when the masking of deleterious alleles favors the evolution of a dominant diploid phase in organisms that alternate between haploid and diploid phases ("alternation of generations"). It is hypothesized that diploidy will be favored whenever masking occurs ("the masking hypothesis"). Using analytical methods, we confirm that this masking hypothesis is essentially correct under free recombination: as long as the heterozygous expression of deleterious alleles is sufficiently masked by the wild-type allele, diploidy is favored over haploidy. When the rate of recombination is lower, however, diploidy is much less likely to be favored over haploidy. In fact, according to our model, the evolution of diploidy is impossible without significant levels of recombination even when masking is fairly strong.