Genome-wide single-cell analysis of recombination activity and de novo mutation rates in human sperm

Next-Generation Genotoxicology: Using Modern Sequencing Technologies to Assess Somatic Mutagenesis and Cancer Risk

Mutations have a profound effect on human health, particularly through an increased risk of carcinogenesis and genetic disease. The strong correlation between mutagenesis and carcinogenesis has been a driving force behind genotoxicity research for more than 50 years. The stochastic and infrequent nature of mutagenesis makes it challenging to observe and to study. Indeed, decades have been spent developing increasingly sophisticated assays and methods to study these low-frequency genetic errors, in hopes of better predicting which chemicals may be carcinogens, understanding their mode of action, and informing guidelines to prevent undue human exposure. While effective, widely used genetic selection-based technologies have a number of limitations that have hampered major advancements in the field of genotoxicity. Emerging new tools, in the form of enhanced next-generation sequencing platforms and methods, are changing this paradigm. In this review, we discuss rapidly evolving sequencing tools and technologies, such as error-corrected sequencing and single cell analysis, which we anticipate will fundamentally reshape the field. In addition, we consider a variety emerging applications for these new technologies, including the detection of DNA adducts, inference of mutational processes based on genomic site and local sequence contexts, and evaluation of genome engineering fidelity, as well as other cutting-edge challenges for the next 50 years of environmental and molecular mutagenesis research. Environ. Mol. Mutagen. 61:135-151, 2020. © 2019 The Authors. Environmental and Molecular Mutagenesis published by Wiley Periodicals, Inc. on behalf of Environmental Mutagen Society.

Recent advances in mammalian reproductive biology

Reproductive biology is a uniquely important topic since it is about germ cells, which are central for transmitting genetic information from generation to generation. In this review, we discuss recent advances in mammalian germ cell development, including preimplantation development, fetal germ cell development and postnatal development of oocytes and sperm. We also discuss the etiologies of female and male infertility and describe the emerging technologies for studying reproductive biology such as gene editing and single-cell technologies.

One-step liquid molding based modular microfluidic circuits

We present an easy-to-follow modular method that combines liquid molding with standard SU-8 lithography to create customized integrated microfluidic devices for the changing needs of users.

Aneuploidy and recombination in the human preimplantation embryo. Copy number variation analysis and genome-wide polymorphism genotyping

RESEARCH QUESTION

What are the incidence and patterns of meiotic trisomies and recombination separately and in relation to each other at the blastocyst stage via single nucleotide polymorphism genotyping combined with array comparative genomic hybridization.

DESIGN

Single nucleotide polymorphism microarrays were carried out on a total of 1442 blastocyst stage embryos derived from 268 fertile couples undergoing preimplantation genetic diagnosis for the purposes of avoiding transmittance of known single gene disorders to their offspring; 24-chromosome aneuploidy screening via array comparative genomic hybridization was carried out in parallel.

RESULTS

One hundred per cent of meiotic trisomies identified in these embryos were of maternal origin and their incidence increased significantly with advancing maternal age (P < 0.0001). A total of 55.8% of meiotic trisomies were meiosis I-type and 44.2% were meiosis II-type. Certain chromosomes were affected more by meiosis I-type errors, whereas others experienced more meiosis II-type errors. A detailed recombination analysis was carried out for 11,476 chromosomes and 17,763 recombination events were recorded. The average number of recombination sites was 24.0 ± 0.3 for male meiosis and 41.2 ± 0.6 for female meiosis (autosomes only). Sex-specific differences were observed in the locations of recombination sites. Comparative analysis conducted between 190 euploid embryos and 69 embryos presenting maternal meiotic trisomies showed similar recombination rates (P = 0.425) and non-recombinant chromatid rates (P = 0.435) between the two categories; differences, however, were observed when analysing embryos affected with specific maternal meiotic trisomies.

CONCLUSIONS

This study yielded unique data concerning recombination and the origin of aneuploidies observed during the first few days of life and provides a novel insight into these important biological processes.

Crossover interference and sex-specific genetic maps shape identical by descent sharing in close relatives

Simulations of close relatives and identical by descent (IBD) segments are common in genetic studies, yet most past efforts have utilized sex averaged genetic maps and ignored crossover interference, thus omitting features known to affect the breakpoints of IBD segments. We developed Ped-sim, a method for simulating relatives that can utilize either sex-specific or sex averaged genetic maps and also either a model of crossover interference or the traditional Poisson model for inter-crossover distances. To characterize the impact of previously ignored mechanisms, we simulated data for all four combinations of these factors. We found that modeling crossover interference decreases the standard deviation of pairwise IBD proportions by 10.4% on average in full siblings through second cousins. By contrast, sex-specific maps increase this standard deviation by 4.2% on average, and also impact the number of segments relatives share. Most notably, using sex-specific maps, the number of segments half-siblings share is bimodal; and when combined with interference modeling, the probability that sixth cousins have non-zero IBD sharing ranges from 9.0 to 13.1%, depending on the sexes of the individuals through which they are related. We present new analytical results for the distributions of IBD segments under these models and show they match results from simulations. Finally, we compared IBD sharing rates between simulated and real relatives and find that the combination of sex-specific maps and interference modeling most accurately captures IBD rates in real data. Ped-sim is open source and available from https://github.com/williamslab/ped-sim.

Novel Techniques of Sperm Selection for Improving IVF and ICSI Outcomes

