Fine Mapping of qAL5.2 Controlling Anther Length in Oryza sativa

Anther length is the critical floral trait determining hybrid rice seed production and is controlled by many quantitative trait loci (QTL). However, the cloning of genes specifically controlling anther size has yet to be reported. Here, we report the fine mapping of qAL5.2 for anther size using backcross inbred lines (BILs) in the genetic background of Oryza sativa indica Huazhan (HZ). Gene chip analysis on the BC4F2 and BC5F1 population identified effective loci on Chr1, Chr5, and Chr8 and two genomic regions on Chr5, named qAL5.1 and qAL5.2. qAL5.2 was identified in both populations with LOD values of 17.54 and 10.19, which explained 35.73% and 25.1% of the phenotypic variances, respectively. Ultimately qAL5.2 was localized to a 73 kb region between HK139 and HK140 on chromosome 5. And we constructed two near-isogenic lines (NILs) for RNA-seq analysis, named NIL-qAL5.2HZ and NIL-qAL5.2KLY, respectively. The result of the GO enrichment analysis revealed that differential genes were significantly enriched in the carbohydrate metabolic process, extracellular region, and nucleic acid binding transcription, and KEGG enrichment analysis revealed that alpha-linolenic acid metabolism was significantly enriched. Meanwhile, candidate genes of qAL5.2 were analyzed in RNA-seq, and it was found that ORF8 is differentially expressed between NIL-qAL5.2HZ and NIL-qAL5.2KLY. The fine mapping of qAL5.2 conferring anther length will promote the breed improvement of the restorer line and understanding of the mechanisms driving crop mating patterns.

QTL Analysis for Rice Quality-Related Traits and Fine Mapping of qWCR3

The quality of rice, evaluated using multiple quality-related traits, is the main determinant of its market competitiveness. In this study, two japonica rice varieties with significant differences in quality-related traits were used as parents to construct two populations, BC3F2 and BC3F2:3, with Kongyu131 (KY131) as the recurrent parent. A genetic linkage map was constructed using the BC3F2 population based on 151 pairs of SSR/InDel polymorphic markers selected between the parents. Grain-shape-related traits (grain length GL, grain width GW, and length-to-width ratio LWR), chalkiness-related traits (white-core rate WCR, white-belly rate WBR, white-back rate BR, and chalkiness rate CR), and amylose content (AC) were investigated in the two populations in 2017 and 2018. Except for BR and CR, the traits showed similar characteristics with a normal distribution in both populations. Genetic linkage analysis was conducted for these quality-related traits, and a total of 37 QTLs were detected in the two populations. Further validation was performed on the newly identified QTLs with larger effects, and three grain shape QTLs and four chalkiness QTLs were successfully validated in different environments. One repeatedly validated QTL, qWCR3, was selected for fine mapping and was successfully narrowed down to a 100 kb region in which only two genes, LOC_0s03g45210 and LOC_0s03g45320, exhibited sequence variations between the parents. Furthermore, the variation of LOC_Os03g45210 leads to a frameshift mutation and premature protein termination. The results of this study provide a theoretical basis for positional cloning of the qWCR3 gene, thus offering new genetic resources for rice quality improvement.

Genome-wide association study identified five quantitative trait loci and two candidate genes for digestive traits in Suhuai pigs

This work aimed to identify markers and candidate genes underlying porcine digestive traits. In total, 331 pigs were genotyped by 80 K Chip data or 50 K Chip data. For apparent neutral detergent fiber digestibility, a total of 19 and 21 candidate single nucleotide polymorphisms (SNP) were respectively identified using a genome-wide efficient mixed-model association algorithm and linkage-disequilibrium adjusted kinship. Among them, three quantitative trait locus (QTL) regions were identified. For apparent acid detergent fiber digestibility, a total of 16 and 17 SNPs were identified by these two methods, respectively. Of these, three QTL regions were also identified. Moreover, two candidate genes (MST1 and LATS1), which are functionally related to intestinal homeostasis and health, were detected near these significant SNPs. Taken together, our results could provide a basis for deeper research on digestive traits in pigs.

Fine mapping and analysis of a candidate gene controlling Phytophthora blight resistance in cucumber

Cucumber blight is a destructive disease. The best way to control this disease is resistance breeding. This study focuses on disease resistance gene mapping and molecular marker development. We used the resistant cucumber, JSH, and susceptible cucumber, B80, as parents to construct F2 populations. Bulked segregant analysis (BSA) combined with specific length amplified fragment sequencing (SLAF-seq) were used, from which we developed cleaved amplified polymorphic sequence (CAPs) markers to map the resistance gene. Resistance in F2 individuals showed a segregation ratio of resistance:susceptibility close to 3:1. The gene in JSH resistant cucumber was mapped to an interval of 9.25 kb, and sequencing results for the three genes in the mapped region revealed three mutations at base sites 225, 302, and 591 in the coding region of Csa5G139130 between JSH and B80, but no mutations in coding regions of Csa5G139140 and Csa5G139150. The mutations caused changes in amino acids 75 and 101 of the protein encoded by Csa5G139130, suggesting that Csa5G139130 is the most likely resistance candidate gene. We developed a molecular marker, CAPs-4, as a closely linked marker for the cucumber blight resistance gene. This is the first report on mapping of a cucumber blight resistance gene and will provideg a useful marker for molecular breeding of cucumber resistance to Phytophthora blight.

Novel breast cancer susceptibility loci under linkage peaks identified in African ancestry consortia

BACKGROUND

Expansion of genome-wide association studies across population groups is needed to improve our understanding of shared and unique genetic contributions to breast cancer. We performed association and replication studies guided by a priori linkage findings from African ancestry (AA) relative pairs.

METHODS

We performed fixed-effect inverse-variance weighted meta-analysis under three significant AA breast cancer linkage peaks (3q26-27, 12q22-23, and 16q21-22) in 9241 AA cases and 10 193 AA controls. We examined associations with overall breast cancer as well as estrogen receptor (ER)-positive and negative subtypes (193,132 SNPs). We replicated associations in the African-ancestry Breast Cancer Genetic Consortium (AABCG).

RESULTS

In AA women, we identified two associations on chr12q for overall breast cancer (rs1420647, OR = 1.15, p = 2.50×10-6; rs12322371, OR = 1.14, p = 3.15×10-6), and one for ER-negative breast cancer (rs77006600, OR = 1.67, p = 3.51×10-6). On chr3, we identified two associations with ER-negative disease (rs184090918, OR = 3.70, p = 1.23×10-5; rs76959804, OR = 3.57, p = 1.77×10-5) and on chr16q we identified an association with ER-negative disease (rs34147411, OR = 1.62, p = 8.82×10-6). In the replication study, the chr3 associations were significant and effect sizes were larger (rs184090918, OR: 6.66, 95% CI: 1.43, 31.01; rs76959804, OR: 5.24, 95% CI: 1.70, 16.16).

CONCLUSION

The two chr3 SNPs are upstream to open chromatin ENSR00000710716, a regulatory feature that is actively regulated in mammary tissues, providing evidence that variants in this chr3 region may have a regulatory role in our target organ. Our study provides support for breast cancer variant discovery using prioritization based on linkage evidence.

A Bayesian fine-mapping model using a continuous global-local shrinkage prior with applications in prostate cancer analysis

The aim of fine mapping is to identify genetic variants causally contributing to complex traits or diseases. Existing fine-mapping methods employ Bayesian discrete mixture priors and depend on a pre-specified maximum number of causal variants, which may lead to sub-optimal solutions. In this work, we propose a Bayesian fine-mapping method called h2-D2, utilizing a continuous global-local shrinkage prior. We also present an approach to define credible sets of causal variants in continuous prior settings. Simulation studies demonstrate that h2-D2 outperforms current state-of-the-art fine-mapping methods such as SuSiE and FINEMAP in accurately identifying causal variants and estimating their effect sizes. We further applied h2-D2 to prostate cancer analysis and discovered some previously unknown causal variants. In addition, we inferred 369 target genes associated with the detected causal variants and several pathways that were significantly over-represented by these genes, shedding light on their potential roles in prostate cancer development and progression.

Unravelling novel and pleiotropic genes for cannon bone circumference and bone mineral density in Yorkshire pigs

