Genetic variation in a tepary bean (Phaseolus acutifolius A. Gray) diversity panel reveals loci associated with biotic stress resistance.

Tepary bean (Phaseolus acutifolius A. Gray), indigenous to the arid climates of northern Mexico and the Southwest United States, diverged from common bean (Phaseolus vulgaris L.), approximately 2 million years ago and exhibits a wide range of resistance to biotic stressors. The tepary genome is highly syntenic to the common bean genome providing a foundation for discovery and breeding of agronomic traits between these two crop species. Although a limited number of adaptive traits from tepary bean have been introgressed into common bean, hybridization barriers between these two species required the development of bridging lines to alleviate this barrier. Thus, to fully utilize the extant tepary bean germplasm as both a crop and as a donor of adaptive traits, we developed a diversity panel of 422 cultivated, weedy, and wild tepary bean accessions which were then genotyped and phenotyped to enable population genetic analyses and genome-wide association studies for their response to a range of biotic stressors. Population structure analyses of the panel revealed eight subpopulations and the differentiation of botanical varieties within P. acutifolius. Genome-wide association studies revealed loci and candidate genes underlying biotic stress resistance including quantitative trait loci for resistance to weevils, common bacterial blight, Fusarium wilt, and bean common mosaic necrosis virus that can be harnessed not only for tepary bean but also common bean improvement.

Leveraging prior biological knowledge improves prediction of tocochromanols in maize grain.

With an essential role in human health, tocochromanols are mostly obtained by consuming seed oils; however, the vitamin E content of the most abundant tocochromanols in maize (Zea mays L.) grain is low. Several large-effect genes with cis-acting variants affecting messenger RNA (mRNA) expression are mostly responsible for tocochromanol variation in maize grain, with other relevant associated quantitative trait loci (QTL) yet to be fully resolved. Leveraging existing genomic and transcriptomic information for maize inbreds could improve prediction when selecting for higher vitamin E content. Here, we first evaluated a multikernel genomic best linear unbiased prediction (MK-GBLUP) approach for modeling known QTL in the prediction of nine tocochromanol grain phenotypes (12-21 QTL per trait) within and between two panels of 1,462 and 242 maize inbred lines. On average, MK-GBLUP models improved predictive abilities by 7.0-13.6% when compared with GBLUP. In a second approach with a subset of 545 lines from the larger panel, the highest average improvement in predictive ability relative to GBLUP was achieved with a multi-trait GBLUP model (15.4%) that had a tocochromanol phenotype and transcript abundances in developing grain for a few large-effect candidate causal genes (1-3 genes per trait) as multiple response variables. Taken together, our study illustrates the enhancement of prediction models when informed by existing biological knowledge pertaining to QTL and candidate causal genes.

Combining GWAS and TWAS to identify candidate causal genes for tocochromanol levels in maize grain.

Tocochromanols (tocopherols and tocotrienols, collectively vitamin E) are lipid-soluble antioxidants important for both plant fitness and human health. The main dietary sources of vitamin E are seed oils that often accumulate high levels of tocopherol isoforms with lower vitamin E activity. The tocochromanol biosynthetic pathway is conserved across plant species but an integrated view of the genes and mechanisms underlying natural variation of tocochromanol levels in seed of most cereal crops remains limited. To address this issue, we utilized the high mapping resolution of the maize Ames panel of ∼1,500 inbred lines scored with 12.2 million single-nucleotide polymorphisms to generate metabolomic (mature grain tocochromanols) and transcriptomic (developing grain) data sets for genetic mapping. By combining results from genome- and transcriptome-wide association studies, we identified a total of 13 candidate causal gene loci, including 5 that had not been previously associated with maize grain tocochromanols: 4 biosynthetic genes (arodeH2 paralog, dxs1, vte5, and vte7) and a plastid S-adenosyl methionine transporter (samt1). Expression quantitative trait locus (eQTL) mapping of these 13 gene loci revealed that they are predominantly regulated by cis-eQTL. Through a joint statistical analysis, we implicated cis-acting variants as responsible for colocalized eQTL and GWAS association signals. Our multiomics approach provided increased statistical power and mapping resolution to enable a detailed characterization of the genetic and regulatory architecture underlying tocochromanol accumulation in maize grain and provided insights for ongoing biofortification efforts to breed and/or engineer vitamin E and antioxidant levels in maize and other cereals.