Almost 50% of the infertility cases are due to male factors. Assisted reproductive technologies (ARTs) allow to overcome the incapacity of these patients' spermatozoa to fertilize the oocyte and produce a viable and healthy offspring, but the efficiency of the different techniques has still the potential to improve. According to the latest reports of the European Society of Human Reproduction and Embryology (ESHRE) and the Centers for Disease Control and Prevention of the United States (CDC), the percentages of deliveries per ART cycle in 2014 and 2016 were 21 and 22%, respectively. Among the reasons for this relatively low efficiency, the quality of the spermatozoa has been pointed out as critical, and the presence of high percentages of DNA-damaged spermatozoa in patients' ejaculates is possibly one of the main factors reducing the ARTs outcomes. Thus, one of the main challenges in reproductive medicine is to ensure the highest quality of the spermatozoa used in ARTs, and specifically, in terms of genetic integrity. The latest techniques for the preparation and selection of human spermatozoa are herein discussed focusing on those proven to improve one or several of the following parameters: sperm genetic integrity, fertilization capacity, embryo production, and in vitro survival, as well as pregnancy and delivery rates following in vitro fertilization (IVF) and intracytoplasmic sperm injection (ICSI). In addition, we discuss the potential of techniques developed in non-human mammals that could be further transferred to the clinic.

Copy Number Variations in a Cohort of 420 Individuals with Neurodevelopmental Disorders From the South of Brazil

Chromosomal microarray (CMA) is now recommended as first tier for the evaluation in individuals with unexplained neurodevelopmental disorders (ND). However, in developing countries such as Brazil, classical cytogenetic tests are still the most used in clinical practice, as reflected by the scarcity of publications of microarray investigation in larger cohorts. This is a retrospective study which analyses the reading files of CMA and available clinical data from 420 patients from the south of Brazil, mostly children, with neurodevelopmental disorders requested by medical geneticists and neurologists for diagnostic purpose. Previous karyotyping was reported for 138 and includes 17 with abnormal results. The platforms used for CMA were CYTOSCAN 750K (75%) and CYTOSCAN HD (25%). The sex ratio of the patients was 1.625 males :1 female and the mean age was 9.5 years. A total of 96 pathogenic copy number variations (CNVs), 58 deletions and 38 duplications, were found in 18% of the patients and in all chromosomes, except chromosome 11. For 12% of the patients only variants of uncertain clinical significance were found. No clinically relevant CNV was found in 70%. The main referrals for chromosomal microarrays (CMA) were developmental delay (DD), intellectual disability (ID), facial dysmorphism and autism spectrum disorder (ASD). DD/ID were present in 80%, facial dysmorphism in 52% and ASD in 32%. Some phenotypes in this population could be predictive of a higher probability to carry a pathogenic CNV, as follows: dysmorphic facial features (p-value = < 0.0001, OR = 0.32), obesity (p-value = 0.006, OR = 0.20), short stature (p-value = 0.032, OR = 0.44), genitourinary anomalies (p-value = 0.032, OR = 0.63) and ASD (p-value = 0.039, OR = 1.94). The diagnostic rate for CMA in this study was 18%. We present the largest report of CMA data in a cohort with ND in Brazil. We characterize the rare CNVs found together with the main phenotypes presented by each patient, list phenotypes which could predict a higher diagnostic probability by CMA in patients with a neurodevelopmental disorder and show how CMA and classical karyotyping results are complementary.

Copolymerization of single-cell nucleic acids into balls of acrylamide gel

We show the use of 5′-Acrydite oligonucleotides to copolymerize single-cell DNA or RNA into balls of acrylamide gel (BAGs). Combining this step with split-and-pool techniques for creating barcodes yields a method with advantages in cost and scalability, depth of coverage, ease of operation, minimal cross-contamination, and efficient use of samples. We perform DNA copy number profiling on mixtures of cell lines, nuclei from frozen prostate tumors, and biopsy washes. As applied to RNA, the method has high capture efficiency of transcripts and sufficient consistency to clearly distinguish the expression patterns of cell lines and individual nuclei from neurons dissected from the mouse brain. By using varietal tags (UMIs) to achieve sequence error correction, we show extremely low levels of cross-contamination by tracking source-specific SNVs. The method is readily modifiable, and we will discuss its adaptability and diverse applications.

Haploid Selection in “Diploid” Organisms

Evolutionary rates and strength of selection differ markedly between haploid and diploid genomes. Any genes expressed in a haploid state will be directly exposed to selection, whereas alleles in a diploid state may be partially or fully masked by a homologous allele. This difference may shape key evolutionary processes, including rates of adaptation and inbreeding depression, but also the evolution of sex chromosomes, heterochiasmy, and stable sex ratio biases. All diploid organisms carry haploid genomes, most notably the haploid genomes in gametes produced by every sexually reproducing eukaryote. Furthermore, haploid expression occurs in genes with monoallelic expression, in sex chromosomes, and in organelles, such as mitochondria and plastids. A comparison of evolutionary rates among these haploid genomes reveals striking parallels. Evidence suggests that haploid selection has the potential to shape evolution in predominantly diploid organisms, and taking advantage of the rapidly developing technologies, we are now in the position to quantify the importance of such selection on haploid genomes.

National Cancer Institute Alliance for nanotechnology in cancer-Catalyzing research and translation toward novel cancer diagnostics and therapeutics

Nanotechnology has been a burgeoning research field, which is finding compelling applications in several practical areas of everyday life. It has provided novel, paradigm shifting solutions to medical problems and particularly to cancer. In order to accelerate integration of nanotechnology into cancer research and oncology, the National Cancer Institute (NCI) of the National Institutes of Health (NIH) established the NCI Alliance for Nanotechnology in Cancer program in 2005. This effort brought together scientists representing physical sciences, chemistry, and engineering working at the nanoscale with biologists and clinicians working on cancer to form a uniquely multidisciplinary cancer nanotechnology research community. The last 14 years of the program have produced a remarkable body of scientific discovery and demonstrated its utility to the development of practical cancer interventions. This paper takes stock of how the Alliance program influenced melding of disparate research disciplines into the field of nanomedicine and cancer nanotechnology, has been highly productive in the scientific arena, and produced a mechanism of seamless transfer of novel technologies developed in academia to the clinical and commercial space. This article is categorized under: Toxicology and Regulatory Issues in Nanomedicine > Regulatory and Policy Issues in Nanomedicine Therapeutic Approaches and Drug Discovery > Nanomedicine for Oncologic Disease Diagnostic Tools > in vivo Nanodiagnostics and Imaging.

