Major quantitative trait locus qLA3.1 is related to tomato leaf angle by regulating cell length at the petiole base.

KEY MESSAGE: qLA3.1, controlling leaf angle in tomato, was fine-mapped to an interval of 4.45 kb on chromosome A03, and one gene encoding auxin response factor was identified as a candidate gene. Leaf angle is a crucial trait in plant architecture that plays an important role in achieving optimal plant structure. However, there are limited reports on gene localization, cloning, and the function of plant architecture in horticultural crops, particularly regarding leaf angle. In this study, we selected 'Z3' with erect leaves and 'Heinz1706' with horizontal leaves as the phenotype and cytological observation. We combined bulked segregant analysis and fine genetic mapping to identify a candidate gene, known as, i.e., qLA3.1, which was related to tomato leaf angle. Through multiple analyses, we found that Solyc03g113410 was the most probably candidate for qLA3.1, which encoded the auxin response factor SlARF11 in tomato and was homologous to OsARF11 related to leaf angle in rice. We discovered that silencing SlARF11 resulted in upright leaves, while plants with over-expressed SlARF11 exhibited horizontal leaves. We also found that cultivars with erect leaves had a mutation from base G to base A. Moreover, quantitative analysis of plants treated with hormones indicated that SlARF11 might participate in cell elongation and the activation of genes related to auxin and brassinosteroid pathways. Transcriptome analysis further validated that SlARF11 may regulate leaf angle through hormone signaling pathways. These data support the idea that the auxin response factor SlARF11 may have an important function in tomato leaf petiole angles.

QTL analysis on a lemon population provides novel insights on the genetic regulation of the tolerance to the two-spotted spider mite attack.

BACKGROUND: Among the Citrus species, lemon (Citrus limon Burm f.) is one of the most affected by the two-spotted spider mite (Tetranychus urticae Koch). Moreover, chemical control is hampered by the mite's ability to develop genetic resistance against acaricides. In this context, the identification of the genetic basis of the host resistance could represent a sustainable strategy for spider mite control. In the present study, a marker-trait association analysis was performed on a lemon population employing an association mapping approach. An inter-specific full-sib population composed of 109 accessions was phenotyped through a detached-leaf assays performed in modified Huffaker cells. Those individuals, complemented with two inter-specific segregating populations, were genotyped using a target-sequencing approach called SPET (Single Primer Enrichment Technology), the resulting SNPs were employed for the generation of an integrated genetic map. RESULTS: The percentage of damaged area in the full-sib population showed a quantitative distribution with values ranging from 0.36 to 9.67%. A total of 47,298 SNPs were selected for an association mapping study and a significant marker linked with resistance to spider mite was detected on linkage group 5. In silico gene annotation of the QTL interval enabled the detection of 13 genes involved in immune response to biotic and abiotic stress. Gene expression analysis showed an over expression of the gene encoding for the ethylene-responsive transcription factor ERF098-like, already characterized in Arabidopsis and in rice for its involvement in defense response. CONCLUSION: The identification of a molecular marker linked to the resistance to spider mite attack can pave the way for the development of marker-assisted breeding plan for the development of novel selection coupling favorable agronomical traits (e.g. fruit quality, yield) with a higher resistance toward the mite.

Validation and identification of promising gene specific markers governing foliar disease resistance in groundnut (Arachis hypogaea L.).

BACKGROUND: Groundnut is vulnerable to the major foliar fungal disease viz., late leaf spot (LLS) and rust in kharif season, which results in severe yield losses. Until now, LLS and rust resistance linked markers were developed based on GPBD 4 as a major donor source and were validated in its derivatives only, which restricted their use in marker assisted selection (MAS) involving other donors. METHODS AND RESULTS: The current study focused to validate LLS and rust resistance linked markers employing advanced breeding lines of F6 generation, derived from nine different crosses involving nine diverse parents, to identify potential markers for marker-assisted breeding of LLS and rust resistance in groundnut. Out of 28-trait linked markers used for validation, 8 were polymorphic (28.57%). Marker-trait association (MTA) and Single Marker Analysis (SMA) revealed that the SSR marker pPGPseq5D05 is significantly associated with both LLS (15.8% PVE) and rust (17.5% PVE) resistance, whereas, the marker IPAHM103 is tightly linked with rust resistance (26.8% PVE) alone. In silico analysis revealed that the marker gene for IPAHM103 is a zinc finger protein and the marker gene for pPGPseq5D05 is an ADP-ribosylation factor GTPase-activating protein. Both these protein products impart resistance or tolerance to biotic stress in crop plants. Two other markers namely, GMLQ975 and pPGPseq13A10 were also found to be associated with LLS resistance explaining MTA up to 60%. CONCLUSION: These gene specific markers will enable us to screen more number of germplasm lines or newly developed lines in MAS schemes for LLS and rust resistance using a wide range of resistant sources.

