In-depth genome characterization of a Brazilian common bean core collection using DArTseq high-density SNP genotyping

BACKGROUND

Common bean is a legume of social and nutritional importance as a food crop, cultivated worldwide especially in developing countries, accounting for an important source of income for small farmers. The availability of the complete sequences of the two common bean genomes has dramatically accelerated and has enabled new experimental strategies to be applied for genetic research. DArTseq has been widely used as a method of SNP genotyping allowing comprehensive genome coverage with genetic applications in common bean breeding programs.

RESULTS

Using this technology, 6286 SNPs (1 SNP/86.5 Kbp) were genotyped in genic (43.3%) and non-genic regions (56.7%). Genetic subdivision associated to the common bean gene pools (K = 2) and related to grain types (K = 3 and K = 5) were reported. A total of 83% and 91% of all SNPs were polymorphic within the Andean and Mesoamerican gene pools, respectively, and 26% were able to differentiate the gene pools. Genetic diversity analysis revealed an average H E of 0.442 for the whole collection, 0.102 for Andean and 0.168 for Mesoamerican gene pools (F ST  = 0.747 between gene pools), 0.440 for the group of cultivars and lines, and 0.448 for the group of landrace accessions (F ST  = 0.002 between cultivar/line and landrace groups). The SNP effects were predicted with predominance of impact on non-coding regions (77.8%). SNPs under selection were identified within gene pools comparing landrace and cultivar/line germplasm groups (Andean: 18; Mesoamerican: 69) and between the gene pools (59 SNPs), predominantly on chromosomes 1 and 9. The LD extension estimate corrected for population structure and relatedness (r2SV) was ~ 88 kbp, while for the Andean gene pool was ~ 395 kbp, and for the Mesoamerican was ~ 130 kbp.

CONCLUSIONS

For common bean, DArTseq provides an efficient and cost-effective strategy of generating SNPs for large-scale genome-wide studies. The DArTseq resulted in an operational panel of 560 polymorphic SNPs in linkage equilibrium, providing high genome coverage. This SNP set could be used in genotyping platforms with many applications, such as population genetics, phylogeny relation between common bean varieties and support to molecular breeding approaches.

SNP discovery in common bean by restriction-associated DNA (RAD) sequencing for genetic diversity and population structure analysis

Researchers have made great advances into the development and application of genomic approaches for common beans, creating opportunities to driving more real and applicable strategies for sustainable management of the genetic resource towards plant breeding. This work provides useful polymorphic single-nucleotide polymorphisms (SNPs) for high-throughput common bean genotyping developed by RAD (restriction site-associated DNA) sequencing. The RAD tags were generated from DNA pooled from 12 common bean genotypes, including breeding lines of different gene pools and market classes. The aligned sequences identified 23,748 putative RAD-SNPs, of which 3357 were adequate for genotyping; 1032 RAD-SNPs with the highest ADT (assay design tool) score are presented in this article. The RAD-SNPs were structurally annotated in different coding (47.00 %) and non-coding (53.00 %) sequence components of genes. A subset of 384 RAD-SNPs with broad genome distribution was used to genotype a diverse panel of 95 common bean germplasms and revealed a successful amplification rate of 96.6 %, showing 73 % of polymorphic SNPs within the Andean group and 83 % in the Mesoamerican group. A slightly increased He (0.161, n = 21) value was estimated for the Andean gene pool, compared to the Mesoamerican group (0.156, n = 74). For the linkage disequilibrium (LD) analysis, from a group of 580 SNPs (289 RAD-SNPs and 291 BARC-SNPs) genotyped for the same set of genotypes, 70.2 % were in LD, decreasing to 0.10 %in the Andean group and 0.77 % in the Mesoamerican group. Haplotype patterns spanning 310 Mb of the genome (60 %) were characterized in samples from different origins. However, the haplotype frameworks were under-represented for the Andean (7.85 %) and Mesoamerican (5.55 %) gene pools separately. In conclusion, RAD sequencing allowed the discovery of hundreds of useful SNPs for broad genetic analysis of common bean germplasm. From now, this approach provides an excellent panel of molecular tools for whole genome analysis, allowing integrating and better exploring the common bean breeding practices.

