Marie Denis

Genomic selection prediction accuracy in a perennial crop: case study of oil palm (Elaeis guineensis Jacq.).

Key messageGenomic selection empirically appeared valuable for reciprocal recurrent selection in oil palm as it could account for family effects and Mendelian sampling terms, despite small populations and low marker density.AbstractGenomic selection (GS) can increase the genetic gain in plants. In perennial crops, this is expected mainly through shortened breeding cycles and increased selection intensity, which requires sufficient GS accuracy in selection candidates, despite often small training populations. Our objective was to obtain the first empirical estimate of GS accuracy in oil palm (Elaeis guineensis), the major world oil crop. We used two parental populations involved in conventional reciprocal recurrent selection (Deli and Group B) with 131 individuals each, genotyped with 265 SSR. We estimated within-population GS accuracies when predicting breeding values of non-progeny-tested individuals for eight yield traits. We used three methods to sample training sets and five statistical methods to estimate genomic breeding values. The results showed that GS could account for family effects and Mendelian sampling terms in Group B but only for family effects in Deli. Presumably, this difference between populations originated from their contrasting breeding history. The GS accuracy ranged from −0.41 to 0.94 and was positively correlated with the relationship between training and test sets. Training sets optimized with the so-called CDmean criterion gave the highest accuracies, ranging from 0.49 (pulp to fruit ratio in Group B) to 0.94 (fruit weight in Group B). The statistical methods did not affect the accuracy. Finally, Group B could be preselected for progeny tests by applying GS to key yield traits, therefore increasing the selection intensity. Our results should be valuable for breeding programs with small populations, long breeding cycles, or reduced effective size.

Experimental assessment of the accuracy of genomic selection in sugarcane

Sugarcane cultivars are interspecific hybrids with an aneuploid, highly heterozygous polyploid genome. The complexity of the sugarcane genome is the main obstacle to the use of marker-assisted selection in sugarcane breeding. Given the promising results of recent studies of plant genomic selection, we explored the feasibility of genomic selection in this complex polyploid crop. Genetic values were predicted in two independent panels, each composed of 167 accessions representing sugarcane genetic diversity worldwide. Accessions were genotyped with 1,499 DArT markers. One panel was phenotyped in Reunion Island and the other in Guadeloupe. Ten traits concerning sugar and bagasse contents, digestibility and composition of the bagasse, plant morphology, and disease resistance were used. We used four statistical predictive models: bayesian LASSO, ridge regression, reproducing kernel Hilbert space, and partial least square regression. The accuracy of the predictions was assessed through the correlation between observed and predicted genetic values by cross validation within each panel and between the two panels. We observed equivalent accuracy among the four predictive models for a given trait, and marked differences were observed among traits. Depending on the trait concerned, within-panel cross validation yielded median correlations ranging from 0.29 to 0.62 in the Reunion Island panel and from 0.11 to 0.5 in the Guadeloupe panel. Cross validation between panels yielded correlations ranging from 0.13 for smut resistance to 0.55 for brix. This level of correlations is promising for future implementations. Our results provide the first validation of genomic selection in sugarcane.

Maternal Early Pregnancy Serum Metabolites and Risk of Gestational Diabetes Mellitus