High-density mapping of durable and broad-spectrum stripe rust resistance gene Yr30 in wheat.

KEY MESSAGE: The durable stripe rust resistance gene Yr30 was fine-mapped to a 610-kb region in which five candidate genes were identified by expression analysis and sequence polymorphisms. The emergence of genetically diverse and more aggressive races of Puccinia striiformis f. sp. tritici (Pst) in the past twenty years has resulted in global stripe rust outbreaks and the rapid breakdown of resistance genes. Yr30 is an adult plant resistance (APR) gene with broad-spectrum effectiveness and its durability. Here, we fine-mapped the YR30 locus to a 0.52-cM interval using 1629 individuals derived from residual heterozygous F5:6 plants in a Yaco"S"/Mingxian169 recombinant inbred line population. This interval corresponded to a 610-kb region in the International Wheat Genome Sequencing Consortium (IWGSC) RefSeq version 2.1 on chromosome arm 3BS harboring 30 high-confidence genes. Five genes were identified as candidate genes based on functional annotation, expression analysis by RNA-seq and sequence polymorphisms between cultivars with and without Yr30 based on resequencing. Haplotype analysis of the target region identified six haplotypes (YR30_h1-YR30_h6) in a panel of 1215 wheat accessions based on the 660K feature genotyping array. Lines with YR30_h6 displayed more resistance to stripe rust than the other five haplotypes. Near-isogenic lines (NILs) with Yr30 showed a 32.94% higher grain yield than susceptible counterparts when grown in a stripe rust nursery, whereas there was no difference in grain yield under rust-free conditions. These results lay a foundation for map-based cloning Yr30.

A high-density linkage map and sex-determination loci in Pacific white shrimp (Litopenaeus vannamei).

BACKGROUND: Expansion of genomic resources for the Pacific white shrimp (Litopenaeus vannamei), such as the construction of dense genetic linkage maps, is crucial for the application of genomic tools in order to improve economically relevant traits. Sexual dimorphism exists in Pacific white shrimp, and the mapping of the sex-determination region in this species may help in future reproductive applications. We have constructed male, female, and sex-averaged high-density genetic maps using a 50 K single-nucleotide polymorphism (SNP) array, followed by a genome-wide association study (GWAS) to identify genomic regions associated with sex in white shrimp. RESULTS: The genetic map yielded 15,256 SNPs assigned to 44 linkage groups (LG). The lengths of the male, female, and sex-averaged maps were 5,741.36, 5,461.20 and 5,525.26 cM, respectively. LG18 was found to be the largest for both sexes, whereas LG44 was the shortest for males and LG31 for females. A sex-determining region was found in LG31 with 21 statistically significant SNPs. The most important SNP was previously identified as a sex-linked marker and was able to identify 99% of the males and 88% of the females. Although other significant markers had a lower ability to determine sex, putative genes were intercepted or close to them. The oplophorus-luciferin 2-monooxygenase, serine/arginine repetitive matrix protein and spermine oxidase genes were identified as candidates with possible participation in important processes of sexual differentiation in shrimp. CONCLUSIONS: Our results provide novel genomic resources for shrimp, including a high-density linkage map and new insights into the sex-determining region in L. vannamei, which may be usefulfor future genetics and reproduction applications.

Using expression data to fine map QTL associated with fertility in dairy cattle.

