Genome-Wide Association Study of Agronomic and Physiological Traits Related to Drought Tolerance in Potato.

Potato (Solanum tuberosum L.) is often considered a water-sensitive crop and its production can be threatened by drought events, making water stress tolerance a trait of increasing interest. In this study, a panel of 144 tetraploid potato genotypes was evaluated for two consecutive years (2019 and 2020) to observe the variation of several physiological traits such as chlorophyll content and fluorescence, stomatal conductance, NDVI, and leaf area and circumference. In addition, agronomic parameters such as yield, tuber fresh weight, tuber number, starch content, dry matter and reducing sugars were determined. GGP V3 Potato array was used to genotype the population, obtaining a total of 18,259 high-quality SNP markers. Marker-trait association was performed using GWASpoly package in R software and Q + K linear mixed models were considered. This approach allowed us to identify eighteen SNP markers significantly associated with the studied traits in both treatments and years, which were related to genes with known functions. Markers related to chlorophyll content and number of tubers under control and stress conditions, and related to stomatal conductance, NDVI, yield and reducing sugar content under water stress, were identified. Although these markers were distributed throughout the genome, the SNPs associated with the traits under control conditions were found mainly on chromosome 11, while under stress conditions they were detected on chromosome 4. These results contribute to the knowledge of the mechanisms of potato tolerance to water stress and are useful for future marker-assisted selection programs.

Transcriptome Analysis of Two Tetraploid Potato Varieties under Water-Stress Conditions.

Potato (Solanum tuberosum L.) is one of the most important crops worldwide, but due to its sensitivity to drought, its production can be affected by water availability. In this study, the varieties Agria and Zorba were used to determine the expression differences between control and water-stressed plants. For this purpose, they were sequenced by RNAseq, obtaining around 50 million transcripts for each variety and treatment. When comparing the significant transcripts obtained from control and drought-stressed plants of the Agria variety, we detected 931 genes that were upregulated and 2077 genes that were downregulated under stress conditions. When both treatments were compared in Zorba plants, 735 genes were found to be upregulated and 923 genes were found to be downregulated. Significantly more DEGs were found in the Agria variety, indicating a good stress response of this variety. "Abscisic acid and environmental stress-inducible protein TAS14-like" was the most overexpressed gene under drought conditions in both varieties, but expression differences were also found in numerous transcription factors and heat shock proteins. The principal GO term found was "cellular components", more specifically related to the cell membrane and the cell wall, but other metabolic pathways such as carbohydrate metabolism and osmotic adjustment were also identified. These results provide valuable information related to the molecular mechanisms of tolerance to water stress in order to establish the basis for breeding new, more tolerant varieties.

Analysis of genetic diversity and discrimination of Oil Palm DxP populations based on the origins of pisifera elite parents.

A total of 251 Dura cross Pisifera (DxP) hybrid palms from six populations descending from six parental African Pisifera origins and involving 12 progenies were analyzed with 19 selected Simple Sequence Repeats (SSR) markers. A total of 110 alleles were produced, ranging from three to eight per SSR, with a mean of 5.8 alleles per SSR locus. Of these, 68.5% were considered shared alleles by more than one population and the remaining 31.5% were population specific alleles. They generated between six and 21 haplotypes in all populations, and depending on the SSR marker, between one and 10 haplotypes within populations. Various parameters for analyzing genetic variability, differentiation and genetic structure were computed using GenAlEx, Structure and Darwin software. The obtained results confirmed microsatellites as a robust, feasible and trustful method for obtaining DNA fingerprints, tracing the source of oil palm samples. With respect to the authenticity of materials or for solving legitimacy issues, accession belonging to each population by SSR markers could be distinguished, but additional SSR should be screened for improving this process.

Construction of a high density linkage map in Oil Palm using SPET markers.

A high-density genetic linkage map from a controlled cross of two oil palm (Elaeis guineensis) genotypes was constructed based on Single Primer Enrichment Technology (SPET) markers. A 5K panel of hybridization probes were used for this purpose which was derived from previously developed SNP primers in oil palm. Initially, 13,384 SNPs were detected which were reduced to 13,073 SNPs after filtering for only bi-allelic SNP. Around 75% of the markers were found to be monomorphic in the progeny, reducing the markers left for linkage mapping to 3,501. Using Lep-MAP3 software, a linkage map was constructed which contained initially 2,388 markers and had a total length of 1,370 cM. In many cases several adjacent SNP were located on the same locus, due to missing recombination events between them, leading to a total of 1,054 loci on the 16 LG. Nevertheless, the marker density of 1.74 markers per cM (0.57 cM/marker) should allow the detection of QTLs in the future. This study shows that cost efficient SPET markers are suitable for linkage map construction in oil palm and probably, also in other species.