Leg weakness is a prevalent health condition in pig farms. The augmentation of cannon bone circumference and bone mineral density can effectively improve limb strength in pigs and alleviate leg weakness. This study measured forelimb cannon bone circumference (fCBC) and rear limb cannon bone circumference (rCBC) using an inelastic tapeline and rear limb metatarsal area bone mineral density (raBMD) using a dual-energy X-ray absorptiometry bone density scanner. The samples of Yorkshire castrated boars were genotyped using a 50K single-nucleotide polymorphism (SNP) array. The SNP-chip data were imputed to the level of whole-genome sequencing data (iWGS). This study used iWGS data to perform genome-wide association studies and identified novel significant SNPs associated with fCBC on SSC6, SSC12, and SSC13, rCBC on SSC12 and SSC14, and raBMD on SSC7. Based on the high phenotypic and genetic correlations between CBC and raBMD, multi-trait meta-analysis was performed to identify pleiotropic SNPs. A significant potential pleiotropic quantitative trait locus (QTL) regulating both CBC and raBMD was identified on SSC15. Bayes fine mapping was used to establish the confidence intervals for these novel QTLs with the most refined confidence interval narrowed down to 56 kb (15.11 to 15.17 Mb on SSC12 for fCBC). Furthermore, the confidence interval for the potential pleiotropic QTL on SSC15 in the meta-analysis was narrowed down to 7.45 kb (137.55 to137.56 Mb on SSC15). Based on the biological functions of genes, the following genes were identified as novel regulatory candidates for different phenotypes: DDX42, MYSM1, FTSJ3, and MECOM for fCBC; SMURF2, and STC1 for rCBC; RGMA for raBMD. Additionally, RAMP1, which was determined to be located 23.68 kb upstream of the confidence interval of the QTL on SSC15 in the meta-analysis, was identified as a potential pleiotropic candidate gene regulating both CBC and raBMD. These findings offered valuable insights for identifying pathogenic genes and elucidating the genetic mechanisms underlying CBC and BMD.

Fine mapping of Green a, Ga, on chromosome 27 in Bombyx mori

Some strains of silkworms produce green cocoons of varying intensities. This results from quantitative and qualitative differences in flavonoid pigments, which are influenced by the environment and genetic background. We discovered that the appearance of a faint green cocoon is regulated by a gene (G27) located on chromosome 27. Through mating experiments, we found that G27 is identical to an essential flavonoid cocoon gene, Ga. This locus has not been previously described. Furthermore, we narrowed down the Ga region to 438 kbp using molecular markers. Within this region, several predicted genes for sugar transporters form a cluster structure, suggesting that Ga is among them.

Fine mapping of qAHPS07 and functional studies of AhRUVBL2 controlling pod size in peanut (Arachis hypogaea L.)

Cultivated peanut (Arachis hypogaea L.) is an important oil and cash crop. Pod size is one of the major traits determining yield and commodity characteristic of peanut. Fine mapping of quantitative trait locus (QTL) and identification of candidate genes associated with pod size are essential for genetic improvement and molecular breeding of peanut varieties. In this study, a major QTL related to pod size, qAHPS07, was fine mapped to a 36.46 kb interval on chromosome A07 using F2 , recombinant inbred line (RIL) and secondary F2 populations. qAHPS07 explained 38.6%, 23.35%, 37.48%, 25.94% of the phenotypic variation for single pod weight (SPW), pod length (PL), pod width (PW) and pod shell thickness (PST), respectively. Whole genome resequencing and gene expression analysis revealed that a RuvB-like 2 protein coding gene AhRUVBL2 was the most likely candidate for qAHPS07. Overexpression of AhRUVBL2 in Arabidopsis led to larger seeds and plants than the wild type. AhRUVBL2-silenced peanut seedlings represented small leaves and shorter main stems. Three haplotypes were identified according to three SNPs in the promoter of AhRUVBL2 among 119 peanut accessions. Among them, SPW, PW and PST of accessions carrying Hap_ATT represent 17.6%, 11.2% and 26.3% higher than those carrying Hap_GAC，respectively. In addition, a functional marker of AhRUVBL2 was developed. Taken together, our study identified a key functional gene of peanut pod size, which provides new insights into peanut pod size regulation mechanism and offers practicable markers for the genetic improvement of pod size-related traits in peanut breeding.

Functionally distinct regions of the locus Leishmania major response 15 control IgE or IFNγ level in addition to skin lesions

Leishmaniasis, a disease caused by parasites of Leishmania spp., endangers more than 1 billion people living in endemic countries and has three clinical forms: cutaneous, mucocutaneous, and visceral. Understanding of individual differences in susceptibility to infection and heterogeneity of its pathology is largely lacking. Different mouse strains show a broad and heterogeneous range of disease manifestations such as skin lesions, splenomegaly, hepatomegaly, and increased serum levels of immunoglobulin E and several cytokines. Genome-wide mapping of these strain differences detected more than 30 quantitative trait loci (QTLs) that control the response to Leishmania major. Some control different combinations of disease manifestations, but the nature of this heterogeneity is not yet clear. In this study, we analyzed the L. major response locus Lmr15 originally mapped in the strain CcS-9 which carries 12.5% of the genome of the resistant strain STS on the genetic background of the susceptible strain BALB/c. For this analysis, we used the advanced intercross line K3FV between the strains BALB/c and STS. We confirmed the previously detected loci Lmr15, Lmr18, Lmr24, and Lmr27 and performed genetic dissection of the effects of Lmr15 on chromosome 11. We prepared the interval-specific recombinant strains 6232HS1 and 6229FUD, carrying two STS-derived segments comprising the peak linkage of Lmr15 whose lengths were 6.32 and 17.4 Mbp, respectively, and analyzed their response to L. major infection. These experiments revealed at least two linked but functionally distinct chromosomal regions controlling IFNγ response and IgE response, respectively, in addition to the control of skin lesions. Bioinformatics and expression analysis identified the potential candidate gene Top3a. This finding further clarifies the genetic organization of factors relevant to understanding the differences in the individual risk of disease.

Fine Mapping and Candidate Gene Analysis of Rice Grain Length QTL qGL9.1

Grain length (GL) is one of the crucial determinants of rice yield and quality. However, there is still a shortage of knowledge on the major genes controlling the inheritance of GL in japonica rice, which severely limits the improvement of japonica rice yields. Here, we systemically measured the GL of 667 F₂ and 1570 BC₃F₃ individuals derived from two cultivated rice cultivars, Pin20 and Songjing15, in order to identify the major genomic regions associated with GL. A novel major QTL, qGL9.1, was mapped on chromosome 9, which is associated with the GL, using whole-genome re-sequencing with bulked segregant analysis. Local QTL linkage analysis with F₂ and fine mapping with the recombinant plant revealed a 93-kb core region on qGL9.1 encoding 15 protein-coding genes. Only the expression level of LOC_Os09g26970 was significantly different between the two parents at different stages of grain development. Moreover, haplotype analysis revealed that the alleles of Pin20 contribute to the optimal GL (9.36 mm) and GL/W (3.31), suggesting that Pin20 is a cultivated species carrying the optimal GL variation of LOC_Os09g26970. Furthermore, a functional-type mutation (16398989-bp, G>A) located on an exon of LOC_Os09g26970 could be used as a molecular marker to distinguish between long and short grains. Our experiments identified LOC_Os09g26970 as a novel gene associated with GL in japonica rice. This result is expected to further the exploration of the genetic mechanism of rice GL and improve GL in rice japonica varieties by marker-assisted selection.

Genetic Mapping, Candidate Gene Identification and Marker Validation for Host Plant Resistance to the Race 4 of Fusarium oxysporum f. sp. cubense Using Musa acuminata ssp. malaccensis

Fusarium wilt of banana is a devastating disease that has decimated banana production worldwide. Host resistance to Fusarium oxysporum f. sp. Cubense (Foc), the causal agent of this disease, is genetically dissected in this study using two Musa acuminata ssp. Malaccensis segregating populations, segregating for Foc Tropical (TR4) and Subtropical (STR4) race 4 resistance. Marker loci and trait association using 11 SNP-based PCR markers allowed the candidate region to be delimited to a 12.9 cM genetic interval corresponding to a 959 kb region on chromosome 3 of 'DH-Pahang' reference assembly v4. Within this region, there was a cluster of pattern recognition receptors, namely leucine-rich repeat ectodomain containing receptor-like protein kinases, cysteine-rich cell-wall-associated protein kinases, and leaf rust 10 disease-resistance locus receptor-like proteins, positioned in an interspersed arrangement. Their transcript levels were rapidly upregulated in the resistant progenies but not in the susceptible F₂ progenies at the onset of infection. This suggests that one or several of these genes may control resistance at this locus. To confirm the segregation of single-gene resistance, we generated an inter-cross between the resistant parent 'Ma850' and a susceptible line 'Ma848', to show that the STR4 resistance co-segregated with marker '28820' at this locus. Finally, an informative SNP marker 29730 allowed the locus-specific resistance to be assessed in a collection of diploid and polyploid banana plants. Of the 60 lines screened, 22 lines were predicted to carry resistance at this locus, including lines known to be TR4-resistant, such as 'Pahang', 'SH-3362', 'SH-3217', 'Ma-ITC0250', and 'DH-Pahang/CIRAD 930'. Additional screening in the International Institute for Tropical Agriculture's collection suggests that the dominant allele is common among the elite 'Matooke' NARITA hybrids, as well as in other triploid or tetraploid hybrids derived from East African highland bananas. Fine mapping and candidate gene identification will allow characterization of molecular mechanisms underlying the TR4 resistance. The markers developed in this study can now aid the marker-assisted selection of TR4 resistance in breeding programs around the world.