De Novo Mutations Reflect Development and Aging of the Human Germline

Human germline de novo mutations (DNMs) are both a driver of evolution and an important cause of genetic diseases. In the past few years, whole-genome sequencing (WGS) of parent-offspring trios has facilitated the large-scale detection and study of human DNMs, which has led to exciting discoveries. The overarching theme of all of these studies is that the DNMs of an individual are a complex mixture of mutations that arise through different biological processes acting at different times during human development and life.

DNA variants are an unlikely explanation for the changing quality of spermatozoa within the same individual

It has recently been suggested that the human sperm genome is highly unstable, which may be a reasonable explanation as to why men, even fertile men, produce defective spermatozoa. Furthermore, an unstable genome may also explain why the semen profile of the same man changes from one ejaculate to the next. As such, we took multiple ejaculates (between 3 and 6) from 7 individuals over a 6-month period and isolated sperm through density gradients. We then compared the DNA of: (i) good and poor-quality spermatozoa within the same ejaculate; and (ii) from multiple ejaculates from the same individual. Our results suggest that on a global level, DNA present within spermatozoa is actually quite stable and similar between both good and poor sperm. This is important information for the assisted reproductive community when it comes to sperm selection.

Single-Cell Sequencing Reveals the Relationship between Phenotypes and Genotypes of Klinefelter Syndrome

Klinefelter syndrome (KS) is one of the most common congenital disorders of male infertility. Given its high heterogeneity in clinical and genetic presentation, the relationship between transcriptome, clinical phenotype, and associated co-morbidities seen in KS has not been fully clarified. Here, we report a 47,XXY Chinese male with infertility and analyzed the differences in gene expression patterns of peripheral blood mononuclear cells (PBMCs) with regard to a Chinese male and a female control with normal karyotype by single-cell sequencing. A total of 24,439 cells were analyzed and divided into 5 immune cell types (including B cells, T cells, macrophage cells, dendritic cells, and natural killer cells) according to marker genes. Using unsupervised dimensionality reduction and clustering algorithms, we identified molecularly distinct subpopulations of cells between the KS patient and both controls. Gene ontology enrichment analyses yielded terms associated with well-known comorbidities seen in KS as well as an affected immune system and type I diabetes mellitus. Based on our data, we identified several candidate genes which may be implicated in regulating the phenotype of KS. Overall, this analysis provides a comprehensive map of the cell types of PBMCs in a KS patient at the single-cell level, which will contribute to the prevention of comorbidity and improvement of the life quality of KS patients.

Quantification of Recombination Rate and Segregation Distortion by Genotyping and Sequencing of Single Pollen Nuclei

Meiosis is a specialized cell division during which homologous chromosomes can exchange genetic material through recombination. This mechanism generates novel allelic combinations, which can be exploited by plant breeders to achieve crop improvement. Pollen grains are the haploid products of meiosis required in fertilization. Here, we describe two approaches to measure meiotic recombination in single haploid pollen nuclei. Pollen nuclei are first separated by fluorescence-activated cell-sorting. Afterwards, the DNA of single pollen nuclei can be amplified by multiple-displacement-amplification using Phi29 DNA polymerase and meiotic recombination events can be measured using KASP markers. Alternatively, the PicoPLEX DNA-seq kit can be used to amplify the DNA of single pollen nuclei followed by library preparation for whole-genome sequencing and subsequent bioinformatic analysis.

Genome-wide recombination map construction from single individuals using linked-read sequencing

Meiotic recombination rates vary across the genome, often involving localized crossover hotspots and coldspots. Studying the molecular basis and mechanisms underlying this variation has been challenging due to the high cost and effort required to construct individualized genome-wide maps of recombination crossovers. Here we introduce a new method, called ReMIX, to detect crossovers from gamete DNA of a single individual using Illumina sequencing of 10X Genomics linked-read libraries. ReMIX reconstructs haplotypes and identifies the valuable rare molecules spanning crossover breakpoints, allowing quantification of the genomic location and intensity of meiotic recombination. Using a single mouse and stickleback fish, we demonstrate how ReMIX faithfully recovers recombination hotspots and landscapes that have previously been built using hundreds of offspring. ReMIX provides a high-resolution, high-throughput, and low-cost approach to quantify recombination variation across the genome, providing an exciting opportunity to study recombination among multiple individuals in diverse organisms.

Variation of recombination rates within genomes has important implications in genetics and evolution. Here, the authors develop a method for building genome-wide recombination maps from single individuals using linked-read sequencing data, and report its application in mouse and stickleback fish.

High-Throughput Single-Cell Sequencing with Linear Amplification

Conventional methods for single-cell genome sequencing are limited with respect to uniformity and throughput. Here, we describe sci-L3, a single-cell sequencing method that combines combinatorial indexing (sci-) and linear (L) amplification. The sci-L3 method adopts a 3-level (3) indexing scheme that minimizes amplification biases while enabling exponential gains in throughput. We demonstrate the generalizability of sci-L3 with proof-of-concept demonstrations of single-cell whole-genome sequencing (sci-L3-WGS), targeted sequencing (sci-L3-target-seq), and a co-assay of the genome and transcriptome (sci-L3-RNA/DNA). We apply sci-L3-WGS to profile the genomes of >10,000 sperm and sperm precursors from F1 hybrid mice, mapping 86,786 crossovers and characterizing rare chromosome mis-segregation events in meiosis, including instances of whole-genome equational chromosome segregation. We anticipate that sci-L3 assays can be applied to fully characterize recombination landscapes, to couple CRISPR perturbations and measurements of genome stability, and to other goals requiring high-throughput, high-coverage single-cell sequencing.