CONTEXT Significant gaps remain in the understanding of genetic and environmental risk factors, as well as related mechanisms that contribute to gestational diabetes mellitus (GDM). OBJECTIVES This study aimed to investigate early pregnancy maternal serum metabolites and subsequent risk of GDM. DESIGN Information on participant characteristics and GDM diagnosis was collected using in-person interviews and medical record abstraction, respectively. Early pregnancy serum samples were used for nontargeted metabolite profiling using a gas chromatography-mass spectrometry platform. Lasso regression was used to select a set of metabolites that are jointly associated with GDM case-control status. We evaluated the predictive performance of the set of selected metabolites using a receiver operating characteristics curve and area under the curve. PARTICIPANTS A total of 178 GDM cases and 180 controls participated in a pregnancy cohort study. RESULTS A set of 17 metabolites (linoleic acid, oleic acid, myristic acid, d-galactose, d-sorbitol, o-phosphocolamine, l-alanine, l-valine, 5-hydroxy-l-tryptophan, l-serine, sarcosine, l-pyroglutamic acid, l-mimosine, l-lactic acid, glycolic acid, fumaric acid, and urea) differentiated GDM cases from controls. Fold changes of relative abundance of these metabolites among GDM cases compared with controls ranged from 1.47 (linoleic acid) to 0.78 (5-hydroxy-l-tryptophan). Addition of these selected metabolites to a set of well-known GDM risk factors improved the area under the curve significantly from 0.71 to 0.87 (P = 3.97E-07). CONCLUSIONS We identified combinations of metabolites in early pregnancy that are associated with subsequent risk of GDM. Replication of findings may improve understanding of GDM pathogenesis and may have implications for the design of GDM prevention and early diagnosis protocols.

Placental genome and maternal-placental genetic interactions: a genome-wide and candidate gene association study of placental abruption.

While available evidence supports the role of genetics in the pathogenesis of placental abruption (PA), PA-related placental genome variations and maternal-placental genetic interactions have not been investigated. Maternal blood and placental samples collected from participants in the Peruvian Abruptio Placentae Epidemiology study were genotyped using Illumina’s Cardio-Metabochip platform. We examined 118,782 genome-wide SNPs and 333 SNPs in 32 candidate genes from mitochondrial biogenesis and oxidative phosphorylation pathways in placental DNA from 280 PA cases and 244 controls. We assessed maternal-placental interactions in the candidate gene SNPS and two imprinted regions (IGF2/H19 and C19MC). Univariate and penalized logistic regression models were fit to estimate odds ratios. We examined the combined effect of multiple SNPs on PA risk using weighted genetic risk scores (WGRS) with repeated ten-fold cross-validations. A multinomial model was used to investigate maternal-placental genetic interactions. In placental genome-wide and candidate gene analyses, no SNP was significant after false discovery rate correction. The top genome-wide association study (GWAS) hits were rs544201, rs1484464 (CTNNA2), rs4149570 (TNFRSF1A) and rs13055470 (ZNRF3) (p-values: 1.11e-05 to 3.54e-05). The top 200 SNPs of the GWAS overrepresented genes involved in cell cycle, growth and proliferation. The top candidate gene hits were rs16949118 (COX10) and rs7609948 (THRB) (p-values: 6.00e-03 and 8.19e-03). Participants in the highest quartile of WGRS based on cross-validations using SNPs selected from the GWAS and candidate gene analyses had a 8.40-fold (95% CI: 5.8–12.56) and a 4.46-fold (95% CI: 2.94–6.72) higher odds of PA compared to participants in the lowest quartile. We found maternal-placental genetic interactions on PA risk for two SNPs in PPARG (chr3∶12313450 and chr3∶12412978) and maternal imprinting effects for multiple SNPs in the C19MC and IGF2/H19 regions. Variations in the placental genome and interactions between maternal-placental genetic variations may contribute to PA risk. Larger studies may help advance our understanding of PA pathogenesis.

Genome-Wide Prediction Methods in Highly Diverse and Heterozygous Species: Proof-of-Concept through Simulation in Grapevine

Nowadays, genome-wide association studies (GWAS) and genomic selection (GS) methods which use genome-wide marker data for phenotype prediction are of much potential interest in plant breeding. However, to our knowledge, no studies have been performed yet on the predictive ability of these methods for structured traits when using training populations with high levels of genetic diversity. Such an example of a highly heterozygous, perennial species is grapevine. The present study compares the accuracy of models based on GWAS or GS alone, or in combination, for predicting simple or complex traits, linked or not with population structure. In order to explore the relevance of these methods in this context, we performed simulations using approx 90,000 SNPs on a population of 3,000 individuals structured into three groups and corresponding to published diversity grapevine data. To estimate the parameters of the prediction models, we defined four training populations of 1,000 individuals, corresponding to these three groups and a core collection. Finally, to estimate the accuracy of the models, we also simulated four breeding populations of 200 individuals. Although prediction accuracy was low when breeding populations were too distant from the training populations, high accuracy levels were obtained using the sole core-collection as training population. The highest prediction accuracy was obtained (up to 0.9) using the combined GWAS-GS model. We thus recommend using the combined prediction model and a core-collection as training population for grapevine breeding or for other important economic crops with the same characteristics.

