Genetic Diversity and Population Structure of Traditional Chinese Herb “Chai-Hu” Resources Using Genome-Wide SNPs Through Genotyping-By-Sequencing

Bupleurum (named “Chai-hu”) is an important traditional Chinese medicine resource in China. It has been widely used since ancient times and has antipyretic, analgesic and cholagogic functions, but there is little research on its genetic diversity. In this study, genotyping-by-sequencing (GBS) was used to detect SNP loci in 39 Bupleurum germplasm resources from different regions in China and analyse their genetic diversity. A total of 25.1 Gb of data was obtained by sequencing, with an average of 0.64 Gb per sample. After screening, 83898 high-quality SNPs were obtained. The results of genetic research were obtained by phylogenetic tree, principal component analysis and population structure analysis, and the 39 experimental materials were divided into three groups. The average observed heterozygosity and expected heterozygosity of Bupleurum populations were 0.24 and 0.17, respectively, indicating that Bupleurum populations from five different provinces had a low level of genetic diversity. Population nucleotide diversity analysis and analysis of molecular variance showed that the percentage of intrapopulation variation was 120.88%, while the percentage of interpopulation variation was only 2.46%. There was relative aggregation of Bupleurum samples with the same geographical origin, but the division of population structure was not completely correlated with sample origin. The results showed that the genetic diversity of the materials was low and that the genetic variation was narrow. This provides a good basis for the genetic breeding and protection of species diversity of Bupleurum.

Characterization and Mapping of retr04, retr05 and retr06 Broad-Spectrum Resistances to Turnip Mosaic Virus in Brassica juncea, and the Development of Robust Methods for Utilizing Recalcitrant Genotyping Data

Turnip mosaic virus (TuMV) induces disease in susceptible hosts, notably impacting cultivation of important crop species of the Brassica genus. Few effective plant viral disease management strategies exist with the majority of current approaches aiming to mitigate the virus indirectly through control of aphid vector species. Multiple sources of genetic resistance to TuMV have been identified previously, although the majority are strain-specific and have not been exploited commercially. Here, two Brassica juncea lines (TWBJ14 and TWBJ20) with resistance against important TuMV isolates (UK 1, vVIR24, CDN 1, and GBR 6) representing the most prevalent pathotypes of TuMV (1, 3, 4, and 4, respectively) and known to overcome other sources of resistance, have been identified and characterized. Genetic inheritance of both resistances was determined to be based on a recessive two-gene model. Using both single nucleotide polymorphism (SNP) array and genotyping by sequencing (GBS) methods, quantitative trait loci (QTL) analyses were performed using first backcross (BC1) genetic mapping populations segregating for TuMV resistance. Pairs of statistically significant TuMV resistance-associated QTLs with additive interactive effects were identified on chromosomes A03 and A06 for both TWBJ14 and TWBJ20 material. Complementation testing between these B. juncea lines indicated that one resistance-linked locus was shared. Following established resistance gene nomenclature for recessive TuMV resistance genes, these new resistance-associated loci have been termed retr04 (chromosome A06, TWBJ14, and TWBJ20), retr05 (A03, TWBJ14), and retr06 (A03, TWBJ20). Genotyping by sequencing data investigated in parallel to robust SNP array data was highly suboptimal, with informative data not established for key BC1 parental samples. This necessitated careful consideration and the development of new methods for processing compromised data. Using reductive screening of potential markers according to allelic variation and the recombination observed across BC1 samples genotyped, compromised GBS data was rendered functional with near-equivalent QTL outputs to the SNP array data. The reductive screening strategy employed here offers an alternative to methods relying upon imputation or artificial correction of genotypic data and may prove effective for similar biparental QTL mapping studies.

QTL-seq for the identification of candidate genes for days to flowering and leaf shape in pigeonpea

To identify genomic segments associated with days to flowering (DF) and leaf shape in pigeonpea, QTL-seq approach has been used in the present study. Genome-wide SNP profiling of extreme phenotypic bulks was conducted for both the traits from the segregating population (F2) derived from the cross combination- ICP 5529 × ICP 11605. A total of 126.63 million paired-end (PE) whole-genome resequencing data were generated for five samples, including one parent ICP 5529 (obcordate leaf and late-flowering plant), early and late flowering pools (EF and LF) and obcordate and lanceolate leaf shape pools (OLF and LLS). The QTL-seq identified two significant genomic regions, one on CcLG03 (1.58 Mb region spanned from 19.22 to 20.80 Mb interval) for days to flowering (LF and EF pools) and another on CcLG08 (2.19 Mb region spanned from 6.69 to 8.88 Mb interval) for OLF and LLF pools, respectively. Analysis of genomic regions associated SNPs with days to flowering and leaf shape revealed 5 genic SNPs present in the unique regions. The identified genomic regions for days to flowering were also validated with the genotyping-by-sequencing based classical QTL mapping method. A comparative analysis of the identified seven genes associated with days to flowering on 12 Fabaceae genomes, showed synteny with 9 genomes. A total of 153 genes were identified through the synteny analysis ranging from 13 to 36. This study demonstrates the usefulness of QTL-seq approach in precise identification of candidate gene(s) for days to flowering and leaf shape which can be deployed for pigeonpea improvement.