Detection of significant SNP associated with production and oil quality traits in interspecific oil palm hybrids using RARSeq.

A RARSeq based Association mapping study was performed in a population of 104 Elaeis oleifera x E. guineensis hybrids of five origins with the aim of finding functional markers associated to six productive and 19 oil quality traits. For this purpose mRNA of each genotype was isolated and double stranded cDNA was synthesized. Following digestion with two restriction enzymes and adapter ligation, a size selected pool of barcoded amplicons was produced and sequenced using Illumina MiSeq. The obtained sequences were processed with a "snakemake" pipeline, filtered and missing values were imputed. For all traits except two significant effects of the origin was observed. Genetic diversity analyses revealed high variability within origins and an excess of heterozygosity in the population. Two GLM models with Q matrix or PCA matrix as covariates and two MLM models incorporating in addition a Kinship matrix were tested for genotype-phenotype associations using GAPIT software. Using unadjusted p values (< 0.01) 78 potential associations were detected involving 25 SNP and 20 traits. When applying FDR multiple testing with p < 0.05, 25 significant associations remained involving eight SNP and six quality traits. Four SNP were located in genes with a potential relevant biological meaning.

Association Mapping Between Candidate Gene SNP and Production and Oil Quality Traits in Interspecific Oil Palm Hybrids.

Oil palm production is gaining importance in Central and South America. However, the main species Elaeis guineensis (Eg) is suffering severely from bud rod disease, restricting the potential cultivation areas. Therefore, breeding companies have started to work with interspecific Elaeis oleifera × Eg (Eo × Eg) hybrids which are tolerant to this disease. We performed association studies between candidate gene (CG) single nucleotide polymorphisms (SNP) and six production and 19 oil quality traits in 198 accessions of interspecific oil palm hybrids from five different origins. For this purpose, barcoded amplicons of initially 167 CG were produced from each genotype and sequenced with Ion Torrent. After sequence cleaning 115 SNP remained targeting 62 CG. The influence of the origins on the different traits was analyzed and a genetic diversity study was performed. Two generalized linear models (GLM) with principle component analysis (PCA) or structure (Q) matrixes as covariates and two mixed linear models (MLM) which included in addition a Kinship (K) matrix were applied for association mapping using GAPIT. False discovery rate (FDR) multiple testing corrections were applied in order to avoid Type I errors. However, with FDR adjusted p values no significant associations between SNP and traits were detected. If using unadjusted p values below 0.05, seven of the studied CG showed potential associations with production traits, while 23 CG may influence different quality traits. Under these conditions the current approach and the detected candidate genes could be exploited for selecting genotypes with superior CG alleles in Marker Assisted Selection systems.

Physical mapping of QTL for tuber yield, starch content and starch yield in tetraploid potato (Solanum tuberosum L.) by means of genome wide genotyping by sequencing and the 8.3 K SolCAP SNP array.

BACKGROUND: Tuber yield and starch content of the cultivated potato are complex traits of decisive importance for breeding improved varieties. Natural variation of tuber yield and starch content depends on the environment and on multiple, mostly unknown genetic factors. Dissection and molecular identification of the genes and their natural allelic variants controlling these complex traits will lead to the development of diagnostic DNA-based markers, by which precision and efficiency of selection can be increased (precision breeding). RESULTS: Three case-control populations were assembled from tetraploid potato cultivars based on maximizing the differences between high and low tuber yield (TY), starch content (TSC) and starch yield (TSY, arithmetic product of TY and TSC). The case-control populations were genotyped by restriction-site associated DNA sequencing (RADseq) and the 8.3 k SolCAP SNP genotyping array. The allele frequencies of single nucleotide polymorphisms (SNPs) were compared between cases and controls. RADseq identified, depending on data filtering criteria, between 6664 and 450 genes with one or more differential SNPs for one, two or all three traits. Differential SNPs in 275 genes were detected using the SolCAP array. A genome wide association study using the SolCAP array on an independent, unselected population identified SNPs associated with tuber starch content in 117 genes. Physical mapping of the genes containing differential or associated SNPs, and comparisons between the two genome wide genotyping methods and two different populations identified genome segments on all twelve potato chromosomes harboring one or more quantitative trait loci (QTL) for TY, TSC and TSY. CONCLUSIONS: Several hundred genes control tuber yield and starch content in potato. They are unequally distributed on all potato chromosomes, forming clusters between 0.5-4 Mbp width. The largest fraction of these genes had unknown function, followed by genes with putative signalling and regulatory functions. The genetic control of tuber yield and starch content is interlinked. Most differential SNPs affecting both traits had antagonistic effects: The allele increasing TY decreased TSC and vice versa. Exceptions were 89 SNP alleles which had synergistic effects on TY, TSC and TSY. These and the corresponding genes are primary targets for developing diagnostic markers.

