Genetic diversity and population structure of trifoliate yam (Dioscorea dumetorum Kunth) in Cameroon revealed by genotyping-by-sequencing (GBS)

Unravelling the genetic diversity of water yam (Dioscorea alata L.) accessions from Tanzania using simple sequence repeat (SSR) markers

Water yam (Dioscorea alata L.) is among the most cultivated species used as a source of food and income for small-scale farmers in Tanzania. However, little is documented about Dioscorea species available in Tanzania, including their genetic diversity. This study used ten polymorphic microsatellite markers to determine the genetic diversity and relationship of 63 D. alata accessions from six major producing regions. Results revealed a polymorphic information content (PIC) of 0.63, while the number of alleles per locus ranged from 4 to 12 with a mean of 7.60. The expected heterozygosity ranged from 0.20to 0.76, with a mean of 0.53, which suggests moderate genetic diversity of D. alata accessions. Kagera region had the highest mean number of (1.5) private alleles. Analysis of molecular variance revealed that 54% of the variation was attributed to within individual, 39% was attributed to among individual while among population contributed 7% of the total variation. The highest Nei's genetic distance (0.43) was for accessions sampled from Kilimanjaro and Mtwara regions. Principal coordinate analysis and cluster analysis using Unweighted Paired Group Method using Arithmetic (UPGMA) grouped D. alata accessions into two major clusters regardless of geographical origin and local names. The Bayesian structure analysis confirmed the two clusters obtained in UPGMA and revealed an admixture of D. alata accessions in all six regions suggesting farmers' extensive exchange of planting materials. These results are helpful in the selection of D. alata accessions for breeding programs in Tanzania.

Infrared spectroscopy for ploidy estimation: An example in two species of Veronica using fresh and herbarium specimens

PREMISE

Polyploidy has become a central factor in plant evolutionary biological research in recent decades. Methods such as flow cytometry have revealed the widespread occurrence of polyploidy; however, its inference relies on expensive lab equipment and is largely restricted to fresh or recently dried material.

METHODS

Here, we assess the applicability of infrared spectroscopy to infer ploidy in two related species of Veronica (Plantaginaceae). Infrared spectroscopy relies on differences in the absorbance of tissues, which could be affected by primary and secondary metabolites related to polyploidy. We sampled 33 living plants from the greenhouse and 74 herbarium specimens with ploidy known through flow cytometrical measurements and analyzed the resulting spectra using discriminant analysis of principal components (DAPC) and neural network (NNET) classifiers.

RESULTS

Living material of both species combined was classified with 70% (DAPC) to 75% (NNET) accuracy, whereas herbarium material was classified with 84% (DAPC) to 85% (NNET) accuracy. Analyzing both species separately resulted in less clear results.

DISCUSSION

Infrared spectroscopy is quite reliable but is not a certain method for assessing intraspecific ploidy level differences in two species of Veronica. More accurate inferences rely on large training data sets and herbarium material. This study demonstrates an important way to expand the field of polyploid research to herbaria.

Sex Determination in Dioscorea dumetorum: Evidence of Heteromorphic Sex Chromosomes and Sex-Linked NORs

Yams (Dioscorea spp.) are a pantropical genus located worldwide that constitute an important source of nutrients and pharmaceutical substances. Some Dioscorea crop species are widely grown in West Africa. One species that is mainly cultivated in Cameroon is Dioscorea dumetorum. This is a dioecious root crop whose sex-determining system was unknown until now. To address the possible presence of sex chromosomes in D. dumetorum, we performed a karyotype characterization of male and female individuals using classical and molecular cytogenetic approaches. It was determined that 2n = 40 was the most common number of chromosomes in all of the investigated samples. One chromosome pair was longer than the others in the chromosome set and was a heteromorph in male and homomorph in female individuals. This pair corresponded to sex chromosomes, and we also confirmed this with molecular cytogenetic experiments. The results of chromomycin banding revealed the presence of strong positive signals on this chromosome pair. The signals, corresponding to GC-rich DNA regions, were similar in size on the chromosomes of the female individuals, whereas they were different in size in the male individuals. This size difference in the GC-rich heterochromatin regions was also apparent in the interphase nuclei as one small and one large fluorescent spot. The results of the in situ hybridization experiment showed that these chromomycin positive signals on the sex chromosomes also corresponded to the 35S rDNA cluster. The mean 2C DNA value (genome size) obtained for D. dumentorum was 0.71 pg (±0.012), which represents a small genome size. We found no difference in the genome size between the male and female individuals. The results of this study contribute to increasing our knowledge of sex determination in D. dumetorum (standard sex-determining XX/XY system) and may have some agronomic applications.