Estimation of genealogical coancestry in plant species using a pedigree reconstruction algorithm and application to an oil palm breeding population

AbstractKey messageExplicit pedigree reconstruction by simulated annealing gave reliable estimates of genealogical coancestry in plant species, especially when selfing rate was lower than 0.6, using a realistic number of markers. Genealogical coancestry information is crucial in plant breeding to estimate genetic parameters and breeding values. The approach of Fernández and Toro (Mol Ecol 15:1657–1667, 2006) to estimate genealogical coancestries from molecular data through pedigree reconstruction was limited to species with separate sexes. In this study it was extended to plants, allowing hermaphroditism and monoecy, with possible selfing. Moreover, some improvements were made to take previous knowledge on the population demographic history into account. The new method was validated using simulated and real datasets. Simulations showed that accuracy of estimates was high with 30 microsatellites, with the best results obtained for selfing rates below 0.6. In these conditions, the root mean square error (RMSE) between the true and estimated genealogical coancestry was small (<0.07), although the number of ancestors was overestimated and the selfing rate could be biased. Simulations also showed that linkage disequilibrium between markers and departure from the Hardy–Weinberg equilibrium in the founder population did not affect the efficiency of the method. Real oil palm data confirmed the simulation results, with a high correlation between the true and estimated genealogical coancestry (>0.9) and a low RMSE (<0.08) using 38 markers. The method was applied to the Deli oil palm population for which pedigree data were scarce. The estimated genealogical coancestries were highly correlated (>0.9) with the molecular coancestries using 100 markers. Reconstructed pedigrees were used to estimate effective population sizes. In conclusion, this method gave reliable genealogical coancestry estimates. The strategy was implemented in the software MOLCOANC 3.0.

Genomic selection in tree breeding: testing accuracy of prediction models including dominance effect

Background The concept of Marker Assisted Selection (MAS) is rapidly evolving in animal and plant breeding. With the advent of high throughput molecular technology, numerous molecular markers distributed throughout the whole genome can be produced to characterize many genetic entries involving new perspectives in methodology of selection. An important research activity has begun in the animal world given the first theoretical framework for a methodology called genomic selection (GS) [1]. Several statistical approaches have been proposed for the prediction of genomic breeding values and numerous results are available that validates the interest of this method in animal breeding. In plants the GS is still limited to very advanced model species involved in genetic improvement and especially from scenario-based simulation [2,3]. In tree breeding the GS could significantly reduce the cost of genetic improvement schemes by limiting the size and number of field experiments; and facilitating the early selection at the nursery stage [4]. If most of the studies on GS have addressed the prediction of breeding value, taking into account the gene additive effects, there is still a lack of analyses dealing with the total genetic value (genotypic value) including both additive and dominance effects. This aspect is important in plant and especially in tree breeding where the goal of some programs is the production of clones or elite families. The aim of this study is to investigate the performance of GS in the context of tree breeding when the selection is based on genotypic value. The proposed approach allows taking into account both additive and dominance effect [5] for each marker in the statistical model. Six scenarios are simulated to test the reliability of the GS in the frame of recurrent selection scheme.

Identification of Ganoderma Disease Resistance Loci Using Natural Field Infection of an Oil Palm Multiparental Population