Fine Mapping of fw6.3, a Major-Effect Quantitative Trait Locus That Controls Fruit Weight in Tomato

Tomato (Solanum lycopersicum) is a widely consumed vegetable, and the tomato fruit weight is a key yield component. Many quantitative trait loci (QTLs) controlling tomato fruit weight have been identified, and six of them have been fine-mapped and cloned. Here, four loci controlling tomato fruit weight were identified in an F₂ population through QTL seq.; fruit weight 6.3 (fw6.3) was a major-effect QTL and its percentage of variation explanation (R²) was 0.118. This QTL was fine-mapped to a 62.6 kb interval on chromosome 6. According to the annotated tomato genome (version SL4.0, annotation ITAG4.0), this interval contained seven genes, including Solyc06g074350 (the SELF-PRUNING gene), which was likely the candidate gene underlying variation in fruit weight. The SELF-PRUNING gene contained a single-nucleotide polymorphism that resulted in an amino acid substitution in the protein sequence. The large-fruit allele of fw6.3 (fw6.3_HG) was overdominant to the small-fruit allele fw6.3_RG. The soluble solids content was also increased by fw6.3_HG. These findings provide valuable information that will aid the cloning of the FW6.3 gene and ongoing efforts to breed tomato plants with higher yield and quality via molecular marker-assisted selection.

Map-based cloning of the APRR2 gene controlling green stigma in bitter gourd (Momordica charantia)

Bitter gourd is an economically important vegetable and medicinal crop distinguished by its bitter fruits. Its stigma color is widely used to assess the distinctiveness, uniformity, and stability of bitter gourd varieties. However, limited researches have been dedicated to genetic basis of its stigma color. In this study, we employed bulked segregant analysis (BSA) sequencing to identify a single dominant locus McSTC1 located on pseudochromosome 6 through genetic mapping of an F₂ population (n =241) derived from the cross between green and yellow stigma parental lines. An F₂-derived F₃ segregation population (n = 847) was further adopted for fine mapping, which delimited the McSTC1 locus to a 13.87 kb region containing one predicted gene McAPRR2 (Mc06g1638), a homolog of the Arabidopsis two-component response regulator-like gene AtAPRR2. Sequence alignment analysis of McAPRR2 revealed that a 15 bp insertion at exon 9 results in a truncated GLK domain of its encoded protein, which existed in 19 bitter gourd varieties with yellow stigma. A genome-wide synteny search of the bitter gourd McAPRR2 genes in Cucurbitaceae family revealed its close relationship with other cucurbits APRR2 genes that are corresponding to white or light green fruit skin. Our findings provide insights into the molecular marker-assisted breeding of bitter gourd stigma color and the mechanism of gene regulation for stigma color.

Bayesian multivariant fine mapping using the Laplace prior

Currently, the only effect size prior that is routinely implemented in a Bayesian fine-mapping multi-single-nucleotide polymorphism (SNP) analysis is the Gaussian prior. Here, we show how the Laplace prior can be deployed in Bayesian multi-SNP fine mapping studies. We compare the ranking performance of the posterior inclusion probability (PIP) using a Laplace prior with the ranking performance of the corresponding Gaussian prior and FINEMAP. Our results indicate that, for the simulation scenarios we consider here, the Laplace prior can lead to higher PIPs than either the Gaussian prior or FINEMAP, particularly for moderately sized fine-mapping studies. The Laplace prior also appears to have better worst-case scenario properties. We reanalyse the iCOGS case-control data from the CASP8 region on Chromosome 2. Even though this study has a total sample size of nearly 90,000 individuals, there are still some differences in the top few ranked SNPs if the Laplace prior is used rather than the Gaussian prior. R code to implement the Laplace (and Gaussian) prior is available at https://github.com/Kevin-walters/lapmapr.

Identification and fine mapping of a new bacterial blight resistance gene, Xa43(t), in Zhangpu wild rice (Oryza rufipogon)

Bacterial blight (BB) is currently considered one of the most serious rice diseases and is caused by Xanthomonas oryzae pv. oryzae (Xoo). Numerous studies have shown that breeding resistant rice varieties is one of the most effective methods to prevent BB, and it is important to identify and isolate more BB resistance (R) genes from different rice resources. Using a map-based approach, we identified a new QTL/gene, Xa43(t), from ZhangPu wild rice, which was highly resistant to the BB isolate PX099. We performed bulked segregant analysis combined with candidate gene prediction to identify the candidate gene. The Xa43(t) gene was narrowed down to a 29-kb region containing four putative genes. More importantly, the candidate gene Xa43(t) did not affect the main agronomic traits of rice. We also identified a widely applicable molecular marker, namely Inde1-18, which co-segregates with the Xa43(t) gene. The Xa43(t) gene is a new broad-spectrum BB resistance gene without identified alleles and has good application prospects for rice disease resistance breeding.

Identification of CaPs locus involving in purple stripe formation on unripe fruit, reveals allelic variation and alternative splicing of R2R3-MYB transcription factor in pepper (Capsicum annuum L.)

The purple color of unripe pepper fruit is attributed to the accumulation of anthocyanins. Only a few genes controlling the biosynthesis and regulation of anthocyanins have been cloned in Capsicum. In this study, we performed a bulked segregant analysis of the purple striped trait using an F² population derived from a cross between the immature purple striped fruit line Chen12-4-1-1-1-1 and the normal green fruit line Zhongxian101-M-F⁹. We mapped the CaPs locus to an 841.39 kb region between markers M-CA690-Xba and MCA710-03 on chromosome 10. CA10g11690 encodes an R2R3-MYB transcription factor that is involved in the biosynthesis of anthocyanins as the best candidate gene. Overexpression and silencing in transformed tobacco (Nicotiana tabacum) lines indicated that CA10g11690 is involved in the formation of purple stripes in the exocarp. A comparison of parental sequences identified an insertion fragment of 1,926 bp in the second intron region of Chen12-4, and eight SNPs were detected between the two parents. Additionally, there were 49 single nucleotide polymorphic variations, two sequence deletions, and four sequence insertions in the promoter region. We found that CA10g11690 undergoes alternative splicing and generates different transcripts. Thus, the functional transcript of CA10g11690 appeared to be primarily involved in the development of purple phenotype in the exocarp. Our data provide new insight into the mechanism of anthocyanin biosynthesis and a theoretical basis for the future breeding of purple striped pepper varieties.

The lingering effects of Neanderthal introgression on human complex traits

The genetic variants introduced into the ancestors of modern humans from interbreeding with Neanderthals have been suggested to contribute an unexpected extent to complex human traits. However, testing this hypothesis has been challenging due to the idiosyncratic population genetic properties of introgressed variants. We developed rigorous methods to assess the contribution of introgressed Neanderthal variants to heritable trait variation and applied these methods to analyze 235,592 introgressed Neanderthal variants and 96 distinct phenotypes measured in about 300,000 unrelated white British individuals in the UK Biobank. Introgressed Neanderthal variants make a significant contribution to trait variation (explaining 0.12% of trait variation on average). However, the contribution of introgressed variants tends to be significantly depleted relative to modern human variants matched for allele frequency and linkage disequilibrium (about 59% depletion on average), consistent with purifying selection on introgressed variants. Different from previous studies (McArthur et al., 2021), we find no evidence for elevated heritability across the phenotypes examined. We identified 348 independent significant associations of introgressed Neanderthal variants with 64 phenotypes. Previous work (Skov et al., 2020) has suggested that a majority of such associations are likely driven by statistical association with nearby modern human variants that are the true causal variants. Applying a customized fine-mapping led us to identify 112 regions across 47 phenotypes containing 4303 unique genetic variants where introgressed variants are highly likely to have a phenotypic effect. Examination of these variants reveals their substantial impact on genes that are important for the immune system, development, and metabolism.

Upgrading the genome of an elite japonica rice variety Kongyu 131 for lodging resistance improvement

Developing a new rice variety requires tremendous efforts and years of input. To improve the defect traits of the excellent varieties becomes more cost and time efficient than breeding a completely new variety. Kongyu 131 is a high-performing japonica variety with early maturity, high yield, wide adaptability and cold resistance, but the poor-lodging resistance hinders the industrial production of Kongyu 131 in the Northeastern China. In this study, we attempted to improve the lodging resistance of Kongyu 131 from perspectives of both gene and trait. On the one hand, by QTL analysis and fine mapping we discovered the candidate gene loci. The following CRISPR/Cas9 and transgenic complementation study confirmed that Sd1 dominated the lodging resistance and favourable allele was mined for precise introduction and improvement. On the other hand, the Sd1 allelic variant was identified in Kongyu 131 by sequence alignment, then introduced another excellent allelic variation by backcrossing. Then, the two new resulting Kongyu 131 went through the field evaluation under different environments, planting densities and nitrogen fertilizer conditions. The results showed that the plant height of upgraded Kongyu 131 was 17%-26% lower than Kongyu 131 without penalty in yield. This study demonstrated a precise and targeted way to update the rice genome and upgrade the elite rice varieties by improving only a few gene defects from the perspective of breeding.