Genetic Mapping of Resistance to Meloidogyne arenaria in Arachis stenosperma: A New Source of Nematode Resistance for Peanut

Root-knot nematodes (RKN; Meloidogyne sp.) are a major threat to crops in tropical and subtropical regions worldwide. The use of resistant crop varieties is the preferred method of control because nematicides are expensive, and hazardous to humans and the environment. Peanut (Arachis hypogaea) is infected by four species of RKN, the most damaging being M. arenaria, and commercial cultivars rely on a single source of resistance. In this study, we genetically characterize RKN resistance of the wild Arachis species A. stenosperma using a population of 93 recombinant inbred lines developed from a cross between A. duranensis and A. stenosperma. Four quantitative trait loci (QTL) located on linkage groups 02, 04, and 09 strongly influenced nematode root galling and egg production. Drought-related, domestication and agronomically relevant traits were also evaluated, revealing several QTL. Using the newly available Arachis genome sequence, easy-to-use KASP (kompetitive allele specific PCR) markers linked to the newly identified RKN resistance loci were developed and validated in a tetraploid context. Therefore, we consider that A. stenosperma has high potential as a new source of RKN resistance in peanut breeding programs.

Expression of drought tolerance genes in tropical upland rice cultivars (Oryza sativa)

Gene expression related to drought response in the leaf tissues of two Brazilian upland cultivars, the drought-tolerant Douradão and the drought-sensitive Primavera, was analyzed. RNA-seq identified 27,618 transcripts in the Douradão cultivar, with 24,090 (87.2%) homologous to the rice database, and 27,221 transcripts in the Primavera cultivar, with 23,663 (86.9%) homologous to the rice database. Gene-expression analysis between control and water-deficient treatments revealed 493 and 1154 differentially expressed genes in Douradão and Primavera cultivars, respectively. Genes exclusively expressed under drought were identified for Douradão, including two genes of particular interest coding for the protein peroxidase precursor, which is involved in three distinct metabolic pathways. Comparisons between the two drought-exposed cultivars revealed 2314 genes were differentially expressed (978 upregulated, 1336 downregulated in Douradão). Six genes distributed across 4 different transcription factor families (bHLH, MYB, NAC, and WRKY) were identified, all of which were upregulated in Douradão compared to Primavera during drought. Most of the genes identified in Douradão activate metabolic pathways responsible for production of secondary metabolites and genes coding for enzymatically active signaling receptors. Quantitative PCR validation showed that most gene expression was in agreement with computational prediction of these transcripts. The transcripts identified here will define molecular markers for identification of Cis-acting elements to search for allelic variants of these genes through analysis of polymorphic SNPs in GenBank accessions of upland rice, aiming to develop cultivars with the best combination of these alleles, resulting in materials with high yield potential in the event of drought during the reproductive phase.

Discrimination of common bean cultivars using multiplexed microsatellite markers

