Linkage map construction involving a reciprocal translocation

Identification of quantitative trait loci associated with R1-mediated resistance to powdery mildew and sex determination in hop (Humulus lupulus L.)

KEY MESSAGE

Two QTL were identified using linkage mapping approaches, one on hop linkage group 3 (qHl_Chr3.PMR1) associated with powdery mildew resistance and a second on linkage group 10 (cqHl_ChrX.SDR1) associated with sex determination. Hop (Humulus lupulus L.) is a dioecious species cultivated for use in beer. Hop powdery mildew, caused by Podosphaera macularis, is a constraint in many growing regions. Thus, identifying markers associated with powdery mildew resistance and sex provides the opportunity to pyramid R-genes and select female plants as seedlings, respectively. Our objectives were to characterize the genetic basis of R1-mediated resistance in the cultivar Zenith which provides resistance to pathogen races in the US, identify quantitative trait loci (QTL) associated with R1 and sex, and develop markers for molecular breeding-based approaches. Phenotypic evaluation of the population indicated that R1-based resistance and sex are inherited monogenically. We constructed a genetic map using 1339 single nucleotide polymorphisms (SNPs) based upon genotype-by-sequencing of 128 F₁ progeny derived from a Zenith × USDA 21058M biparental population. SNPs were assigned to 10 linkage groups comprising a map length of 1204.97 cM with an average density of 0.94 cM/marker. Quantitative trait locus mapping identified qHl_Chr3.PMR1, associated with R1 on linkage group 3 (LOD = 23.57, R² = 57.2%), and cqHl_ChrX.SDR1, associated with sex on linkage group 10 (LOD = 5.42, R² = 25.0%). Kompetitive allele-specific PCR (KASP) assays were developed for both QTL and assessed against diverse germplasm. Our results indicate that KASP markers associated with R1 may be limited to materials that are pedigree-related to Zenith, whereas markers associated with sex may be transferable across populations. The high-density map, QTL, and associated KASP markers will enable selecting for sex and R1-mediated resistance in hop.

Genetic analysis and QTL mapping of domestication-related traits in chili pepper (Capsicum annuum L.)

Chili pepper (Capsicum annuum L.) is one of the oldest and most phenotypically diverse pre-Columbian crops of the Americas. Despite the abundance of genetic resources, the use of wild germplasm and landraces in chili pepper breeding is limited. A better understanding of the evolutionary history in chili peppers, particularly in the context of traits of agronomic interest, can contribute to future improvement and conservation of genetic resources. In this study, an F₂:₃ mapping population derived from a cross between a C. annuum wild accession (Chiltepin) and a cultivated variety (Puya) was used to identify genomic regions associated with 19 domestication and agronomic traits. A genetic map was constructed consisting of 1023 single nucleotide polymorphism (SNP) markers clustered into 12 linkage groups and spanning a total of 1,263.87 cM. A reciprocal translocation that differentiates the domesticated genome from its wild ancestor and other related species was identified between chromosomes 1 and 8. Quantitative trait locus (QTL) analysis detected 20 marker-trait associations for 13 phenotypes, from which 14 corresponded to previously identified loci, and six were novel genomic regions related to previously unexplored domestication-syndrome traits, including form of unripe fruit, seedlessness, deciduous fruit, and growth habit. Our results revealed that the genetic architecture of Capsicum domestication is similar to other domesticated species with few loci with large effects, the presence of QTLs clusters in different genomic regions, and the predominance of domesticated recessive alleles. Our analysis indicates the domestication process in chili pepper has also had an effect on traits not directly related to the domestication syndrome. The information obtained in this study provides a more complete understanding of the genetic basis of Capsicum domestication that can potentially guide strategies for the exploitation of wild alleles.

High density linkage maps, genetic architecture, and genomic prediction of growth and wood properties in Pinus radiata

Background

The growing availability of genomic resources in radiata pine paves the way for significant advances in fundamental and applied genomic research. We constructed robust high-density linkage maps based on exome-capture genotyping in two F₁ populations, and used these populations to perform quantitative trait locus (QTL) scans, genomic prediction and quantitative analyses of genetic architecture for key traits targeted by tree improvement programmes.

Results