A 163-bp insertion in the Capana10g000198 encoding a MYB transcription factor causes male sterility in pepper (Capsicum annuum L.)

Male sterility provides an efficient approach for commercial exploitation of heterosis. Despite more than 20 genic male sterile (GMS) mutants documented in pepper (Capsicum annuum L.), only two causal genes have been successfully identified. Here, a novel spontaneous recessive GMS mutant, designated msc-3, is identified and characterized at both phenotypic and histological levels. Pollen abortion of msc-3 mutant may be due to the delayed tapetum degradation, leading to the non-degeneration of tetrads callosic wall. Then, a modified MutMap method and molecular marker linkage analysis were employed to fine mapping the msc-3 locus, which was delimited to the ~139.91-kb region harboring 10 annotated genes. Gene expression and structure variation analyses indicate the Capana10g000198, encoding a R2R3-MYB transcription factor, is the best candidate gene for the msc-3 locus. Expression profiling analysis shows the Capana10g000198 is an anther-specific gene, and a 163-bp insertion in the Capana10g000198 is highly correlated with the male sterile (MS) phenotype. Additionally, downregulation of Capana10g000198 in male fertile plants through virus-induced gene silencing resulted in male sterility. Finally, possible regulatory relationships of the msc-3 gene with the other two reported pepper GMS genes, msc-1 and msc-2, have been studied, and comparative transcriptome analysis reveals the expression of 16 GMS homologs are significantly downregulated in the MS anthers. Overall, our results reveal that Capana10g000198 is the causal gene underlying the msc-3 locus, providing important theoretical clues and basis for further in-depth study on the regulatory mechanisms of pollen development in pepper.

Mapping a leaf rust resistance gene LrOft in durum wheat Ofanto and its suppressor SuLrOft in common wheat

Epidemics of leaf rust (caused by the fungal pathogen Puccinia triticina Erikss., Pt) raise concerns regarding sustainability of wheat production. Deployment of resistant cultivars is the most effective and economic strategy for combating this disease. Ofanto is a durum wheat cultivar that exhibits high resistance to Pt race PHT throughout its entire growing period. In the present study, we identified a leaf rust resistance gene in Ofanto and temporarily designated it as LrOft. LrOft was mapped to a 2.5 cM genetic interval in chromosome arm 6BL between Indel markers 6B6941 and 6B50L24. During introgression of LrOft from Ofanto to common wheat it was observed that F₁ plants of Ofanto crossed with Shi4185 exhibited leaf rust resistance whereas the F₁ of Ofanto crossed with ND4503 was susceptible. In order to map the presumed suppressor locus, a Shi4185/ND4503//Ofanto three-way pentaploid population was generated and SuLrOft was mapped on chromosome arm 2AS. SuLrOft was mapped within a 2.6 cM genetic interval flanked by 2AS50L14 and 2AS50L6. Fine mapping using 2,268 plants of the three-way cross narrowed the suppressor locus to a 68.2-kbp physical interval according to IWGSC RefSeq v1.1. Sequence analysis of genes in the physical interval revealed that TraesCS2A02G110800 encoding an RPP-13-like protein with an NB-ARC domain was a potential candidate for SuLrOft.

Fine mapping of the Hairy glume (Hg) gene in a chromosome variation region at the distal terminus of 1AS

Trichomes are differentiated epidermal cells and exist on above-ground organs of nearly all land plants with important roles in resistance to a wide range of biotic and abiotic stresses. We attempted to obtain candidate gene (s) for Hairy glume (Hg), responsible for the trichome on wheat glume, by using bulked segregant exome capture sequencing (BSE-Seq), while Hg was only mapped in 0.52-3.26 Mb of 1AS. To further fine map this gene and identify candidate genes in this region, a near isogenic line-derived population consisting of 2,050 F₂ lines was generated in the present study. By analyzing this population, Hg was fine mapped into a 0.90 cM region covering a physical distance of ~825.03 Kb encompassing 6 high- and 23 low-confidence genes in the reference genome of Chinese Spring. A presence-absence variation was identified in the fine mapping region through analyses of sequence-tagged sites markers and genome sequences of the hairy glume parent of the near isogenic lines. The results presented here will be useful for further cloning Hg in wheat.

Characterization and Fine Mapping of qRPR1-3 and qRPR3-1, Two Major QTLs for Rind Penetrometer Resistance in Maize

Stalk strength is one of the most important traits in maize, which affects stalk lodging resistance and, consequently, maize harvestable yield. Rind penetrometer resistance (RPR) as an effective and reliable measurement for evaluating maize stalk strength is positively correlated with stalk lodging resistance. In this study, one F₂ and three F_2:3 populations derived from the cross of inbred lines 3705I (the low RPR line) and LH277 (the high RPR line) were constructed for mapping quantitative trait loci (QTL), conferring RPR in maize. Fourteen RPR QTLs were identified in four environments and explained the phenotypic variation of RPR from 4.14 to 15.89%. By using a sequential fine-mapping strategy based on the progeny test, two major QTLs, qRPR1-3 and qRPR3-1, were narrowed down to 4-Mb and 550-kb genomic interval, respectively. The quantitative real-time PCR (qRT-PCR) assay was adopted to identify 12 candidate genes responsible for QTL qRPR3-1. These findings should facilitate the identification of the polymorphism loci underlying QTL qRPR3-1 and molecular breeding for RPR in maize.

Fine Mapping and Characterization of an Aphid-Resistance Gene in the Soybean Landrace Fangzheng Moshidou

The soybean aphid poses a severe threat to soybean quality and yield by sucking phloem sap and transmitting plant viruses. An early-maturing and highly resistant soybean landrace, Fangzheng Moshidou, with markedly reduced aphid colonization has been identified by screening of aphid-resistant soybean accessions. In a population derived from the cross of Fangzheng Moshidou with the susceptible cultivar Beifeng 9, resistance was conferred by a single dominant gene. Three linked markers, Satt114, Satt334, and Sct_033, on chromosome 13 were identified by bulked-segregant analysis. Additional simple-sequence repeat and single-nucleotide polymorphism (SNP) markers were developed for gene mapping. The resistance of Fangzheng Moshidou was fine-mapped to the interval between the SNP markers YCSNP20 and YCSNP80, corresponding to 152.8 kb in the Williams 82 assembly 2 genome. This region was near the reported loci Rag2 and Rag5 but did not overlap the interval containing them. A unique haplotype is described for Fangzheng Moshidou that distinguishes it from soybean accessions PI 587972, PI 594879, and PI 567301B in the interval containing Rag2 and Rag5. These results indicate that Fangzheng Moshidou harbors a novel gene at a tightly linked resistance locus, designated as RagFMD. Fourteen candidate genes were annotated in the fine-mapping region, including seven NBS-LRR genes, which are usually considered resistance genes in plant defense. Most of these candidate genes showed variations distinguishing the resistant and susceptible parents and some genes also showed differences in expression between the two parental lines and at several times after aphid infestation. Isolation of RagFMD would advance the study of molecular mechanisms of soybean aphid resistance and contribute to precise selection of resistant soybeans.

Identification of qGL3.5, a Novel Locus Controlling Grain Length in Rice Through Bulked Segregant Analysis and Fine Mapping

Grain length (GL) directly affects the yield and quality of rice. Very few cloned GL-related genes are applied in production because their yield-increasing effects are not obvious, and the overall regulatory networks underlying the associated processes remain poorly understood. DNA samples from two bulk DNA pools (L-pool and S-pool) and their parents (KJ01 and Huaye 3) were subjected to high-throughput sequencing. Using bulked segregant analysis (BSA), qGL3.5 was mapped to a 0.34-Mb "hotspot" region on chromosome 3 that contains 37 genes related to various traits. Then, qGL3.5 was mapped to the genomic interval between the flanking markers M2 and M3 using 2786 BC₄F₂ individuals. Because the region from B5 to B6 was not the associated region under BSA-seq analysis, qGL3.5 was narrowed down to the interval between B6 and M3, which spanned 24.0-kb. Of all 37 genes with non-synonymous single-nucleotide polymorphisms (SNPs) between KJ01 and Huaye 3 based on BSA-seq analysis, only one complete annotated gene, ORF18 (Gene ID: LOC_Os03g42790.1) was found. ORF18 encodes an IBR-RING zinc-finger-related protein, with one really interesting new gene (RING) and two in between ring finger (IBR) domains. The knockout of ORF18 derived from Huaye 3 using clustered, regularly interspaced, short palindromic repeat (CRISPR)/CRISPR-associated 9 (Cas9) editing technology increased the GL of the mutant by approximately 2.2 mm. The novel locus qGL3.5 negatively regulated GL by promoting elongation of the longitudinal cell of the grain outer glume. These results provide a new genetic resource for rice grain shape breeding and a starting point for the functional characterization of the wild rice GL gene.