Genome-wide association identifies a missing hydrolase for tocopherol synthesis in plants.

The tocopherol biosynthetic pathway, encoded by VTE genes 1 through 6, is highly conserved in plants but most large effect quantitative trait loci for seed total tocopherols (totalT) lack VTE genes, indicating other activities are involved. A genome-wide association study of Arabidopsis seed tocopherols showed five of seven significant intervals lacked VTE genes, including the most significant, which mapped to an uncharacterized, seed-specific, envelope-localized, alpha/beta hydrolase with esterase activity, designated AtVTE7. Atvte7 null mutants decreased seed totalT 55% while a leaky allele of the maize ortholog, ZmVTE7, decreased kernel and leaf totalT 38% and 49%, respectively. Overexpressing AtVTE7 or ZmVTE7 partially or fully complemented the Atvte7 seed phenotype and increased leaf totalT by 3.6- and 6.9-fold, respectively. VTE7 has the characteristics of an esterase postulated to provide phytol from chlorophyll degradation for tocopherol synthesis, but bulk chlorophyll levels were unaffected in vte7 mutants and overexpressing lines. Instead, levels of specific chlorophyll biosynthetic intermediates containing partially reduced side chains were impacted and strongly correlated with totalT. These intermediates are generated by a membrane-associated biosynthetic complex containing protochlorophyllide reductase, chlorophyll synthase, geranylgeranyl reductase (GGR) and light harvesting-like 3 protein, all of which are required for both chlorophyll and tocopherol biosynthesis. We propose a model where VTE7 releases prenyl alcohols from chlorophyll biosynthetic intermediates, which are then converted to the corresponding diphosphates for tocopherol biosynthesis.

Phased, chromosome-scale genome assemblies of tetraploid potato reveal a complex genome, transcriptome, and predicted proteome landscape underpinning genetic diversity.

Cultivated potato is a clonally propagated autotetraploid species with a highly heterogeneous genome. Phased assemblies of six cultivars including two chromosome-scale phased genome assemblies revealed extensive allelic diversity, including altered coding and transcript sequences, preferential allele expression, and structural variation that collectively result in a highly complex transcriptome and predicted proteome, which are distributed across the homologous chromosomes. Wild species contribute to the extensive allelic diversity in tetraploid cultivars, demonstrating ancestral introgressions predating modern breeding efforts. As a clonally propagated autotetraploid that undergoes limited meiosis, dysfunctional and deleterious alleles are not purged in tetraploid potato. Nearly a quarter of the loci bore mutations are predicted to have a high negative impact on protein function, complicating breeder's efforts to reduce genetic load. The StCDF1 locus controls maturity, and analysis of six tetraploid genomes revealed that 12 allelic variants of StCDF1 are correlated with maturity in a dosage-dependent manner. Knowledge of the complexity of the tetraploid potato genome with its rampant structural variation and embedded deleterious and dysfunctional alleles will be key not only to implementing precision breeding of tetraploid cultivars but also to the construction of homozygous, diploid potato germplasm containing favorable alleles to capitalize on heterosis in F1 hybrids.

Genomic variation within the maize stiff-stalk heterotic germplasm pool.

The stiff-stalk heterotic group in Maize (Zea mays L.) is an important source of inbreds used in U.S. commercial hybrid production. Founder inbreds B14, B37, B73, and, to a lesser extent, B84, are found in the pedigrees of a majority of commercial seed parent inbred lines. We created high-quality genome assemblies of B84 and four expired Plant Variety Protection (ex-PVP) lines LH145 representing B14, NKH8431 of mixed descent, PHB47 representing B37, and PHJ40, which is a Pioneer Hi-Bred International (PHI) early stiff-stalk type. Sequence was generated using long-read sequencing achieving highly contiguous assemblies of 2.13-2.18 Gbp with N50 scaffold lengths >200 Mbp. Inbred-specific gene annotations were generated using a core five-tissue gene expression atlas, whereas transposable element (TE) annotation was conducted using de novo and homology-directed methodologies. Compared with the reference inbred B73, synteny analyses revealed extensive collinearity across the five stiff-stalk genomes, although unique components of the maize pangenome were detected. Comparison of this set of stiff-stalk inbreds with the original Iowa Stiff Stalk Synthetic breeding population revealed that these inbreds represent only a proportion of variation in the original stiff-stalk pool and there are highly conserved haplotypes in released public and ex-Plant Variety Protection inbreds. Despite the reduction in variation from the original stiff-stalk population, substantial genetic and genomic variation was identified supporting the potential for continued breeding success in this pool. The assemblies described here represent stiff-stalk inbreds that have historical and commercial relevance and provide further insight into the emerging maize pangenome.