Genome Analysis Revives a Forgotten Hybrid Crop Edo-dokoro in the Genus Dioscorea

A rhizomatous Dioscorea crop 'Edo-dokoro' was described in old records of Japan, but its botanical identity has not been characterized. We found that Edo-dokoro is still produced by four farmers in Tohoku-machi of the Aomori prefecture, Japan. The rhizomes of Edo-dokoro are a delicacy to the local people and are sold in the markets. Morphological characters of Edo-dokoro suggest its hybrid origin between the two species, Dioscorea tokoro and Dioscorea tenuipes. Genome analysis revealed that Edo-dokoro likely originated by hybridization of a male D. tokoro to a female D. tenuipes, followed by a backcross with a male plant of D. tokoro. Edo-dokoro is a typical minor crop possibly maintained for more than 300 years but now almost forgotten by the public. We hypothesize that there are many such uncharacterized genetic heritages passed over generations by small-scale farmers that await serious scientific investigation for future use and improvement by using modern genomics information.

Automating microsatellite screening and primer design from multi-individual libraries using Micro-Primers

Analysis of intra- and inter-population diversity has become important for defining the genetic status and distribution patterns of a species and a powerful tool for conservation programs, as high levels of inbreeding could lead into whole population extinction in few generations. Microsatellites (SSR) are commonly used in population studies but discovering highly variable regions across species' genomes requires demanding computation and laboratorial optimization. In this work, we combine next generation sequencing (NGS) with automatic computing to develop a genomic-oriented tool for characterizing SSRs at the population level. Herein, we describe a new Python pipeline, named Micro-Primers, designed to identify, and design PCR primers for amplification of SSR loci from a multi-individual microsatellite library. By combining commonly used programs for data cleaning and microsatellite mining, this pipeline easily generates, from a fastq file produced by high-throughput sequencing, standard information about the selected microsatellite loci, including the number of alleles in the population subset, and the melting temperature and respective PCR product of each primer set. Additionally, potential polymorphic loci can be identified based on the allele ranges observed in the population, to easily guide the selection of optimal markers for the species. Experimental results show that Micro-Primers significantly reduces processing time in comparison to manual analysis while keeping the same quality of the results. The elapsed times at each step can be longer depending on the number of sequences to analyze and, if not assisted, the selection of polymorphic loci from multiple individuals can represent a major bottleneck in population studies.

Transcriptome Sequence Reveals Candidate Genes Involving in the Post-Harvest Hardening of Trifoliate Yam Dioscorea dumetorum

Storage ability of trifoliate yam (Dioscorea dumetorum) is restricted by a severe post-harvest hardening (PHH) phenomenon, which starts within the first 24 h after harvest and renders tubers inedible. Previous work has only focused on the biochemical changes affecting PHH in D. dumetorum. To the best of our knowledge, the candidate genes responsible for the hardening of D. dumetorum have not been identified. Here, transcriptome analyses of D. dumetorum tubers were performed in yam tubers of four developmental stages: 4 months after emergence (4MAE), immediately after harvest (AH), 3 days after harvest (3DAH) and 14 days after harvest (14DAH) of four accessions (Bangou 1, Bayangam 2, Fonkouankem 1, and Ibo sweet 3) using RNA-Seq. In total, between AH and 3DAH, 165, 199, 128 and 61 differentially expressed genes (DEGs) were detected in Bayangam 2, Fonkouankem 1, Bangou 1 and Ibo sweet 3, respectively. Functional analysis of DEGs revealed that genes encoding for CELLULOSE SYNTHASE A (CESA), XYLAN O-ACETYLTRANSFERASE (XOAT), CHLOROPHYLL A/B BINDING PROTEIN1, 2, 3, 4 (LHCB1, LHCB2, LHCB3, and LCH4) and an MYB transcription factor were predominantly and significantly up-regulated 3DAH, implying that these genes were potentially involved in the PHH as confirmed by qRT-PCR. A hypothetical mechanism of this phenomenon and its regulation has been proposed. These findings provide the first comprehensive insights into gene expression in yam tubers after harvest and valuable information for molecular breeding against the PHH.

Genetic parameters, prediction, and selection in a white Guinea yam early-generation breeding population using pedigree information

