Legume-wide comparative analysis of pod shatter locus PDH1 reveals phaseoloid specificity, high cowpea expression, and stress responsive genomic context.

Pod dehiscence is a major source of yield loss in legumes, which is exacerbated by aridity. Disruptive mutations in "Pod indehiscent 1" (PDH1), a pod sclerenchyma-specific lignin biosynthesis gene, has been linked to significant reductions in dehiscence in several legume species. We compared syntenic PDH1 regions across 12 legumes and two outgroups to uncover key historical evolutionary trends at this important locus. Our results clarified the extent to which PDH1 orthologs are present in legumes, showing the typical genomic context surrounding PDH1 has only arisen relatively recently in certain phaseoloid species (Vigna, Phaseolus, Glycine). The notable absence of PDH1 in Cajanus cajan may be a major contributor to its indehiscent phenotype compared with other phaseoloids. In addition, we identified a novel PDH1 ortholog in Vigna angularis and detected remarkable increases in PDH1 transcript abundance during Vigna unguiculata pod development. Investigation of the shared genomic context of PDH1 revealed it lies in a hotspot of transcription factors and signaling gene families that respond to abscisic acid and drought stress, which we hypothesize may be an additional factor influencing expression of PDH1 under specific environmental conditions. Our findings provide key insights into the evolutionary history of PDH1 and lay the foundation for optimizing the pod dehiscence role of PDH1 in major and understudied legume species.

The complex relationship between disease resistance and yield in crops.

In plants, growth and defence are controlled by many molecular pathways that are antagonistic to one another. This results in a 'growth-defence trade-off', where plants temporarily reduce growth in response to pests or diseases. Due to this antagonism, genetic variants that improve resistance often reduce growth and vice versa. Therefore, in natural populations, the most disease resistant individuals are often the slowest growing. In crops, slow growth may translate into a yield penalty, but resistance is essential for protecting yield in the presence of disease. Therefore, plant breeders must balance these traits to ensure optimal yield potential and yield stability. In crops, both qualitative and quantitative disease resistance are often linked with genetic variants that cause yield penalties, but this is not always the case. Furthermore, both crop yield and disease resistance are complex traits influenced by many aspects of the plant's physiology, morphology and environment, and the relationship between the molecular growth-defence trade-off and disease resistance-yield antagonism is not well-understood. In this article, we highlight research from the last 2 years on the molecular mechanistic basis of the antagonism between defence and growth. We then discuss the interaction between disease resistance and crop yield from a breeding perspective, outlining the complexity and nuances of this relationship and where research can aid practical methods for simultaneous improvement of yield potential and disease resistance.

crosshap: R package for local haplotype visualization for trait association analysis.