BACKGROUND: Female fertility is an important trait in dairy cattle. Identifying putative causal variants associated with fertility may help to improve the accuracy of genomic prediction of fertility. Combining expression data (eQTL) of genes, exons, gene splicing and allele specific expression is a promising approach to fine map QTL to get closer to the causal mutations. Another approach is to identify genomic differences between cows selected for high and low fertility and a selection experiment in New Zealand has created exactly this resource. Our objective was to combine multiple types of expression data, fertility traits and allele frequency in high- (POS) and low-fertility (NEG) cows with a genome-wide association study (GWAS) on calving interval in Australian cows to fine-map QTL associated with fertility in both Australia and New Zealand dairy cattle populations. RESULTS: Variants that were significantly associated with calving interval (CI) were strongly enriched for variants associated with gene, exon, gene splicing and allele-specific expression, indicating that there is substantial overlap between QTL associated with CI and eQTL. We identified 671 genes with significant differential expression between POS and NEG cows, with the largest fold change detected for the CCDC196 gene on chromosome 10. Our results provide numerous candidate genes associated with female fertility in dairy cattle, including GYS2 and TIGAR on chromosome 5 and SYT3 and HSD17B14 on chromosome 18. Multiple QTL regions were located in regions with large numbers of copy number variants (CNV). To identify the causal mutations for these variants, long read sequencing may be useful. CONCLUSIONS: Variants that were significantly associated with CI were highly enriched for eQTL. We detected 671 genes that were differentially expressed between POS and NEG cows. Several QTL detected for CI overlapped with eQTL, providing candidate genes for fertility in dairy cattle.

Meta-QTL analysis and identification of candidate genes for multiple-traits associated with spot blotch resistance in bread wheat.

In bread wheat, a literature search gave 228 QTLs for six traits, including resistance against spot blotch and the following five other related traits: (i) stay green; (ii) flag leaf senescence; (iii) green leaf area duration; (iv) green leaf area of the main stem; and (v) black point resistance. These QTLs were used for metaQTL (MQTL) analysis. For this purpose, a consensus map with 72,788 markers was prepared; 69 of the above 228 QTLs, which were suitable for MQTL analysis, were projected on the consensus map. This exercise resulted in the identification of 16 meta-QTLs (MQTLs) located on 11 chromosomes, with the PVE ranging from 5.4% (MQTL7) to 21.8% (MQTL5), and the confidence intervals ranging from 1.5 to 20.7 cM (except five MQTLs with a range of 36.1-57.8 cM). The number of QTLs associated with individual MQTLs ranged from a maximum of 17 in MQTL3 to 8 each in MQTL5 and MQTL8 and 5 each in MQTL7 and MQTL14. The 16 MQTLs, included 12 multi-trait MQTLs; one of the MQTL also overlapped a genomic region carrying the major spot blotch resistance gene Sb1. Of the total 16 MQTLs, 12 MQTLs were also validated through marker-trait associations that were available from earlier genome-wide association studies. The genomic regions associated with MQTLs were also used for the identification of candidate genes (CGs) and led to the identification of 516 CGs encoding 508 proteins; 411 of these proteins are known to be associated with resistance against several biotic stresses. In silico expression analysis of CGs using transcriptome data allowed the identification of 71 differentially expressed CGs, which were examined for further possible studies. The findings of the present study should facilitate fine-mapping and cloning of genes, enabling Marker Assisted Selection.

Cytological and molecular characterization of wheat lines carrying leaf rust and stem rust resistance genes Lr24 and Sr24.

Previous studies showed that Australian wheat cultivars Janz and Sunco carry leaf rust and stem rust resistance genes Lr24 and Sr24 derived from Thinopyrum ponticum chromosome arm 3AgL. However, the size of the alien segments carrying Lr24 and Sr24 in the lines were not determined. In this study, we used non-denaturing fluorescence in situ hybridization (ND-FISH), genomic in situ hybridization (GISH), and PCR-based landmark unique gene (PLUG) markers to visualize the alien segments in Janz and Sunco, and further compared them with the segments in US cultivars Agent and Amigo. The fraction length (FL) of the alien translocation in Agent was 0.70-1.00, whereas those in Janz, Sunco, and Amigo were smaller, at FL 0.85-1.00. It was deduced that the alien gene RAg encoding for red grain color and rust resistance genes Lr24 and Sr24 on chromosome arm 3AgL were in bins of FL 0.70-0.85 and 0.85-1.00, respectively. We retrieved and extracted nucleotide-binding site-leucine-rich repeat (NBS-LRR) receptor genes corresponding to the region of Lr24 and Sr24 on chromosomes 3E, and 3J, 3Js and 3St from the reference genome sequences of Th. elongatum and Th. intermedium, respectively. A set of molecular markers developed for Lr24 and Sr24 from those extracted NBS-LRR genes will provide valuable information for fine mapping and cloning of these genes.

