Vincent Le Guen

De novo transcriptome analysis of Hevea brasiliensis tissues by RNA-seq and screening for molecular markers

BackgroundThe rubber tree, Hevea brasiliensis, is a species native to the Brazilian Amazon region and it supplies almost all the world’s natural rubber, a strategic raw material for a variety of products. One of the major challenges for developing rubber tree plantations is adapting the plant to biotic and abiotic stress. Transcriptome analysis is one of the main approaches for identifying the complete set of active genes in a cell or tissue for a specific developmental stage or physiological condition.ResultsHere, we report on the sequencing, assembling, annotation and screening for molecular markers from a pool of H. brasiliensis tissues. A total of 17,166 contigs were successfully annotated. Then, 2,191 Single Nucleotide Variation (SNV) and 1.397 Simple Sequence Repeat (SSR) loci were discriminated from the sequences. From 306 putative, mainly non-synonymous SNVs located in CDS sequences, 191 were checked for their ability to characterize 23 Hevea genotypes by an allele-specific amplification technology. For 172 (90%), the nucleotide variation at the predicted genomic location was confirmed, thus validating the different steps from sequencing to the in silico detection of the SNVs.ConclusionsThis is the first study of the H. brasiliensis transcriptome, covering a wide range of tissues and organs, leading to the production of the first developed SNP markers. This process could be amplified to a larger set of in silico detected SNVs in expressed genes in order to increase the marker density in available and future genetic maps. The results obtained in this study will contribute to the H. brasiliensis genetic breeding program focused on improving of disease resistance and latex yield.

QTL Mapping of Growth-Related Traits in a Full-Sib Family of Rubber Tree ( Hevea brasiliensis ) Evaluated in a Sub-Tropical Climate

The rubber tree (Hevea spp.), cultivated in equatorial and tropical countries, is the primary plant used in natural rubber production. Due to genetic and physiological constraints, inbred lines of this species are not available. Therefore, alternative approaches are required for the characterization of this species, such as the genetic mapping of full-sib crosses derived from outbred parents. In the present study, an integrated genetic map was obtained for a full-sib cross family with simple sequence repeats (SSRs) and expressed sequence tag (EST-SSR) markers, which can display different segregation patterns. To study the genetic architecture of the traits related to growth in two different conditions (winter and summer), quantitative trait loci (QTL) mapping was also performed using the integrated map. Traits evaluated were height and girth growth, and the statistical model was based in an extension of composite interval mapping. The obtained molecular genetic map has 284 markers distributed among 23 linkage groups with a total length of 2688.8 cM. A total of 18 QTLs for growth traits during the summer and winter seasons were detected. A comparison between the different seasons was also conducted. For height, QTLs detected during the summer season were different from the ones detected during winter season. This type of difference was also observed for girth. Integrated maps are important for genetics studies in outbred species because they represent more accurately the polymorphisms observed in the genitors. QTL mapping revealed several interesting findings, such as a dominance effect and unique segregation patterns that each QTL could exhibit, which were independent of the flanking markers. The QTLs identified in this study, especially those related to phenotypic variation associated with winter could help studies of marker-assisted selection that are particularly important when the objective of a breeding program is to obtain phenotypes that are adapted to sub-optimal regions.

A rubber tree’s durable resistance to Microcyclus ulei is conferred by a qualitative gene and a major quantitative resistance factor

The components of genetic resistance from the Hevea brasiliensis cultivar MDF 180 against South American Leaf Blight (SALB) caused by Microcyclus ulei were investigated by QTL mapping. MDF 180 has already been described as a cultivar with a high level of partial and long-lasting resistance. The resistance of progeny individuals from a cross between a susceptible cultivar and MDF 180 was assessed both under controlled conditions of inoculation by three M. ulei isolates and under natural infection in a field trial. Genetic maps of the two parents of this progeny were mainly established based on microsatellites and AFLP markers. No resistance QTL were found in the susceptible parent. In the resistant parent, we identified a qualitative gene responsible for the resistance against isolates from French Guiana and a major quantitative resistance factor determining the resistance against isolates from the state of Bahia (Brazil). The qualitative resistance gene was denominated M15md and was located in the linkage group g15. Four minor resistance QTLs were also identified, two of which showed an epistatic interaction with M15md. The durability of the resistance of MDF 180 is discussed in light of these data.

