An identical-by-descent segment harbors a 12-bp insertion determining fruit softening during domestication and speciation in Pyrus
BMC Biology volume 20, Article number: 215 (2022)
Although the wild relatives of pear originated in southwest China, this fruit crop was independently domesticated and improved in Asia and Europe, and there are major phenotypic differences (e.g., maturity and fruit firmness) between Asian and European pears.
In this study, we examined the genomes of 113 diverse pear accessions using an identity-by-descent (IBD) approach to investigate how historical gene flow has shaped fruit firmness traits in Asian and European pears. We found a 3-Mbp IBD-enriched region (IBD-ER) that has undergone “convergent domestication” in both the Asian and European pear lineages, and a genome-wide association study (GWAS) of fruit firmness phenotypes strongly implicated the TRANSLOCON AT THE INNER CHLOROPLAST ENVELOPE55 (TIC55) locus within this 3-Mbp IBD-ER. Furthermore, we identified a tandem duplication that includes a 12-bp insertion located in the first exon of TIC55 that is uniquely present in Asian pears, and expression analysis showed that the pear TIC55 gene is highly expressed in Asian pear, while it is weakly or not expressed in European pear; this could contribute to the differences in fruit firmness between Asian and European pear fruits.
Our findings provide insights into how pear fruit softening has been impacted during domestication, and we identified candidate genes associated with fruit softening that can contribute to the breeding and improvement of pear and other fruit crops.
Pears (species in the genus Pyrus L.), which belong to the subfamily Pomoideae in the botanical family Rosaceae, comprise the third most important fruit crop in the world, ranking just after apple with a yield of ~ 16.1 million tons (2020, FAOSTAT) per year . There are at least 22 widely recognized pear species , and more than 5000 pear accessions are currently distributed worldwide [3, 4]. In 2018, a pear domestication model proposed that the pear originated in southwest China, was disseminated through central Asia, and eventually underwent independent domestication in the Asian and European pear populations .
Beneficial reciprocal gene flow via interspecific hybridization is an important way to improve phenotypic traits during domestication [5, 6]. With the release of more genomic DNA resequencing data, gene flow has been widely investigated in many plant species. In rice, introgression analysis has determined that most rice types exhibit distinct degrees of admixture, with many accessions within a population conserving particular identity-by-descent (IBD) tracts due to artificial selection . In maize, long conserved IBD tracts have been detected between tropical and temperate inbred lines, which may be attributed to gene flow between tropical and temperate maize that must have occurred during domestication and genetic improvement . Pear, a typical gametophytic self-incompatible plant, is an obligate out-crosser . Therefore, it is expected that there is a more likely gene flow between wild and cultivated pear populations. During domestication, pears underwent dramatic morphological and physiological changes, and these differences are considered to be viable candidate targets for investigating pear domestication and improvement.
Firmness is one of the most important fruit traits because it determines the shelf-life and commercial value of the fruit. Compared with wild pears, the firmness of cultivated pear fruits is markedly lower; moreover, the firmness of cultivated Asian pear fruits is higher than that of cultivated European pears (Fig. 1A, B; Additional file 1: Table S1). This suggests that fruit firmness is likely to be one of the typical phenotypic traits that must have been involved in the independent domestication of Asian and European pears. Fruit softening is closely related to the degradation of cell walls that are mainly composed of pectin, cellulose, and hemicellulose . The degradation of cell walls involves the interaction of many cell wall-related enzymes, including polygalacturonase (PG) , pectin methylesterase (PME) , β-galactosidase (β-GAL) , and pectate lyases (PLs) . Ethylene (ET) is the main hormone that promotes fruit ripening and senescence, and it plays an important role in the fruit softening process. Therefore, ET biosynthesis and signal transduction are important factors that affect fruit softening. For example, the ACS1 and ACO1 genes have been shown to regulate ET biosynthesis in apple fruit during ripening [15, 16]. During the signal transduction process, ethylene response factor (ERF) family members are widely reported to be involved in regulating downstream ET-responsive genes, such as SR1 , LeERF2 , MaERF9 and MaERF11 , and ERF4 . SR1, a member of the tomato CALMODULIN BINDING TRANSCRIPTION ACTIVATOR (CAMTA) family, is involved in ET signal transduction by regulating the expression level of ETHYLENE INSENSITIVE3 (EIN3) . In addition, fruit softening is one of the important main characteristics of senescence, and genes involved in plant senescence have been widely studied. In Arabidopsis, knocking out the expression of CAMTA3, a gene that encodes a calmodulin-binding protein, contributes to early senescence . Moreover, it has been reported that the TRANSLOCON AT THE INNER CHLOROPLAST ENVELOPE55 (TIC55) gene plays a vital role in chlorophyll breakdown and is associated with a translocon for proteins via an inner chloroplast membrane .
In this study, resequencing data from 113 pear accessions was analyzed to explore the potential role(s) of gene flow and the phenotypic effect(s) on fruit firmness during the domestication of pear (Additional file 2: Table S2). Our findings included (1) a landscape of the incidence of gene flow in the genomes of Asian and European populations and (2) fruit firmness-related phenotypic effects resulting from gene flow during domestication must have accompanied the divergence of the Asian and European pear lineages.
Genome-wide identification of putative identity-by-descent (IBD) tracts in 113 Pyrus accessions
In our study, we used genome sequencing data from 113 pear accessions collected from previous research; the accessions came from four pear populations and included 31 cultivated Asian, 32 wild Asian, 25 cultivated European, and 25 wild European accessions. We cataloged IBD tracts at the whole-genome level to explore the gene flow within/across these four populations. A total of 3,320,796 IBD tracts were identified in pairwise comparisons of the 113 pear accessions. We divided the identified IBD tracts from all 10 pairwise comparisons into two categories: intra-population tracts and inter-population tracts. Among the intra-population tracts, the cultivated Asian-cultivated Asian (CA-CA) comparison showed the highest extent of IBD sharing (2411), and the wild European-wild European (WE-WE) comparison showed the lowest extent of IBD sharing (1367). These results suggest that the cultivated Asian pear population has experienced the highest extent of intra-population genetic exchange among the four populations. For the inter-population comparisons, we found that cultivated Asian and wild Asian showed slightly higher levels of IBD sharing compared to that between cultivated European and wild European pears (Fig. 2A), indicating a relatively higher level of gene flow between wild Asian and cultivated Asian pears. Notably, the IBD tract data corresponds with our finding from an FST analysis showing that the degree of divergence in cultivated Asian and wild Asian pear (FST = 0.021) is reduced compared to cultivated European and wild European pear (FST = 0.035) (Fig. 2B).
We also estimated the length distribution of the IBD tracts. Most IBD tracts had lengths between 1500 and 3500 bp in all 10 pairwise comparisons (Additional file 3: Table S3 and Additional file 14: Fig. S1). The number of IBD tracts decreased sharply with increasing tract length for both the intra-population and inter-population comparisons; this can likely be explained by the fact that longer IBD tracts are more likely to be disrupted by genetic recombination [6, 22]. A heatmap and principal component analysis (PCA) based on the IBD tract data showed that the 113 pear accessions were divided into two groups that corresponded to Asian and European populations (Fig. 2C, D). Except for three Pyrus sinkiangensis and a single Pyrus armeniacaefolia accessions, Asian pears in this PCA model were distinct from the European pears; this separation likely reflects the long-term geographical separation of Asian and European pears.
To investigate which accessions have contributed the highest proportion of genetic background to the pear gene pool, we generated an IBD sharing network and used it to explore the relationships among the pear accessions. This is shown in Fig. 2E, where the dot sizes correspond to the number(s) of shared IBD tracts. We found that the Pyc_py10 (“Wasekozo”) accession of Pyrus pyrifolia contained the highest level of IBD sharing (140,546), reflecting the known fact that “Wasekozo” has been widely used as a hybrid parent in Asian pear breeding programs . For example, the Japanese pear variety “Kosui,” one of the most widely grown Asian pear cultivars, was derived from a cross between “Wasekozo” and “Kikusui” (MAFF, https://www.maff.go.jp/j/tokei/index.html).
Identification of IBD-enriched chromosomal regions
We also measured the number of IBD tracts in 500-kb windows across the 17 pear chromosomes (Additional file 4: Table S4). Extensive but uneven distributions of IBD tracts were identified in all the 10 pairwise comparisons (Fig. 3). At the chromosome level, the fewest IBD tracts were detected for chromosome 4 (2.87%) in the comparison between cultivated Asian and wild Asian (CA-WA) pears, and the fewest IBD tracts from the cultivated European-wild European (CE-WE) pear comparison were for chromosome 7 (3.09%). It was also interesting to note that chromosome 15 had the most significant number of IBD tracts for both the CA-WA (12.52%) and CE-WE (11.99%) comparisons (Additional file 15: Fig. S2A and Fig. S2B).
We next defined the top 5% of windows (39) with the most significant number of IBD segments as IBD-enriched regions (IBD-ERs) (Additional file 5: Table S5). In the CA-WA comparison, a total of 1725 genes were identified in IBD-ERs (Asian IBD-ERs), while a slightly lower number of genes (1354) was identified within IBD-ERs in the CE-WE comparisons (European IBD-ERs) (Additional file 6: Table S6 and Additional file 7: Table S7). A Gene Ontology (GO) term analysis revealed that genes with functional annotations related to multiple aspects of carbohydrate metabolism were enriched in the Asian pear IBD-ERs (p-value < 0.05, Additional file 15: Fig. S2C and Fig. S2D). Of particular note is a 3-Mbp region on chromosome 15 containing 446 protein-encoding genes that was detected as a common IBD-ER in both the Asian and European pear populations (Fig. 3). Notably, a previous study reported a QTL (MEST050) related to harvest time that overlaps this 3-Mbp region , and a KEGG pathway analysis we conducted indicated that genes in this region were enriched for annotated functions related to fruit ripening (e.g., “Zeatin biosynthesis”) (Additional file 16: Fig. S3). These findings suggest that fruit maturity-related traits have been of particular interest to humans and were selected during the domestication processes for both Asian and European pears.
A 3-Mbp IBD-ER shared between Asian and European populations experienced positive selection during pear domestication
We used two methods to identify selective sweeps among the IBD-ERs to explore which, if any, have undergone artificial selection. First, we performed a cross-population likelihood ratio analysis of the whole genomes between wild and cultivated populations using the XP-CLR program . We considered regions harboring XP-CLR values in the top 20% (Asian: XP-CLR value ≥ 5.17; European, XP-CLR value ≥ 2.33) as putative domestication regions exhibiting selection signals (Fig. 4A–C) . Approximately 17.9% and 20.5% of the IBD-ERs showed signals of positive selection among Asian and European pears, respectively. It is notable that only two IBD-ERs (500-kb windows; Chr15: 6,000,001–6,500,000, Chr15: 7,500,001–8,000,000) displayed selection signatures in both the Asian and European populations, and these two selection signals were both within the aforementioned 3-Mbp IBD-ER on chromosome 15 (Chr15: 5,000,000–8,000,000). This finding supports our earlier conclusion that the 3-Mbp IBD-ER did indeed undergo positive selection during the domestication processes for both Asian and European pears.
Second, we identified the population differentiation level (FST) at the whole-genome level between wild and cultivated populations using the VCFtools program. Briefly, FST values for a single genomic region between 0.05 and 0.15 indicate that two populations have a moderate degree of genetic differentiation; an FST value > 0.15 suggests that two populations are separated. Our analysis showed that 17.9% and 15.4% of the IBD-ERs had FST values > 0.15 in the Asian and European populations, respectively (Fig. 4D, E). Conspicuously, the only IBD-ER (Chr15: 7,500,001–8,000,000) that had an FST value > 0.15 in both the Asian and European populations was located within the 3-Mbp IBD-ER on chromosome 15. Thus, both of the methods we used to evaluate selective sweeps support the conclusion that the 3-Mbp IBD-ER on chromosome 15 has experienced strong selection. In addition, cross-population extended haplotype homozygosity XP-EHH analysis also supported this conclusion. There were a total of 448 and 91 selection signals in the 3-Mb IBD-ER in the wild Asian-vs-cultivated Asian (WA-CA) and wild European-vs-cultivated European (WE-CE) comparisons, respectively (Additional file 17: Fig. S4). The detection of positive selection for this common region in the two independently domesticated pear populations suggests that we have detected an example of “convergent domestication” in our study [27, 28]; this is apparently similar to the previously reported case in Brassica where 15 convergent selection loci were identified during the independent domestications of Brassica rapa and Brassica oleracea .
GWAS for fruit firmness identifies candidate causal genes within the 3-Mbp IBD-ER shared by Asian and European pears
We also used the 200 K AXIOM® PyrSNP array to genotype a germplasm diversity panel specifically assembled for cultivated pears (comprising 164 pear cultivars) that we used in a GWAS examining fruit traits (Additional file 8: Table S8). Very briefly, 29,269 high-quality SNPs were retained after filtering (Fig. 5A and Additional file 9: Table S9 and Additional file 10: Table S10), among which 3566 were both positioned within protein-coding regions and resulted in predicted changes in amino acids or altered start/stop properties for transcription (Fig. 5B and Additional file 11: Table S11). Population structure, PCA, and neighbor-joining phylogenetic tree analyses all supported dividing these 164 cultivated pear accessions into two groups, with group II comprising mostly Asian pears and group I comprising mostly European pears (Fig. 5C–E). Of note, and consistent with both results from our analysis of the four pear populations and conclusions from a previous study , nucleotide diversity (π) and linkage disequilibrium (LD) analysis strongly support the notion that Asian pears are substantially more genetically diverse than European pears.
Before conducting GWAS, phenotypic data for fruit firmness was collected for all 164 accessions from the cultivated pear diversity panel analyzed with the 200 K AXIOM® PyrSNP array. The fruit firmness data demonstrated a normal distribution (Shapiro–Wilk normality test, W = 0.988, p-value = 0.180) (Fig. 6A, B) and was used for genome-wide association analysis based on a mixed linear model (MLM). The GWAS detected two loci, AX-179248435 on chromosome 13 and AX-179263231 on chromosome 15, that were significantly associated with fruit firmness (Fig. 6C and Additional file 18: Fig. S5). Based on a previous genome-wide association study in pear, we used 100-kb windows to detect candidate genes . A total of 32 candidate genes were identified from the flanking regions of the two fruit firmness-associated sites: 12 genes from the site on chromosome 13 and 20 genes on chromosome 15 (Additional file 12: Table S12).
Two of these genes are particularly notable: first, the pear ortholog of the previously characterized calmodulin-binding protein gene CAMTA3 (Pbr035644.2) had haplotypes that were divergent between Asian and European pears (Figs. 6C and 7A, B); second, the cultivated Asian pears carry the AA allele and have firmer fruits than cultivated European pears, which have the GG allele. The SNP results in a predicted change from asparagine to aspartic in the CAMTA3 protein (Fig. 7C). We also found a 9-bp insertion in the CDS of CAMTA3 in the Asian pear population that is absent in the European pear population. This indel caused a 3-amino acid insertion (A-G-L: alanine, glycine, and lysine) in the CAMTA3 protein in Asian pears (Additional file 19: Fig. S6). Previous studies in Arabidopsis have demonstrated that disruption of CAMTA3 results in early senescence [17, 31]. In general, fruit firmness slowly decreases during fruit ripening (Additional file 20: Fig. S7). These results indicate that CAMTA3 might regulate fruit firmness by mediating maturity and senescence in pears.
The second notable gene is Pbr019967.1 (Fig. 6C), which is located within the aforementioned 3-Mbp IBD-ER that is common to both Asian and European pears. A BLAST search identified Pbr019967.1 as the likely pear ortholog of Arabidopsis TIC55. To further determine whether TIC55 contributes to the differences in firmness between Asian and European pears, we investigated the gene structure of TIC55 and its genetic divergence between Asian and European pears. The PbTIC55 gene contains three exons and two introns (Fig. 7D) and encodes a predicted protein of 575 amino acid residues with a calculated molecular mass of 64.21 kDa and an isoelectric point of 8.89. Alignments of the coding region sequences (CDS) of the PbTIC55 gene from European and Asian pear accessions revealed the presence of a nucleotide polymorphism: exon 1 of PbTIC55 from the cultivated European accessions is 12 bp shorter compared to the Asian accessions, and this was confirmed by Sanger sequencing of the gene from seven Asian and six European pears with significant firmness differences (Fig. 7E). The length difference is due to the duplication of a 12-bp sequence within exon 1, and all seven of the Asian pear accessions were found to carry this 12-bp tandem repeat (Fig. 7F) (Additional file 21: Fig. S8).
To further investigate the effect(s) of this 12-bp tandem repeat on TIC55, we compared the expression level of TIC55 between the Asian pear group (nine cultivated varieties) and the European pear group (six cultivated varieties) using quantitative real-time PCR (qRT-PCR) analysis. The results indicated that the expression of TIC55 is significantly higher in the Asian pear group than in the European pear group (T-test, p < 0.001) (Fig. 7G). Asian pears are known to have a relatively slower rate of fruit softening compared to European pears. The gene expression results further supported the idea that 12-bp tandem repeat is associated with fruit softening in Asian pear fruits.
In general, fruit firmness slowly decreases during the fruit development process (Additional file 20: Fig. S7). To verify the function of the TIC55 gene, we investigated the expression pattern of TIC55 during pear fruit development in both Asian and European pears using transcriptome data. An increase in the expression level of TIC55 was found during the development of Asian pear fruits, while TIC55 was expressed at very low levels or was not expressed at all during the development of European pear fruits (Fig. 7H). We next selected two pear accessions to quantify the relative expression level of TIC55 during fruit development using qRT-PCR analysis (Additional file 22: Fig. S9). The qRT-PCR results were consistent with the results of transcriptome data analysis. In general, fruit firmness slowly decreases during fruit development, and European pears soften faster than Asian pears. These results indicated that TIC55 expression acts to repress fruit softening.
We examined the genomes of 113 pear accessions to estimate the demographic history and identify the landscape of gene flow during the independent domestications of Asian and European pears. A 3-Mbp IBD-ER showed strong selective signals that occurred during the domestication of both Asian and European pears, and this IBD-ER harbors the TIC55 locus, which is significantly associated with fruit softening. We identified a 12-bp insertion mutation in the first exon of TIC55 in Asian pear accessions that potentially explains the differences in fruit softening between Asian and European pears. This study provides novel insights into our understanding of the genetic basis of fruit softening during pear domestication. Identification of genes associated with fruit firmness will facilitate future molecular studies and will also aid in the genetic improvement and breeding efforts for pears.
In maize, it has been reported that the lengths of IBD tracts range between 3 and 4 Mbp [8, 32]. However, in this study, we found that most of the IBD tracts were between 1500 and 3500 bp in length, which is substantially shorter than the IBD tracts detected in maize. This finding may be attributed to self-incompatibility in pears, because a higher incidence of genetic recombination can break up long IBD tracts . Compared to IBD-ERs, chromosomal segments lacking obvious IBD regions exhibit higher genetic diversity and typically undergo relatively weaker selective pressure . In pear, the S-locus for self-incompatibility is located at the end of chromosome 17 and shows a high level of diversity and a rapid evolution rate . In our study, we found no shared IBD tracts at the S-locus between wild and cultivated pears (X2, p-value < 2.2e − 16). This finding is consistent with the observations in maize showing that the genomic regions near the known cross-incompatibility loci are likely to be particularly resistant to gene flow .
Hormone levels are an important factor that is associated with fruit firmness and fruit size. In our study, the KEGG pathway analysis of genes located in the 3-Mb IBD-ER indicated that these genes are related to zeatin biosynthesis. Zeatin is a natural cytokinin hormone that is produced in plants . Cytokinin levels are decreased during fruit ripening, and elevated levels of cytokinin contribute to delayed ripening [35, 36]. During the fruit ripening process, ethylene is involved in a range of physiological and biochemical changes including chlorophyll degradation, fruit softening, and the production of volatiles . The effect of ethylene on fruit firmness is mainly reflected in ethylene’s effect on pectin metabolism enzymes in the cell wall [38, 39]. In addition, Xin et al. reported that ethylene is required for cucumber fruit elongation, which further affects cucumber fruit size . ABA regulates fruit development by interacting with other hormones. For example, ABA induces cell enlargement during tomato fruit growth by decreasing the synthesis of ethylene, which further influences fruit size in tomatoes . A decade of studies have shown that ABA is an important promotor of fruit softening that mediates transcriptional changes in genes that encode enzymes related to cell wall degradation [42, 43]. For example, the expression level of FaRGlyase1, a gene that encodes the enzyme rhamnogalacturonate lyase, was shown to be positively regulated by ABA and negatively regulated by auxins .
Fruit firmness is a critical quality trait because it determines the quality grading of fruit and strongly affects shelf life . Fruits of cultivated pears are less firm than fruits of wild pear relatives, although it is important to note that excessive fruit softening is deleterious for long-term storage. While firmness has been widely investigated in many fruits, including sweet cherry, strawberry, and peach [46,47,48,49], knowledge of the genetic basis of the differential fruit firmness in pears is rather limited. The process of fruit ripening is usually accompanied by fruit softening (Additional file 20: Fig. S7) [50, 51]. It has been determined that CAMTA3 plays a critical role in ethylene signaling by directly regulating the expression of both NON-RACE-SPECIFIC DISEASE RESISTANCE1 (NDR1) and EIN3 in Arabidopsis [17, 31]. In this study, we have identified another gene, TIC55, that is found within a 3-Mbp IBD-ER in the pear genome that has undergone strong positive selection during domestication. Previous studies reported that TIC55 functions in chlorophyll breakdown during leaf senescence and fruit ripening [52, 53].
Asian and European pears display extensive phenotypic differences, such as their stone cell content, fruit firmness, and flavor , and these differences reflect their independent domestications. Of these traits, probably the most apparent trait difference between these two populations is that of fruit firmness. Compared to European pears, the fruits of Asian pears are significantly firmer, yet the underlying genetic mechanisms contributing to this trait are unclear. Based on genotypic divergence and PCR amplification verification, we have found a 12-bp insertion mutation in the first exon of the TIC55 gene that is unique to Asian pears. This mutation is likely to be involved in the observed differences in fruit softening phenotypes and could potentially explain the differences in the trait between Asian and European pears (Fig. 8).
The objectives of our study were to assess the identity-by-descent (IBD) landscape of 113 diverse pear accessions to investigate how historical gene flow has shaped fruit traits in Asian and European pears. A total of 3,320,796 IBD tracts were identified in pairwise comparisons between the 113 pear accessions. Notably, a 3-Mbp region that is shared between Asian IBD-enriched regions (IBD-ERs) and European IBD-ERs was identified on chromosome 15, and both XP-CLR and FST analysis strongly supported the hypothesis that the 3-Mbp IBD-ERs had experienced selection during pear domestication. One fruit firmness-associated gene, TIC55, was found to be located within the 3-Mbp IBD-ER that is shared between Asian and European pears. Furthermore, we identified a tandem repeat that results from a 12-bp duplication mutation found in the first exon of TIC55 that is unique to Asian pears, which have a relatively slower rate of fruit softening than European pears. This study investigated the historical gene flow and highlighted its phenotypic effects on pear fruit softening, with important implications for understanding the independent domestications of Asian and European pears.
Acquisition of genome sequences
We used genomic resequencing data for 113 pear accessions that represent the five major cultivated pear groups and the most important wild Pyrus species in the world that was generated in our previous study. The resequenced pear accessions included 63 Asian (31 cultivated and 32 wild) and 50 European (25 cultivated and 25 wild) accessions .
SNP source and pairwise IBD detection
SNP calling was performed following previously described protocols . A total of 18,302,883 SNPs were identified in all 113 accession genomes, and these were used as inputs for IBD identification. First, all 113 accessions were phased in Beagle (version 4.0). Second, the IBD tracts were extracted with the Beagle IBD function, as described by Browning and Browning , with the LOD score set to > 3.0. To profile the frequency of IBD tracts along each chromosome, the chromosomes were divided into 500-kb windows, and the numbers of detected IBD tracts for each group in each window were counted. These results were illustrated using the Circos tool (version 0.69) .
FST analysis, XP-CLR analysis, and XP-EHH analysis
The package VCFtools (version 0.1.17) was used to compute pairwise FST values between the wild and cultivated pear populations . The XP-CLR software was used to estimate XP-CLR values between the wild and cultivated populations . The regions harboring XP-CLR values in the top 20% (Asian: XP-CLR value ≥ 5.17; European, XP-CLR value ≥ 2.33) were deemed to be the selection domestication region. The FST values and XP-CLR values were calculated using a sliding window of 500 kb and a step length of 50 kb. The regions harboring FST values > 1.5 were considered to be selection domestication regions. We also performed comparisons between the wild and cultivated populations in both Asian and European pears using cross-population extended haplotype homozygosity (XP-EHH). For XP-EHH analysis, we first phased and imputed genotyping data using the Beagle software. We then used the Selscan software to measure the XP-EHH scores with the default parameters. The top 1% of XP-EHH signals were considered to be selection signals.
Genotyping using a 200 K AXIOM® PyrSNP array
A total of 164 diverse pear accessions were selected and genotyped using the 200 K AXIOM® PyrSNP array as previously described . Twelve of the pear accessions were shared between the 113 resequenced pear accessions and the 164 accessions used in the GWAS in this study. A kit (Tiangen Biotech Co., Ltd., Beijing, China) was used to extract genomic DNA from frozen fruit tissues from each of the 164 pear accessions. All raw hybridization intensity data were analyzed using the Axiom Analysis Suite (v.1.1.1) with a diploid threshold configuration (Affymetrix, Santa Clara, CA). Samples with call rates (CR) > 97 and dish quality control (DQC) values > 0.82 were retained and subjected to further analysis.
Based on the reliability of the SNPs, SNPs from the results obtained from the Axiom Analysis Suite were classified into the following six categories using the “SNPolisher” R package: “PolyHighResolution” (PHR), “NoMinorHom” (NMH), “OTV,” “MonoHighResolution,” “CallRateBelowThreshold,” and “other.” Those SNPs classified as “other” had the lowest quality and were excluded from further analysis.
A genome-wide association study (GWAS) for pear fruit firmness
The pear fruit firmness trait was measured using a Brookfield CT3 Texture Analyzer (kg/cm2) (AMETEK Brookfield, Middleboro, MA). Three pear accessions were randomly selected as three biological replicates. Three surface points along each pear fruit were randomly selected to measure fruit firmness, and the average value of these measurements was designated as the fruit firmness of each pear fruit analyzed.
The Genomic Association and Prediction Integrated Tool (GAPIT) package was used to conduct GWAS using the mixed linear model (MLM) . The GWAS significance threshold was set to 3.41 × 10−5 (1/n, where n is the number of SNPs used in GWAS). QQman (R package) was used to create Manhattan plots .
Linkage disequilibrium (LD) analyses were performed to search for candidate genes in regions 50 kb upstream and 50 kb downstream of significant SNPs using plink (v1.07) . LDheatmap (R package) was used to display linkage disequilibrium values of SNPs within a 100-kb candidate region .
Quantitative real-time PCR analysis (qRT-PCR)
We selected nine Asian and seven European pear accessions in which to assay the expression pattern of PbTIC55 using qRT-PCR. The Asian pear group comprised the cultivars “Hongxiao,” “Oushuu,” “Saigyoku,” “Qiuxiangshui,” “Hehuali,” “Mansoo,” “Lelingshali,” “Balengli,” and “Meirensu.” The European pear group included the cultivars “La France,” “Jana,” “Bartlett,” “Bosco,” “Clarke,” and “Comice.”
Total RNA was extracted using the Plant Total RNA isolation Plus kit (Foregene; http://www.foregene.com), and first-strand cDNA was synthesized using the EasyScript® First-Strand cDNA Synthesis SuperMix Kit (Transgene; https://www.transgen.com.cn/rt_pcr.html). Primers were designed to amplify genes using online software available from NCBI (National Center for Biotechnology Information) (https://www.ncbi.nlm.nih.gov/tools/primerblast/) (Additional file 13: Table S13). As described previously , LightCycler 480 SYBR GREEN I Master (Roche) was used for the qRT-PCR analyses. Each 10 μl reaction contained 150 ng of template cDNA, 0.5 μM of each primer, and 5 μl of LightCycler 480 SYBR GREEN I Master. All reactions were performed in 96-well plates with three replicates for each cDNA sample. The qRT-PCR amplification conditions were as follows: 3 min at 95 °C followed by 45 cycles of 95 °C for 3 s, 60 °C for 10 s, and 72 °C for 30 s. Fluorescence signal data was collected at 60 °C. The Pyrus TUB gene was used as the internal control. The average threshold cycle (Ct) of each cDNA sample was calculated. Relative gene expression levels were calculated using the 2−ΔΔCt method as described .
A Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis was performed using KOBAS 3.0 .
All statistical tests were performed using R language basic packages (version 3.5.1) (https://www.R-project.org/).
Genome-wide association analysis
Mixed linear model
Kyoto Encyclopedia of Genes and Genomes
Single nucleotide polymorphism
Minor allele frequency
Genomic Association and Prediction Integrated Tool
- CAMTA3 :
Calmodulin Binding Transcription Activator3
- TIC55 :
Translocon At The Inner Chloroplast Envelope55
Mirabdulbaghi M. Investigations on determination of nutritional status of pear trees according to a new index - Deviation From Optimum Percentage (DOP). Cercet Agron Mold. 2015;47(4):83–98.
Janick J, Moore JN. Fruit breeding. New York: Wiley; 1996.
Wu J, Wang Y, Xu J, Korban SS, Fei Z, Tao S, Ming R, Tai S, Khan AM, Postman JD, et al. Diversification and independent domestication of Asian and European pears. Genome Biol. 2018;19(1):77.
Li J, Zhang M, Li X, Khan A, Kumar S, Allan AC, Lin-Wang K, Espley RV, Wang C, Wang R, et al. Pear genetics: recent advances, new prospects, and a roadmap for the future. Hortic Res. 2022;9:uhab040.
Ma Y, Wang J, Hu Q, Li J, Sun Y, Zhang L, Abbott RJ, Liu J, Mao K. Ancient introgression drives adaptation to cooler and drier mountain habitats in a cypress species complex. Commun Biol. 2019;2(1):213.
Martin SH, Jiggins CD. Interpreting the genomic landscape of introgression. Curr Opin Genet Dev. 2017;47:69–74.
Zhao K, Wright M, Kimball J, Eizenga G, McClung A, Kovach M, Tyagi W, Ali ML, Tung CW, Reynolds A, et al. Genomic diversity and introgression in O. sativa reveal the impact of domestication and breeding on the rice genome. PloS One. 2010;5(5):e10780.
Li X, Jian Y, Xie C, Wu J, Xu Y, Zou C. Fast diffusion of domesticated maize to temperate zones. Sci Rep. 2017;7(1):2077.
Wu J, Wang Z, Shi Z, Shu Z, Zhang S. The genome of the pear (Pyrus bretschneideri Rehd.). Genome Res. 2013;23(2):396–408.
Julian JD, Zabotina OA. Xyloglucan biosynthesis: from genes to proteins and their functions. Front Plant Sci. 2022;13:920494.
Bird CR, Smith CJS, Ray JA, Moureau P, Bevan MW, Bird AS, Hughes S, Morris PC, Grierson D, Schuch W. The tomato polygalacturonase gene and ripening-specific expression in transgenic plants. Plant Mol Biol. 1988;11(5):651–62.
Sanchez N, Gutiérrez-López GF, Cáez-Ramírez G. Correlation among PME activity, viscoelastic, and structural parameters for Carica papaya edible tissue along ripening. J Food Sci. 2020;85(6):1805–14.
Smith DL, Abbott JA, Gross KC. Down-regulation of tomato β-galactosidase 4 results in decreased fruit softening. Plant Physiol. 2002;129(4):1755–62.
Uluisik S, Chapman NH, Smith R, Poole M, Adams G, Gillis RB, Besong TMD, Sheldon J, Stiegelmeyer S, Perez L, et al. Genetic improvement of tomato by targeted control of fruit softening. Nat Biotechnol. 2016;34(9):950–2.
Harada T, Sunako T, Wakasa Y, Soejima J, Satoh T, Niizeki M. An allele of the 1-aminocyclopropane-1-carboxylate synthase gene (Md-ACS1) accounts for the low level of ethylene production in climacteric fruits of some apple cultivars. Theor Appl Genet. 2000;101(5):742–6.
Costa F, Stella S, Van de Weg WE, Guerra W, Cecchinel M, Dallavia J, Koller B, Sansavini S. Role of the genes Md-ACO1 and Md-ACS1 in ethylene production and shelf life of apple (Malus domestica Borkh). Euphytica. 2005;141(1):181–90.
Nie H, Zhao C, Wu G, Wu Y, Chen Y, Tang D. SR1, a calmodulin-binding transcription factor, modulates plant defense and ethylene-induced senescence by directly regulating NDR1 and EIN3. Plant Physiol. 2012;158(4):1847–59.
Zhang Z, Zhang H, Quan R, Wang X-C, Huang R. Transcriptional regulation of the ethylene response factor LeERF2 in the expression of ethylene biosynthesis genes controls ethylene production in tomato and tobacco. Plant Physiol. 2009;150(1):365–77.
Xiao Y-y, Chen J-y, Kuang J-f. Shan W, Xie H, Jiang Y-m, Lu W-j: Banana ethylene response factors are involved in fruit ripening through their interactions with ethylene biosynthesis genes. J Exp Bot. 2013;64(8):2499–510.
Hu Y, Han Z, Sun Y, Wang S, Wang T, Wang Y, Xu K, Zhang X, Xu X, Han Z, et al. ERF4 affects fruit firmness through TPL4 by reducing ethylene production. Plant J. 2020;103(3):937–50.
Kuai B, Chen J, Hörtensteiner S. The biochemistry and molecular biology of chlorophyll breakdown. J Exp Bot. 2017;69(4):751–67.
Wang X, Chen L, Ma J. Genomic introgression through interspecific hybridization counteracts genetic bottleneck during soybean domestication. Genome Biol. 2019;20(1):22.
Saito T. Advances in Japanese pear breeding in Japan. Breeding Sci. 2016;66(1):46–59.
Yamamoto T, Terakami S, Moriya S, Hosaka F, Kurita K, Kanamori H, Katayose Y, Saito T, Nishitani C. DNA markers developed from genome sequencing analysis in Japanese pear (Pyrus pyrifolia). Acta Hortic. 2013;976(976):477–83.
Chen H, Patterson N, Reich D. Population differentiation as a test for selective sweeps. Genome Res. 2010;20(3):393–402.
Wang L, He F, Huang Y, He J, Yang S, Zeng J, Deng C, Jiang X, Fang Y, Wen S, et al. Genome of wild mandarin and domestication history of mandarin. Mol Plant. 2018;11(8):1024–37.
Paterson AH, Lin Y-R, Li Z, Schertz KF, Doebley JF, Pinson SRM, Liu S-C, Stansel JW, Irvine JE. Convergent domestication of cereal crops by independent mutations at corresponding genetic loci. Science. 1995;269(5231):1714–8.
Lenser T, Theißen G. Molecular mechanisms involved in convergent crop domestication. Trends Plant Sci. 2013;18(12):704–14.
Cheng F, Sun R, Hou X, Zheng H, Zhang F, Zhang Y, Liu B, Liang J, Zhuang M, Liu Y, et al. Subgenome parallel selection is associated with morphotype diversification and convergent crop domestication in Brassica rapa and Brassica oleracea. Nat Genet. 2016;48(10):1218–24.
Li X, Singh J, Qin M, Li S, Zhang X, Zhang M, Khan A, Zhang S, Wu J. Development of an integrated 200K SNP genotyping array and application for genetic mapping, genome assembly improvement and genome wide association studies in pear (Pyrus). Plant Biotechnol J. 2019;17(8):1582–94.
Consonni C, Humphry ME, Hartmann HA, Livaja M, Durner J, Westphal L, Vogel J, Lipka V, Kemmerling B, Schulze-Lefert P, et al. Conserved requirement for a plant host cell protein in powdery mildew pathogenesis. Nat Genet. 2006;38(6):716–20.
Zhou Z, Zhang C, Zhou Y, Hao Z, Wang Z, Zeng X, Di H, Li M, Zhang D, Yong H, et al. Genetic dissection of maize plant architecture with an ultra-high density bin map based on recombinant inbred lines. BMC Genomics. 2016;17(1):178.
Hufford MB, Lubinksy P, Pyhäjärvi T, Devengenzo MT, Ellstrand NC, Ross-Ibarra J. The genomic signature of crop-wild introgression in maize. PLoS Genet. 2013;9(5):e1003477.
Martin RC, Mok MC, Mok DWS. A gene encoding the cytokinin enzyme zeatinO-xylosyltransferase of Phaseolus vulgaris1. Plant Physiol. 1999;120(2):553–8.
McAtee P, Karim S, Schaffer R, David K. A dynamic interplay between phytohormones is required for fruit development, maturation, and ripening. Front Plant Sci. 2013;4:79.
Kumar R, Khurana A, Sharma AK. Role of plant hormones and their interplay in development and ripening of fleshy fruits. J Exp Bot. 2014;65(16):4561–75.
Aubert C, Chanforan C. Postharvest changes in physicochemical properties and volatile constituents of apricot (Prunus armeniaca L.). Characterization of 28 cultivars. J Agr Food Chem. 2007;55(8):3074–82.
Merchante C, Vallarino JG, Osorio S, Aragüez I, Villarreal N, Ariza MT, Martínez GA, Medina-Escobar N, Civello MP, Fernie AR, et al. Ethylene is involved in strawberry fruit ripening in an organ-specific manner. J Agr Food Chem. 2013;64(14):4421–39.
Villarreal NM, Marina M, Nardi CF, Civello PM, Martínez GA. Novel insights of ethylene role in strawberry cell wall metabolism. Plant Sci. 2016;252:1–11.
Xin T, Zhang Z, Li S, Zhang S, Li Q, Zhang Z-H, Huang S, Yang X. Genetic regulation of ethylene dosage for cucumber fruit elongation. Plant Cell. 2019;31(5):1063–76.
Nitsch L, Kohlen W, Oplaat C, Charnikhova T, Cristescu S, Michieli P, Wolters-Arts M, Bouwmeester H, Mariani C, Vriezen WH, et al. ABA-deficiency results in reduced plant and fruit size in tomato. J Plant Physiol. 2012;169(9):878–83.
Moya-León MA, Mattus-Araya E, Herrera R. Molecular events occurring during softening of strawberry fruit. Front Plant Sci. 2019;10:615.
Lindo-García V, Muñoz P, Larrigaudière C, Munné-Bosch S, Giné-Bordonaba J. Interplay between hormones and assimilates during pear development and ripening and its relationship with the fruit postharvest behaviour. Plant Sci. 2020;291:110339.
Molina-Hidalgo FJ, Franco AR, Villatoro C, Medina-Puche L, Mercado JA, Hidalgo MA, Monfort A, Caballero JL, Muñoz-Blanco J, Blanco-Portales R. The strawberry (Fragaria×ananassa) fruit-specific rhamnogalacturonate lyase 1 (FaRGLyase1) gene encodes an enzyme involved in the degradation of cell-wall middle lamellae. J Exp Bot. 2013;64(6):1471–83.
Yang L, Huang W, Xiong F, Xian Z, Su D, Ren M, Li Z. Silencing of SlPL, which encodes a pectate lyase in tomato, confers enhanced fruit firmness, prolonged shelf-life and reduced susceptibility to grey mould. Plant Biotechnol J. 2017;15(12):1544–55.
Calle A, Balas F, Cai L, Iezzoni A, López-Corrales M, Serradilla MJ, Wünsch A. Fruit size and firmness QTL alleles of breeding interest identified in a sweet cherry ‘Ambrunés’ × ‘Sweetheart’ population. Mol Breeding. 2020;40(9):86.
Ornelas-Paz JDJ, Quintana-Gallegos BM, Escalante-Minakata P, Reyes-Hernández J, Ruiz-Cruz S. Relationship between the firmness of Golden Delicious apples and the physicochemical characteristics of the fruits and their pectin during development and ripening. J Food Sci Tech. 2017;55(1):1–9.
Küçükönder H, Vursavuş KK, Üçkardeş F. Determining the factors affecting fruit hardness of different peach types with meta analysis. Turkish J Agric Food Sci Technol. 2014;3(2):98.
Døving A, Måge F. Methods of testing strawberry fruit firmness. Acta Agric Scand B Soil Plant Sci. 2002;52(1):43–51.
Valero C, Crisosto CH, Slaughter D. Relationship between nondestructive firmness measurements and commercially important ripening fruit stages for peaches, nectarines and plums. Postharvest Biol Tec. 2007;44(3):248–53.
Brummell DA. Cell wall disassembly in ripening fruit. Funct Plant Biol. 2006;33(2):103–19.
Hauenstein M, Christ B, Das A, Aubry S, Hörtensteiner S. A role for TIC55 as a hydroxylase of phyllobilins, the products of chlorophyll breakdown during plant senescence. Plant Cell. 2016;28(10):2510–27.
Chou M-L, Liao W-Y, Wei W-C. Li AY-S, Chu C-Y, Wu C-L, Liu C-L, Fu T-H, Lin L-F: The direct involvement of dark-induced Tic55 protein in chlorophyll catabolism and its indirect role in the MYB108-NAC signaling pathway during leaf senescence in Arabidopsis thaliana. Int J Mol Sci. 2018;19(7):1854.
Browning SR, Browning BL. Rapid and accurate haplotype phasing and missing-data inference for whole-genome association studies by use of localized haplotype clustering. Am J Hum Genet. 2007;81(5):1084–97.
Krzywinski M, Schein J, Birol I, Connors J, Gascoyne R, Horsman D, Jones SJ, Marra MA. Circos: an information aesthetic for comparative genomics. Genome Res. 2009;19(9):1639–45.
Danecek P, Auton A, Abecasis G, Albers CA, Banks E, DePristo MA, Handsaker RE, Lunter G, Marth GT, Sherry ST, et al. The variant call format and VCFtools. Bioinformatics. 2011;27(15):2156–8.
Lipka AE, Tian F, Wang Q, Peiffer J, Li M, Bradbury PJ, Gore MA, Buckler ES, Zhang Z. GAPIT: genome association and prediction integrated tool. Bioinformatics. 2012;28(18):2397–9.
Turner SD. qqman: an R package for visualizing GWAS results using Q-Q and manhattan plots. J Open Source Softw. 2018;3(25):731.
Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MAR, Bender D, Maller J, Sklar P, Bakker PIWD, Daly MJ. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007;81(3):559–75.
Shin J-H, Blay S, McNeney B, Graham J. LDheatmap: an R function for graphical display of pairwise linkage disequilibria between single nucleotide polymorphisms. J Stat Softw. 2006;16(3):1–9.
Li X, Xue C, Li J, Qiao X, Li L, Yu L, Huang Y, Wu J. Genome-wide identification, evolution and functional divergence of MYB transcription factors in Chinese white pear (Pyrus bretschneideri). Plant Cell Physiol. 2016;57(4):824.
Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2−ΔΔCT method. Methods. 2001;25(4):402–8.
Bu D, Luo H, Huo P, Wang Z, Zhang S, He Z, Wu Y, Zhao L, Liu J, Guo J et al. KOBAS-i: intelligent prioritization and exploratory visualization of biological functions for gene enrichment analysis. Nucleic Acids Res. 2021;49(W1):W317-W325.
Wu J, et al. The genome of the pear (Pyrus bretschneideri Rehd.). NCBI accession: PRJNA157875. 2012. https://www.ncbi.nlm.nih.gov/bioproject/PRJNA157875.
Wu J, et al. Diversification and independent domestication of Asian and European pears. NCBI accession: PRJNA381668. 2018. https://www.ncbi.nlm.nih.gov/bioproject/PRJNA381668.
The authors wish to acknowledge the help of Dr. Barbara Gilmore (USDA-ARS, National Clonal Germplasm Repository at Corvallis, OR) in providing us with pear fruit samples and the assistance of Dr. Youngsoo Lee (Department of Food Science and Human Nutrition, University of Illinois at Urbana-Champaign) who provided help with the texture analyzer for firmness evaluation. We thank the high-performance computing platforms at the Bioinformatics Center of Nanjing Agricultural University and the Centre of Pear Engineering Technology Research for supporting this project.
This work was funded by the National Key Research and Development Program (2018YFD1000200), the National Natural Science Foundation of China (31725024 and 31801835), the Earmarked Fund for China Agriculture Research System (CARS-28), and the Earmarked Fund for Jiangsu Agricultural Industry Technology System (JATS453).
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. Fruit firmness in wild and cultivated pears.
Table S2. Whole genome sequencing (WGS) data for 113 pear accessions used in our study.
Table S3. Length distribution of putative IBD tracts between populations and within populations. WE, wild European; CA, cultivated Asian; WE, wild European; CE, cultivated European. Bin size = 500 bp.
Table S4. Numbers of IBD tracts in pear genomes (500 kb window) in 10 comparisons.
Table S5. The IBD-enriched regions (IBD-ERs) in Asian and European pear populations.
Table S6. Genes in IBD-enriched regions (IBD-ERs) identified between wild Asian and cultivated Asian pears.
Table S7. Genes in IBD-enriched regions (IBD-ERs) identified between wild European and cultivated European pears.
Table S8. The collection sites of 164 pear accessions used in the GWAS in this study.
Table S9. Quality control of 200,481 SNPs from the SNP array.
Table S10. SNPs, SNP density, and gene density of all 29,269 high-quality SNPs on all 17 pear chromosomes.
Table S11. A summary of the categorized SNPs.
Table S12. Two significant SNPs associated with fruit firmness.
Table S13. Names and sequences (5’ to 3’) of all oligonucleotide primers used in this study.
Fig. S1. Distribution of IBD tract lengths within populations (cultivated Asian-cultivated Asian, wild Asian-wild Asian, cultivated European-cultivated European, and wild European-wild European) and between populations (cultivated Asian-wild Asian, cultivated Asian-cultivated European, cultivated Asian-wild European, wild Asian-cultivated European, wild Asian-wild European, and cultivated European-wild European). WE, wild European; CA, cultivated Asian; WE, wild European; and CE, cultivated European.
Fig. S2. Distribution of IBD tracts along the 17 pear chromosomes. (A) Distribution of IBD tracts in the cultivated Asian-wild Asian (CA-WA; red bars) and cultivated European-wild European (CE-WE; blue bars) comparisons. (B) Lengths of all 17 pear chromosomes. (C and D) KEGG functional enrichment of genes from IBD-ERs of Asian (C) and European (D) pears. Pathways related to carbohydrate metabolism are enclosed in red boxes. The sizes of the dots correspond to the number of genes.
Fig. S3. KEGG analysis of genes within the 3 Mbp IBD-ER.
Fig. S4. The distribution of XP-EHH scores on the 17 pear chromosomes in the WA-CA comparison (A) and WE-CE comparison (B). The cutoff lines represent the top 1% of selection signals (Asian pear population: XP-EHH > 0.828041; European pear population: XP-EHH > 0.680248).
Fig. S5. GWAS of fruit firmness. Quantile-Quantile (Q-Q) plots of a genome-wide association analysis for fruit firmness.
Fig. S6. Sanger DNA sequencing of the CAMTA3 gene in Asian and European pears showing the 9-bp insertion mutation unique to Asian pears. (A) Sequence alignment of the coding sequences of the PbCAMTA3 gene from eight European and nine Asian pear accessions. (B) Gene structure showing the introns and exons of PbCAMTA3 (upper) and the predicted domains (IPT (transcription factor immunoglobulin), IQ motifs (calmodulin-binding), CG-1 (a DNA-binding domain specific to sequence), and ankyrin (ANK) repeats) in the PbCAMTA3 protein structure.
Fig. S7. Fruit firmness of four pear accessions, including ‘Hosui’ (cultivated Asian), ‘Nanguoli’ (cultivated Asian), ‘Starkrimson’ (cultivated European), and ‘Yali’ (cultivated Asian) during five developmental stages of fruit growth. S1 corresponds to the physiological fruit drop stage at 30 days after flowering (DAF), S2 corresponds to the rapid fruit enlargement stage at 55 DAF, S3 corresponds to the fruiting stage at 85 DAF, S4 corresponds to the pre-maturity fruit stage at 115 DAF, and S5 corresponds to the mature fruit stage.
Fig. S8. Distribution of sequencing reads that mapped to the promoter region (upstream 2,000 bp) of the PbTIC55 gene from 113 pear accessions derived from four populations (wild Asian, cultivated Asian, wild European, and cultivated European populations). The results of this analysis did not identify any potential indels or SVs in the promoter region of TIC55 between the Asian and European pear populations.
Fig. S9. The relative expression levels of TIC55 during pear fruit development in ‘Dangshansuli’ (Asian pear) and ‘Starkrimson’ (European pear) based on qRT-PCR. DAF: days after flowering; DAH: days after harvest.
About this article
Cite this article
Song, B., Li, X., Cao, B. et al. An identical-by-descent segment harbors a 12-bp insertion determining fruit softening during domestication and speciation in Pyrus. BMC Biol 20, 215 (2022). https://doi.org/10.1186/s12915-022-01409-w