Analysis of DNA polymorphisms allows for the genetic identification and precise discrimination of species with a narrow genetic base such as common bean. The primary objectives of the present study were to molecularly characterize commercial common bean varieties developed at various research institutions using microsatellite markers and to determine the degree of genetic diversity among the bean varieties analyzed. Fifty cultivars representing 12 grain classes and 64 genitors, i.e., accessions used to develop these cultivars, were characterized. Based on an analysis of 24 simple sequence repeats, the estimates for the average number of alleles and genetic diversity were 8.29 and 0.646, respectively. The combined probability of identity was estimated at 7.05 x 10(-17), indicating a high individual discriminatory power. Thirty-two percent of the cultivars exhibited heterogeneity for multiple loci that reflected either homozygosity for different alleles of a given locus in different individuals or heterozygosity for the locus. The average genetic diversity for the groups of cultivars and genitors was 0.605 and 0.660, respectively, with no genetic differentiation (FST) between these groups. Although similar estimates of expected heterozygosity were observed when the cultivars were grouped by release date, a greater number of private alleles was observed in the most recent cultivars. The genetic differentiation among cultivars originating from different institutions was not different from zero (FST = 0.01). The molecular profile database derived from these analyses may increase the statistical power of genetic estimates and may be incorporated into breeding programs for common bean. Furthermore, the profiles obtained for the different cultivars may be used as molecular descriptors to complement traditional descriptors used in distinctiveness, uniformity and stability tests, thereby improving the traceability of samples and their derivatives and helping to protect the intellectual property rights of breeders.

Molecular characterization of high performance inbred lines of Brazilian common beans

The identification of germplasm genetic variability in breeding programs of the common bean (Phaseolus vulgaris) is essential for determining the potential of each combination of parent plants to obtain superior genotypes. The present study aimed to estimated the extent of genetic diversity in 172 lineages and cultivars of the common bean by integrating five tests of value for cultivation and use (VCU) that were conducted over the last eight years by the breeding program of Embrapa Arroz e Feijão in Brazil. Nine multilocus genotyping systems composed of 36 fluorescent microsatellite markers distributed across 11 different chromosomes of the common bean were used, of which 24 were polymorphic in all trials. One hundred and eighty-seven alleles were identified, with an average of 7.79 alleles per locus and an average gene diversity of 0.65. The combined probability of identity for all loci was 1.32 x 10(-16). Lineages that are more genetically divergent between the selection cycles were identified, allowing the breeding program to develop a crossbreed between elite genotypes with a low degree of genetic relatedness. HE values ranged from 0.31 to 0.63, with a large reduction in the genetic base over successive selection cycles. The test showed a significant degree of differentiation (FST = 0.159). Private alleles (26%) were identified and can be directly incorporated into the gene pool of cultivated germplasm, thereby contributing effectively to the expansion of genetic diversity in this bean-breeding program.

Development and characterization of new microsatellites for Eugenia dysenterica DC (Myrtaceae)

Microsatellite markers were developed for population genetic analyses of the Neotropical tree Eugenia dysenterica DC (Myrtaceae), after construction of a shotgun genomic library for microsatellite discovery. Nine primers were designed, of which 5 yielded amplified product. These primers were polymorphic for 97 individuals collected in 3 distinct localities. The number of alleles per locus (primer) ranged from 3 to 11 and expected heterozygosities varied from 0.309 to 0.884. The probability of locus identity was ~1.88 x 10(-4) and the probability of paternity exclusion was ~0.9367. The 5 microsatellite primer pairs may be suitable for population genetic studies such as parentage and fine-scale genetic analyses of this species.

Development of microsatellite markers for the endangered Neotropical tree species Tibouchina papyrus (Melastomataceae)

We isolated and characterized 12 microsatellite loci for Tibouchina papyrus (Melastomataceae), an endangered species with narrow and disjunct range, endemics to a few localities in "cerrado rupestre" from Central Brazil. These microsatellites were obtained by sequencing of a genomic shotgun library for primer design. Leaves from 96 individuals collected in the three known local populations were genotyped using the 12 primers designed to analyze the polymorphisms at each locus. The number of alleles per locus ranged from one to six; two loci were monomorphic. Among the polymorphic loci, expected heterozygosities ranged from 0.161 to 0.714. Combined paternity exclusion probability was 0.957 and combined genetic identity (0.051) was high for studies on parentage. Tibouchina papyrus is a rare and endemic tree species of outcrop quartzite and sandstone soils, with highly isolated populations, which may have lead to the low degree of polymorphism that we detected. Also, motifs of most loci are larger than dinucleotide, which typically display lower levels of polymorphism.