Better understanding of the genetic control of traits in breeding populations is crucial for the selection of superior varieties and parents. This study aimed to assess genetic parameters and breeding values for six essential traits in a white Guinea yam (Dioscorea rotundata Poir.) breeding population. For this, pedigree-based best linear unbiased prediction (P-BLUP) was used. The results revealed significant nonadditive genetic variances and medium to high (.45-.79) broad-sense heritability estimates for the traits studied. The pattern of associations among the genetic values of the traits suggests that selection based on a multiple-trait selection index has potential for identifying superior breeding lines. Parental breeding values predicted using progeny performance identified 13 clones with high genetic potential for simultaneous improvement of the measured traits in the yam breeding program. Subsets of progeny were identified for intermating or further variety testing based on additive genetic and total genetic values. Selection of the top 5% progenies based on the multi-trait index revealed positive genetic gains for fresh tuber yield (t ha-1), tuber yield (kg plant-1), and average tuber weight (kg). However, genetic gain was negative for tuber dry matter content and Yam mosaic virus resistance in comparison with standard varieties. Our results show the relevance of P-BLUP for the selection of superior parental clones and progenies with higher breeding values for interbreeding and higher genotypic value for variety development in yam.

Population structure, genetic diversity and genomic selection signatures among a Brazilian common bean germplasm

Brazil is the world's largest producer of common bean. Knowledge of the genetic diversity and relatedness of accessions adapted to Brazilian conditions is of great importance for the conservation of germplasm and for directing breeding programs aimed at the development of new cultivars. In this context, the objective of this study was to analyze the genetic diversity, population structure, and linkage disequilibrium (LD) of a diversity panel consisting of 219 common bean accessions, most of which belonging to the Mesoamerican gene pool. Genotyping by sequencing (GBS) of these accessions allowed the identification of 49,817 SNPs with minor allele frequency > 0.05. Of these, 17,149 and 12,876 were exclusive to the Mesoamerican and Andean pools, respectively, and 11,805 SNPs could differentiate the two gene pools. Further the separation according to the gene pool, bayesian analysis of the population structure showed a subdivision of the Mesoamerican accessions based on the origin and color of the seed tegument. LD analysis revealed the occurrence of long linkage blocks and low LD decay with physical distance between SNPs (LD half decay in 249 kb, corrected for population structure and relatedness). The GBS technique could effectively characterize the Brazilian common bean germplasms, and the diversity panel used in this study may be of great use in future genome-wide association studies.

Comparative Phylogeography of Veronica spicata and V. longifolia (Plantaginaceae) Across Europe: Integrating Hybridization and Polyploidy in Phylogeography

Climatic fluctuations in the Pleistocene caused glacial expansion-contraction cycles in Eurasia and other parts of the world. Consequences of these cycles, such as population expansion and subsequent subdivision, have been studied in many taxa at intraspecific population level across much of the Northern Hemisphere. However, the consequences for the potential of hybridization and polyploidization are poorly understood. Here, we investigated the phylogeographic structure of two widespread, closely related species, Veronica spicata and Veronica longifolia, across their European distribution ranges. We assessed the extent and the geographic pattern of polyploidization in both species and hybridization between them. We used genome-scale SNP data to clarify phylogenetic relationships and detect possible hybridization/introgression events. In addition, crossing experiments were performed in different combination between V. spicata and V. longifolia individuals of two ploidy levels and of different geographic origins. Finally, we employed ecological niche modeling to infer macroclimatic differences between both species and both ploidy levels. We found a clear genetic structure reflecting the geographical distribution patterns in both species, with V. spicata showing higher genetic differentiation than V. longifolia. We retrieved significant signals of hybridization and introgression in natural populations from the genetic data and corroborated this with crossing experiments. However, there were no clear phylogeographic patterns and unequivocal macroclimatic niche differences between diploid and tetraploid lineages. This favors the hypothesis, that autopolyploidization has happened frequently and in different regions. The crossing experiments produced viable hybrids when the crosses were made between plants of the same ploidy levels but not in the interploidy crosses. The results suggest that hybridization occurs across the overlapping areas of natural distribution ranges of both species, with apparently directional introgression from V. spicata to V. longifolia. Nevertheless, the two species maintain their species-level separation due to their adaptation to different habitats and spatial isolation rather than reproductive isolation.

Genetic diversity and population structure of Quercus fabri Hance in China revealed by genotyping-by-sequencing

