Expanded functionality, increased accuracy, and enhanced speed in the de novo genotyping-by-sequencing pipeline GBS-SNP-CROP

A first look at sea-lavenders genomics - can genome wide SNP information tip the scales of controversy in the Limonium vulgare species complex?

BACKGROUND

Sea-lavenders (Limonium Mill., Plumbaginaceae) are a cosmopolitan group of diploid and polyploid plants often adapted to extreme saline environments, with a mostly Tethyan distribution, occurring in the Mediterranean, Irano-Turanian, Euro-Siberian and in the New World. The halophylic Limonium vulgare polyploid complex in particular, presents a large distribution throughout extreme salt-marsh habitats and shows little morphological but high taximetric variation, frequently blurring species delimitation. In this work we pursue three main goals: assert whether SNP data from polyploid individuals has the resolution to distinguish the seven sampled species, to better understand how genetically structured Limonium vulgare is, and attempt to identify specific molecular mechanisms for the differentiation between L. maritimum and L. vulgare. For this purpose, 95 individuals were genotyped using Genotyping by Sequencing (GBS), which were assembled as two independent datasets using IPYRAD. All analyses performed downstream of assembly were fully automated. Phylogenetic inference, PCA, and admixture plots were used to infer answers to the study's main goals.

RESULTS

Close to 10,000 SNPs were obtained for each dataset. Phylogenetic analyses reveal that polyploid data can be used to infer species relationships. Population structure analyses suggest a genetically structured L. vulgare. A set of 34 SNPs were found to be fully segregated between L. vulgare and L. maritimum, two of which are potentially linked to proteins that might be involved in the speciation process.

CONCLUSION

Despite polyploid data analyses shortcomings, GBS generated SNPs have the resolution to discern all seven included species. Limonium vulgare revealed pronounced genetic structure along a geographical north-south cline. L. maritimum always appears as a distinct genetic entity. Segregated SNPs between L. vulgare and L. maritimum indicate salinity response and morphological trait control genes as potentially interesting to follow up for studying these species' divergence process.

Evidence of Asexual Overwintering of Melampsora paradoxa and Mapping of Stem Rust Host Resistance in Salix

Melampsora rust is a devastating disease of shrub willow in North America. Previous work has identified Melampsora paradoxa as one of two identified rust species in New York State that infect Salix purpurea and other important Salix host species, however little is known about the population of this rust species in this region. Genotyping-by-sequencing was used to identify single nucleotide polymorphisms (SNPs) and assess population diversity of M. paradoxa isolates collected from three Salix breeding populations in Geneva, NY between 2015 and 2020. Statistical analyses of SNP revealed that all isolates collected were clonally derived even though they were collected across years. In 2020, isolates were collected from stem infections where uredospore pustules were observed, and these isolates were also identical to M. paradoxa collected in previous seasons. These data suggest that M. paradoxa sampled across multiple years overwintered and reproduced asexually and that stem infection is a possible mechanism for overwintering, both of which are novel findings for this rust species. Additionally, field disease ratings were conducted on a S. purpurea × S. suchowensis F₁ breeding population with high disease severity, enabling the discovery of QTL for resistance on chromosomes 1 and 19. Lastly, Colletotrichum salicis was frequently associated with stem rust and may play a role in M. paradoxa stem infection. Together, this work is the first substantial exploration into M. paradoxa population biology, stem infection, and host resistance in Salix.

On genetic diversity in caraway: Genotyping of a large germplasm collection

Caraway (Carum carvi) is a widespread and frequently used spice and medicinal plant with a long history of cultivation. However, due to ongoing climatic changes, the cultivation is becoming increasingly risky. To secure caraway cultivation in future, timely breeding efforts to develop adapted material are necessary. Analysis of genetic diversity can accompany this process, for instance, by revealing untapped gene pools. Here, we analyzed 137 accessions using genotyping by sequencing (GBS). Hence, we can report a broad overview of population structure and genetic diversity of caraway. Population structure was determined using a principal coordinate analysis, a Bayesian clustering analysis, phylogenetic trees and a neighbor network based on 13,155 SNPs. Genotypic data indicate a clear separation of accessions into two subpopulations, which correlates with the flowering type (annual vs. biennial). Four winter-annual accessions were closer related to biennial accessions. In an analysis of molecular variance, genetic variation between the two subpopulations was 7.84%. In addition, we estimated the genome size for 35 accessions by flow cytometry. An average genome size of 4.282 pg/2C (± 0.0096 S.E.) was estimated. Therefore, we suggest a significantly smaller genome size than stated in literature.