Our mapping approach used probabilistic error correction of the marker data, followed by an iterative approach based on stringent parameters. This approach proved highly effective in producing high-density maps with robust marker orders and realistic map lengths (1285–4674 markers per map, with sizes ranging from c. 1643–2292 cM, and mean marker intervals of 0.7–2.1 cM). Colinearity was high between parental linkage maps, although there was evidence for a large chromosomal rearrangement (affecting ~ 90 cM) in one of the parental maps. In total, 28 QTL were detected for growth (stem diameter) and wood properties (wood density and fibre properties measured by Silviscan) in the QTL discovery population, with 1–3 QTL of small to moderate effect size detected per trait in each parental map. Four of these QTL were validated in a second, unrelated F₁ population. Results from genomic prediction and analyses of genetic architecture were consistent with those from QTL scans, with wood properties generally having moderate to high genomic heritabilities and predictive abilities, as well as somewhat less complex genetic architectures, compared to growth traits.

Conclusions

Despite the economic importance of radiata pine as a plantation forest tree, robust high-density linkage maps constructed from reproducible, sequence-anchored markers have not been published to date. The maps produced in this study will be a valuable resource for several applications, including the selection of marker panels for genomic prediction and anchoring a recently completed de novo whole genome assembly. We also provide the first map-based evidence for a large genomic rearrangement in radiata pine. Finally, results from our QTL scans, genomic prediction, and genetic architecture analyses are informative about the genomic basis of variation in important phenotypic traits.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08950-6.

A 2.09 Mb fragment translocation on chromosome 6 causes abnormalities during meiosis and leads to less seed watermelon

Seedlessness is a valuable agronomic trait in watermelon (Citrullus lanatus) breeding. Conventional less seed watermelons are mainly triploid, which has many disadvantages due to unbalanced genome content. Less seed watermelon can be achieved at the diploid level when certain reproductive genes are mutated or by chromosome translocation, which leads to defects during meiosis. However, the formation mechanism of diploid less seed watermelons remains largely unknown. Here, we identified a spontaneous mutant line, watermelon line "148", which can set seeds normally when self-pollinated. A total of 148 × JM F₁ hybrid plants exhibited seed number reductions to 50.3% and 47.3% of those of the two parental lines, respectively, which are considered to be less seed. Examination of pollen viability and hybridization experiments revealed that F₁ hybrids produce semisterile pollen and ovules. Further cytological observations indicated that semisterility was a result of a reciprocal translocation of chromosomes, which exhibited one quadrivalent ring of four chromosomes at prometaphase I during meiosis. RT-qPCR analysis indirectly confirmed that the semisterile phenotype is caused by chromosome translocation rather than disruption of specific meiotic gene expression. F₂ population genetic analysis indicated that the "148" watermelon line is a homozygous translocation and that the less seed phenotype of the F₁ hybrid is prompted by one chromosome fragment translocation. The translocated fragment was further fine mapped to a 2.09 Mb region on chromosome 6 by whole-genome resequencing and genetic map cloning procedures. Our work revealed that a 2.09 Mb chromosome fragment translocation on chromosome 6, causing meiotic defects at metaphase I during meiosis, leads to diploid less seed watermelon. Our findings provide a new promising method for less seed watermelon breeding at the diploid level, as well as a fragment size reference for breeding less seed watermelon through artificially induced chromosome translocation.

Analysis of barley mutants ert-c.1 and ert-d.7 reveals two loci with additive effect on plant architecture

Both mutant ert-c.1 and ert-d.7 carry T2-T3 translocations in the Ert-c gene. Principal coordinate analyses revealed the translocation types and translocation breakpoints. Mutant ert-d.7 is an Ert-c Ert-d double mutant. Mutations in the Ert-c and Ert-d loci are among the most common barley mutations affecting plant architecture. The mutants have various degrees of erect and compact spikes, often accompanied with short and stiff culms. In the current study, complementation tests, linkage mapping, principal coordinate analyses and fine mapping were conducted. We conclude that the original ert-d.7 mutant does not only carry an ert-d mutation but also an ert-c mutation. Combined, mutations in Ert-c and Ert-d cause a pyramid-dense spike phenotype, whereas mutations in only Ert-c or Ert-d give a pyramid and dense phenotype, respectively. Associations between the Ert-c gene and T2-T3 translocations were detected in both mutant ert-c.1 and ert-d.7. Different genetic association patterns indicate different translocation breakpoints in these two mutants. Principal coordinate analysis based on genetic distance and screening of recombinants from all four ends of polymorphic regions was an efficient way to narrow down the region of interest in translocation-involved populations. The Ert-c gene was mapped to the marker interval of 2_0801to1_0224 on 3HL near the centromere. The results illuminate a complex connection between two single genes having additive effects on barley spike architecture and will facilitate the identification of the Ert-c and Ert-d genes.

The origin and widespread occurrence of Sli-based self-compatibility in potato