Clustering and classification methods for single-cell RNA-sequencing data

Appropriate ways to measure the similarity between single-cell RNA-sequencing (scRNA-seq) data are ubiquitous in bioinformatics, but using single clustering or classification methods to process scRNA-seq data is generally difficult. This has led to the emergence of integrated methods and tools that aim to automatically process specific problems associated with scRNA-seq data. These approaches have attracted a lot of interest in bioinformatics and related fields. In this paper, we systematically review the integrated methods and tools, highlighting the pros and cons of each approach. We not only pay particular attention to clustering and classification methods but also discuss methods that have emerged recently as powerful alternatives, including nonlinear and linear methods and descending dimension methods. Finally, we focus on clustering and classification methods for scRNA-seq data, in particular, integrated methods, and provide a comprehensive description of scRNA-seq data and download URLs.

Gene Expression Noise Produces Cell-to-Cell Heterogeneity in Eukaryotic Homologous Recombination Rate

Variation in gene expression among genetically identical individual cells (called gene expression noise) directly contributes to phenotypic diversity. Whether such variation can impact genome stability and lead to variation in genotype remains poorly explored. We addressed this question by investigating whether noise in the expression of genes affecting homologous recombination (HR) activity either directly (RAD52) or indirectly (RAD27) confers cell-to-cell heterogeneity in HR rate in Saccharomyces cerevisiae. Using cell sorting to isolate subpopulations with various expression levels, we show that spontaneous HR rate is highly heterogeneous from cell-to-cell in clonal populations depending on the cellular amount of proteins affecting HR activity. Phleomycin-induced HR is even more heterogeneous, showing that RAD27 expression variation strongly affects the rate of recombination from cell-to-cell. Strong variations in HR rate between subpopulations are not correlated to strong changes in cell cycle stage. Moreover, this heterogeneity occurs even when simultaneously sorting cells at equal expression level of another gene involved in DNA damage response (BMH1) that is upregulated by DNA damage, showing that the initiating DNA damage is not responsible for the observed heterogeneity in HR rate. Thus gene expression noise seems mainly responsible for this phenomenon. Finally, HR rate non-linearly scales with Rad27 levels showing that total amount of HR cannot be explained solely by the time- or population-averaged Rad27 expression. Altogether, our data reveal interplay between heterogeneity at the gene expression and genetic levels in the production of phenotypic diversity with evolutionary consequences from microbial to cancer cell populations.

Per-Nucleus Crossover Covariation and Implications for Evolution

Crossing over is a nearly universal feature of sexual reproduction. Here, analysis of crossover numbers on a per-chromosome and per-nucleus basis reveals a fundamental, evolutionarily conserved feature of meiosis: within individual nuclei, crossover frequencies covary across different chromosomes. This effect results from per-nucleus covariation of chromosome axis lengths. Crossovers can promote evolutionary adaptation. However, the benefit of creating favorable new allelic combinations must outweigh the cost of disrupting existing favorable combinations. Covariation concomitantly increases the frequencies of gametes with especially high, or especially low, numbers of crossovers, and thus might concomitantly enhance the benefits of crossing over while reducing its costs. A four-locus population genetic model suggests that such an effect can pertain in situations where the environment fluctuates: hyper-crossover gametes are advantageous when the environment changes while hypo-crossover gametes are advantageous in periods of environmental stasis. These findings reveal a new feature of the basic meiotic program and suggest a possible adaptive advantage.

Factors influencing meiotic recombination revealed by whole-genome sequencing of single sperm

Recombination is critical to meiosis and evolution, yet many aspects of the physical exchange of DNA via crossovers remain poorly understood. We report an approach for single-cell whole-genome DNA sequencing by which we sequenced 217 individual hybrid mouse sperm, providing a kilobase-resolution genome-wide map of crossovers. Combining this map with molecular assays measuring stages of recombination, we identified factors that affect crossover probability, including PRDM9 binding on the non-initiating template homolog and telomere proximity. These factors also influence the time for sites of recombination-initiating DNA double-strand breaks to find and engage their homologs, with rapidly engaging sites more likely to form crossovers. We show that chromatin environment on the template homolog affects positioning of crossover breakpoints. Our results also offer insights into recombination in the pseudoautosomal region.

Live sperm trap microarray for high throughput imaging and analysis

There is a growing appreciation and understanding of cell-to-cell variability in biological samples. However, research and clinical practice in male fertility has relied on population, or sample-based characteristics. Single-cell resolution is particularly important given the winner-takes-all nature of both natural and in vitro fertilization: it is the properties of a single cell, not the population, that are passed to the next generation. While there are a range of methods for single cell analysis, arraying a larger number of live sperm has not been possible due to the strong locomotion of the cells. Here we present a 103-trap microarray that traps, aligns and arrays individual live sperm. The method enables high-resolution imaging of the aligned cell head, the application of dye-based DNA and mitochondrial analyses, and the quantification of motility characteristics, such as tail beat. In testing, a 2400-post array trapped ∼400 sperm for individual analyses of tail beating frequency and amplitude, DNA integrity via acridine orange staining, and mitochondrial activity via staining. While literature results are mixed regarding a possible correlation between motility and DNA integrity of sperm at sample-level, results here find no statistical correlation between tail beat characteristics and DNA integrity at the cell-level. The trap array uniquely enables the high-throughput study of individual live sperm in semen samples - assessing the inherently single-cell selection process of fertilization, with single-cell resolution.

Selection for longer lived sperm within ejaculate reduces reproductive ageing in offspring

