Homozygosity Mapping Reveals Population History and Trait Architecture in Self-Incompatible Pear (Pyrus spp.)

Runs of homozygosity reveal contrasting histories of inbreeding across global lineages of the edible porcini mushroom, Boletus edulis

Inbreeding, the mating of individuals that are related through common ancestry, is of central importance in evolutionary and conservation biology due to its impacts on individual fitness and population dynamics. However, while advanced genomic approaches have revolutionised the study of inbreeding in animals, genomic studies of inbreeding are rare in plants and lacking in fungi. We investigated global patterns of inbreeding in the prized edible porcini mushroom Boletus edulis using 225 whole genomes from seven lineages distributed across the northern hemisphere. Genomic inbreeding was quantified using runs of homozygosity (ROHs). We found appreciable variation both among and within lineages, with some individuals having over 20% of their genomes in ROHs. Much of this variation could be explained by a combination of elevation and latitude, and to a lesser extent by predicted habitat suitability during the last glacial maximum. In line with this, the majority of ROHs were short, reflecting ancient common ancestry dating back approximately 200-1700 generations ago, while longer ROHs indicative of recent common ancestry (less than approximately 50 generations ago) were infrequent. Our study reveals the inbreeding legacy of major climatic events in a widely distributed forest mutualist, aligning with prevailing theories and empirical studies of the impacts of historical glaciation events on the dominant forest tree species of the northern hemisphere.

Genomic prediction and genome-wide association study using combined genotypic data from different genotyping systems: application to apple fruit quality traits

With advances in next-generation sequencing technologies, various marker genotyping systems have been developed for genomics-based approaches such as genomic selection (GS) and genome-wide association study (GWAS). As new genotyping platforms are developed, data from different genotyping platforms must be combined. However, the potential use of combined data for GS and GWAS has not yet been clarified. In this study, the accuracy of genomic prediction (GP) and the detection power of GWAS increased for most fruit quality traits of apples when using combined data from different genotyping systems, Illumina Infinium single-nucleotide polymorphism array and genotyping by random amplicon sequencing-direct (GRAS-Di) systems. In addition, the GP model, which considered the inbreeding effect, further improved the accuracy of the seven fruit traits. Runs of homozygosity (ROH) islands overlapped with the significantly associated regions detected by the GWAS for several fruit traits. Breeders may have exploited these regions to select promising apples by breeders, increasing homozygosity. These results suggest that combining genotypic data from different genotyping platforms benefits the GS and GWAS of fruit quality traits in apples. Information on inbreeding could be beneficial for improving the accuracy of GS for fruit traits of apples; however, further analysis is required to elucidate the relationship between the fruit traits and inbreeding depression (e.g. decreased vigor).

Deleterious phenotypes in wild Arabidopsis arenosa populations are common and linked to runs of homozygosity

In this study, we aimed to systematically assess the frequency at which potentially deleterious phenotypes appear in natural populations of the outcrossing model plant Arabidopsis arenosa, and to establish their underlying genetics. For this purpose, we collected seeds from wild A. arenosa populations and screened over 2,500 plants for unusual phenotypes in the greenhouse. We repeatedly found plants with obvious phenotypic defects, such as small stature and necrotic or chlorotic leaves, among first-generation progeny of wild A. arenosa plants. Such abnormal plants were present in about 10% of maternal sibships, with multiple plants with similar phenotypes in each of these sibships, pointing to a genetic basis of the observed defects. A combination of transcriptome profiling, linkage mapping and genome-wide runs of homozygosity patterns using a newly assembled reference genome indicated a range of underlying genetic architectures associated with phenotypic abnormalities. This included evidence for homozygosity of certain genomic regions, consistent with alleles that are identical by descent being responsible for these defects. Our observations suggest that deleterious alleles with different genetic architectures are segregating at appreciable frequencies in wild A. arenosa populations.

Genetic diversity, population structure, and genome-wide association analysis of ginkgo cultivars

Ginkgo biloba is an economically valuable tree worldwide. The species has nearly become extinct during the Quaternary, which has likely resulted in reduction of its genetic variability. The genetic variability is now conserved in few natural populations in China and a number of cultivars that are, however, derived from a few ancient trees, helping the species survive in China through medieval times. Despite the recent interest in ginkgo, however, detailed knowledge of its genetic diversity, conserved in cultivated trees and cultivars, has remained poor. This limits efficient conservation of its diversity as well as efficient use of the existing germplasm resources. Here we performed genotyping-by-sequencing (GBS) on 102 cultivated germplasms of ginkgo collected to explore their genetic structure, kinship, and inbreeding prediction. For the first time in ginkgo, a genome-wide association analysis study (GWAS) was used to attempt gene mapping of seed traits. The results showed that most of the germplasms did not show any obvious genetic relationship. The size of the ginkgo germplasm population expanded significantly around 1500 years ago during the Sui and Tang dynasties. Classification of seed cultivars based on a phylogenetic perspective does not support the current classification criteria based on phenotype. Twenty-four candidate genes were localized after performing GWAS on the seed traits. Overall, this study reveals the genetic basis of ginkgo seed traits and provides insights into its cultivation history. These findings will facilitate the conservation and utilization of the domesticated germplasms of this living fossil plant.