Analysis of genetic diversity and population structure among Quercus fabri populations is essential for the conservation and utilization of Q. fabri resources. Here, the genetic diversity and structure of 158 individuals from 13 natural populations of Quercus fabri in China were analyzed using genotyping-by-sequencing (GBS). A total of 459,564 high-quality single nucleotide polymorphisms (SNPs) were obtained after filtration for subsequent analysis. Genetic structure analysis revealed that these individuals can be clustered into two groups and the structure can be explained mainly by the geographic barrier, showed gene introgression from coastal to inland areas and high mountains could significantly hinder the mutual introgression of genes. Genetic diversity analysis indicated that the individual differences within groups are greater than the differences between the two groups. These results will help us better understand the genetic backgrounds of Q. fabri.

Mycelial biomass and concentration of loline alkaloids driven by complex population structure in Epichloë uncinata and meadow fescue (Schedonorus pratensis)

Many efforts have been made to select and isolate naturally occurring animal-friendly Epichloë strains for later reinfection into elite cultivars. Often this process involves large-scale screening of Epichloë-infected wild grass populations where strains are characterized and alkaloids measured. Here, we describe for the first time the use of genotyping-by-sequencing (GBS) on a collection of 217 Epichloë-infected grasses (7 S. arundinaceum, 4 L. perenne, and 206 S. pratensis). This genotyping strategy is cheaper than complete genome sequencing, is suitable for a large number of individuals, and, when applied to endophyte-infected grasses, conveniently genotypes both organisms. In total, 6273 single nucleotide polymorphisms (SNPs) in the endophyte data set and 38 323 SNPs in the host data set were obtained. Our findings reveal a composite structure with three distinct endophyte clusters unrelated to the three main S. pratensis gene pools that have most likely spread from different glacial refugia in Eurasia. All three gene pools can establish symbiosis with E. uncinata. A comparison of the endophyte clusters with microsatellite-based fingerprinting of the same samples allows a quick test to discriminate between these clusters using two simple sequence repeats (SSRs). Concentrations of loline alkaloids and mycelial biomass are correlated and differ significantly among the plant and endophyte subpopulations; one endophyte strain has higher levels of lolines than others, and one specific host genotype is particularly suitable to host E. uncinata. These findings pave the way for targeted artificial inoculations of specific host-endophyte combinations to boost loline production in the symbiota and for genome association studies with the aim of isolating genes involved in the compatibility between meadow fescue and E. uncinata.

High Contiguity De Novo Genome Sequence Assembly of Trifoliate Yam (Dioscorea dumetorum) Using Long Read Sequencing

Trifoliate yam (Dioscorea dumetorum) is one example of an orphan crop, not traded internationally. Post-harvest hardening of the tubers of this species starts within 24 h after harvesting and renders the tubers inedible. Genomic resources are required for D. dumetorum to improve breeding for non-hardening varieties as well as for other traits. We sequenced the D. dumetorum genome and generated the corresponding annotation. The two haplophases of this highly heterozygous genome were separated to a large extent. The assembly represents 485 Mbp of the genome with an N50 of over 3.2 Mbp. A total of 35,269 protein-encoding gene models as well as 9941 non-coding RNA genes were predicted, and functional annotations were assigned.

Identification of Loci Controlling the Dwarfism Trait in the White Sailfin Molly (Poecilia latipinna) Using Genome-Wide Association Studies Based on Genotyping-By-Sequencing

Dwarfism is a condition defined by low harvest weight in fish, but also results in strange body figures which may have potential for the selective breeding of new ornamental fish strains. The objectives of this study are to reveal the physiological causes of dwarfism and identify the genetic loci controlling this trait in the white sailfin molly. Skeletons of dwarf and normal sailfin mollies were observed by X-ray radioscopy and skeletal staining. Genome-wide association studies based on genotyping-by-sequencing (n = 184) were used to map candidate genomic regions associated with the dwarfism trait. Quantitative real-time PCR was performed to determine the expression level of candidate genes in normal (n = 8) and dwarf (n = 8) sailfin mollies. We found that the dwarf sailfin molly has a short and dysplastic spine in comparison to the normal fish. Two regions, located at NW_015112742.1 and NW_015113621.1, were significantly associated with the dwarfism trait. The expression level of three candidate genes, ADAMTS like 1, Larp7 and PPP3CA, were significantly different between the dwarf and normal sailfin mollies in the hepatopancreas, with PPP3CA also showing significant differences in the vertebrae and Larp7 showing significant differences in the muscle. This study identified genomic regions and candidate genes associated with the dwarfism trait in the white sailfin molly and would provide a reference to determine dwarf-causing variations.