Genetic Determinism of Sensitivity to Corynespora cassiicola Exudates in Rubber Tree (Hevea brasiliensis)

An indirect phenotyping method was developed in order to estimate the susceptibility of rubber tree clonal varieties to Corynespora Leaf Fall (CLF) disease caused by the ascomycete Corynespora cassiicola. This method consists in quantifying the impact of fungal exudates on detached leaves by measuring the induced electrolyte leakage (EL%). The tested exudates were either crude culture filtrates from diverse C. cassiicola isolates or the purified cassiicolin (Cas1), a small secreted effector protein produced by the aggressive isolate CCP. The test was found to be quantitative, with the EL% response proportional to toxin concentration. For eight clones tested with two aggressive isolates, the EL% response to the filtrates positively correlated to the response induced by conidial inoculation. The toxicity test applied to 18 clones using 13 toxinic treatments evidenced an important variability among clones and treatments, with a significant additional clone x treatment interaction effect. A genetic linkage map was built using 306 microsatellite markers, from the F1 population of the PB260 x RRIM600 family. Phenotyping of the population for sensitivity to the purified Cas1 effector and to culture filtrates from seven C. cassiicola isolates revealed a polygenic determinism, with six QTL detected on five chromosomes and percentages of explained phenotypic variance varying from 11 to 17%. Two common QTL were identified for the CCP filtrate and the purified cassiicolin, suggesting that Cas1 may be the main effector of CCP filtrate toxicity. The CCP filtrate clearly contrasted with all other filtrates. The toxicity test based on Electrolyte Leakage Measurement offers the opportunity to assess the sensitivity of rubber genotypes to C. cassiicola exudates or purified effectors for genetic investigations and early selection, without risk of spreading the fungus in plantations. However, the power of this test for predicting field susceptibility of rubber clones to CLF will have to be further investigated.

Genetic Diversity Strategy for the Management and Use of Rubber Genetic Resources: More than 1,000 Wild and Cultivated Accessions in a 100-Genotype Core Collection.

The rubber tree [Hevea brasiliensis (Willd. ex Adr. de Juss.) Muell. Arg.] is the only plant species worldwide that is cultivated for the commercial production of natural rubber. This study describes the genetic diversity of the Hevea spp. complex that is available in the main ex situ collections of South America, including Amazonian populations that have never been previously described. Genetic data were analyzed to determine the genetic structure of the wild populations, quantify the allelic diversity and suggest the composition of a core collection to capture the maximum genetic diversity within a minimal sample size. A total of 1,117 accessions were genotyped with 13 microsatellite markers. We identified a total of 408 alleles, 319 of which were shared between groups and 89 that were private in different groups of accessions. In a population structure and principal component analysis, the level of clustering reflected a primary division into the following two subgroups: cluster 1, which consisted of varieties from the advanced breeding germplasm that originated from the Wickham and Mato Grosso accessions; and cluster 2, which consisted of the wild germplasm from the Acre, Amazonas, Pará and Rondônia populations and Hevea spp. The analyses revealed a high frequency of gene flow between the groups, with the genetic differentiation coefficient (GST) estimated to be 0.018. Additionally, no distinct separation among the H. brasiliensis accessions and the other species from Amazonas was observed. A core collection of 99 accessions was identified that captured the maximum genetic diversity. Rubber tree breeders can effectively utilize this core collection for cultivar improvement. Furthermore, such a core collection could provide resources for forming an association panel to evaluate traits with agronomic and commercial importance. Our study generated a molecular database that should facilitate the management of the Hevea germplasm and its use for subsequent genetic and genomic breeding.

A newly identified locus controls complete resistance to Microcyclus ulei in the Fx2784 rubber clone