First genetic linkage map of Taraxacum koksaghyz Rodin based on AFLP, SSR, COS and EST-SSR markers.

Taraxacum koksaghyz Rodin (TKS) has been studied in many occasions as a possible alternative source for natural rubber production of good quality and for inulin production. Some tire companies are already testing TKS tire prototypes. There are also many investigations on the production of bio-fuels from inulin and inulin applications for health improvement and in the food industry. A limited amount of genomic resources exist for TKS and particularly no genetic linkage map is available in this species. We have constructed the first TKS genetic linkage map based on AFLP, COS, SSR and EST-SSR markers. The integrated linkage map with eight linkage groups (LG), representing the eight chromosomes of Russian dandelion, has 185 individual AFLP markers from parent 1, 188 individual AFLP markers from parent 2, 75 common AFLP markers and 6 COS, 1 SSR and 63 EST-SSR loci. Blasting the EST-SSR sequences against known sequences from lettuce allowed a partial alignment of our TKS map with a lettuce map. Blast searches against plant gene databases revealed some homologies with useful genes for downstream applications in the future.

Genetic architecture of palm oil fatty acid composition in cultivated oil palm (Elaeis guineensis Jacq.) compared to its wild relative E. oleifera (H.B.K) Cortés.

We searched for quantitative trait loci (QTL) associated with the palm oil fatty acid composition of mature fruits of the oil palm E. guineensis Jacq. in comparison with its wild relative E. oleifera (H.B.K) Cortés. The oil palm cross LM2T x DA10D between two heterozygous parents was considered in our experiment as an intraspecific representative of E. guineensis. Its QTLs were compared to QTLs published for the same traits in an interspecific Elaeis pseudo-backcross used as an indirect representative of E. oleifera. Few correlations were found in E. guineensis between pulp fatty acid proportions and yield traits, allowing for the rather independent selection of both types of traits. Sixteen QTLs affecting palm oil fatty acid proportions and iodine value were identified in oil palm. The phenotypic variation explained by the detected QTLs was low to medium in E. guineensis, ranging between 10% and 36%. The explained cumulative variation was 29% for palmitic acid C16:0 (one QTL), 68% for stearic acid C18:0 (two QTLs), 50% for oleic acid C18:1 (three QTLs), 25% for linoleic acid C18:2 (one QTL), and 40% (two QTLs) for the iodine value. Good marker co-linearity was observed between the intraspecific and interspecific Simple Sequence Repeat (SSR) linkage maps. Specific QTL regions for several traits were found in each mapping population. Our comparative QTL results in both E. guineensis and interspecific materials strongly suggest that, apart from two common QTL zones, there are two specific QTL regions with major effects, which might be one in E. guineensis, the other in E. oleifera, which are independent of each other and harbor QTLs for several traits, indicating either pleiotropic effects or linkage. Using QTL maps connected by highly transferable SSR markers, our study established a good basis to decipher in the future such hypothesis at the Elaeis genus level.

Construction of a 10,000-marker ultradense genetic recombination map of potato: providing a framework for accelerated gene isolation and a genomewide physical map.