Identification of a Candidate restorer-of-fertility Gene Rf3 Encoding a Pentatricopeptide Repeat Protein for the Cytoplasmic Male Sterility in Soybean

The cytoplasmic male sterility/restorer-of-fertility (CMS/Rf) system plays a vital role in high-efficiency hybrid seed production in crops, including soybean (Glycine max (L.) Merr.). The markers linked to fertility restoration and the restorer-of-fertility (Rf) genes are essential because they can facilitate the breeding of new CMS lines and production of commercial hybrid soybean seeds. To date, several soybean Rf genes have been mapped to various genetic loci in diverse genetic populations. However, the mapping range of restorer genes remains narrow, with relatively limited practical applicability. Therefore, in the present study, F2 and F3 segregating populations derived from the CMS line JLCMS5A crossed with the restorer line JLR2 were developed and used for Rf3 gene fine mapping. Genetic investigation indicated that the restorer line JLR2 was controlled by a single dominant gene, Rf3. By integrating bulk-segregant analysis and next-generation sequencing, a 4 Mb region on chromosome 9 was identified, which was most likely the target region harboring the candidate gene responsible for fertility restoration. This region was further narrowed down to 86.44 Kb via fine mapping in F2 and F3 populations using SSR, InDel, and dCAPS markers. This region contained 10 putative genes (Glyma.09G171100-Glyma.09G172000). Finally, Glyma.09G171200, which encodes a mitochondria-targeted pentatricopeptide repeat protein, was proposed as the potential candidate for Rf3 using sequence alignment and expression analysis in restorer and CMS lines. Based on single-nucleotide polymorphisms in Glyma.09G171200, a CAPS marker co-segregated with Rf3 named CAPS1712 was developed. Our results will be fundamental in the assisted selection and creation of potent lines for the production and rapid selection of novel restorer lines.

Discovery of Genomic Regions and Candidate Genes Controlling Root Development Using a Recombinant Inbred Line Population in Rapeseed (Brassica napus L.)

Marker-assisted selection enables breeders to quickly select excellent root architectural variations, which play an essential role in plant productivity. Here, ten root-related and shoot biomass traits of a new F₆ recombinant inbred line (RIL) population were investigated under hydroponics and resulted in high heritabilities from 0.61 to 0.83. A high-density linkage map of the RIL population was constructed using a Brassica napus 50k Illumina single nucleotide polymorphism (SNP) array. A total of 86 quantitative trait loci (QTLs) explaining 4.16–14.1% of the phenotypic variances were detected and integrated into eight stable QTL clusters, which were repeatedly detected in different experiments. The codominant markers were developed to be tightly linked with three major QTL clusters, qcA09-2, qcC08-2, and qcC08-3, which controlled both root-related and shoot biomass traits and had phenotypic contributions greater than 10%. Among these, qcA09-2, renamed RT.A09, was further fine-mapped to a 129-kb interval with 19 annotated genes in the B. napus reference genome. By integrating the results of real-time PCR and comparative sequencing, five genes with expression differences and/or amino acid differences were identified as important candidate genes for RT.A09. Our findings laid the foundation for revealing the molecular mechanism of root development and developed valuable markers for root genetic improvement in rapeseed.

Fine mapping and cloning of the major seed protein quantitative trait loci on soybean chromosome 20

Soybean is the most important source of protein meal worldwide and the quantitative trait loci (QTL) cqSeed protein‐003 on chromosome 20 exerts the greatest additive effect of any protein QTL mapped in the crop. Through genetic mapping and candidate gene downregulation, we identified that an insertion/deletion variant in Glyma.20G85100 is the likely gene that underlies this important QTL.

Fine mapping of Rf5 region for a sorghum fertility restorer gene and microsynteny analysis across grass species

Cytoplasmic male sterility (CMS) is widely used to control pollination in the production of commercial F₁ hybrid seed in sorghum. So far, 6 major fertility restorer genes, Rf1 to Rf6, have been reported in sorghum. Here, we fine-mapped the Rf5 locus on sorghum chromosome 5 using descendant populations of a 'Nakei MS-3A' × 'JN43' cross. The Rf5 locus was narrowed to a 140-kb region in BTx623 genome (161-kb in JN43) with 16 predicted genes, including 6 homologous to the rice fertility restorer Rf1 (PPR.1 to PPR.6). These 6 homologs have tandem pentatricopeptide repeat (PPR) motifs. Many Rf genes encode PPR proteins, which bind RNA transcripts and modulate gene expression at the RNA level. No PPR genes were detected at the Rf5 locus on the corresponding homologous chromosome of rice, foxtail millet, or maize, so this gene cluster may have originated by chromosome translocation and duplication after the divergence of sorghum from these species. Comparison of the sequences of these genes between fertile and CMS lines identified PPR.4 as the most plausible candidate gene for Rf5.

Bulk segregation analysis in the NGS era: a review of its teenage years

Bulk segregation analysis (BSA) utilizes a strategy of pooling individuals with extreme phenotypes to conduct economical and rapidly linked marker screening or quantitative trait locus (QTL) mapping. With the development of next-generation sequencing (NGS) technology in the past 10 years, BSA methods and technical systems have been gradually developed and improved. At the same time, the ever-decreasing costs of sequencing accelerate NGS-based BSA application in different species, including eukaryotic yeast, grain crops, economic crops, horticultural crops, trees, aquatic animals, and insects. This paper provides a landscape of BSA methods and reviews the BSA development process in the past decade, including the sequencing method for BSA, different populations, different mapping algorithms, associated region threshold determination, and factors affecting BSA mapping. Finally, we summarize related strategies in QTL fine mapping combining BSA.

Fine Mapping Analysis of the MHC Region to Identify Variants Associated With Chinese Vitiligo and SLE and Association Across These Diseases

The important role of MHC in the pathogenesis of vitiligo and SLE has been confirmed in various populations. To map the most significant MHC variants associated with the risk of vitiligo and SLE, we conducted fine mapping analysis using 1117 vitiligo cases, 1046 SLE cases and 1693 healthy control subjects in the Han-MHC reference panel and 1000 Genomes Project phase 3. rs113465897 (P=1.03×10^-13, OR=1.64, 95%CI =1.44–1.87) and rs3129898 (P=4.21×10^-17, OR=1.93, 95%CI=1.66–2.25) were identified as being most strongly associated with vitiligo and SLE, respectively. Stepwise conditional analysis revealed additional independent signals at rs3130969(p=1.48×10^-7, OR=0.69, 95%CI=0.60–0.79), HLA-DPB1*03:01 (p=1.07×10^-6, OR=1.94, 95%CI=1.49–2.53) being linked to vitiligo and HLA-DQB1*0301 (P=4.53×10^-7, OR=0.62, 95%CI=0.52-0.75) to SLE. Considering that epidemiological studies have confirmed comorbidities of vitiligo and SLE, we used the GCTA tool to analyse the genetic correlation between these two diseases in the HLA region, the correlation coefficient was 0.79 (P=5.99×10^-10, SE=0.07), confirming their similar genetic backgrounds. Our findings highlight the value of the MHC region in vitiligo and SLE and provide a new perspective for comorbidities among autoimmune diseases.

Candidate resistance genes to foliar phylloxera identified at Rdv3 of hybrid grape

The foliage of the native grape species Vitis riparia and certain cold-hardy hybrid grapes are particularly susceptible to the insect pest phylloxera, Daktulosphaira vitifoliae Fitch. A previous study using a cold-hardy hybrid grape biparental F1 population (N~125) detected the first quantitative trait locus (QTL) for foliar resistance on chromosome 14, designated as resistance to Daktulosphaira vitifoliae 3 (Rdv3). This locus spans a ~7-Mbp (10-20 cM) region and is too wide for effective marker-assisted selection or identification of candidate genes. Therefore, we fine mapped the QTL using a larger F1 population, GE1783 (N~1023), and genome-wide rhAmpSeq haplotype markers. Through three selective phenotyping experiments replicated in the greenhouse, we screened 184 potential recombinants of GE1783 using a 0 to 7 severity rating scale among other phylloxera severity traits. A 500-kb fine mapped region at 4.8 Mbp on chromosome 14 was identified. The tightly linked rhAmpSeq marker 14_4805213 and flanking markers can be used for future marker-assisted breeding. This region contains 36 candidate genes with predicted functions in disease resistance (R genes and Bonzai genes) and gall formation (bifunctional 3-dehydroquinate dehydratase/shikimate dehydrogenase). Disease resistance genes suggest a traditional R-gene-mediated resistance mechanism often accompanied by a hypersensitive response, which has been widely studied in the plant pathology field. A novel resistance mechanism, non-responsiveness to phylloxera gall formation is proposed as a function of the bifunctional dehydratase gene, which plays a role in gallic acid biosynthesis and is important in gall formation. This study has implications for improvement of foliar phylloxera resistance in cold-hardy hybrid germplasm and is a starting place to understand the mechanism of resistance in crops to gall-forming insects.