Selection, recombination and population history effects on runs of homozygosity (ROH) in wild red deer (Cervus elaphus)

The distribution of runs of homozygosity (ROH) may be shaped by a number of interacting processes such as selection, recombination and population history, but little is known about the importance of these mechanisms in shaping ROH in wild populations. We combined an empirical dataset of >3000 red deer genotyped at >35,000 genome-wide autosomal SNPs and evolutionary simulations to investigate the influence of each of these factors on ROH. We assessed ROH in a focal and comparison population to investigate the effect of population history. We investigated the role of recombination using both a physical map and a genetic linkage map to search for ROH. We found differences in ROH distribution between both populations and map types indicating that population history and local recombination rate have an effect on ROH. Finally, we ran forward genetic simulations with varying population histories, recombination rates and levels of selection, allowing us to further interpret our empirical data. These simulations showed that population history has a greater effect on ROH distribution than either recombination or selection. We further show that selection can cause genomic regions where ROH is common, only when the effective population size (Ne) is large or selection is particularly strong. In populations having undergone a population bottleneck, genetic drift can outweigh the effect of selection. Overall, we conclude that in this population, genetic drift resulting from a historical population bottleneck is most likely to have resulted in the observed ROH distribution, with selection possibly playing a minor role.

Homozygosity mapping identified loci and candidate genes responsible for freezing tolerance in Camelina sativa

Homozygosity mapping is an effective tool for detecting genomic regions responsible for a given trait when the phenotype is controlled by a limited number of dominant or co-dominant loci. Freezing tolerance is a major attribute in agricultural crops such as camelina. Previous studies indicated that freezing tolerance differences between a tolerant (Joelle) and susceptible (CO46) variety of camelina were controlled by a small number of dominant or co-dominant genes. We performed whole genome homozygosity mapping to identify markers and candidate genes responsible for freezing tolerance difference between these two genotypes. A total of 28 F3 RILs were sequenced to ∼30× coverage, and parental lines were sequenced to >30-40× coverage with Pacific Biosciences high fidelity technology and 60× coverage using Illumina whole genome sequencing. Overall, about 126k homozygous single nucleotide polymorphism markers were identified that differentiate both parents. Moreover, 617 markers were also homozygous in F3 families fixed for freezing tolerance/susceptibility. All these markers mapped to two contigs forming a contiguous stretch of chromosome 11. The homozygosity mapping detected 9 homozygous blocks among the selected markers and 22 candidate genes with strong similarity to regions in or near the homozygous blocks. Two such genes were differentially expressed during cold acclimation in camelina. The largest block contained a cold-regulated plant thionin and a putative rotamase cyclophilin 2 gene previously associated with freezing resistance in arabidopsis (Arabidopsis thaliana). The second largest block contains several cysteine-rich RLK genes and a cold-regulated receptor serine/threonine kinase gene. We hypothesize that one or more of these genes may be primarily responsible for freezing tolerance differences in camelina varieties.

Combined transcriptomic and proteomic analysis reveals multiple pathways involved in self-pollen tube development and the potential roles of FviYABBY1 in self-incompatibility in Fragaria viridis

Fragaria viridis exhibits S-RNase-based gametophytic self-incompatibility, in which S-RNase is the major factor inhibiting pollen tube growth. However, the pathways involved in and the immediate causes of the inhibition of pollen tube growth remain unknown. Here, interactive RNA sequencing and proteome analysis revealed changes in the transcriptomic and proteomic profiles of F. viridis styles harvested at 0 and 24 h after self-pollination. A total of 2,181 differentially expressed genes and 200 differentially abundant proteins were identified during the pollen development stage of self-pollination. Differentially expressed genes and differentially abundant proteins associated with self-incompatible pollination were further mined, and multiple pathways were found to be involved. Interestingly, the expression pattern of the transcription factor FviYABBY1, which is linked to polar growth, differed from those of other genes within the same family. Specifically, FviYABBY1 expression was extremely high in pollen, and its expression trend in self-pollinated styles was consistent with that of S-RNase. Furthermore, FviYABBY1 interacted with S-RNase in a non-S haplotype way. Therefore, FviYABBY1 affects the expression of polar growth-related genes in self-pollen tubes and is positively regulated by S-RNase.