Using runs of homozygosity and machine learning to disentangle sources of inbreeding and infer self-fertilization rates.

Runs of homozygosity (ROHs) are indicative of elevated homozygosity and inbreeding due to mating of closely related individuals. Self-fertilization can be a major source of inbreeding which elevates genome-wide homozygosity and thus should also create long ROHs. While ROHs are frequently used to understand inbreeding in the context of conservation and selective breeding, as well as for consanguinity of populations and their demographic history, it remains unclear how ROH characteristics are altered by selfing and if this confounds expected signatures of inbreeding due to demographic change. Using simulations, we study the impact of the mode of reproduction and demographic history on ROHs. We apply random forests to identify unique characteristics of ROHs, indicative of different sources of inbreeding. We pinpoint distinct features of ROHs that can be used to better characterize the type of inbreeding the population was subjected to and to predict outcrossing rates and complex demographic histories. Using additional simulations and four empirical datasets, two from highly selfing species and two from mixed-maters, we predict the selfing rate and validate our estimations. We find that self-fertilization rates are successfully identified even with complex demography. Population genetic summary statistics improve algorithm accuracy particularly in the presence of additional inbreeding, e.g., from population bottlenecks. Our findings highlight the importance of ROHs in disentangling confounding factors related to various sources of inbreeding and demonstrate situations where such sources cannot be differentiated. Additionally, our random forest models provide a novel tool to the community for inferring selfing rates using genomic data.

A genome-wide comprehensive analysis of nucleosome positioning in yeast.

In eukaryotic cells, the one-dimensional DNA molecules need to be tightly packaged into the spatially constraining nucleus. Folding is achieved on its lowest level by wrapping the DNA around nucleosomes. Their arrangement regulates other nuclear processes, such as transcription and DNA repair. Despite strong efforts to study nucleosome positioning using Next Generation Sequencing (NGS) data, the mechanism of their collective arrangement along the gene body remains poorly understood. Here, we classify nucleosome distributions of protein-coding genes in Saccharomyces cerevisiae according to their profile similarity and analyse their differences using functional Principal Component Analysis. By decomposing the NGS signals into their main descriptive functions, we compared wild type and chromatin remodeler-deficient strains, keeping position-specific details preserved whilst considering the nucleosome arrangement as a whole. A correlation analysis with other genomic properties, such as gene size and length of the upstream Nucleosome Depleted Region (NDR), identified key factors that influence the nucleosome distribution. We reveal that the RSC chromatin remodeler-which is responsible for NDR maintenance-is indispensable for decoupling nucleosome arrangement within the gene from positioning outside, which interfere in rsc8-depleted conditions. Moreover, nucleosome profiles in chd1Δ strains displayed a clear correlation with RNA polymerase II presence, whereas wild type cells did not indicate a noticeable interdependence. We propose that RSC is pivotal for global nucleosome organisation, whilst Chd1 plays a key role for maintaining local arrangement.

Purging due to self-fertilization does not prevent accumulation of expansion load.

As species expand their geographic ranges, colonizing populations face novel ecological conditions, such as new environments and limited mates, and suffer from evolutionary consequences of demographic change through bottlenecks and mutation load accumulation. Self-fertilization is often observed at species range edges and, in addition to countering the lack of mates, is hypothesized as an evolutionary advantage against load accumulation through increased homozygosity and purging. We study how selfing impacts the accumulation of genetic load during range expansion via purging and/or speed of colonization. Using simulations, we disentangle inbreeding effects due to demography versus due to selfing and find that selfers expand faster, but still accumulate load, regardless of mating system. The severity of variants contributing to this load, however, differs across mating system: higher selfing rates purge large-effect recessive variants leaving a burden of smaller-effect alleles. We compare these predictions to the mixed-mating plant Arabis alpina, using whole-genome sequences from refugial outcrossing populations versus expanded selfing populations. Empirical results indicate accumulation of expansion load along with evidence of purging in selfing populations, concordant with our simulations, suggesting that while purging is a benefit of selfing evolving during range expansions, it is not sufficient to prevent load accumulation due to range expansion.

Functional interplay between Mediator and RSC chromatin remodeling complex controls nucleosome-depleted region maintenance at promoters.

Chromatin organization is crucial for transcriptional regulation in eukaryotes. Mediator is an essential and conserved co-activator thought to act in concert with chromatin regulators. However, it remains largely unknown how their functions are coordinated. Here, we provide evidence in the yeast Saccharomyces cerevisiae that Mediator establishes physical contact with RSC (Remodels the Structure of Chromatin), a conserved and essential chromatin remodeling complex that is crucial for nucleosome-depleted region (NDR) formation. We determine the role of Mediator-RSC interaction in their chromatin binding, nucleosome occupancy, and transcription on a genomic scale. Mediator and RSC co-localize on wide NDRs of promoter regions, and specific Mediator mutations affect nucleosome eviction and TSS-associated +1 nucleosome stability. This work shows that Mediator contributes to RSC remodeling function to shape NDRs and maintain chromatin organization on promoter regions. It will help in our understanding of transcriptional regulation in the chromatin context relevant for severe diseases.

A quantitative modelling approach for DNA repair on a population scale.