Combining Fine Mapping, Whole-Genome Re-Sequencing, and RNA-Seq Unravels Candidate Genes for a Soybean Mutant with Short Petioles and Weakened Pulvini

A short petiole is an important agronomic trait for the development of plant ideotypes with high yields. However, the genetic basis underlying this trait remains unclear. Here, we identified and characterized a novel soybean mutant with short petioles and weakened pulvini, designated as short petioles and weakened pulvini (spwp). Compared with the wild type (WT), the spwp mutant displayed shortened petioles, owing to the longitudinally decreased cell length, and exhibited a smaller pulvinus structure due to a reduction in motor cell proliferation and expansion. Genetic analysis showed that the phenotype of the spwp mutant was controlled by two recessive nuclear genes, named as spwp1 and spwp2. Using a map-based cloning strategy, the spwp1 locus was mapped in a 183 kb genomic region on chromosome 14 between markers S1413 and S1418, containing 15 annotated genes, whereas the spwp2 locus was mapped in a 195 kb genomic region on chromosome 11 between markers S1373 and S1385, containing 18 annotated genes. Based on the whole-genome re-sequencing and RNA-seq data, we identified two homologous genes, Glyma.11g230300 and Glyma.11g230600, as the most promising candidate genes for the spwp2 locus. In addition, the RNA-seq analysis revealed that the expression levels of genes involved in the cytokinin and auxin signaling transduction networks were altered in the spwp mutant compared with the WT. Our findings provide new gene resources for insights into the genetic mechanisms of petiole development and pulvinus establishment, as well as soybean ideotype breeding.

A Major Quantitative Trait Loci Cluster Controlling Three Components of Yield and Plant Height Identified on Chromosome 4B of Common Wheat

Wheat yield is not only affected by three components of yield, but also affected by plant height (PH). Identification and utilization of the quantitative trait loci (QTL) controlling these four traits is vitally important for breeding high-yielding wheat varieties. In this work, we conducted a QTL analysis using the recombinant inbred lines (RILs) derived from a cross between two winter wheat varieties of China, "Nongda981" (ND981) and "Nongda3097" (ND3097), exhibiting significant differences in spike number per unit area (SN), grain number per spike (GNS), thousand grain weight (TGW), and PH. A total of 11 environmentally stable QTL for these four traits were detected. Among them, four major and stable QTLs (QSn.cau-4B-1.1, QGns.cau-4B-1, QTgw.cau-4B-1.1, and QPh.cau-4B-1.2) explaining the highest phenotypic variance for SN, GNS, TGW, and PH, respectively, were mapped on the same genomic region of chromosome 4B and were considered a QTL cluster. The QTL cluster spanned a genetic distance of about 12.3 cM, corresponding to a physical distance of about 8.7 Mb. Then, the residual heterozygous line (RHL) was used for fine mapping of the QTL cluster. Finally, QSn.cau-4B-1.1, QGns.cau-4B-1, and QPh.cau-4B-1.2 were colocated to the physical interval of about 1.4 Mb containing 31 annotated high confidence genes. QTgw.cau-4B-1.1 was divided into two linked QTL with opposite effects. The elite NILs of the QTL cluster increased SN and PH by 55.71-74.82% and 14.73-23.54%, respectively, and increased GNS and TGW by 29.72-37.26% and 5.81-11.24% in two environments. Collectively, the QTL cluster for SN, GNS, TGW, and PH provides a theoretical basis for improving wheat yield, and the fine-mapping result will be beneficial for marker-assisted selection and candidate genes cloning.

Fine Mapping of a Vigor QTL in Chickpea (Cicer arietinum L.) Reveals a Potential Role for Ca4_TIFY4B in Regulating Leaf and Seed Size

Plant vigor is a complex trait for which the underlying molecular control mechanisms remain unclear. Vigorous plants tend to derive from larger seeds and have greater early canopy cover, often with bigger leaves. In this study, we delimited the size of a major vigor quantitative trait locus (QTL) on chickpea chromosome 4-104.4 kb, using recombinant association analysis in 15 different heterogeneous inbred families, derived from a Rupali/Genesis836 recombinant inbred line population. The phenotypic and molecular genetic analysis provided evidence for a role of the gene Ca4_TIFY4B, in determining leaf and seed size in chickpea. A non-synonymous single-nucleotide polymorphism (SNP) in the high-vigor parent was located inside the core motif TIFYCG, resulting in a residue change T[I/S]FYCG. Complexes formed by orthologs of Ca4_TIFY4B (PEAPOD in Arabidopsis), Novel Interactor of JAZ (CaNINJA), and other protein partners are reported to act as repressors regulating the transcription of downstream genes that control plant organ size. When tested in a yeast 2-hybrid (Y2H) assay, this residue change suppressed the interaction between Ca4_TIFY4B and CaNINJA. This is the first report of a naturally occurring variant of the TIFY family in plants. A robust gene-derived molecular marker is available for selection in chickpea for seed and plant organ size, i.e., key component traits of vigor.

Integrating physiology, genetics, and transcriptome to decipher a new thermo-sensitive and light-sensitive virescent leaf gene mutant in cucumber

Chloroplasts are the material basis of photosynthesis, and temperature and light severely affect chloroplast development and thus influence photosynthetic efficiency. This study identified a spontaneous virescent leaf mutant, SC311Y, whose cotyledons and true leaves were yellow and gradually turned green. However, temperature and light affected the process of turning green. In addition, this mutant (except at the seedling stage) had ruffled leaves with white stripes, sterile males, and poorly fertile female flowers. Genetic characteristics analysis revealed that the recessive gene controlled the virescent leaf. Two F₂ populations mapped v-3 to the interval of 33.54-35.66 Mb on chromosome 3. In this interval, BSA-Seq, RNA-Seq, and cDNA sequence analyses revealed only one nonsynonymous mutation in the Csa3G042730 gene, which encoded the RNA exosome supercomplex subunit resurrection1 (RST1). Csa3G042730 was predicted to be the candidate gene controlling the virescent leaf, and the candidate gene may regulate chloroplast development by regulating plastid division2 (PDV2). A transcriptome analysis showed that different factors caused the reduced chlorophyll and carotenoid content in the mutants. To our knowledge, this study is the first report of map-based cloning related to virescent leaf, male-sterile, and chloroplast RNA regulation in cucumber. The results could accelerate the study of the RNA exosome supercomplex for the dynamic regulation of chloroplast RNA.

Fine Mapping of Clubroot Resistance Loci CRA8.1 and Candidate Gene Analysis in Chinese Cabbage (Brassica rapa L.)

Clubroot is caused by Plasmodiophora brassicae, which threatens Brassicaceae crop production worldwide. In recent years, there has been an outbreak and rapid spread of clubroot in many major cruciferous crop-producing areas of China. In this study, we identified a cabbage material DingWen (DW) with different resistant capabilities from Huashuang5R (H5R) and Huayouza62R of Brassica napus, which are currently used as the main resistant cultivars for clubroot management in China. We used a next-generation sequencing-based bulked segregant analysis approach, combined with genetic mapping to identify clubroot-resistant (CR) genes from F₁ population generated from a cross between the DW (CR) and HZSX (clubroot susceptible). The CR locus of DW (named CRA8.1) was mapped to a region between markers A08-4346 and A08-4853, which contains two different loci CRA8.1a and CRA8.1b after fine mapping. The CRA8.1b loci contain a fragment of 395 kb between markers A08-4624 and A08-4853 on A08 chromosome, and it is responsible for the resistance to PbZj and PbXm isolates. However, together with CRA8.1a, corresponding to a 765-kb region between markers A08-4346 and A08-4624, then it can confer resistance to PbXm ⁺. Finally, through expression analysis between resistant and susceptible materials, two genes encoding TIR-NBS-LRR proteins (BraA08g039211E and BraA08g039212E) and one gene encoding an RLP protein (BraA08g039193E) were identified to be the most likely CR candidates for the peculiar resistance in DW.

Fine Mapping of qd1, a Dominant Gene that Regulates Stem Elongation in Bread Wheat