Using cultivars which are genetically resistant to South American leaf blight (SALB) caused by the fungus Microcyclus ulei is the only way to plant rubber trees in disease-affected areas. Numerous field observations led to the hypothesis that the resistance of the cultivar Fx2784 to SALB is likely to be monogenic. In this study, we investigated this hypothesis by examining the distribution of the trait in a cross between the resistant cultivar and a susceptible one. The individuals resulting from this cross were planted in field trials in French Guiana and Brazil. The resistance of all the trees was assessed by field observations. Bulk segregant analysis (BSA) using microsatellite markers was performed in French Guiana to determine which markers were genetically linked to resistance, and the results were validated by field observations in Brazil. In both locations, a 1:1 segregation of the resistance trait was observed, thus reinforcing the monogenic hypothesis. BSA showed tight linkage between resistance and the microsatellite markers located in linkage group 2 in the Hevea genome and enabled to pinpoint the resistance locus. The location was confirmed by observations on the trees planted in Brazil. This result should facilitate the use of Fx2784 resistance in future breeding programs for SALB resistance. This is the third major locus conferring resistance to SALB identified in rubber tree (Hevea spp.). These three loci are genetically independent, a favorable situation for genetic improvement of SALB resistance.

High-resolution genetic map and QTL analysis of growth-related traits of Hevea brasiliensis cultivated under suboptimal temperature and humidity conditions

Rubber tree (Hevea brasiliensis) cultivation is the main source of natural rubber worldwide and has been extended to areas with suboptimal climates and lengthy drought periods; this transition affects growth and latex production. High-density genetic maps with reliable markers support precise mapping of quantitative trait loci (QTL), which can help reveal the complex genome of the species, provide tools to enhance molecular breeding, and shorten the breeding cycle. In this study, QTL mapping of the stem diameter, tree height, and number of whorls was performed for a full-sibling population derived from a GT1 and RRIM701 cross. A total of 225 simple sequence repeats (SSRs) and 186 single-nucleotide polymorphism (SNP) markers were used to construct a base map with 18 linkage groups and to anchor 671 SNPs from genotyping by sequencing (GBS) to produce a very dense linkage map with small intervals between loci. The final map was composed of 1,079 markers, spanned 3,779.7 cM with an average marker density of 3.5 cM, and showed collinearity between markers from previous studies. Significant variation in phenotypic characteristics was found over a 59-month evaluation period with a total of 38 QTLs being identified through a composite interval mapping method. Linkage group 4 showed the greatest number of QTLs (7), with phenotypic explained values varying from 7.67 to 14.07%. Additionally, we estimated segregation patterns, dominance, and additive effects for each QTL. A total of 53 significant effects for stem diameter were observed, and these effects were mostly related to additivity in the GT1 clone. Associating accurate genome assemblies and genetic maps represents a promising strategy for identifying the genetic basis of phenotypic traits in rubber trees. Then, further research can benefit from the QTLs identified herein, providing a better understanding of the key determinant genes associated with growth of Hevea brasiliensis under limiting water conditions.

Linkage disequilibrium and population structure in wild and cultivated populations of rubber tree (Hevea brasiliensis).

Among rubber tree species, which belong to the Hevea genus of the Euphorbiaceae family, Hevea brasiliensis (Willd. ex Adr.de Juss.) Muell. Arg. is the main commercial source of natural rubber production worldwide. Knowledge of the population structure and linkage disequilibrium (LD) of this species is essential for the efficient organization and exploitation of genetic resources. Here, we obtained single-nucleotide polymorphisms (SNPs) using a genotyping-by-sequencing (GBS) approach and then employed the SNPs for the following objectives: (i) to identify the positions of SNPs on a genetic map of a segregating mapping population, (ii) to evaluate the population structure of a germplasm collection, and (iii) to detect patterns of LD decay among chromosomes for future genetic association studies in rubber tree. A total of 626 genotypes, including both germplasm accessions (368) and individuals from a genetic mapping population (254), were genotyped. A total of 77,660 and 21,283 SNPs were detected by GBS in the germplasm and mapping populations, respectively. The mapping population, which was previously mapped, was constructed with 1,062 markers, among which only 576 SNPs came from GBS, reducing the average interval between two adjacent markers to 4.4 cM. SNPs from GBS genotyping were used for the analysis of genetic structure and LD estimation in the germplasm accessions. Two groups, which largely corresponded to the cultivated and wild populations, were detected using STRUCTURE and via principal coordinate analysis. LD analysis, also using the mapped SNPs, revealed that non-random associations varied along chromosomes, with regions of high LD interspersed with regions of low LD. Considering the length of the genetic map (4,693 cM) and the mean LD (0.49 for cultivated and 0.02 for wild populations), a large number of evenly spaced SNPs would be needed to perform genome-wide association studies in rubber tree, and the wilder the genotypes used, the more difficult the mapping saturation.