Molecular phylogeny and taxonomy of the Hydrangea serrata complex (Hydrangeaceae) in western Japan, including a new subspecies of H. acuminata from Yakushima

According to the contemporary classification of Hydrangea native to Japan, H. serrata is a polymorphic species including six varieties. We discovered a plant identified as H. serrata, but morphologically distinct from previously known varieties, in Yakushima island where approximately 50 endemic species are known. To determine the relationship of this plant with previously known varieties, we examined morphology and constructed a highly resolved phylogeny of H. serrata and its relatives using three chloroplast genomic regions, rbcL, trnL intron, psbA-trnH, and two nuclear genomic regions, ITS1 and ITS2, and Multiplex ISSR genotyping by sequencing (MIG-seq). Based on these morphological and phylogenetic observations, we describe Hydrangea acuminata subsp. yakushimensissubsp. nov. as a newly discovered lineage in Yakushima, Japan and propose Hydrangea minamitaniistat. nov. and Hydrangea acuminata subsp. australisstat. nov. which were previously treated as varieties of H. serrata.

The diversity of Passalora fulva Isolates Collected from Tomato Plants in US high tunnels

High tunnels extend the growing season of high value crops, including tomatoes, but the environmental conditions within high tunnels favor the spread of the tomato leaf mold pathogen, Passalora fulva (syn. Cladosporium fulvum). Tomato leaf mold results in defoliation, and if severe, losses in yield. Despite substantial research, little is known regarding the genetic structure and diversity of populations of P. fulva associated with high tunnel tomato production in the United States. From 2016 to 2019, a total of 50 P. fulva isolates were collected from tomato leaf samples in high tunnels in the Northeast and Minnesota. Other Cladosporium species were also isolated from the leaf surfaces. Koch’s postulates were conducted to confirm that P. fulva was the cause of the disease symptoms observed. Race determination experiments revealed that the isolates belonged to either race 0 (six isolates) or race 2 (44 isolates). Polymorphisms were identified within four previously characterized effector genes Avr2, Avr4, Avr4e, and Avr9. The largest number of polymorphisms were observed for Avr2. Both mating type genes, MAT1-1-1 and MAT1-2-1, were present in the isolate collection. For further insights into the pathogen diversity, the 50 isolates were genotyped at 7,514 single-nucleotide polymorphism loci using genotyping-by-sequencing: differentiation by region but not by year was observed. Within the collection of 50 isolates, there were 18 distinct genotypes. Information regarding P. fulva population diversity will enable better management recommendations for growers, as high tunnel production of tomatoes expands.

Multi-GBS: A massively multiplexed GBS-based protocol optimized for large, repetitive conifer genomes

Reduced representation library approaches are still a valuable tool for breeding and population and ecological genomics, even with impressive increases in sequencing capacity in recent years. Unfortunately, current approaches only allow for multiplexing up to 384 samples. To take advantage of increased sequencing capacity, we present Multi-GBS, a massively multiplexable extension to Genotyping-by-Sequencing that is also optimized for large conifer genomes. In Norway Spruce, a highly repetitive 20Gbp diploid genome with high population genetic variation, we call over a million variants in 32 genotypes from three populations, two natural forest in the Alps and Bohemian Alps, and a managed population from southeastern Austria using the existing TASSEL GBSv2 pipeline. Metric MDS analysis of replicated genotypes shows that technical bias in resulting genotype calling is minimal and that populations cluster in biologically meaningful ways.

Genomic regions associated with herbicide tolerance in a worldwide faba bean (Vicia faba L.) collection