Self-compatible (SC) diploid potatoes allow innovative potato breeding. Therefore, the Sli gene, originally described in S. chacoense, has received much attention. In elite S. tuberosum diploids, spontaneous berry set is occasionally observed. We aimed to map SC from S. tuberosum origin. Two full-sib mapping populations from non-inbred diploids were used. Bulks were composed based on both pollen tube growth and berry set upon selfing. After DNA sequencing of the parents and bulks, we generated k-mer tables. Set algebra and depth filtering were used to identify bulk-specific k-mers. Coupling and repulsion phase k-mers, transmitted from the SC parent, mapped in both populations to the distal end of chromosome 12. Intersection between the k-mers from both populations, in coupling phase with SC, exposed a shared haplotype of approximately 1.5 Mb. Subsequently, we screened read archives of potatoes and wild relatives for k-mers specific to this haplotype. The well-known SC clones US-W4 and RH89-039-16, but surprisingly, also S. chacoense clone M6 were positives. Hence, the S. tuberosum source of SC seems identical to Sli. Furthermore, the candidate region drastically reduced to 333 kb. Haplotype-specific KASP markers were designed and validated on a panel of diploid clones including another renown SC dihaploid G254. Interestingly, k-mers specific to the SC haplotype were common in tetraploid varieties. Pedigree information suggests that the SC haplotype was introduced into tetraploid varieties via the founder "Rough Purple Chili". We show that Sli is surprisingly widespread and indigenous to the cultivated gene pool of potato.

Identification and characterization of a female gametophyte defect in sdk1-7 +/- abi3-6 +/- heterozygotes of Arabidopsis thaliana

Successful reproduction in angiosperms is dependent on the highly synchronous development of their male and female gametophytes and the ensuing fusion of the gametes from these reproductive tissue types. When crossing a T-DNA insertion line sdk1-7-/-(Salk_024564), one of the S-domain receptor kinases involved in ABA responses with a fast neutron deletion line abi3-6-/-, the F1 heterozygotes (sdk1-7+/-abi3-6 +/-) displayed 50% ovule abortion suggesting a likely gametophytic defects. We identified and characterized an early stage female gametophyte developmental defect in the heterozygous mutant ovules. Recombination frequency analysis of the F2 progenies from selfed heterozygotes revealed a possible pseudo-linkage of sdk1-7 and abi3-6 suggesting a reciprocal translocation event in the heterozygote. Our study emphasizes the importance of robust analysis to distinguish gametophytic defect phenotypes caused by genetic interactions and that resulting from possible chromosomal translocation events.

QTL mapping of mycelial growth and aggressiveness to distinct hosts in Ceratocystis pathogens

Some species of Ceratocystis display strong host specificity, such as C. fimbriata sensu stricto that is restricted to sweet potato (Ipomoea batatas) as host. In contrast, the closely related C. manginecans, infects Acacia mangium and Mangifera indica but is not pathogenic to I. batatas. Despite the economic importance of these fungi, knowledge regarding the genetic factors that influence their pathogenicity and host specificity is limited. A recent inheritance study, based on an interspecific cross between C. fimbriata and C. manginecans and the resultant 70 F₁ progeny, confirmed that traits such as mycelial growth rate, spore production and aggressiveness on A. mangium and I. batatas are regulated by multiple genes. In the present study, a quantitative trait locus (QTL) analysis was performed to determine the genomic loci associated with these traits. All 70 progeny isolates were genotyped with SNP markers and a linkage map was constructed. The map contained 467 SNPs, distributed across nine linkage groups, with a total length of 1203 cm. Using the progeny genotypes and phenotypes, one QTL was identified on the linkage map for mycelial growth rate, one for aggressiveness to A. mangium and two for aggressiveness to I. batatas (P < 0.05). Two candidate genes, likely associated with mycelial growth rate, were identified in the QTL region. The three QTLs associated with aggressiveness to different hosts contained candidate genes involved in protein processing, detoxification and regions with effector genes and high transposable element density. The results provide a foundation for studies considering the function of genes regulating various quantitative traits in Ceratocystis.

High-density genetic map using whole-genome resequencing for fine mapping and candidate gene discovery for disease resistance in peanut