Multi-parental populations are promising tools for identifying quantitative disease resistance loci. Stem rot caused by Ganoderma boninense is a major threat to palm oil production, with yield losses of up to 80% prompting premature replantation of palms. There is evidence of genetic resistance sources, but the genetic architecture of Ganoderma resistance has not yet been investigated. This study aimed to identify Ganoderma resistance loci using an oil palm multi-parental population derived from nine major founders of ongoing breeding programs. A total of 1200 palm trees of the multi-parental population was planted in plots naturally infected by Ganoderma, and their health status was assessed biannually over 25 yr. The data were treated as survival data, and modeled using the Cox regression model, including a spatial effect to take the spatial component in the spread of Ganoderma into account. Based on the genotypes of 757 palm trees out of the 1200 planted, and on pedigree information, resistance loci were identified using a random effect with identity-by-descent kinship matrices as covariance matrices in the Cox model. Four Ganoderma resistance loci were identified, two controlling the occurrence of the first Ganoderma symptoms, and two the death of palm trees, while favorable haplotypes were identified among a major gene pool for ongoing breeding programs. This study implemented an efficient and flexible QTL mapping approach, and generated unique valuable information for the selection of oil palm varieties resistant to Ganoderma disease.

Mixed model approach for IBD-based QTL mapping in a complex oil palm pedigree

BackgroundElaeis guineensis is the world’s leading source of vegetable oil, and the demand is still increasing. Oil palm breeding would benefit from marker-assisted selection but genetic studies are scarce and inconclusive. This study aims to identify genetic bases of oil palm production using a pedigree-based approach that is innovative in plant genetics.ResultsA quantitative trait locus (QTL) mapping approach involving two-step variance component analysis was employed using phenotypic data on 30852 palms from crosses between more than 300 genotyped parents of two heterotic groups. Genome scans were performed at parental level by modeling QTL effects as random terms in linear mixed models with identity-by-descent (IBD) kinship matrices. Eighteen QTL regions controlling production traits were identified among a large genetically diversified sample from breeding program. QTL patterns depended on the genetic origin, with only one region shared between heterotic groups. Contrasting effects of QTLs on bunch number and weights reflected the close negative correlation between the two traits.ConclusionsThe pedigree-based approach using data from ongoing breeding programs is a powerful, relevant and economic approach to map QTLs. Genetic determinisms contributing to heterotic effects have been identified and provide valuable information for orienting oil palm breeding strategies.

Centrosomal targeting of Syk kinase is controlled by its catalytic activity and depends on microtubules and the dynein motor

The nonreceptor Syk kinase is detected in epithelial cells, where it acts as a tumor suppressor, in addition to its well‐established role in immunoreceptor‐based signal transduction in hematopoietic cells. Thus, several carcinomas and melanomas have subnormal concentrations of Syk. Although Syk is mainly localized at the plasma membrane, it is also present in centrosomes, where it is involved in the control of cell division. The mechanisms responsible for its centrosomal localization and action are unknown. We used wild‐type and mutant fluorescent Syk fusion proteins in live‐cell imaging (fluorescence recovery after photobleaching, total internal reflection fluorescence, and photoactivation) combined with mathematical modeling to demonstrate that Syk is actively transported to the centrosomes via the microtubules and that this transport depends on the dynein/dynactin molecular motor. Syk can only target the centrosomes if its kinase activity is intact and it is catalytically active at the centrosomes. We showed that the autophosphorylated Y130 Syk residue helps to uncouple Syk from the plasma membrane and to promote its translocation to the centrosome, suggesting that the subcellular location of Syk depends on its autophosphorylation on specific tyrosine residues. We have thus established the details of how Syk is trafficked intracellularly and found evidence that its targeting to the centrosomes is controlled by autophosphorylation.—Fargier, G., Favard, C., Parmeggiani, A., Sahuquet, A., Mérezègue, F., Morel, A., Denis, M., Molinari, N., Mangeat, P. H., Coopman, P. J., Montcourrier, P. Centrosomal targeting of Syk kinase is controlled by its catalytic activity and depends on microtubules and the dynein motor. FASEB J. 27, 109–122 (2013). www.fasebj.org