Weeds represent one of the major constraints for faba bean crop. The identification of molecular markers associated with key genes imparting tolerance to herbicides can facilitate and fasten the efficient and effective development of herbicide tolerant cultivars. We phenotyped 140 faba bean genotypes in three open field experiments at two locations in Lebanon and Morocco against three herbicide treatments (T1 metribuzin 250 g ai/ha; T2 imazethapyr 75 g ai/ha; T3 untreated) and one in greenhouse where T1 and T3 were applied. The same set was genotyped using genotyping by sequencing (GBS) which yield 10,794 high quality single nucleotide polymorphisms (SNPs). ADMIXTURE software was used to infer the population structure which revealed two ancestral subpopulations. To identify SNPs associated with phenological and yield related traits under herbicide treatments, Single-trait (ST) and Multi-trait (MT) Genome Wide Association Studies (GWAS) were fitted using GEMMA software, showing 10 and 14 highly significant associations, respectively. Genomic sequences containing herbicide tolerance associated SNPs were aligned against the NCBI database using BLASTX tool using default parameters to annotate candidate genes underlying the causal variants. SNPs from acidic endochitinase, LRR receptor-like serine/threonine-protein kinase RCH1, probable serine/threonine-protein kinase NAK, malate dehydrogenase, photosystem I core protein PsaA and MYB-related protein P-like were significantly associated with herbicide tolerance traits.

Insights into opium poppy (Papaver spp.) genetic diversity from genotyping-by-sequencing analysis

Opium poppy (Papaver somniferum) is one of the world’s oldest medicinal plants and a versatile model system to study secondary metabolism. However, our knowledge of its genetic diversity is limited, restricting utilization of the available germplasm for research and crop improvement. We used genotyping-by-sequencing to investigate the extent of genetic diversity and population structure in a collection of poppy germplasm consisting of 91 accessions originating in 30 countries of Europe, North Africa, America, and Asia. We identified five genetically distinct subpopulations using discriminate analysis of principal components and STRUCTURE analysis. Most accessions obtained from the same country were grouped together within subpopulations, likely a consequence of the restriction on movement of poppy germplasm. Alkaloid profiles of accessions were highly diverse, with morphine being dominant. Phylogenetic analysis identified genetic groups that were largely consistent with the subpopulations detected and that could be differentiated broadly based on traits such as number of branches and seed weight. These accessions and the associated genotypic data are valuable resources for further genetic diversity analysis, which could include definition of poppy core sets to facilitate genebank management and use of the diversity for genetic improvement of this valuable crop.

A High-Density Genetic Map Enables Genome Synteny and QTL Mapping of Vegetative Growth and Leaf Traits in Gardenia

The gardenia is a traditional medicinal horticultural plant in China, but its molecular genetic research has been largely hysteretic. Here, we constructed an F1 population with 200 true hybrid individuals. Using the genotyping-by-sequencing method, a high-density sex-average genetic map was generated that contained 4,249 SNPs with a total length of 1956.28 cM and an average genetic distance of 0.46 cM. We developed 17 SNP-based Kompetitive Allele-Specific PCR markers and found that 15 SNPs were successfully genotyped, of which 13 single-nucleotide polymorphism genotypings of 96 F1 individuals showed genotypes consistent with GBS-mined genotypes. A genomic collinearity analysis between gardenia and the Rubiaceae species Coffea arabica, Coffea canephora and Ophiorrhiza pumila showed the relativity strong conservation of LG11 with NC_039,919.1, HG974438.1 and Bliw01000011.1, respectively. Lastly, a quantitative trait loci analysis at three phenotyping time points (2019, 2020, and 2021) yielded 18 QTLs for growth-related traits and 31 QTLs for leaf-related traits, of which qBSBN7-1, qCD8 and qLNP2-1 could be repeatably detected. Five QTL regions (qCD8 and qSBD8, qBSBN7 and qSI7, qCD4-1 and qLLLS4, qLNP10 and qSLWS10-2, qSBD10 and qLLLS10) with potential pleiotropic effects were also observed. This study provides novel insight into molecular genetic research and could be helpful for further gene cloning and marker-assisted selection for early growth and development traits in the gardenia.

Mapping of dwarfing QTL of Ari1327, a semi-dwarf mutant of upland cotton

Background Upland Cotton (Gossypium hirsutum L.) has few cotton varieties suitable for mechanical harvesting. The plant height of the cultivar is one of the key features that need to modify. Hence, this study was planned to locate the QTL for plant height in a 60Co γ treated upland cotton semi-dwarf mutant Ari1327. Results Interestingly, bulk segregant analysis (BSA) and genotyping by sequencing (GBS) methods exhibited that candidate QTL was co-located in the region of 5.80–9.66 Mb at D01 chromosome in two F2 populations. Using three InDel markers to genotype a population of 1241 individuals confirmed that the offspring’s phenotype is consistent with the genotype. Comparative analysis of RNA-seq between the mutant and wild variety exhibited that Gh_D01G0592 was identified as the source of dwarfness from 200 genes. In addition, it was also revealed that the appropriate use of partial separation markers in QTL mapping can escalate linkage information. Conclusions Overwhelmingly, the results will provide the basis to reveal the function of candidate genes and the utilization of excellent dwarf genetic resources in the future.