Whole-genome resequencing (WGRS) of mapping populations has facilitated development of high-density genetic maps essential for fine mapping and candidate gene discovery for traits of interest in crop species. Leaf spots, including early leaf spot (ELS) and late leaf spot (LLS), and Tomato spotted wilt virus (TSWV) are devastating diseases in peanut causing significant yield loss. We generated WGRS data on a recombinant inbred line population, developed a SNP-based high-density genetic map, and conducted fine mapping, candidate gene discovery and marker validation for ELS, LLS and TSWV. The first sequence-based high-density map was constructed with 8869 SNPs assigned to 20 linkage groups, representing 20 chromosomes, for the 'T' population (Tifrunner × GT-C20) with a map length of 3120 cM and an average distance of 1.45 cM. The quantitative trait locus (QTL) analysis using high-density genetic map and multiple season phenotyping data identified 35 main-effect QTLs with phenotypic variation explained (PVE) from 6.32% to 47.63%. Among major-effect QTLs mapped, there were two QTLs for ELS on B05 with 47.42% PVE and B03 with 47.38% PVE, two QTLs for LLS on A05 with 47.63% and B03 with 34.03% PVE and one QTL for TSWV on B09 with 40.71% PVE. The epistasis and environment interaction analyses identified significant environmental effects on these traits. The identified QTL regions had disease resistance genes including R-genes and transcription factors. KASP markers were developed for major QTLs and validated in the population and are ready for further deployment in genomics-assisted breeding in peanut.

High density SNP mapping and QTL analysis for time of leaf budburst in Corylus avellana L

The growing area of European hazelnut (Corylus avellana L.) is increasing, as well as the number of producing countries, and there is a pressing need for new improved cultivars. Hazelnut conventional breeding process is slow, due to the length of juvenile phase and the high heterozygosity level. The development of genetic linkage maps and the identification of molecular markers tightly linked to QTL (quantitative trait loci) of agronomic interest are essential tools for speeding up the selection of seedlings carrying desired traits through marker-assisted selection. The objectives of this study were to enrich a previous linkage map and confirm QTL related to time of leaf budburst, using an F₁ population obtained by crossing Tonda Gentile delle Langhe with Merveille de Bollwiller. Genotyping-by-Sequencing was used to identify a total of 9,999 single nucleotide polymorphism markers. Well saturated linkage maps were constructed for each parent using the double pseudo-testcross mapping strategy. A reciprocal translocation was detected in Tonda Gentile delle Langhe between two non-homologous chromosomes. Applying a bioinformatic approach, we were able to disentangle ‘pseudo-linkage’ between markers, removing markers around the translocation breakpoints and obtain a linear order of the markers for the two chromosomes arms, for each linkage group involved in the translocation. Twenty-nine QTL for time of leaf budburst were identified, including a stably expressed region on LG_02 of the Tonda Gentile delle Langhe map. The stability of these QTL and their coding sequence content indicates promise for the identification of specific chromosomal regions carrying key genes involved in leaf budburst.

The potential of probabilistic graphical models in linkage map construction

Probabilistic graphical models show great potential for robust and reliable construction of linkage maps. We show how to use probabilistic graphical models to construct high-quality linkage maps in the face of data perturbations caused by genotyping errors and reciprocal translocations. It has been shown that linkage map construction can be hampered by the presence of genotyping errors and chromosomal rearrangements such as inversions and translocations. Here, we report a novel method for linkage map construction using probabilistic graphical models. The method is proven, both theoretically and practically, to be effective in filtering out markers that contain genotyping errors. In particular, it carries out marker filtering and ordering simultaneously, and is therefore superior to the standard post hoc filtering using nearest-neighbour stress. Furthermore, we demonstrate empirically that the proposed method offers a promising solution to linkage map construction in the case of a reciprocal translocation.

Computationally efficient map construction in the presence of segregation distortion

We present a novel estimator for map construction in the presence of segregation distortion which is highly computationally efficient. For multi-parental designs this estimator outperforms methods that do not account for segregation distortion, at no extra computational cost. Inclusion of genetic markers exhibiting segregation distortion in a linkage map can result in biased estimates of genetic distance and distortion of map positions. Removal of distorted markers is hence a typical filtering criterion; however, this may result in exclusion of biologically interesting regions of the genome such as introgressions and translocations. Estimation of additional parameters characterizing the distortion is computationally slow, as it relies on estimation via the Expectation Maximization algorithm or a higher dimensional numerical optimisation. We propose a robust M-estimator (RM) capable of handling tens of thousands of distorted markers from a single linkage group. We show via simulation that for multi-parental designs the RM estimator can perform much better than uncorrected estimation, at no extra computational cost. We then apply the RM estimator to chromosome 2B in wheat in a multi-parent population segregating for the Sr36 introgression, a known transmission distorter. The resulting map contains over 700 markers, and is consistent with maps constructed from crosses which do not exhibit segregation distortion.