SUMMARY: Genome-wide association studies (GWAS) excels at harnessing dense genomic variant datasets to identify candidate regions responsible for producing a given phenotype. However, GWAS and traditional fine-mapping methods do not provide insight into the complex local landscape of linkage that contains and has been shaped by the causal variant(s). Here, we present crosshap, an R package that performs robust density-based clustering of variants based on their linkage profiles to capture haplotype structures in a local genomic region of interest. Following this, crosshap is equipped with visualization tools for choosing optimal clustering parameters (ɛ) before producing an intuitive figure that provides an overview of the complex relationships between linked variants, haplotype combinations, phenotype, and metadata traits. AVAILABILITY AND IMPLEMENTATION: The crosshap package is freely available under the MIT license and can be downloaded directly from CRAN with R >4.0.0. The development version is available on GitHub alongside issue support (https://github.com/jacobimarsh/crosshap). Tutorial vignettes and documentation are available (https://jacobimarsh.github.io/crosshap/).

Genome composition in Brassica interspecific hybrids affects chromosome inheritance and viability of progeny.

Interspecific hybridization is widespread in nature and can result in the formation of new hybrid species as well as the transfer of traits between species. However, the fate of newly formed hybrid lineages is relatively understudied. We undertook pairwise crossing between multiple genotypes of three Brassica allotetraploid species Brassica juncea (2n = AABB), Brassica carinata (2n = BBCC), and Brassica napus (2n = AACC) to generate AABC, BBAC, and CCAB interspecific hybrids and investigated chromosome inheritance and fertility in these hybrids and their self-pollinated progeny. Surprisingly, despite the presence of a complete diploid genome in all hybrids, hybrid fertility was very low. AABC and BBAC first generation (F1) hybrids both averaged ~16% pollen viability compared to 3.5% in CCAB hybrids: most CCAB hybrid flowers were male-sterile. AABC and CCAB F1 hybrid plants averaged 5.5 and 0.5 seeds per plant, respectively, and BBAC F1 hybrids ~56 seeds/plant. In the second generation (S1), all confirmed self-pollinated progeny resulting from CCAB hybrids were sterile, producing no self-pollinated seeds. Three AABC S1 hybrids putatively resulting from unreduced gametes produced 3, 14, and 182 seeds each, while other AABC S1 hybrids averaged 1.5 seeds/plant (0-8). BBAC S1 hybrids averaged 44 seeds/plant (range 0-403). We also observed strong bias towards retention rather than loss of the haploid genomes, suggesting that the subgenomes in the Brassica allotetraploids are already highly interdependent, such that loss of one subgenome is detrimental to fertility and viability. Our results suggest that relationships between subgenomes determine hybridization outcomes in these species.

Characterisation of common hypothetical surface peptides between protozoan parasites (Perkinsus olseni) originating from different geographical locations.

Perkinsus olseni and P. marinus are classified as notifiable pathogens by the World Organisation for Animal Health and are known to cause perkinsosis in a variety of molluscs globally. Mass mortalities due to these parasites in farms and in the wild have been a recurrent issue. Diagnosis for these protozoans is currently done using Ray's fluid thioglycollate medium method followed by optical microscopy or molecular assays. Both require a high level of skill and are time-consuming. An immunoassay method would make the diagnosis of perkinsosis quicker and cheaper. The present study used mass spectrometry-based proteomics to investigate common hypothetical surface peptides between different geographical isolates of P. olseni, which could be used to develop immunoassays in the future. Two peptides were identified: POLS_08089, which is a 42.7 kDa peptide corresponding to the 60S ribosomal subunit protein L4; and POLS_15916, which is a conserved hypothetical protein of 55.6 kDa. The identification of peptides may allow the development of immunoassays through a more targeted approach.

Allele segregation analysis of F1 hybrids between independent Brassica allohexaploid lineages.

In the Brassica genus, we find both diploid species (one genome) and allotetraploid species (two different genomes) but no naturally occurring hexaploid species (three different genomes, AABBCC). Although hexaploids can be produced via human intervention, these neo-polyploids have quite unstable genomes and usually suffer from severe genome reshuffling. Whether these genome rearrangements continue in later generations and whether genomic arrangements follow similar, reproducible patterns between different lineages is still unknown. We crossed Brassica hexaploids resulting from different species combinations to produce five F1 hybrids and analyzed the karyotypes of the parents and the F1 hybrids, as well as allele segregation in a resulting test-cross population via molecular karyotyping using SNP array genotyping. Although some genomic regions were found to be more likely to be duplicated, deleted, or rearranged, a consensus pattern was not shared between genotypes. Brassica hexaploids had a high tolerance for fixed structural rearrangements, but which rearrangements occur and become fixed over many generations does not seem to show either strong reproducibility or to indicate selection for stability. On average, we observed 10 de novo chromosome rearrangements contributed almost equally from both parents to the F1 hybrids. At the same time, the F1 hybrid meiosis produced on average 8.6 new rearrangements. Hence, the increased heterozygosity in the F1 hybrid did not significantly improve genome stability in our hexaploid hybrids and might have had the opposite effect. However, hybridization between lineages was readily achieved and may be exploited for future genetics and breeding purposes.

Pangenomics Comes of Age: From Bacteria to Plant and Animal Applications.

The pangenome refers to a collection of genomic sequence found in the entire species or population rather than in a single individual; the sequence can be core, present in all individuals, or accessory (variable or dispensable), found in a subset of individuals only. While pangenomic studies were first undertaken in bacterial species, developments in genome sequencing and assembly approaches have allowed construction of pangenomes for eukaryotic organisms, fungi, plants, and animals, including two large-scale human pangenome projects. Analysis of the these pangenomes revealed key differences, most likely stemming from divergent evolutionary histories, but also surprising similarities.

Comparative pangenome analyses provide insights into the evolution of Brassica rapa resistance gene analogues (RGAs).

Brassica rapa is grown worldwide as economically important vegetable and oilseed crop. However, its production is challenged by yield-limiting pathogens. The sustainable control of these pathogens mainly relies on the deployment of genetic resistance primarily driven by resistance gene analogues (RGAs). While several studies have identified RGAs in B. rapa, these were mainly based on a single genome reference and do not represent the full range of RGA diversity in B. rapa. In this study, we utilized the B. rapa pangenome, constructed from 71 lines encompassing 12 morphotypes, to describe a comprehensive repertoire of RGAs in B. rapa. We show that 309 RGAs were affected by presence-absence variation (PAV) and 223 RGAs were missing from the reference genome. The transmembrane leucine-rich repeat (TM-LRR) RGA class had more core gene types than variable genes, while the opposite was observed for nucleotide-binding site leucine-rich repeats (NLRs). Comparative analysis with the B. napus pangenome revealed significant RGA conservation (93%) between the two species. We identified 138 candidate RGAs located within known B. rapa disease resistance QTL, of which the majority were under negative selection. Using blackleg gene homologues, we demonstrated how these genes in B. napus were derived from B. rapa. This further clarifies the genetic relationship of these loci, which may be useful in narrowing-down candidate blackleg resistance genes. This study provides a novel genomic resource towards the identification of candidate genes for breeding disease resistance in B. rapa and its relatives.

Key Advances in the New Era of Genomics-Assisted Disease Resistance Improvement of Brassica Species.

Disease resistance improvement remains a major focus in breeding programs as diseases continue to devastate Brassica production systems due to intensive cultivation and climate change. Genomics has paved the way to understand the complex genomes of Brassicas, which has been pivotal in the dissection of the genetic underpinnings of agronomic traits driving the development of superior cultivars. The new era of genomics-assisted disease resistance breeding has been marked by the development of high-quality genome references, accelerating the identification of disease resistance genes controlling both qualitative (major) gene and quantitative resistance. This facilitates the development of molecular markers for marker assisted selection and enables genome editing approaches for targeted gene manipulation to enhance the genetic value of disease resistance traits. This review summarizes the key advances in the development of genomic resources for Brassica species, focusing on improved genome references, based on long-read sequencing technologies and pangenome assemblies. This is further supported by the advances in pathogen genomics, which have resulted in the discovery of pathogenicity factors, complementing the mining of disease resistance genes in the host. Recognizing the co-evolutionary arms race between the host and pathogen, it is critical to identify novel resistance genes using crop wild relatives and synthetic cultivars or through genetic manipulation via genome-editing to sustain the development of superior cultivars. Integrating these key advances with new breeding techniques and improved phenotyping using advanced data analysis platforms will make disease resistance improvement in Brassica species more efficient and responsive to current and future demands.

CRISPR-Based Genome Editing Tools: An Accelerator in Crop Breeding for a Changing Future.

Genome editing is an important strategy to maintain global food security and achieve sustainable agricultural development. Among all genome editing tools, CRISPR-Cas is currently the most prevalent and offers the most promise. In this review, we summarize the development of CRISPR-Cas systems, outline their classification and distinctive features, delineate their natural mechanisms in plant genome editing and exemplify the applications in plant research. Both classical and recently discovered CRISPR-Cas systems are included, detailing the class, type, structures and functions of each. We conclude by highlighting the challenges that come with CRISPR-Cas and offer suggestions on how to tackle them. We believe the gene editing toolbox will be greatly enriched, providing new avenues for a more efficient and precise breeding of climate-resilient crops.

Machine learning models outperform deep learning models, provide interpretation and facilitate feature selection for soybean trait prediction.

Recent growth in crop genomic and trait data have opened opportunities for the application of novel approaches to accelerate crop improvement. Machine learning and deep learning are at the forefront of prediction-based data analysis. However, few approaches for genotype to phenotype prediction compare machine learning with deep learning and further interpret the models that support the predictions. This study uses genome wide molecular markers and traits across 1110 soybean individuals to develop accurate prediction models. For 13/14 sets of predictions, XGBoost or random forest outperformed deep learning models in prediction performance. Top ranked SNPs by F-score were identified from XGBoost, and with further investigation found overlap with significantly associated loci identified from GWAS and previous literature. Feature importance rankings were used to reduce marker input by up to 90%, and subsequent models maintained or improved their prediction performance. These findings support interpretable machine learning as an approach for genomic based prediction of traits in soybean and other crops.

Extensive polyploid clonality was a successful strategy for seagrass to expand into a newly submerged environment.

Polyploidy has the potential to allow organisms to outcompete their diploid progenitor(s) and occupy new environments. Shark Bay, Western Australia, is a World Heritage Area dominated by temperate seagrass meadows including Poseidon's ribbon weed, Posidonia australis. This seagrass is at the northern extent of its natural geographic range and experiences extremes in temperature and salinity. Our genomic and cytogenetic assessments of 10 meadows identified geographically restricted, diploid clones (2n = 20) in a single location, and a single widespread, high-heterozygosity, polyploid clone (2n = 40) in all other locations. The polyploid clone spanned at least 180 km, making it the largest known example of a clone in any environment on earth. Whole-genome duplication through polyploidy, combined with clonality, may have provided the mechanism for P. australis to expand into new habitats and adapt to new environments that became increasingly stressful for its diploid progenitor(s). The new polyploid clone probably formed in shallow waters after the inundation of Shark Bay less than 8500 years ago and subsequently expanded via vegetative growth into newly submerged habitats.

Haplotype mapping uncovers unexplored variation in wild and domesticated soybean at the major protein locus cqProt-003.

KEY MESSAGE: The major soy protein QTL, cqProt-003, was analysed for haplotype diversity and global distribution, and results indicate 304 bp deletion and variable tandem repeats in protein coding regions are likely causal candidates. Here, we present association and linkage analysis of 985 wild, landrace and cultivar soybean accessions in a pan genomic dataset to characterize the major high-protein/low-oil associated locus cqProt-003 located on chromosome 20. A significant trait-associated region within a 173 kb linkage block was identified, and variants in the region were characterized, identifying 34 high confidence SNPs, 4 insertions, 1 deletion and a larger 304 bp structural variant in the high-protein haplotype. Trinucleotide tandem repeats of variable length present in the second exon of gene Glyma.20G085100 are strongly correlated with the high-protein phenotype and likely represent causal variation. Structural variation has previously been found in the same gene, for which we report the global distribution of the 304 bp deletion and have identified additional nested variation present in high-protein individuals. Mapping variation at the cqProt-003 locus across demographic groups suggests that the high-protein haplotype is common in wild accessions (94.7%), rare in landraces (10.6%) and near absent in cultivated breeding pools (4.1%), suggesting its decrease in frequency primarily correlates with domestication and continued during subsequent improvement. However, the variation that has persisted in under-utilized wild and landrace populations holds high breeding potential for breeders willing to forego seed oil to maximize protein content. The results of this study include the identification of distinct haplotype structures within the high-protein population, and a broad characterization of the genomic context and linkage patterns of cqProt-003 across global populations, supporting future functional characterization and modification.

In silico integration of disease resistance QTL, genes and markers with the Brassica juncea physical map.

Brassica juncea (AABB), Indian mustard, is a source of disease resistance genes for a wide range of pathogens. The availability of reference genome sequences for B. juncea has made it possible to characterise the genomic structure and distribution of these disease resistance genes. Potentially functional disease resistance genes can be identified by co-localization with genetically mapped disease resistance quantitative trait loci (QTL). Here we identify and characterise disease resistance gene analogs (RGAs), including nucleotide-binding site-leucine-rich repeat (NLR), receptor-like kinase (RLK) and receptor-like protein (RLP) classes, and investigate their association with disease resistance QTL intervals. The molecular genetic marker sequences for four white rust (Albugo candida) disease resistance QTL, six blackleg (Leptosphaeria maculans) disease resistance QTL and BjCHI1, a gene cloned from B. juncea for hypocotyl rot disease, were extracted from previously published studies and used to compare with candidate RGAs. Our results highlight the complications for the identification of functional resistance genes, including the duplicated appearance of genetic markers for several resistance loci, including Ac2(t), AcB1-A4.1, AcB1-A5.1, Rlm6 and PhR2 in both the A and B genomes, due to the presence of homoeologous regions. Furthermore, the white rust loci, Ac2(t) and AcB1-A4.1, mapped to the same position on chromosome A04 and may be different alleles of the same gene. Despite these challenges, a total of nine candidate genomic regions hosting 14 RLPs, 28 NLRs and 115 RLKs were identified. This study facilitates the mapping and cloning of functional resistance genes for applications in crop improvement programs. Supplementary Information: The online version contains supplementary material available at 10.1007/s11032-022-01309-5.

The miR393-Target Module Regulates Plant Development and Responses to Biotic and Abiotic Stresses.

MicroRNAs (miRNAs), a class of endogenous small RNAs, are broadly involved in plant development, morphogenesis and responses to various environmental stresses, through manipulating the cleavage, translational expression, or DNA methylation of target mRNAs. miR393 is a conserved miRNA family present in many plants, which mainly targets genes encoding the transport inhibitor response1 (TIR1)/auxin signaling F-box (AFB) auxin receptors, and thus greatly affects the auxin signal perception, Aux/IAA degradation, and related gene expression. This review introduces the advances made on the miR393/target module regulating plant development and the plant's responses to biotic and abiotic stresses. This module is valuable for genetic manipulation of optimized conditions for crop growth and development and would also be helpful in improving crop yield through molecular breeding.

Expanding Gene-Editing Potential in Crop Improvement with Pangenomes.

Pangenomes aim to represent the complete repertoire of the genome diversity present within a species or cohort of species, capturing the genomic structural variance between individuals. This genomic information coupled with phenotypic data can be applied to identify genes and alleles involved with abiotic stress tolerance, disease resistance, and other desirable traits. The characterisation of novel structural variants from pangenomes can support genome editing approaches such as Clustered Regularly Interspaced Short Palindromic Repeats and CRISPR associated protein Cas (CRISPR-Cas), providing functional information on gene sequences and new target sites in variant-specific genes with increased efficiency. This review discusses the application of pangenomes in genome editing and crop improvement, focusing on the potential of pangenomes to accurately identify target genes for CRISPR-Cas editing of plant genomes while avoiding adverse off-target effects. We consider the limitations of applying CRISPR-Cas editing with pangenome references and potential solutions to overcome these limitations.

Pangenomes as a Resource to Accelerate Breeding of Under-Utilised Crop Species.

Pangenomes are a rich resource to examine the genomic variation observed within a species or genera, supporting population genetics studies, with applications for the improvement of crop traits. Major crop species such as maize (Zea mays), rice (Oryza sativa), Brassica (Brassica spp.), and soybean (Glycine max) have had pangenomes constructed and released, and this has led to the discovery of valuable genes associated with disease resistance and yield components. However, pangenome data are not available for many less prominent crop species that are currently under-utilised. Despite many under-utilised species being important food sources in regional populations, the scarcity of genomic data for these species hinders their improvement. Here, we assess several under-utilised crops and review the pangenome approaches that could be used to build resources for their improvement. Many of these under-utilised crops are cultivated in arid or semi-arid environments, suggesting that novel genes related to drought tolerance may be identified and used for introgression into related major crop species. In addition, we discuss how previously collected data could be used to enrich pangenome functional analysis in genome-wide association studies (GWAS) based on studies in major crops. Considering the technological advances in genome sequencing, pangenome references for under-utilised species are becoming more obtainable, offering the opportunity to identify novel genes related to agro-morphological traits in these species.

Isolation and Sequencing of Chromosome Arm 7RS of Rye, Secale cereale.

Rye (Secale cereale) is a climate-resilient cereal grown extensively as grain or forage crop in Northern and Eastern Europe. In addition to being an important crop, it has been used to improve wheat through introgression of genomic regions for improved yield and disease resistance. Understanding the genomic diversity of rye will assist both the improvement of this crop and facilitate the introgression of more valuable traits into wheat. Here, we isolated and sequenced the short arm of rye chromosome 7 (7RS) from Triticale 380SD using flow cytometry and compared it to the public Lo7 rye whole genome reference assembly. We identify 2747 Lo7 genes present on the isolated chromosome arm and two clusters containing seven and sixty-five genes that are present on Triticale 380SD 7RS, but absent from Lo7 7RS. We identified 29 genes that are not assigned to chromosomal locations in the Lo7 assembly but are present on Triticale 380SD 7RS, suggesting a chromosome arm location for these genes. Our study supports the Lo7 reference assembly and provides a repertoire of genes on Triticale 7RS.

Current status of structural variation studies in plants.

Structural variations (SVs) including gene presence/absence variations and copy number variations are a common feature of genomes in plants and, together with single nucleotide polymorphisms and epigenetic differences, are responsible for the heritable phenotypic diversity observed within and between species. Understanding the contribution of SVs to plant phenotypic variation is important for plant breeders to assist in producing improved varieties. The low resolution of early genetic technologies and inefficient methods have previously limited our understanding of SVs in plants. However, with the rapid expansion in genomic technologies, it is possible to assess SVs with an ever-greater resolution and accuracy. Here, we review the current status of SV studies in plants, examine the roles that SVs play in phenotypic traits, compare current technologies and assess future challenges for SV studies.

Modelling of gene loss propensity in the pangenomes of three Brassica species suggests different mechanisms between polyploids and diploids.

Plant genomes demonstrate significant presence/absence variation (PAV) within a species; however, the factors that lead to this variation have not been studied systematically in Brassica across diploids and polyploids. Here, we developed pangenomes of polyploid Brassica napus and its two diploid progenitor genomes B. rapa and B. oleracea to infer how PAV may differ between diploids and polyploids. Modelling of gene loss suggests that loss propensity is primarily associated with transposable elements in the diploids while in B. napus, gene loss propensity is associated with homoeologous recombination. We use these results to gain insights into the different causes of gene loss, both in diploids and following polyploidization, and pave the way for the application of machine learning methods to understanding the underlying biological and physical causes of gene presence/absence.