Stem elongation is a critical phase for yield determination and, as a major trait, is targeted for manipulation for improvement in bread wheat (Triticum aestivum L.). In a previous study, we characterized a mutant showing rapid stem elongation but with no effect on plant height at maturity. The present study aimed to finely map the underlying mutated gene, qd1, in this mutant. By analyzing an F₂ segregating population consisting of 606 individuals, we found that the qd1 gene behaved in a dominant manner. Moreover, by using the bulked segregant RNA sequencing (BSR-seq)-based linkage analysis method, we initially mapped the qd1 gene to a 13.55 Mb region on chromosome 4B (from 15.41 to 28.96 Mb). This result was further confirmed in F₂ and BC₃F₂ segregating populations. Furthermore, by using transcriptome sequencing data, we developed 14 Kompetitive Allele-Specific PCR (KASP) markers and then mapped the qd1 gene to a smaller and more precise 5.08 Mb interval from 26.80 to 31.88 Mb. To develop additional markers to finely map the qd1 gene, a total of 4,481 single-nucleotide polymorphisms (SNPs) within the 5.08 Mb interval were screened, and 25 KASP markers were developed based on 10x-depth genome resequencing data from both wild-type (WT) and mutant plants. The qd1 gene was finally mapped to a 1.33 Mb interval from 28.86 to 30.19 Mb on chromosome 4B. Four candidate genes were identified in this region. Among them, the expression pattern of only TraesCS4B02G042300 in the stems was concurrent with the stem development of the mutant and WT. The qd1 gene could be used in conjunction with molecular markers to manipulate stem development in the future.

Genome-Wide Association and Mendelian Randomization Analysis Reveal the Causal Relationship Between White Blood Cell Subtypes and Asthma in Africans

Background: White blood cell (WBC) traits and their subtypes such as basophil count (Bas), eosinophil count (Eos), lymphocyte count (Lym), monocyte count (Mon), and neutrophil counts (Neu) are known to be associated with diseases such as stroke, peripheral arterial disease, and coronary heart disease. Methods: We meta-analyze summary statistics from genome-wide association studies in 17,802 participants from the African Partnership for Chronic Disease Research (APCDR) and African ancestry individuals from the Blood Cell Consortium (BCX2) using GWAMA. We further carried out a Bayesian fine mapping to identify causal variants driving the association with WBC subtypes. To access the causal relationship between WBC subtypes and asthma, we conducted a two-sample Mendelian randomization (MR) analysis using summary statistics of the Consortium on Asthma among African Ancestry Populations (CAAPA: n cases = 7,009, n control = 7,645) as our outcome phenotype. Results: Our metanalysis identified 269 loci at a genome-wide significant value of (p = 5 × 10-9) in a composite of the WBC subtypes while the Bayesian fine-mapping analysis identified genetic variants that are more causal than the sentinel single-nucleotide polymorphism (SNP). We found for the first time five novel genes (LOC126987/MTCO3P14, LINC01525, GAPDHP32/HSD3BP3, FLG-AS1/HMGN3P1, and TRK-CTT13-1/MGST3) not previously reported to be associated with any WBC subtype. Our MR analysis showed that Mon (IVW estimate = 0.38, CI: 0.221, 0.539, p < 0.001), Neu (IVW estimate = 0.189, CI: 0.133, 0.245, p < 0.001), and WBCc (IVW estimate = 0.185, CI: 0.108, 0.262, p < 0.001) are associated with increased risk of asthma. However, there was no evidence of causal relationship between Lym and asthma risk. Conclusion: This study provides insight into the relationship between some WBC subtypes and asthma and potential route in the treatment of asthma and may further inform a new therapeutic approach.

Identification of Candidate Genes for Self-Compatibility in Perennial Ryegrass (Lolium perenne L.)

Self-incompatibility (SI) is a genetic mechanism preventing self-pollination in ~40% of plant species. Two multiallelic loci, called S and Z, control the gametophytic SI system of the grass family (Poaceae), which contains all major forage grasses. Loci independent from S and Z have been reported to disrupt SI and lead to self-compatibility (SC). A locus causing SC in perennial ryegrass (Lolium perenne L.) was previously mapped on linkage group (LG) 5 in an F2 population segregating for SC. Using a subset of the same population (n = 68), we first performed low-resolution quantitative trait locus (QTL) mapping to exclude the presence of additional, previously undetected contributors to SC. The previously reported QTL on LG 5 explained 38.4% of the phenotypic variation, and no significant contribution from other genomic regions was found. This was verified by the presence of significantly distorted markers in the region overlapping with the QTL. Second, we fine mapped the QTL to 0.26 centimorgan (cM) using additional 2,056 plants and 23 novel sequence-based markers. Using Italian ryegrass (Lolium multiflorum Lam.) genome assembly as a reference, the markers flanking SC were estimated to span a ~3 Mb region encoding for 57 predicted genes. Among these, seven genes were proposed as relevant candidate genes based on their annotation and function described in previous studies. Our study is a step forward to identify SC genes in forage grasses and provides diagnostic markers for marker-assisted introgression of SC into elite germplasm.

Fine Mapping of a New Major QTL-qGLS8 for Gray Leaf Spot Resistance in Maize

Gray leaf spot (GLS), caused by different species of Cercospora, is a fungal, non-soil-borne disease that causes serious reductions in maize yield worldwide. The identification of major quantitative trait loci (QTLs) for GLS resistance in maize is essential for developing marker-assisted selection strategies in maize breeding. Previous research found a significant difference (P < 0.01) in GLS resistance between T32 (highly resistant) and J51 (highly susceptible) genotypes of maize. Initial QTL analysis was conducted in an F_{2 : 3} population of 189 individuals utilizing genetic maps that were constructed using 181 simple sequence repeat (SSR) markers. One QTL (qGLS8) was detected, defined by the markers umc1130 and umc2354 in three environments. The qGLS8 QTL detected in the initial analysis was located in a 51.96-Mb genomic region of chromosome 8 and explained 7.89-14.71% of the phenotypic variation in GLS resistance in different environments. We also developed a near isogenic line (NIL) BC₃F₂ population with 1,468 individuals and a BC₃F₂-Micro population with 180 individuals for fine mapping. High-resolution genetic and physical maps were constructed using six newly developed SSRs. The QTL-qGLS8 was narrowed down to a 124-kb region flanked by the markers ym20 and ym51 and explained up to 17.46% of the phenotypic variation in GLS resistance. The QTL-qGLS8 contained seven candidate genes, such as an MYB-related transcription factor 24 and a C ₃ H transcription factor 347), and long intergenic non-coding RNAs (lincRNAs). The present study aimed to provide a foundation for the identification of candidate genes for GLS resistance in maize.

Mapping and Analysis of a Novel Genic Male Sterility Gene in Watermelon (Citrullus lanatus)

Seed production is critical for watermelon production, which mostly involves first-generation hybrid varieties. However, watermelon hybrid seed production currently requires complex procedures, including artificial isolation and pollination. Therefore, the development and use of a male-sterile system to generate watermelon hybrids can simplify the process. The scarcity of male-sterile watermelon germplasm resources necessitates the use of molecular breeding methods. Unfortunately, the genes responsible for male sterility in watermelon have not been cloned. Thus, the genetic basis of the male sterility remains unknown. In this study, two DNA pools derived from male-sterile and normal plants in the F2 population were used for whole-genome resequencing. The Illumina high-throughput sequencing resulted in 62.99 Gbp clean reads, with a Q30 of 80% after filtering. On the basis of the SNP index association algorithm, eight candidate regions (0.32 Mb) related to specific traits were detected on chromosome 6. Expression pattern analyses and watermelon transformation studies generated preliminary evidence that Cla006625 encodes a pollen-specific leucine-rich repeat protein (ClaPEX1) influencing the male sterility of watermelon. The identification and use of genic male sterility genes will promote watermelon male sterility research and lay the foundation for the efficient application of seed production technology.

Fine Mapping of the MHC Region Identifies Novel Variants Associated with HBV-Related Hepatocellular Carcinoma in Han Chinese

Introduction

Genome-wide association studies identified susceptibility loci in the major histocompatibility complex region for hepatitis B virus (HBV)-related hepatocellular carcinoma (HCC). However, the causal variants underlying HBV-related HCC pathogenesis remain elusive.

Methods

With a total of 1,161 HBV-related HCC cases and 1,353 chronic HBV carriers without HCC, we imputed human leukocyte antigen (HLA) variants based on a Chinese HLA reference panel and evaluated the associations of these variants with the risk of HBV-related HCC. Conditional analyses were used to identify independent signals associated with the risk of HBV-related HCC (P false-discovery rate (FDR) <0.20). A total of 14,930 variants within the MHC region were genotyped or imputed.

Results