Identification and validation of two major QTL for grain number per spike on chromosomes 2B and 2D in bread wheat (Triticum aestivum L.).

KEY MESSAGE: Major QTL for grain number per spike were identified on chromosomes 2B and 2D. Haplotypes and candidate genes of QGns.cib-2B.1 were analyzed. Grain number per spike (GNS) is one of the main components of wheat yield. Genetic dissection of their regulatory factors is essential to improve the yield potential. In present study, a recombinant inbred line population comprising 180 lines developed from the cross between a high GNS line W7268 and a cultivar Chuanyu12 was employed to identify quantitative trait loci (QTL) associated with GNS across six environments. Two major QTL, QGns.cib-2B.1 and QGns.cib-2D.1, were detected in at least four environments with the phenotypic variations of 12.99-27.07% and 8.50-13.79%, respectively. And significant interactions were observed between the two major QTL. In addition, QGns.cib-2B.1 is a QTL cluster for GNS, grain number per spikelet and fertile tiller number, and they were validated in different genetic backgrounds using Kompetitive Allele Specific PCR (KASP) markers. QGns.cib-2B.1 showed pleotropic effects on other yield-related traits including plant height, spike length, and spikelet number per spike, but did not significantly affect thousand grain weight which suggested that it might be potentially applicable in breeding program. Comparison analysis suggested that QGns.cib-2B.1 might be a novel QTL. Furthermore, haplotype analysis of QGns.cib-2B.1 indicated that it is a hot spot of artificial selection during wheat improvement. Based on the expression patterns, gene annotation, orthologs analysis and sequence variations, the candidate genes of QGns.cib-2B.1 were predicted. Collectively, the major QTL and KASP markers reported here provided a wealth of information for the genetic basis of GNS and grain yield improvement.

Molecular characterization of QTL for grain zinc and iron concentrations in wheat landrace Chinese Spring.

KEY MESSAGE: Three stable QTL for grain zinc concentration were identified in wheat landrace Chinese Spring. Favorable alleles were more frequent in landraces than in modern wheat cultivars. Wheat is a major source of dietary energy for the growing world population. Developing cultivars with enriched zinc and iron can potentially alleviate human micronutrient deficiency. In this study, a recombinant inbred line (RIL) population with 245 lines derived from cross Zhou 8425B/Chinese Spring was used to detect quantitative trait loci (QTL) for grain zinc concentration (GZnC) and grain iron concentration (GFeC) across four environments. Three stable QTL for GZnC with all favorable alleles from Chinese Spring were identified on chromosomes 3BL, 5AL, and 5BL. These QTL explaining maxima of 8.7%, 5.8%, and 7.1% of phenotypic variances were validated in 125 resequenced wheat accessions encompassing both landraces and modern cultivars using six kompetitive allele specific PCR (KASP) assays. The frequencies of favorable alleles for QGZnCzc.caas-3BL, QGZnCzc.caas-5AL and QGZnCzc.caas-5BL were higher in landraces (90.4%, 68.0%, and 100.0%, respectively) compared to modern cultivars (45.9%, 35.4%, and 40.9%), suggesting they were not selected in breeding programs. Candidate gene association studies on GZnC in the cultivar panel further delimited the QTL into 8.5 Mb, 4.1 Mb, and 47.8 Mb regions containing 46, 4, and 199 candidate genes, respectively. The 5BL QTL located in a region where recombination was suppressed. Two stable and three less stable QTL for GFeC with favorable alleles also from Chinese Spring were identified on chromosomes 4BS (Rht-B1a), 4DS (Rht-D1a), 1DS, 3AS, and 6DS. This study sheds light on the genetic basis of GZnC and GFeC in Chinese Spring and provides useful molecular markers for wheat biofortification.

The wall-associated kinase GWN1 controls grain weight and grain number in rice.

Grain size is a crucial agronomic trait that determines grain weight and final yield. Although several genes have been reported to regulate grain size in rice (Oryza sativa), the function of Wall-Associated Kinase family genes affecting grain size is still largely unknown. In this study, we identified GRAIN WEIGHT AND NUMBER 1 (GWN1) using map-based cloning. GWN1 encodes the OsWAK74 protein kinase, which is conserved in plants. GWN1 negatively regulates grain length and weight by regulating cell proliferation in spikelet hulls. We also found that GWN1 negatively influenced grain number by influencing secondary branch numbers and finally increased plant grain yield. The GWN1 gene was highly expressed in inflorescences and its encoded protein is located at the cell membrane and cell wall. Moreover, we identified three haplotypes of GWN1 in the germplasm. GWN1hap1 showing longer grain, has not been widely utilized in modern rice varieties. In summary, GWN1 played a very important role in regulating grain length, weight and number, thereby exhibiting application potential in molecular breeding for longer grain and higher yield.