Males produce numerous sperm in a single ejaculate that greatly outnumber their potential egg targets. Recent studies found that phenotypic and genotypic variation among sperm in a single ejaculate of a male affects the fitness and performance of the resulting offspring. Specifically, within-ejaculate sperm selection for sperm longevity increased the performance of the resulting offspring in several key life-history traits in early life. Because increased early-life reproductive performance often correlates with rapid ageing, it is possible that within-ejaculate sperm selection increases early-life fitness at the cost of accelerated senescence. Alternatively, within-ejaculate sperm selection could improve offspring quality throughout the life cycle, including reduced age-specific deterioration. We tested the two alternative hypotheses in an experimental setup using zebrafish Danio rerio. We found that within-ejaculate sperm selection for sperm longevity reduced age-specific deterioration of fecundity and offspring survival but had no effect on fertilization success in males. Remarkably, we found an opposing effect of within-ejaculate sperm selection on female fecundity, where selection for sperm longevity resulted in increased early-life performance followed by a slow decline, while females sired by unselected sperm started low but increased their fecundity with age. Intriguingly, within-ejaculate sperm selection also reduced the age-specific decline in fertilization success in females, suggesting that selection for sperm longevity improves at least some aspects of female reproductive ageing. These results demonstrate that within-ejaculate variation in sperm phenotype contributes to individual variation in animal life histories in the two sexes and may have important implications for assisted fertilization programs in livestock and humans.

A rigorous measure of genome-wide genetic shuffling that takes into account crossover positions and Mendel's second law