The great advances of sequencing technologies allow the in vivo measurement of nuclear processes-such as DNA repair after UV exposure-over entire cell populations. However, data sets usually contain only a few samples over several hours, missing possibly important information in between time points. We developed a data-driven approach to analyse CPD repair kinetics over time in Saccharomyces cerevisiae. In contrast to other studies that consider sequencing signals as an average behaviour, we understand them as the superposition of signals from independent cells. By motivating repair as a stochastic process, we derive a minimal model for which the parameters can be conveniently estimated. We correlate repair parameters to a variety of genomic features that are assumed to influence repair, including transcription rate and nucleosome density. The clearest link was found for the transcription unit length, which has been unreported for budding yeast to our knowledge. The framework hence allows a comprehensive analysis of nuclear processes on a population scale.

Male meiotic recombination rate varies with seasonal temperature fluctuations in wild populations of autotetraploid Arabidopsis arenosa.

Meiosis, the cell division by which eukaryotes produce haploid gametes, is essential for fertility in sexually reproducing species. This process is sensitive to temperature, and can fail outright at temperature extremes. At less extreme values, temperature affects the genome-wide rate of homologous recombination, which has important implications for evolution and population genetics. Numerous studies in laboratory conditions have shown that recombination rate plasticity is common, perhaps nearly universal, among eukaryotes. These studies have also shown that variation in the length or timing of stresses can strongly affect results, raising the important question whether these findings translate to more variable field conditions. Moreover, lower or higher recombination rate could cause certain kinds of meiotic aberrations, especially in polyploid species-raising the additional question whether temperature fluctuations in field conditions cause problems. Here, we tested whether (1) recombination rate varies across a season in the wild in two natural populations of autotetraploid Arabidopsis arenosa, (2) whether recombination rate correlates with temperature fluctuations in nature, and (3) whether natural temperature fluctuations might cause meiotic aberrations. We found that plants in two genetically distinct populations showed a similar plastic response with recombination rate increases correlated with both high and low temperatures. In addition, increased recombination rate correlated with increased multivalent formation, especially at lower temperature, hinting that polyploids in particular may suffer meiotic problems in conditions they encounter in nature. Our results show that studies of recombination rate plasticity done in laboratory settings inform our understanding of what happens in nature.

Selective Loss of Diversity in Doubled-Haploid Lines from European Maize Landraces.

Maize landraces are well adapted to their local environments and present valuable sources of genetic diversity for breeding and conservation. But the maintenance of open-pollinated landraces in ex-situ programs is challenging, as regeneration of seed can often lead to inbreeding depression and the loss of diversity due to genetic drift. Recent reports suggest that the production of doubled-haploid (DH) lines from landraces may serve as a convenient means to preserve genetic diversity in a homozygous form that is immediately useful for modern breeding. The production of doubled-haploid (DH) lines presents an extreme case of inbreeding which results in instantaneous homozygosity genome-wide. Here, we analyzed the effect of DH production on genetic diversity, using genome-wide SNP data from hundreds of individuals of five European landraces and their related DH lines. In contrast to previous findings, we observe a dramatic loss of diversity at both the haplotype level and that of individual SNPs. We identify thousands of SNPs that exhibit allele frequency differences larger than expected under models of neutral genetic drift and document losses of shared haplotypes. We find evidence consistent with selection at functional sites that are potentially involved in the diversity differences between landrace and DH populations. Although we were unable to uncover more details about the mode of selection, we conclude that landrace DH lines may be a valuable tool for the introduction of variation into maize breeding programs but come at the cost of decreased genetic diversity.

Crossing Methods and Cultivation Conditions for Rapid Production of Segregating Populations in Three Grain Amaranth Species.

Grain amaranths (Amaranthus spp.) have been cultivated for thousands of years in Central and South America. Their grains are of high nutritional value, but the low yield needs to be increased by selection of superior genotypes from genetically diverse breeding populations. Amaranths are adapted to harsh conditions and can be cultivated on marginal lands although little is known about their physiology. The development of controlled growing conditions and efficient crossing methods is important for research on and improvement of this ancient crop. Grain amaranth was domesticated in the Americas and is highly self-fertilizing with a large inflorescence consisting of thousands of very small flowers. We evaluated three different crossing methods (open pollination, hot water emasculation and hand emasculation) for their efficiency in amaranth and validated them with genetic markers. We identified cultivation conditions that allow an easy control of flowering time by day length manipulation and achieved flowering times of 4 weeks and generation times of 2 months. All three different crossing methods successfully produced hybrid F1 offspring, but with different success rates. Open pollination had the lowest (10%) and hand emasculation the highest success rate (74%). Hot water emasculation showed an intermediate success rate (26%) with a maximum of 94% success. It is simple to perform and suitable for a more large-scale production of hybrids. We further evaluated 11 single nucleotide polymorphism (SNP) markers and found that they were sufficient to validate all crosses of the genotypes used in this study for intra- and interspecific hybridizations. Despite its very small flowers, crosses in amaranth can be carried out efficiently and evaluated with inexpensive SNP markers. Suitable growth conditions strongly reduce the generation time and allow the control of plant height, flowering time, and seed production. In combination, this enables the rapid production of segregating populations which makes amaranth an attractive model for basic plant research but also facilitates further the improvement of this ancient crop by plant breeding.