RicePilaf: a post-GWAS/QTL dashboard to integrate pangenomic, coexpression, regulatory, epigenomic, ontology, pathway, and text-mining information to provide functional insights into rice QTLs and GWAS loci.

BACKGROUND: As the number of genome-wide association study (GWAS) and quantitative trait locus (QTL) mappings in rice continues to grow, so does the already long list of genomic loci associated with important agronomic traits. Typically, loci implicated by GWAS/QTL analysis contain tens to hundreds to thousands of single-nucleotide polmorphisms (SNPs)/genes, not all of which are causal and many of which are in noncoding regions. Unraveling the biological mechanisms that tie the GWAS regions and QTLs to the trait of interest is challenging, especially since it requires collating functional genomics information about the loci from multiple, disparate data sources. RESULTS: We present RicePilaf, a web app for post-GWAS/QTL analysis, that performs a slew of novel bioinformatics analyses to cross-reference GWAS results and QTL mappings with a host of publicly available rice databases. In particular, it integrates (i) pangenomic information from high-quality genome builds of multiple rice varieties, (ii) coexpression information from genome-scale coexpression networks, (iii) ontology and pathway information, (iv) regulatory information from rice transcription factor databases, (v) epigenomic information from multiple high-throughput epigenetic experiments, and (vi) text-mining information extracted from scientific abstracts linking genes and traits. We demonstrate the utility of RicePilaf by applying it to analyze GWAS peaks of preharvest sprouting and genes underlying yield-under-drought QTLs. CONCLUSIONS: RicePilaf enables rice scientists and breeders to shed functional light on their GWAS regions and QTLs, and it provides them with a means to prioritize SNPs/genes for further experiments. The source code, a Docker image, and a demo version of RicePilaf are publicly available at https://github.com/bioinfodlsu/rice-pilaf.

Combining genome-wide association study and linkage mapping in the genetic dissection of amylose content in maize (Zea mays L.).

KEY MESSAGE: Integrated linkage and association analysis revealed genetic basis across multiple environments. The genes Zm00001d003102 and Zm00001d015905 were further verified to influence amylose content using gene-based association study. Maize kernel amylose is an important source of human food and industrial raw material. However, the genetic basis underlying maize amylose content is still obscure. Herein, we used an intermated B73 × Mo17 (IBM) Syn10 doubled haploid population composed of 222 lines and a germplasm set including 305 inbred lines to uncover the genetic control for amylose content under four environments. Linkage mapping detected 16 unique QTL, among which four were individually repeatedly identified across multiple environments. Genome-wide association study revealed 17 significant (P = 2.24E-06) single-nucleotide polymorphisms, of which two (SYN19568 and PZE-105090500) were located in the intervals of the mapped QTL (qAC2 and qAC5-3), respectively. According to the two population co-localized loci, 20 genes were confirmed as the candidate genes for amylose content. Gene-based association analysis indicated that the variants in Zm00001d003102 (Beta-16-galactosyltransferase GALT29A) and Zm00001d015905 (Sugar transporter 4a) affected amylose content across multi-environment. Tissue expression analysis showed that the two genes were specifically highly expressed in the ear and stem, respectively, suggesting that they might participate in sugar transport from source to sink organs. Our study provides valuable genetic information for breeding maize varieties with high amylose.

Improved functional mapping of complex trait heritability with GSA-MiXeR implicates biologically specific gene sets.

While genome-wide association studies are increasingly successful in discovering genomic loci associated with complex human traits and disorders, the biological interpretation of these findings remains challenging. Here we developed the GSA-MiXeR analytical tool for gene set analysis (GSA), which fits a model for the heritability of individual genes, accounting for linkage disequilibrium across variants and allowing the quantification of partitioned heritability and fold enrichment for small gene sets. We validated the method using extensive simulations and sensitivity analyses. When applied to a diverse selection of complex traits and disorders, including schizophrenia, GSA-MiXeR prioritizes gene sets with greater biological specificity compared to standard GSA approaches, implicating voltage-gated calcium channel function and dopaminergic signaling for schizophrenia. Such biologically relevant gene sets, often with fewer than ten genes, are more likely to provide insights into the pathobiology of complex diseases and highlight potential drug targets.