Composite modeling of leaf shape along shoots discriminates Vitis species better than individual leaves.

PREMISE: Leaf morphology is dynamic, continuously deforming during leaf expansion and among leaves within a shoot. Here, we measured the leaf morphology of more than 200 grapevines (Vitis spp.) over four years and modeled changes in leaf shape along the shoot to determine whether a composite leaf shape comprising all the leaves from a single shoot can better capture the variation and predict species identity compared with individual leaves. METHODS: Using homologous universal landmarks found in grapevine leaves, we modeled various morphological features as polynomial functions of leaf nodes. The resulting functions were used to reconstruct modeled leaf shapes across the shoots, generating composite leaves that comprehensively capture the spectrum of leaf morphologies present. RESULTS: We found that composite leaves are better predictors of species identity than individual leaves from the same plant. We were able to use composite leaves to predict the species identity of previously unassigned grapevines, which were verified with genotyping. DISCUSSION: Observations of individual leaf shape fail to capture the true diversity between species. Composite leaf shape-an assemblage of modeled leaf snapshots across the shoot-is a better representation of the dynamic and essential shapes of leaves, in addition to serving as a better predictor of species identity than individual leaves.

Genome sequencing of four culinary herbs reveals terpenoid genes underlying chemodiversity in the Nepetoideae.

Species within the mint family, Lamiaceae, are widely used for their culinary, cultural, and medicinal properties due to production of a wide variety of specialized metabolites, especially terpenoids. To further our understanding of genome diversity in the Lamiaceae and to provide a resource for mining biochemical pathways, we generated high-quality genome assemblies of four economically important culinary herbs, namely, sweet basil (Ocimum basilicum L.), sweet marjoram (Origanum majorana L.), oregano (Origanum vulgare L.), and rosemary (Rosmarinus officinalis L.), and characterized their terpenoid diversity through metabolite profiling and genomic analyses. A total 25 monoterpenes and 11 sesquiterpenes were identified in leaf tissue from the 4 species. Genes encoding enzymes responsible for the biosynthesis of precursors for mono- and sesqui-terpene synthases were identified in all four species. Across all 4 species, a total of 235 terpene synthases were identified, ranging from 27 in O. majorana to 137 in the tetraploid O. basilicum. This study provides valuable resources for further investigation of the genetic basis of chemodiversity in these important culinary herbs.

QTL mapping of post-processing color retention in two black bean populations.

KEY MESSAGE: Several QTL governing color retention in processed black beans were identified by traditional and novel phenotyping methods applied to two black bean mapping populations. When black beans are hydrothermally processed prior to consumption, water-soluble anthocyanins are released from the seed coat, resulting in an undesirable faded brown color in the cooked product. The aim of this research was to develop mapping populations with different genetic sources of color retention in order to identify regions of the bean genome associated with canning quality traits. Two half-sibling black bean recombinant inbred line (RIL) populations segregating for post-processing color retention were developed. These RIL populations were phenotyped for canning quality traits over two years and genotyped using the BARCBean6k_3 BeadChip. In addition to traditional phenotyping by trained panelists, cooked beans were also phenotyped using a novel digital image analysis pipeline. Measurements of post-processing seed coat color from both phenotyping methods were compared, and the digital image analysis was shown to outperform the trained panelists. Quantitative trait loci (QTL) for post-processing color retention were detected on six chromosomes, with QTL on Pv08 and Pv11 consistently detected across phenotyping methods, populations, and years. Color retention QTL on Pv08 explained up to 32% of phenotypic variation but were significant over a large physical interval due to low SNP marker coverage. However, color retention QTL on Pv11 also explained a substantial amount of phenotypic variation (r2 ≈ 25%) and mapped to a small genomic region near 52.5 Mbp. The QTL and methods described in this study will be useful for dry bean breeders and food scientists to produce high quality black beans that meet consumer needs.