Chromosomal rearrangements and the genetics of reproductive barriers in mimulus (monkey flowers)

Chromosomal rearrangements may directly cause hybrid sterility and can facilitate speciation by preserving local adaptation in the face of gene flow. We used comparative linkage mapping with shared gene-based markers to identify potential chromosomal rearrangements between the sister monkeyflowers Mimulus lewisii and Mimulus cardinalis, which are textbook examples of ecological speciation. We then remapped quantitative trait loci (QTLs) for floral traits and flowering time (premating isolation) and hybrid sterility (postzygotic isolation). We identified three major regions of recombination suppression in the M. lewisii × M. cardinalis hybrid map compared to a relatively collinear Mimulus parishii × M. lewisii map, consistent with a reciprocal translocation and two inversions specific to M. cardinalis. These inferences were supported by targeted intraspecific mapping, which also implied a M. lewisii-specific reciprocal translocation causing chromosomal pseudo-linkage in both hybrid mapping populations. Floral QTLs mapped in this study, along with previously mapped adaptive QTLs, were clustered in putatively rearranged regions. All QTLs for male sterility, including two underdominant loci, mapped to regions of recombination suppression. We argue that chromosomal rearrangements may have played an important role in generating and consolidating barriers to gene flow as natural selection drove the dramatic ecological and morphological divergence of these species.

A framework genetic map of Muscadinia rotundifolia

This study presents a framework linkage map based on microsatellite markers for Muscadinia rotundifolia (1n = 20). The mapping population consisted of 206 progeny generated from a cross of two M. rotundifolia varieties, 'Fry' and 'Trayshed'. A total of 884 primers were tested for their ability to amplify markers: 686 amplified and 312 simple sequence repeat (SSR) primer pairs generated 322 polymorphic markers for either one or both parents. The map for the female parent 'Fry' consisted of 212 markers and covered 879 cM on 18 chromosomes. The average distance between the markers was 4.1 cM and chromosome 6 was not represented due to a lack of polymorphic markers. The map for the male parent 'Trayshed' consisted of 191 markers and covered 841 cM on 19 chromosomes. The consensus map consisted of 314 markers on 19 chromosomes with a total distance of 1,088 cM, which represented 66 % of the distance covered by the Vitis vinifera reference linkage map. Marker density varied greatly among chromosomes from 5 to 35 mapped markers. Relatively good synteny was observed across 19 chromosomes based on markers in common with the V. vinifera reference map. Extreme segregation distortion was observed for chromosome 8 and 14 on the female parent map, and 4 on the male parent map. The lack of mapping coverage for the 20th M. rotundifolia chromosome is discussed in relation to possible evolutionary events that led to the reduction in chromosome number from 21 to 19 in the ancestral genome.

Genetic characterization of a reciprocal translocation present in a widely grown barley variety

Artificially induced translocation stocks have been used to physically map the barley genome; however, natural translocations are extremely uncommon in cultivated genotypes. Albacete is a barley variety widely grown in recent decades in Spain and carrying a reciprocal translocation which obviously does not affect its agronomical fitness. This translocation has been characterized by a combination of cytological and molecular genetic approaches. Firstly, recombination frequencies between markers on chromosomes 1H and 3H were estimated to determine the boundaries of the reciprocal interchange. Secondly, 1H-3H wheat barley telosome addition lines were used to assign selected markers to chromosome arms. Thirdly, fluorescence in situ hybridization (FISH) with rDNA probes (5S and 18S-5.8S-26S) and microsatellite probes [(ACT)(5), (AAG)(5) and (CAG)(5)] was used to determine the locations of the translocation breakpoints more precisely. Fourthly, fine-mapping of the regions around the translocation breakpoints was used to increase the marker density for comparative genomics. The results obtained in this study indicate that the translocation is quite large with breakpoints located on the long arms of chromosomes 1H and 3H, between the pericentromeric (AAG)(5) bands and above the (ACT)(5) interstitial distal bands, resulting in the reciprocal translocation 1HS.1HL-3HL and 3HS.3HL-1HL. The gene content around the translocation breakpoints could be inferred from syntenic relationships observed among different species from the grass family Poaceae (rice, Sorghum and Brachypodium) and was estimated at approximately 1,100 and 710 gene models for 1H and 3H, respectively. Duplicated segments between chromosomes Os01 and Os05 in rice derived from ancestral duplications within the grass family overlap with the translocation breakpoints on chromosomes 1H and 3H in the barley variety Albacete.