Multi-trait association mapping for phosphorous efficiency reveals flexible root architectures in sorghum.

BACKGROUND: On tropical regions, phosphorus (P) fixation onto aluminum and iron oxides in soil clays restricts P diffusion from the soil to the root surface, limiting crop yields. While increased root surface area favors P uptake under low-P availability, the relationship between the three-dimensional arrangement of the root system and P efficiency remains elusive. Here, we simultaneously assessed allelic effects of loci associated with a variety of root and P efficiency traits, in addition to grain yield under low-P availability, using multi-trait genome-wide association. We also set out to establish the relationship between root architectural traits assessed in hydroponics and in a low-P soil. Our goal was to better understand the influence of root morphology and architecture in sorghum performance under low-P availability. RESULT: In general, the same alleles of associated SNPs increased root and P efficiency traits including grain yield in a low-P soil. We found that sorghum P efficiency relies on pleiotropic loci affecting root traits, which enhance grain yield under low-P availability. Root systems with enhanced surface area stemming from lateral root proliferation mostly up to 40 cm soil depth are important for sorghum adaptation to low-P soils, indicating that differences in root morphology leading to enhanced P uptake occur exactly in the soil layer where P is found at the highest concentration. CONCLUSION: Integrated QTLs detected in different mapping populations now provide a comprehensive molecular genetic framework for P efficiency studies in sorghum. This indicated extensive conservation of P efficiency QTL across populations and emphasized the terminal portion of chromosome 3 as an important region for P efficiency in sorghum. Increases in root surface area via enhancement of lateral root development is a relevant trait for sorghum low-P soil adaptation, impacting the overall architecture of the sorghum root system. In turn, particularly concerning the critical trait for water and nutrient uptake, root surface area, root system development in deeper soil layers does not occur at the expense of shallow rooting, which may be a key reason leading to the distinctive sorghum adaptation to tropical soils with multiple abiotic stresses including low P availability and drought.

Two genes encoding caffeoyl coenzyme A O-methyltransferase 1 (CCoAOMT1) are candidate genes for physical seed dormancy in cowpea (Vigna unguiculata (L.) Walp.).

KEY MESSAGE: The major QTL Sdp1.1+ controlling seed dormancy in cowpea was finely mapped, and two CCoAOMT1 genes were identified as candidate genes for the dormancy. Seed dormancy in wild cowpea may be useful in breeding cultivated cowpea with pre-harvest sprouting resistance. A previous study identified a major quantitative trait locus (QTL) for seed dormancy, Sdp1.1+ , using the population of the cross between cultivated cowpea 'JP81610' and wild cowpea 'JP89083.' However, the molecular basis of seed dormancy in cowpea is not yet known. In this study, we aimed to finely map the locus Sdp1.1+ and identify candidate gene(s) for it. Germination tests demonstrated that the seed coat is the major factor controlling seed dormancy in the wild cowpea JP89083. Microscopic observations revealed that wild cowpea seeds, unlike cultivated cowpea seeds, possessed a palisade cuticle layer. Fine mapping using a large F2 population of the cross JP81610 × JP89083 grown in Thailand revealed a single QTL, Sdp1.1+ , controlling seed dormancy. The Sdp1.1+ was confirmed using a small F2 population of the same cross grown in Japan. The Sdp1.1+ was mapped to a 37.34-Kb region containing three genes. Two closely linked genes, Vigun03g278900 (VuCCoAOMT1a) and Vigun03g290000 (VuCCoAOMT1b), located 4.844 Kb apart were considered as candidate genes for seed dormancy. The two genes encoded caffeoyl coenzyme A O-methyltransferase 1 (CCoAOMT1). DNA sequencing and alignment of VuCCoAOMT1a and VuCCoAOMT1b between JP89083 and JP81610 revealed a single nucleotide polymorphism (SNP) causing an amino acid change in VuCCoAOMT1a and several SNPs leading to six amino acid changes in VuCCoAOMT1b. Altogether, these results indicate that VuCCoAOMT1a and VuCCoAOMT1b are candidate genes controlling physical seed dormancy in the wild cowpea JP89083.