An ultradense genetic linkage map with >10,000 AFLP loci was constructed from a heterozygous diploid potato population. To our knowledge, this is the densest meiotic recombination map ever constructed. A fast marker-ordering algorithm was used, based on the minimization of the total number of recombination events within a given marker order in combination with genotyping error-detection software. This resulted in "skeleton bin maps," which can be viewed as the most parsimonious marker order. The unit of distance is not expressed in centimorgans but in "bins." A bin is a position on the genetic map with a unique segregation pattern that is separated from adjacent bins by a single recombination event. Putative centromeres were identified by a strong clustering of markers, probably due to cold spots for recombination. Conversely, recombination hot spots resulted in large intervals of up to 15 cM without markers. The current level of marker saturation suggests that marker density is proportional to physical distance and independent of recombination frequency. Most chromatids (92%) recombined once or never, suggesting strong chiasma interference. Absolute chiasma interference within a chromosome arm could not be demonstrated. Two examples of contig construction and map-based cloning have demonstrated that the marker spacing was in accordance with the expected physical distance: approximately one marker per BAC length. Currently, the markers are used for genetic anchoring of a physical map of potato to deliver a sequence-ready minimal tiling path of BAC contigs of specific chromosomal regions for the potato genome sequencing consortium (http://www.potatogenome.net).

A linkage map of an F2 hybrid population of Antirrhinum majus and A. molle.

To increase the utility of Antirrhinum for genetic and evolutionary studies, we constructed a molecular linkage map for an interspecific hybrid A. majus x A. molle. An F(2) population (n = 92) was genotyped at a minimum of 243 individual loci. Although distorted transmission ratios were observed at marker loci throughout the genome, a mapping strategy based on a fixed framework of codominant markers allowed the loci to be placed into eight robust linkage groups consistent with the haploid chromosome number of Antirrhinum. The mapped loci included 164 protein-coding genes and a similar number of unknown sequences mapped as AFLP, RFLP, ISTR, and ISSR markers. Inclusion of sequences from mutant loci allowed provisional alignment of classical and molecular linkage groups. The total map length was 613 cM with an average interval of 2.5 cM, but most of the loci were aggregated into clusters reducing the effective distance between markers. Potential causes of transmission ratio distortion and its effects on map construction were investigated. This first molecular linkage map for Antirrhinum should facilitate further mapping of mutations, major QTL, and other coding sequences in this model genus.

Toward a marker-dense meiotic map of the potato genome: lessons from linkage group I.

Segregation data were obtained for 1260 potato linkage group I-specific AFLP loci from a heterozygous diploid potato population. Analytical tools that identified potential typing errors and/or inconsistencies in the data and that assembled cosegregating markers into bins were applied. Bins contain multiple-marker data sets with an identical segregation pattern, which is defined as the bin signature. The bin signatures were used to construct a skeleton bin map that was based solely on observed recombination events. Markers that did not match any of the bin signatures exactly (and that were excluded from the calculation of the skeleton bin map) were placed on the map by maximum likelihood. The resulting maternal and paternal maps consisted of 95 and 101 bins, respectively. Markers derived from EcoRI/MseI, PstI/MseI, and SacI/MseI primer combinations showed different genetic distributions. Approximately three-fourths of the markers placed into a bin were considered to fit well on the basis of an estimated residual "error rate" of 0-3%. However, twice as many PstI-based markers fit badly, suggesting that parental PstI-site methylation patterns had changed in the population. Recombination frequencies were highly variable across the map. Inert, presumably centromeric, regions caused extensive marker clustering while recombination hotspots (or regions identical by descent) resulted in empty bins, despite the level of marker saturation.

Cold sweetening in diploid potato: mapping quantitative trait loci and candidate genes.

A candidate gene approach has been used as a first step to identify the molecular basis of quantitative trait variation in potato. Sugar content of tubers upon cold storage was the model trait chosen because the metabolic pathways involved in starch and sugar metabolism are well known and many of the genes have been cloned. Tubers of two F(1) populations of diploid potato grown in six environments were evaluated for sugar content after cold storage. The populations were genotyped with RFLP, AFLP, and candidate gene markers. QTL analysis revealed that QTL for glucose, fructose, and sucrose content were located on all potato chromosomes. Most QTL for glucose content mapped to the same positions as QTL for fructose content. QTL explaining >10% of the variability for reducing sugars were located on linkage groups I, III, VII, VIII, IX, and XI. QTL consistent across populations and/or environments were identified. QTL were linked to genes encoding invertase, sucrose synthase 3, sucrose phosphate synthase, ADP-glucose pyrophosphorylase, sucrose transporter 1, and a putative sucrose sensor. The results suggest that allelic variants of enzymes operating in carbohydrate metabolic pathways contribute to the genetic variation in cold sweetening.

Quantitative trait loci controlling vegetative growth rate in the edible basidiomycete Pleurotus ostreatus.

Mycelium growth rate is a quantitative characteristic that exhibits continuous variation. This trait has applied interest, as growth rate is correlated with production yield and increased advantage against competitors. In this work, we studied growth rate variation in the edible basidiomycete Pleurotus ostreatus growing as monokaryotic or dikaryotic mycelium on Eger medium or on wheat straw. Our analysis resulted in identification of several genomic regions (quantitative trait loci [QTLs]) involved in the control of growth rate that can be mapped on the genetic linkage map of this fungus. In some cases monokaryotic and dikaryotic QTLs clustered at the same map position, indicating that there are principal genomic areas responsible for growth rate control. The availability of this linkage map of growth rate QTLs can help in the design of rational strain breeding programs based on genomic information.