Comparative studies in evolutionary genetics rely critically on evaluation of the total amount of genetic shuffling that occurs during gamete production. Such studies have been hampered by the absence of a direct measure of this quantity. Existing measures consider crossing-over by simply counting the average number of crossovers per meiosis. This is qualitatively inadequate, because the positions of crossovers along a chromosome are also critical: a crossover toward the middle of a chromosome causes more shuffling than a crossover toward the tip. Moreover, traditional measures fail to consider shuffling from independent assortment of homologous chromosomes (Mendel's second law). Here, we present a rigorous measure of genome-wide shuffling that does not suffer from these limitations. We define the parameter [Formula: see text] as the probability that the alleles at two randomly chosen loci are shuffled during gamete production. This measure can be decomposed into separate contributions from crossover number and position and from independent assortment. Intrinsic implications of this metric include the fact that [Formula: see text] is larger when crossovers are more evenly spaced, which suggests a selective advantage of crossover interference. Utilization of [Formula: see text] is enabled by powerful emergent methods for determining crossover positions either cytologically or by DNA sequencing. Application of our analysis to such data from human male and female reveals that (i) [Formula: see text] in humans is close to its maximum possible value of 1/2 and that (ii) this high level of shuffling is due almost entirely to independent assortment, the contribution of which is ∼30 times greater than that of crossovers.

Mapping the Mouse Cell Atlas by Microwell-Seq

Single-cell RNA sequencing (scRNA-seq) technologies are poised to reshape the current cell-type classification system. However, a transcriptome-based single-cell atlas has not been achieved for complex mammalian systems. Here, we developed Microwell-seq, a high-throughput and low-cost scRNA-seq platform using simple, inexpensive devices. Using Microwell-seq, we analyzed more than 400,000 single cells covering all of the major mouse organs and constructed a basic scheme for a mouse cell atlas (MCA). We reveal a single-cell hierarchy for many tissues that have not been well characterized previously. We built a web-based "single-cell MCA analysis" pipeline that accurately defines cell types based on single-cell digital expression. Our study demonstrates the wide applicability of the Microwell-seq technology and MCA resource.

Chromosomal scan of single sperm cells by combining fluorescence-activated cell sorting and next-generation sequencing

PURPOSE

The purpose of this study was to develop a feasible approach for single sperm isolation and chromosome analysis by next-generation sequencing (NGS).

METHODS

Single sperm cells were isolated from semen samples of normozoospermic male and an infertile reciprocal translocation (RcT) carrier with the 46,XY,t(7;13)(p12;q12.1) karyotype using the optimized fluorescence-activated cell sorting (FACS) technique. Genome profiling was performed using NGS.

RESULTS

Following whole-genome amplification, NGS, and quality control, the final chromosome analysis was performed on 31 and 6 single cell samples derived from the RcT carrier and normozoospermic male, respectively. All sperm cells from normozoospermic male showed a normal haploid 23-chromosome profile. For the RcT carrier, the sequencing data revealed that 64.5% of sperm cells harbored different variants of chromosome aberrations, involving deletion of 7p or 7q, duplication of 7p, and duplication of 13q, which is concordant with the expected chromosome segregation patterns observed in balanced translocation carriers. In one sample, a duplication of 9q was also detected.

CONCLUSIONS

We optimized FACS protocol for simple and efficient isolation of single human sperm cells that subsequently enabled a successful genome-wide chromosome profiling and identification of segmental aneuploidies from these individual cells, following NGS analysis. This approach may be useful for analyzing semen samples of infertile men or chromosomal aberration carriers to facilitate the reproductive risk assessment.

Zygotic gene activation in the chicken occurs in two waves, the first involving only maternally derived genes

The first wave of transcriptional activation occurs after fertilisation in a species-specific pattern. Despite its importance to initial embryonic development, the characteristics of transcription following fertilisation are poorly understood in Aves. Here, we report detailed insights into the onset of genome activation in chickens. We established that two waves of transcriptional activation occurred, one shortly after fertilisation and another at Eyal-Giladi and Kochav Stage V. We found 1544 single nucleotide polymorphisms across 424 transcripts derived from parents that were expressed in offspring during the early embryonic stages. Surprisingly, only the maternal genome was activated in the zygote, and the paternal genome remained silent until the second-wave, regardless of the presence of a paternal pronucleus or supernumerary sperm in the egg. The identified maternal genes involved in cleavage that were replaced by bi-allelic expression. The results demonstrate that only maternal alleles are activated in the chicken zygote upon fertilisation, which could be essential for early embryogenesis and evolutionary outcomes in birds.

The early stages of animal development involve a handover of genetic control. Initially, the egg cell is maintained by genetic information inherited from the mother, but soon after fertilization it starts to depend on its own genes instead. Activating genes inside the fertilized egg cell (zygote) so that they can take control of development is known as zygotic genome activation.

Despite the fact that birds are often used to study how embryos develop, zygotic genome activation in birds is not well understood. Fertilization in birds, including chickens, is different to mammals in that it requires multiple sperm to fertilize an egg cell. As such, zygotic genome activation in birds is likely to differ from that in mammals.

By examining gene expression in embryos from mixed-breed chickens, Hwang, Seo et al. showed that there are two stages of zygotic genome activation in chickens. The genes derived from the mother become active in the first stage, while genes from the father become active in the second stage. Genome activation in birds is therefore very different to the same process in mammals, which involves genome activation of both parents from the first stage. This extra level of control may help to prevent genetic complications resulting from the presence of multiple sperm, each of which carries a different set of genes from the father.

Elucidation of the developmental mechanism of ovarian mature cystic teratomas using B allele-frequency plots of single nucleotide polymorphism array data

Ovarian mature cystic teratomas (MCTs) originate from post-meiotic germ cells. Conventional methods such as karyotyping or short tandem repeat-polymorphism analysis may be used to better classify MCTs, although such data would be insufficient. The aim of this study was to elucidate the origin of ovarian MCTs using B allele-frequency (BAF) plots of single nucleotide polymorphism array data. MCTs can be classified in terms of the zygosity of the centromeres and distal chromosome regions. We evaluated the zygosity of all chromosomes from 38 MCT specimens using BAF plot data. BAF plots were used to determine the homozygous and heterozygous regions over the whole genome. Theoretically, MCTs originated from the fusion of two ova (previously referred to as type V MCTs) should have a mixed pattern of centromeric zygosity, that is, a combination of heterozygous and homozygous regions in the centromeric regions. However, no MCTs in this study met this criterion. We identified 13 type I MCTs, 14 type II MCTs, and 11 type III MCTs. In addition, BAF plots facilitated the construction of recombination maps at the whole-genome level for type I and II MCTs. No crossover, especially in the short arms, contributed to the failure of meiosis I, resulting in type I MCTs. Crossover in all arms might assure the normal progress of meiosis in human oocytes. In conclusion, our findings indicate that BAF plots can elucidate the developmental mechanism of MCTs, and further serve as useful analytical tools for analyzing human oocyte meiosis, and related aberrations.

Somatic Mutagenesis in Mammals and Its Implications for Human Disease and Aging

DNA mutations as a consequence of errors during DNA damage repair, replication, or mitosis are the substrate for evolution. In multicellular organisms, mutations can occur in the germline and also in somatic tissues, where they are associated with cancer and other chronic diseases and possibly with aging. Recent advances in high-throughput sequencing have made it relatively easy to study germline de novo mutations, but in somatic cells, the vast majority of mutations are low-abundant and can be detected only in clonal lineages, such as tumors, or single cells. Here we review recent results on somatic mutations in normal human and animal tissues with a focus on their possible functional consequences.

Genome wide analysis of meiotic recombination in yeast: For a few SNPs more

Diploid organisms undergo meiosis to produce haploid germ cells. Crossover events during meiosis promote genetic diversity and facilitate accurate chromosome segregation. The baker's yeast Saccharomyces cerevisiae is extensively used as a model for analysis of meiotic recombination. Conventional methods for measuring recombination events in S. cerevisiae have been limited by the number and density of genetic markers. Next generation sequencing (NGS)-based analysis of hybrid yeast genomes bearing thousands of heterozygous single nucleotide polymorphism (SNP) markers has revolutionized analysis of meiotic recombination. By facilitating analysis of marker segregation in the whole genome with unprecedented resolution, this method has resulted in the generation of high-resolution recombination maps in wild-type and meiotic mutants. These studies have provided novel insights into the mechanism of meiotic recombination. In this review, we discuss the methodology, challenges, insights and future prospects of using NGS-based methods for whole genome analysis of meiotic recombination. The objective is to facilitate the use of these high through-put sequencing methods for the analysis of meiotic recombination given their power to provide significant new insights into the process. © 2018 The Authors. IUBMB Life published by Wiley Periodicals, Inc. on behalf of International Union of Biochemistry and Molecular Biology, 70(8):743-752, 2018.

Reproductive genetics and the aging male

PURPOSE

To examine current evidence of the known effects of advanced paternal age on sperm genetic and epigenetic changes and associated birth defects and diseases in offspring.

METHODS

Review of published PubMed literature.

RESULTS

Advanced paternal age (> 40 years) is associated with accumulated damage to sperm DNA and mitotic and meiotic quality control mechanisms (mismatch repair) during spermatogenesis. This in turn causes well-delineated abnormalities in sperm chromosomes, both numerical and structural, and increased sperm DNA fragmentation (3%/year of age) and single gene mutations (relative risk, RR 10). An increase in related abnormalities in offspring has also been described, including miscarriage (RR 2) and fetal loss (RR 2). There is also a significant increase in rare, single gene disorders (RR 1.3 to 12) and congenital anomalies (RR 1.2) in offspring. Current research also suggests that autism, schizophrenia, and other forms of "psychiatric morbidity" are more likely in offspring (RR 1.5 to 5.7) with advanced paternal age. Genetic defects related to faulty sperm quality control leading to single gene mutations and epigenetic alterations in several genetic pathways have been implicated as root causes.

CONCLUSIONS

Advanced paternal age is associated with increased genetic and epigenetic risk to offspring. However, the precise age at which risk develops and the magnitude of the risk are poorly understood or may have gradual effects. Currently, there are no clinical screenings or diagnostic panels that target disorders associated with advanced paternal age. Concerned couples and care providers should pursue or recommend genetic counseling and prenatal testing regarding specific disorders.

Genotyping single-sperm cells by universal MARSALA enables the acquisition of linkage information for combined pre-implantation genetic diagnosis and genome screening

PURPOSE

This paper aims to investigate the feasibility of performing pre-implantation genetic diagnosis (PGD) and pre-implantation genetic screening (PGS) simultaneously by a universal strategy without the requirement of genotyping relevant affected family members or lengthy preliminary work on linkage analysis.

METHODS

By utilizing a universal Mutated Allele Revealed by Sequencing with Aneuploidy and Linkage Analyses (MARSALA) strategy based on low depth whole genome sequencing (~3x), not involving specific primers' design nor the enrichment of SNP markers for haplotype construction. Single-sperm cells and trephectoderm cells from in vitro fertilized embryos from a couple carrying HBB mutations were genotyped. Haplotypes of paternal alleles were constructed and investigated in embryos, and the chromosome copy number profiles were simultaneously analyzed.

RESULTS

The universal MARSALA strategy allows the selection of a euploid embryo free of disease mutations for in uterus transfer and successful pregnancy. A follow-up amniocentesis was performed at 17 weeks of gestation to confirm the PGD/PGS results.

CONCLUSION

We present the first successful PGD procedure based on genotyping multiple single-sperm cells to obtain SNP linkage information. Our improved PGD/PGS procedure does not require genotyping the proband or relevant family members and therefore can be applicable to a wider population of patients when conducting PGD for monogenic disorders.

Distinctive types of postzygotic single-nucleotide mosaicisms in healthy individuals revealed by genome-wide profiling of multiple organs

Postzygotic single-nucleotide mosaicisms (pSNMs) have been extensively studied in tumors and are known to play critical roles in tumorigenesis. However, the patterns and origin of pSNMs in normal organs of healthy humans remain largely unknown. Using whole-genome sequencing and ultra-deep amplicon re-sequencing, we identified and validated 164 pSNMs from 27 postmortem organ samples obtained from five healthy donors. The mutant allele fractions ranged from 1.0% to 29.7%. Inter- and intra-organ comparison revealed two distinctive types of pSNMs, with about half originating during early embryogenesis (embryonic pSNMs) and the remaining more likely to result from clonal expansion events that had occurred more recently (clonal expansion pSNMs). Compared to clonal expansion pSNMs, embryonic pSNMs had higher proportion of C>T mutations with elevated mutation rate at CpG sites. We observed differences in replication timing between these two types of pSNMs, with embryonic and clonal expansion pSNMs enriched in early- and late-replicating regions, respectively. An increased number of embryonic pSNMs were located in open chromatin states and topologically associating domains that transcribed embryonically. Our findings provide new insights into the origin and spatial distribution of postzygotic mosaicism during normal human development.

Genomic mosaicism led by postzygotic mutation is the major cause of cancers and many non-cancer developmental disorders. Theoretically, postzygotic mutations should be accumulated during the developmental process of healthy individuals, but the genome-wide characterization of postzygotic mosaicisms across many organ types of the same individual remained limited. In this study, we identified and validated two types of postzygotic mosaicism from the whole-genomes of 27 organs obtained from five healthy donors. We further found that the postzygotic mosaicisms arising during early embryogenesis and later clonal expansion events show distinct genomic patterns in mutation spectrum, replication timing, and chromatin status.

Blending genomes: distributive conjugal transfer in mycobacteria, a sexier form of HGT

This review discusses a novel form of horizontal gene transfer (HGT) found in mycobacteria called Distributive Conjugal Transfer (DCT). While satisfying the criteria for conjugation, DCT occurs by a mechanism so distinct from oriT-mediated conjugation that it could be considered a fourth category of HGT. DCT involves the transfer of chromosomal DNA between mycobacteria and, most significantly, generates transconjugants with mosaic genomes of the parental strains. Multiple segments of donor chromosomal DNA can be co-transferred regardless of their location or the genetic selection and, as a result, the transconjugant genome contains many donor-derived segments; hence the name DCT. This distinguishing feature of DCT separates it from the other known mechanisms of HGT, which generally result in the introduction of a single, defined segment of DNA into the recipient chromosome (Fig. ). Moreover, these mosaic progeny are generated from a single conjugal event, which provides enormous capacity for rapid adaptation and evolution, again distinguishing it from the three classical modes of HGT. Unsurprisingly, the unusual mosaic products of DCT are generated by a conjugal mechanism that is also unusual. Here, we will describe the unique features of DCT and contrast those to other mechanisms of HGT, both from a mechanistic and an evolutionary perspective. Our focus will be on transfer of chromosomal DNA, as opposed to plasmid mobilization, because DCT mediates transfer of chromosomal DNA and is a chromosomally encoded process.

Single cell genomics to study DNA and chromosome changes in human gametes and embryos

Genomic and chromosomal changes occur with a high rate in the germline and preimplantation embryos. To study such changes directly in the germline of mammals requires access to material as well as single cell genomics. Recent improvements in embryology and single-cell DNA amplification make it possible to study the genomic changes directly in human oocytes, sperm, and preimplantation embryos. This is particularly important for the study of chromosome segregation directly in human oocytes and preimplantation embryos. Here, we present a practical approach how to obtain high quality DNA sequences and genotypes from single cells, using manual handling of the material that makes it possible to detect genomic changes in meiosis and mitosis spanning the entire range from single nucleotide changes to whole chromosome aneuploidies.

Discovery of Lineage-Specific Genome Change in Rice Through Analysis of Resequencing Data

New mutations are rare, which makes their discovery laborious and time-consuming. Arthur and Bennetzen describe an approach for enriching recent mutations that relies only on a reference genome sequence and resequencing data for other...

Genome comparisons provide information on the nature of genetic change, but such comparisons are challenged to differentiate the importance of the actual sequence change processes relative to the role of selection. This problem can be overcome by identifying changes that have not yet had the time to undergo millions of years of natural selection. We describe a strategy to discover accession-specific changes in the rice genome using an abundant resource routinely provided for many genome analyses, resequencing data. The sequence of the fully sequenced rice genome from variety Nipponbare was compared to the pooled (∼114×) resequencing data from 126 japonica rice accessions to discover “Nipponbare-specific” sequences. Analyzing nonrepetitive sequences, 8504 “candidate” Nipponbare-specific changes were detected, of which around two-thirds are true novel sequence changes and the rest are predicted genome sequencing errors. Base substitutions outnumbered indels in this data set by > 28:1, with ∼8:5 bias toward transversions over transitions, and no transposable element insertions or excisions were observed. These results indicate that the strategy employed is effective for finding recent sequence changes, sequencing errors, and rare alleles in any organism that has both a reference genome sequence and a wealth of resequencing data.

Construction of PRDM9 allele-specific recombination maps in cattle using large-scale pedigree analysis and genome-wide single sperm genomics

PRDM9 contributes to hybrid sterility and species evolution. However, its role is to be confirmed in cattle, a major domesticated livestock species. We previously found an association near PRDM9 with cattle recombination features, but the causative variants are still unknown. Using millions of genotyped cattle with pedigree information, we characterized five PRDM9 alleles and generated allele-specific recombination maps. By examining allele-specific recombination patterns, we observed the impact of PRDM9 on global distribution of recombination, especially in the two ends of chromosomes. We also showed strong associations between recombination hotspot regions and functional mutations within PRDM9 zinc finger domain. More importantly, we found one allele of PRDM9 to be very different from others in both protein composition and recombination landscape, indicating the causative role of this allele on the association between PRDM9 and cattle recombination. When comparing recombination maps from sperm and pedigree data, we observed similar genome-wide recombination patterns, validating the quality of pedigree-based results. Collectively, these evidence supported PRDM9 alleles as causal variants for the reported association with cattle recombination. Our study comprehensively surveyed the bovine PRDM9 alleles, generated allele-specific recombination maps, and expanded our understanding of the role of PRDM9 on genome distribution of recombination.

Hydrogel droplet single-cell processing: DNA purification, handling, release, and on-chip linearization

The preparation and handling of mammalian single-cell genomic DNA is limited by the complexity bottleneck inherent to performing multi-step, multi-reagent operations in a microfluidic environment. We have developed a method for benchtop preparation of high-molecular weight, intact, single-cell genomes and demonstrate the extraction of long nucleic acid molecules in a microfluidic system. Lymphoblasts are encapsulated inside of alginate microparticles using a droplet microfluidics, and cells are lysed in bulk. The purified genomes are then delivered to and imaged on a dedicated microfluidic device. High-molecular weight DNA is protected from shear and retains its original cellular identity. Using this encapsulation protocol, we were able to extract individual nucleic acid strands on the millimeter scale inside of a microfluidic channel.

Sperm selection by thermotaxis improves ICSI outcome in mice

The ejaculate is a heterogeneous pool of spermatozoa containing only a small physiologically adequate subpopulation for fertilization. As there is no method to isolate this subpopulation, its specific characteristics are unknown. This is one of the main reasons why we lack effective tools to identify male infertility and for the low efficiency of assisted reproductive technologies. The aim of this study was to improve ICSI outcome by sperm selection through thermotaxis. Here we show that a specific subpopulation of mouse and human spermatozoa can be selected in vitro by thermotaxis and that this subpopulation is the one that enters the fallopian tube in mice. Further, we confirm that these selected spermatozoa in mice and humans show a much higher DNA integrity and lower chromatin compaction than unselected sperm, and in mice, they give rise to more and better embryos through intracytoplasmic sperm injection, doubling the number of successful pregnancies. Collectively, our results indicate that a high quality sperm subpopulation is selected in vitro by thermotaxis and that this subpopulation is also selected in vivo within the fallopian tube possibly by thermotaxis.

Aging and neurodegeneration are associated with increased mutations in single human neurons

It has long been hypothesized that aging and neurodegeneration are associated with somatic mutation in neurons; however, methodological hurdles have prevented testing this hypothesis directly. We used single-cell whole-genome sequencing to perform genome-wide somatic single-nucleotide variant (sSNV) identification on DNA from 161 single neurons from the prefrontal cortex and hippocampus of 15 normal individuals (aged 4 months to 82 years), as well as 9 individuals affected by early-onset neurodegeneration due to genetic disorders of DNA repair (Cockayne syndrome and xeroderma pigmentosum). sSNVs increased approximately linearly with age in both areas (with a higher rate in hippocampus) and were more abundant in neurodegenerative disease. The accumulation of somatic mutations with age-which we term genosenium-shows age-related, region-related, and disease-related molecular signatures and may be important in other human age-associated conditions.

Single-Cell Genomic Analysis in Plants

Individual cells in an organism are variable, which strongly impacts cellular processes. Advances in sequencing technologies have enabled single-cell genomic analysis to become widespread, addressing shortcomings of analyses conducted on populations of bulk cells. While the field of single-cell plant genomics is in its infancy, there is great potential to gain insights into cell lineage and functional cell types to help understand complex cellular interactions in plants. In this review, we discuss current approaches for single-cell plant genomic analysis, with a focus on single-cell isolation, DNA amplification, next-generation sequencing, and bioinformatics analysis. We outline the technical challenges of analysing material from a single plant cell, and then examine applications of single-cell genomics and the integration of this approach with genome editing. Finally, we indicate future directions we expect in the rapidly developing field of plant single-cell genomic analysis.