Mapping and annotating genomic loci to prioritize genes and implicate distinct polygenic adaptations for skin color.

Evidence for adaptation of human skin color to regional ultraviolet radiation suggests shared and distinct genetic variants across populations. However, skin color evolution and genetics in East Asians are understudied. We quantified skin color in 48,433 East Asians using image analysis and identified associated genetic variants and potential causal genes for skin color as well as their polygenic interplay with sun exposure. This genome-wide association study (GWAS) identified 12 known and 11 previously unreported loci and SNP-based heritability was 23-24%. Potential causal genes were determined through the identification of nonsynonymous variants, colocalization with gene expression in skin tissues, and expression levels in melanocytes. Genomic loci associated with pigmentation in East Asians substantially diverged from European populations, and we detected signatures of polygenic adaptation. This large GWAS for objectively quantified skin color in an East Asian population improves understanding of the genetic architecture and polygenic adaptation of skin color and prioritizes potential causal genes.

Diversity in bread and durum wheat stigma morphology and linkage of increased stigma length to dwarfing gene Rht14.

KEY MESSAGE: The dwarfing allele Rht14 of durum wheat associates with greater stigma length, an important trait for hybrid breeding, whilst major dwarfing alleles Rht-B1b and Rht-D1b showed little to no effect. Although much understudied in wheat, the stigma is a crucial component for attaining grain set, the fundamental basis for yield, particularly in hybrid production systems where successful grain set relies on wind-driven pollen dispersal by the male parent and effective pollen capture by the female parent. Females with long stigma that exsert early are thought to be advantageous. Using glasshouse-grown lines, we examined variation in Total Stigma Length (TSL) across diverse panels comprising 27 durum and 116 bread wheat genotypes. Contrasting genotypes were selected for population development and genetic analysis. Quantitative trait loci (QTL) analysis was performed on a durum F2 population and a bread wheat recombinant inbred line (RIL) population. Contrasting with studies of anther length, we found no large effect on TSL of the GA-insensitive semi-dwarfing genes Rht-B1 and Rht-D1 in either durum or bread wheat. However, in durum cultivar Italo, we identified a region on chromosome 6A which is robustly associated with larger TSL and contains the Rht14 allele for reduced plant height, a trait that is favourable for female line development in hybrid systems. This dual effect locus explained 25.2 and 19.2% of TSL phenotypic variation in experiments across two growing seasons, with preliminary results suggesting this locus may increase TSL when transferred to bread wheat. In a bread wheat, RIL population minor QTL on 1A and 2A was indicated, but the strongest association was with Ppd-B1. Methods developed here, and the identification of a TSL-enhancing locus provides advances and further opportunities in the study of wheat stigma.

Identification of CaPCR1, an OFP gene likely involved in pointed versus concave fruit tip regulation in pepper (Capsicum annuum L.) using recombinant inbred lines.

KEY MESSAGE: CaPCR1 (Capana12g002165) was a candidate gene regulating fruit concave/pointed tip shape in pepper. The concave shape of the fruit tip in pepper plants is highly susceptible to drought and low temperature stresses, resulting in the appearance of a pointed tip fruit, which affects its commercial value. However, few studies on the process of fruit tip development and regulatory genes in pepper have been reported. Herein, the developmental process of the ovary before anthesis, especially changes in the shape of the ovary tip, was studied in detail. The results showed that the final fruit tip shape was consistent with the ovary tip shape before anthesis, and a concave tip shape gradually developed. F4 recombinant inbred lines (RILs) were constructed to map the genes regulating fruit tip shape through hybridization of the LRS and SBS pepper inbred lines. CaPCR1 (Capana12g002165), an OFP (OVATE Family Protein) family gene, was located in the candidate region on chr12. Three SNPs were found in the protein coding sequence of CaPCR1 between SBS and LRS, but only one SNP led to amino acid variation. Sequence variations, including base replacements, deletions and insertions, were also detected in the gene promoter region. The relative expression level of the CaPCR1 gene was significantly greater in the concave tip ovary than in the pointed tip ovary. qRT‒PCR analysis revealed that the CaPCR1 gene was expressed mainly in the gynoecium, placenta and green fruit pericarp, which was consistent with its function in ovary and fruit development. Taken together, these results suggested that CaPCR1 is a candidate gene involved in fruit tip shape determination in pepper.