Simulation with RADinitio improves RADseq experimental design and sheds light on sources of missing data

Restriction-site associated DNA sequencing (RADseq) has become a powerful and versatile tool in modern population genomics, enabling large-scale evolutionary and genomic analyses in otherwise inaccessible biological systems. With its widespread use, different variants on the protocol have been developed to suit specific experimental needs. Researchers face the challenge of choosing the optimal molecular and sequencing protocols for their reduced representation experimental design, an often-complicated process. Strategic errors can lead to biased data generation that has reduced power to answer biological questions. Here, we present RADinitio, simulation software for the selection and optimization of RADseq experiments via the generation of sequencing data that behave similarly to empirical sources. RADinitio provides an evolutionary simulation of populations, implementation of various RADseq protocols with customizable parameters, and thorough assessment of missing data. We test the efficacy of the software using different RAD protocols across several organisms, highlighting the importance of protocol selection on the magnitude and quality of data acquired. Additionally, we test the effects of RAD library preparation and sequencing on allelic dropout, observing that library preparation and sequencing often contributes more to missing alleles than population-level variation.

Development of species-specific molecular markers in Vanilla for seedling selection of hybrids

Vanilla planifolia is the primary botanical source of vanilla extract used globally in various foods and beverages. V. planifolia has a global distribution based on a few foundational clones and therefore has limited genetic diversity. Many Vanilla species easily hybridize with V. planifolia and could be a source of valuable genetic traits like increased vanillin content, disease resistance, or early flowering. While breeding Vanilla hybrids may improve plant performance, basic molecular tools for this species are lacking. DNA-based molecular markers are the most efficient method to validate hybrid progeny, detect hybrids in commercial plantings, and identify unknown accessions. This study used publicly available sequence data to develop species-specific, qRT-PCR-based molecular markers for Vanilla. Over 580,000 assembled sequence fragments were filtered for species specificity and twenty-two targets were selected for qRT-PCR screening. Ten targets differentially amplified among V. planifolia, V. pompona, V. phaeantha, and V. palmarum with ΔCT values as high as 17.58 between species. The ten targets were used to validate the parentage of hybrid progeny from controlled crosses with most hybrid progeny showing amplification patterns similar to both parents. The ten targets were also used to screen sixteen Vanilla species for specificity, and supported species assignments for unknown accessions including the detection of putative hybrids. This is the first report using species-specific, qRT-PCR-based molecular markers in Vanilla. These markers are inexpensive, simple to develop, and can rapidly screen large populations. These methods will enable the further development of species-specific molecular markers when creating Vanilla interspecific hybrid populations.

Mapping non-host resistance to the stem rust pathogen in an interspecific barberry hybrid

BACKGROUND

Non-host resistance (NHR) presents a compelling long-term plant protection strategy for global food security, yet the genetic basis of NHR remains poorly understood. For many diseases, including stem rust of wheat [causal organism Puccinia graminis (Pg)], NHR is largely unexplored due to the inherent challenge of developing a genetically tractable system within which the resistance segregates. The present study turns to the pathogen's alternate host, barberry (Berberis spp.), to overcome this challenge.

RESULTS

In this study, an interspecific mapping population derived from a cross between Pg-resistant Berberis thunbergii (Bt) and Pg-susceptible B. vulgaris was developed to investigate the Pg-NHR exhibited by Bt. To facilitate QTL analysis and subsequent trait dissection, the first genetic linkage maps for the two parental species were constructed and a chromosome-scale reference genome for Bt was assembled (PacBio + Hi-C). QTL analysis resulted in the identification of a single 13 cM region (~ 5.1 Mbp spanning 13 physical contigs) on the short arm of Bt chromosome 3. Differential gene expression analysis, combined with sequence variation analysis between the two parental species, led to the prioritization of several candidate genes within the QTL region, some of which belong to gene families previously implicated in disease resistance.

CONCLUSIONS

Foundational genetic and genomic resources developed for Berberis spp. enabled the identification and annotation of a QTL associated with Pg-NHR. Although subsequent validation and fine mapping studies are needed, this study demonstrates the feasibility of and lays the groundwork for dissecting Pg-NHR in the alternate host of one of agriculture's most devastating pathogens.