We identified two variants, rs114401688 (P = 1.05 × 10⁻⁶, P_FDR = 2.43 × 10⁻³) and rs115126566 (P = 9.04 × 10⁻⁵, P_FDR = 1.77 × 10⁻¹), that are independently associated with the risk of HBV-related HCC. Single nucleotide polymorphism (SNP) rs114401688 is in linkage disequilibrium with a previously reported SNP rs9275319. In the current study, we found that its association with HCC could be explained by HLA-DQB1*04 and HLA-DRB1*04. SNP rs115126566 is a novel risk variant and may function by regulating transcriptions of HLA-DPA1/DPB1 through enhancer-mediated mechanisms. HLA zygosity analysis showed that homozygosity at HLA-DQB1 gene is suggestively associated with a higher risk of HCC (P = 0.10) and the risk was more pronounced in the older age group (age ≥50, P = 0.03).

Discussion

Our findings further the understanding of the genetic basis for HBV-related HCC predisposition in chronic HBV carriers.

Identification of Gray Leaf Spot Disease Candidate Gene in Narrow-Leafed Lupin (Lupinus angustifolius L.)

Selection for resistance against gray leaf spot (GLS) is a major objective in the lupin breeding programs. A segregation ratio of 1:1 (resistant:susceptible) in F8 recombinant inbred lines (RIL8) derived from a cross between a breeding line 83A:476 (resistant to GLS) and a wild accession P27255 (susceptible to GLS) indicated that GLS was controlled by a single major gene. To develop molecular markers linked to GLS, in the beginning, only 11 resistant lines and six susceptible lines from the 83A:476 and P27255 population were genotyped with MFLP markers, and three MFLP markers were identified to be co-segregated with GLS. This method was very efficient, but the markers were located outside of the gene, and could not be used in other germplasms. Then QTL analysis and fine mapping were conducted to identify the gene. Finally, the gene was narrowed down to a 241-kb region containing two disease resistance genes. To further identify the candidate gene, DNA variants between accessions Tanjil (resistant to GLS) and Unicrop (susceptible to GLS) were analyzed. The results indicated that only one SNP was detected in the 241 kb region. This SNP was located in the TMV resistance protein N-like gene region and also identified between 83A:476 and P27255. Genotyping the Tanjil/Unicrop RIL8 population showed that this SNP co-segregated with GLS resistance. The phylogenetic tree analysis of this gene among 18 lupin accessions indicates that Australian resistant breeding line/varieties were clustered into one group and carry two resistant alleles, while susceptible accessions were clustered into different groups.

GmST1, which encodes a sulfotransferase, confers resistance to soybean mosaic virus strains G2 and G3

Soybean mosaic virus (SMV) is one of the most widespread and devastating viral diseases worldwide. The genetic architecture of qualitative resistance to SMV in soybean remains unclear. Here, the Rsvg2 locus was identified as underlying soybean resistance to SMV by genome-wide association and linkage analyses. Fine mapping results showed that soybean resistance to SMV strains G2 and G3 was controlled by a single dominant gene, GmST1, on chromosome 13, encoding a sulfotransferase (SOT). A key variation at position 506 in the coding region of GmST1 associated with the structure of the encoded SOT and changed SOT activity levels between RSVG2-S and RSVG2-R alleles. In RSVG2-S allele carrier "Hefeng25", the overexpression of GmST1 carrying the RSVG2-R allele from the SMV-resistant line "Dongnong93-046" conferred resistance to SMV strains G2 and G3. Compared to Hefeng25, the accumulation of SMV was decreased in transgenic plants carrying the RSVG2-R allele. SMV infection differentiated both the accumulation of jasmonates and expression patterns of genes involved in jasmonic acid (JA) signalling, biosynthesis and catabolism in RSVG2-R and RSVG2-S allele carriers. This characterization of GmST1 suggests a new scenario explaining soybean resistance to SMV.

Fine Mapping of the Gene Controlling the Fruit Skin Hairiness of Prunus persica and Its Uses for MAS in Progenies

The fruit skin pubescence of Prunus persica is an economically important characteristic and comprises the classification criteria. The mapping and identification of a complete linkage marker to the fruit skin trichome trait locus of peach fruit are critical for the molecular marker-assisted selection for peach/nectarine. In this study, the BC₁ population was constructed from the parents “Zhongyou No. 4”, the recurrent parent, and “Baihuashanbitao”, the non-recurrent parent. Based on the 38 BC₁ individuals’ phenotypes and their genotyping using next-generation sequencing, the G (Glabrous skin) locus of the gene was first identified between 14.099 and 16.721 Mb on chromosome 5. Using other individuals of this population, the gene was fine-mapped in the range of 481 kb with SNP markers. Based on the resequencing data of other cultivars (lines), the candidate SNP in the gene Prupe.5G196400 was obtained. Subsequently, the SNP marker was designed and applied to natural and hybrid peach populations. Via genotyping analysis, we confirmed co-segregation between the peach/nectarine phenotype, which was used in the identification of peach or nectarine with 100% accuracy.

Fine Mapping without Phenotyping: Identification of Selection Targets in Secondary Evolve and Resequence Experiments

Evolve and Resequence (E&R) studies investigate the genomic selection response of populations in an Experimental Evolution setup. Despite the popularity of E&R, empirical studies in sexually reproducing organisms typically suffer from an excess of candidate loci due to linkage disequilibrium, and single gene or SNP resolution is the exception rather than the rule. Recently, so-called "secondary E&R" has been suggested as promising experimental follow-up procedure to confirm putatively selected regions from a primary E&R study. Secondary E&R provides also the opportunity to increase mapping resolution by allowing for additional recombination events, which separate the selection target from neutral hitchhikers. Here, we use computer simulations to assess the effect of different crossing schemes, population size, experimental duration, and number of replicates on the power and resolution of secondary E&R. We find that the crossing scheme and population size are crucial factors determining power and resolution of secondary E&R: A simple crossing scheme with few founder lines consistently outcompetes crossing schemes where evolved populations from a primary E&R experiment are mixed with a complex ancestral founder population. Regardless of the experimental design tested, a population size of at least 4,800 individuals, which is roughly five times larger than population sizes in typical E&R studies, is required to achieve a power of at least 75%. Our study provides an important step toward improved experimental designs aiming to characterize causative SNPs in Experimental Evolution studies.

Fine Mapping and Candidate Gene Identification for the CapUp Locus Controlling Fruit Orientation in Pepper (Capsicum spp.)

The orientation of fruits is a distinguishing morphological feature of pepper (Capsicum spp.) varieties. The pendent (downward curved) growth of the fruit stalks, known as pedicels, is highly correlated with fruit weight and pedicel length. A previous genetic analysis revealed that the pendent fruit orientation is governed by a dominant gene, and incomplete inheritance is also observed in some Capsicum accessions. To identify and localize this gene, a single quantitative trait locus (QTL) analysis was performed on one F₂ and two recombinant inbred line (RIL) populations, and a genome-wide association study (GWAS) was performed using a core collection. Common QTL regions associated with fruit orientation were detected on chromosome 12. A total of 187,966 SNPs were identified in a genotyping-by-sequencing (GBS) for GWAS analysis of 196 Capsicum annuum, 25 Capsicum baccatum, 21 Capsicum chinense, and 14 Capsicum frutescens accessions, representing the germplasm collection of South Korea. The results of these analyses enabled us to narrow down the CapUp region of interest to 200-250 Mbp on chromosome 12. Seven candidate genes were found to be located between two markers that were completely cosegregated with the fruit orientation phenotype. The findings and markers developed in this study will be helpful for additional understanding of pepper fruit development and breeding for fruit orientation.

Fine Mapping of the Leaf Rust Resistance Gene Lr65 in Spelt Wheat 'Altgold'

Wheat leaf rust (also known as brown rust), caused by the fungal pathogen Puccinia triticina Erikss. (Pt), is one by far the most troublesome wheat disease worldwide. The exploitation of resistance genes has long been considered as the most effective and sustainable method to control leaf rust in wheat production. Previously the leaf rust resistance gene Lr65 has been mapped to the distal end of chromosome arm 2AS linked to molecular marker Xbarc212. In this study, Lr65 was delimited to a 0.8 cM interval between flanking markers Alt-64 and AltID-11, by employing two larger segregating populations obtained from crosses of the resistant parent Altgold Rotkorn (ARK) with the susceptible parents Xuezao and Chinese Spring (CS), respectively. 24 individuals from 622 F₂ plants of crosses between ARK and CS were obtained that showed the recombination between Lr65 gene and the flanking markers Alt-64 and AltID-11. With the aid of the CS reference genome sequence (IWGSC RefSeq v1.0), one SSR marker was developed between the interval matched to the Lr65-flanking marker and a high-resolution genetic linkage map was constructed. The Lr65 was finally located to a region corresponding to 60.11 Kb of the CS reference genome. The high-resolution genetic linkage map founded a solid foundation for the map-based cloning of Lr65 and the co-segregating marker will facilitate the marker-assisted selection (MAS) of the target gene.