New resources for genetic studies in Populus nigra: genome-wide SNP discovery and development of a 12k Infinium array.

Whole genome resequencing of 51 Populus nigra (L.) individuals from across Western Europe was performed using Illumina platforms. A total number of 1 878 727 SNPs distributed along the P. nigra reference sequence were identified. The SNP calling accuracy was validated with Sanger sequencing. SNPs were selected within 14 previously identified QTL regions, 2916 expressional candidate genes related to rust resistance, wood properties, water-use efficiency and bud phenology and 1732 genes randomly spread across the genome. Over 10 000 SNPs were selected for the construction of a 12k Infinium Bead-Chip array dedicated to association mapping. The SNP genotyping assay was performed with 888 P. nigra individuals. The genotyping success rate was 91%. Our high success rate was due to the discovery panel design and the stringent parameters applied for SNP calling and selection. In the same set of P. nigra genotypes, linkage disequilibrium throughout the genome decayed on average within 5-7 kb to half of its maximum value. As an application test, ADMIXTURE analysis was performed with a selection of 600 SNPs spread throughout the genome and 706 individuals collected along 12 river basins. The admixture pattern was consistent with genetic diversity revealed by neutral markers and the geographical distribution of the populations. These newly developed SNP resources and genotyping array provide a valuable tool for population genetic studies and identification of QTLs through natural-population based genetic association studies in P. nigra.

First insights into the transcriptome and development of new genomic tools of a widespread circum-Mediterranean tree species, Pinus halepensis Mill.

Aleppo pine (Pinus halepensis Mill.) is a relevant conifer species for studying adaptive responses to drought and fire regimes in the Mediterranean region. In this study, we performed Illumina next-generation sequencing of two phenotypically divergent Aleppo pine accessions with the aims of (i) characterizing the transcriptome through Illumina RNA-Seq on trees phenotypically divergent for adaptive traits linked to fire adaptation and drought, (ii) performing a functional annotation of the assembled transcriptome, (iii) identifying genes with accelerated evolutionary rates, (iv) studying the expression levels of the annotated genes and (v) developing gene-based markers for population genomic and association genetic studies. The assembled transcriptome consisted of 48,629 contigs and covered about 54.6 Mbp. The comparison of Aleppo pine transcripts to Picea sitchensis protein-coding sequences resulted in the detection of 34,014 SNPs across species, with a Ka /Ks average value of 0.216, suggesting that the majority of the assembled genes are under negative selection. Several genes were differentially expressed across the two pine accessions with contrasted phenotypes, including a glutathione-s-transferase, a cellulose synthase and a cobra-like protein. A large number of new markers (3334 amplifiable SSRs and 28,236 SNPs) have been identified which should facilitate future population genomics and association genetics in this species. A 384-SNP Oligo Pool Assay for genotyping with the Illumina VeraCode technology has been designed which showed an high overall SNP conversion rate (76.6%). Our results showed that Illumina next-generation sequencing is a valuable technology to obtain an extensive overview on whole transcriptomes of nonmodel species with large genomes.

Temporal dynamics in the evolution of the sunflower genome as revealed by sequencing and annotation of three large genomic regions.

Improved knowledge of genome composition, especially of its repetitive component, generates important informations in both theoretical and applied research. In this study, we provide the first insight into the local organization of the sunflower genome by sequencing and annotating 349,380 bp from 3 BAC clones, each including one single-copy gene. These analyses resulted in the identification of 11 putative gene sequences, 18 full-length LTR retrotransposons, 6 incomplete LTR retrotransposons, 2 non-autonomous LTR-retroelements (LINEs), 2 putative DNA transposons fragments and one putative helitron. Among LTR-retrotransposons, non-autonomous elements (the so-called LARDs), which do not carry any protein-encoding sequence, were discovered for the first time in the sunflower. The insertion time of intact retroelements was measured, based on sister LTRs divergence. All isolated elements were inserted relatively recently, especially those belonging to the Gypsy superfamily. Retrotransposon families related to those identified in the BAC clones are present also in other species of Helianthus, both annual and perennial, and even in other Asteraceae. In one of the three BAC clones, we found five copies of a lipid transfer protein (LTP) encoding gene within less than 100,000 bp, four of which are potentially functional. Two of these are interrupted by LTR retrotransposons, in the intron and in the coding sequence, respectively. The divergence between sister LTRs of the retrotransposons inserted within the genes indicates that LTP gene duplication started earlier than 1.749 MYRS ago. On the whole, the results reported in this study confirm that the sunflower is an excellent system to study transposons dynamics and evolution.

An analysis of sequence variability in eight genes putatively involved in drought response in sunflower (Helianthus annuus L.).

With the aim to study variability in genes involved in ecological adaptations, we have analysed sequence polymorphisms of eight unique genes putatively involved in drought response by isolation and analysis of allelic sequences in eight inbred lines of sunflower of different origin and phenotypic characters and showing different drought response in terms of leaf relative water content (RWC). First, gene sequences were amplified by PCR on genomic DNA from a highly inbred line and their products were directly sequenced. In the absence of single nucleotide polymorphisms, the gene was considered as unique. Then, the same PCR reaction was performed on genomic DNAs of eight inbred lines to isolate allelic variants to be compared. The eight selected genes encode a dehydrin, a heat shock protein, a non-specific lipid transfer protein, a z-carotene desaturase, a drought-responsive-element-binding protein, a NAC-domain transcription regulator, an auxin-binding protein, and an ABA responsive-C5 protein. Nucleotide diversity per synonymous and non-synonymous sites was calculated for each gene sequence. The π (a)/π (s) ratio range was usually very low, indicating strong purifying selection, though with locus-to-locus differences. As far as non-coding regions, the intron showed a larger variability than the other regions only in the case of the dehydrin gene. In the other genes tested, in which one or more introns occur, variability in the introns was similar or even lower than in the other regions. On the contrary, 3'-UTRs were usually more variable than the coding regions. Linkage disequilibrium in the selected genes decayed on average within 1,000 bp, with large variation among genes. A pairwise comparison between genetic distances calculated on the eight genes and the difference in RWC showed a significant correlation in the first phases of drought stress. The results are discussed in relation to the function of analysed genes, i.e. involved in gene regulation and signal transduction, or encoding enzymes and defence proteins.

HACRE1, a recently inserted copia-like retrotransposon of sunflower (Helianthus annuus L.).

In this paper we report on the isolation and characterization, for the first time, of a complete 6511 bp retrotransposon of sunflower. Considering its protein domain order and sequence similarity to other copia elements of dicotyledons, this retrotransposon was assigned to the copia retrotransposon superfamily and named HACRE1 (Helianthus annuus copia-like retroelement 1). HACRE1 carries 5' and 3' long terminal repeats (LTRs) flanking an internal region of 4661 bp. The LTRs are identical in their sequence except for two deletions of 7 and 5 nucleotides in the 5' LTR. Based on the sequence identity of the LTRs, HACRE1 was estimated to have inserted within the last approximately 84 000 years. The isolated sequence contains a complete open reading frame with only one complete reading frame. The absence of nonsense mutations agrees with the very high sequence identity between LTRs, confirming that HACRE1 insertion is recent. The haploid genome of sunflower (inbred line HCM) contains about 160 copies of HACRE1. This retrotransposon is expressed in leaflets from 7-day-old plantlets under different light conditions, probably in relation to the occurrence of many putative light-related regulatory cis-elements in the LTRs. However, sequenced cDNAs show less variability than HACRE1 genomic sequences, indicating that only a subset of this family is expressed under these conditions.