Abdelali Barakat

The high-quality draft genome of peach (Prunus persica) identifies unique patterns of genetic diversity, domestication and genome evolution

Rosaceae is the most important fruit-producing clade, and its key commercially relevant genera (Fragaria, Rosa, Rubus and Prunus) show broadly diverse growth habits, fruit types and compact diploid genomes. Peach, a diploid Prunus species, is one of the best genetically characterized deciduous trees. Here we describe the high-quality genome sequence of peach obtained from a completely homozygous genotype. We obtained a complete chromosome-scale assembly using Sanger whole-genome shotgun methods. We predicted 27,852 protein-coding genes, as well as noncoding RNAs. We investigated the path of peach domestication through whole-genome resequencing of 14 Prunus accessions. The analyses suggest major genetic bottlenecks that have substantially shaped peach genome diversity. Furthermore, comparative analyses showed that peach has not undergone recent whole-genome duplication, and even though the ancestral triplicated blocks in peach are fragmentary compared to those in grape, all seven paleosets of paralogs from the putative paleoancestor are detectable.

FLOWERING LOCUS T duplication coordinates reproductive and vegetative growth in perennial poplar

Annual plants grow vegetatively at early developmental stages and then transition to the reproductive stage, followed by senescence in the same year. In contrast, after successive years of vegetative growth at early ages, woody perennial shoot meristems begin repeated transitions between vegetative and reproductive growth at sexual maturity. However, it is unknown how these repeated transitions occur without a developmental conflict between vegetative and reproductive growth. We report that functionally diverged paralogs FLOWERING LOCUS T1 (FT1) and FLOWERING LOCUS T2 (FT2), products of whole-genome duplication and homologs of Arabidopsis thaliana gene FLOWERING LOCUS T (FT), coordinate the repeated cycles of vegetative and reproductive growth in woody perennial poplar (Populus spp.). Our manipulative physiological and genetic experiments coupled with field studies, expression profiling, and network analysis reveal that reproductive onset is determined by FT1 in response to winter temperatures, whereas vegetative growth and inhibition of bud set are promoted by FT2 in response to warm temperatures and long days in the growing season. The basis for functional differentiation between FT1 and FT2 appears to be expression pattern shifts, changes in proteins, and divergence in gene regulatory networks. Thus, temporal separation of reproductive onset and vegetative growth into different seasons via FT1 and FT2 provides seasonality and demonstrates the evolution of a complex perennial adaptive trait after genome duplication.

Conservation and divergence of microRNAs in Populus

BackgroundMicroRNAs (miRNAs) are small RNAs (sRNA) ~21 nucleotides in length that negatively control gene expression by cleaving or inhibiting the translation of target gene transcripts. miRNAs have been extensively analyzed in Arabidopsis and rice and partially investigated in other non-model plant species. To date, 109 and 62 miRNA families have been identified in Arabidopsis and rice respectively. However, only 33 miRNAs have been identified from the genome of the model tree species (Populus trichocarpa), of which 11 are Populus specific. The low number of miRNA families previously identified in Populus, compared with the number of families identified in Arabidopsis and rice, suggests that many miRNAs still remain to be discovered in Populus. In this study, we analyzed expressed small RNAs from leaves and vegetative buds of Populus using high throughput pyrosequencing.ResultsAnalysis of almost eighty thousand small RNA reads allowed us to identify 123 new sequences belonging to previously identified miRNA families as well as 48 new miRNA families that could be Populus-specific. Comparison of the organization of miRNA families in Populus, Arabidopsis and rice showed that miRNA family sizes were generally expanded in Populus. The putative targets of non-conserved miRNA include both previously identified targets as well as several new putative target genes involved in development, resistance to stress, and other cellular processes. Moreover, almost half of the genes predicted to be targeted by non-conserved miRNAs appear to be Populus-specific. Comparative analyses showed that genes targeted by conserved and non-conserved miRNAs are biased mainly towards development, electron transport and signal transduction processes. Similar results were found for non-conserved miRNAs from Arabidopsis.ConclusionOur results suggest that while there is a conserved set of miRNAs among plant species, a large fraction of miRNAs vary among species. The non-conserved miRNAs may regulate cellular, physiological or developmental processes specific to the taxa that produce them, as appears likely to be the case for those miRNAs that have only been observed in Populus. Non-conserved and conserved miRNAs seem to target genes with similar biological functions indicating that similar selection pressures are acting on both types of miRNAs. The expansion in the number of most conserved miRNAs in Populus relative to Arabidopsis, may be linked to the recent genome duplication in Populus, the slow evolution of the Populus genome, or to differences in the selection pressure on duplicated miRNAs in these species.

The cinnamyl alcohol dehydrogenase gene family in Populus: phylogeny, organization, and expression

BackgroundLignin is a phenolic heteropolymer in secondary cell walls that plays a major role in the development of plants and their defense against pathogens. The biosynthesis of monolignols, which represent the main component of lignin involves many enzymes. The cinnamyl alcohol dehydrogenase (CAD) is a key enzyme in lignin biosynthesis as it catalyzes the final step in the synthesis of monolignols. The CAD gene family has been studied in Arabidopsis thaliana, Oryza sativa and partially in Populus. This is the first comprehensive study on the CAD gene family in woody plants including genome organization, gene structure, phylogeny across land plant lineages, and expression profiling in Populus.ResultsThe phylogenetic analyses showed that CAD genes fall into three main classes (clades), one of which is represented by CAD sequences from gymnosperms and angiosperms. The other two clades are represented by sequences only from angiosperms. All Populus CAD genes, except PoptrCAD 4 are distributed in Class II and Class III. CAD genes associated with xylem development (PoptrCAD 4 and PoptrCAD 10) belong to Class I and Class II. Most of the CAD genes are physically distributed on duplicated blocks and are still in conserved locations on the homeologous duplicated blocks. Promoter analysis of CAD genes revealed several motifs involved in gene expression modulation under various biological and physiological processes. The CAD genes showed different expression patterns in poplar with only two genes preferentially expressed in xylem tissues during lignin biosynthesis.ConclusionThe phylogeny of CAD genes suggests that the radiation of this gene family may have occurred in the early ancestry of angiosperms. Gene distribution on the chromosomes of Populus showed that both large scale and tandem duplications contributed significantly to the CAD gene family expansion. The duplication of several CAD genes seems to be associated with a genome duplication event that happened in the ancestor of Salicaceae. Phylogenetic analyses associated with expression profiling and results from previous studies suggest that CAD genes involved in wood development belong to Class I and Class II. The other CAD genes from Class II and Class III may function in plant tissues under biotic stresses. The conservation of most duplicated CAD genes, the differential distribution of motifs in their promoter regions, and the divergence of their expression profiles in various tissues of Populus plants indicate that genes in the CAD family have evolved tissue-specialized expression profiles and may have divergent functions.

The distribution of T-DNA in the genomes of transgenic Arabidopsis and rice

Almost all the nuclear genes of four Gramineae (maize, wheat, barley, rice) and pea are located in DNA fractions covering only a 1–2% GC range and representing between 10 and 25% of the different genomes. These DNA fractions comprise large gene‐rich regions (collectively called the ‘gene space’) separated by vast gene‐empty, repeated sequences. In contrast, in Arabidopsis thaliana, genes are distributed in DNA fractions covering an 8% GC range and representing 85% of the genome. Here, we investigated the integration of a transferred DNA (T‐DNA) in the genomes of Arabidopsis and rice and found different patterns of integration, which are correlated with the different gene distributions. While T‐DNA integrates essentially everywhere in the Arabidopsis genome, integration was detected only in the gene space, namely in the gene‐rich, transcriptionally active, regions of the rice genome. The implications of these results for the integration of foreign DNA are discussed.

Genome wide identification of chilling responsive microRNAs in Prunus persica.

BackgroundMicroRNAs (miRNAs) are small RNAs (sRNAs) approximately 21 nucleotides in length that negatively control gene expression by cleaving or inhibiting the translation of target gene transcripts. Within this context, miRNAs and siRNAs are coming to the forefront as molecular mediators of gene regulation in plant responses to annual temperature cycling and cold stress. For this reason, we chose to identify and characterize the conserved and non-conserved miRNA component of peach (Prunus persica (L.) Batsch) focusing our efforts on both the recently released whole genome sequence of peach and sRNA transcriptome sequences from two tissues representing non-dormant leaves and dormant leaf buds. Conserved and non-conserved miRNAs, and their targets were identified. These sRNA resources were used to identify cold-responsive miRNAs whose gene targets co-localize with previously described QTLs for chilling requirement (CR).ResultsAnalysis of 21 million peach sRNA reads allowed us to identify 157 and 230 conserved and non-conserved miRNA sequences. Among the non-conserved miRNAs, we identified 205 that seem to be specific to peach. Comparative genome analysis between peach and Arabidopsis showed that conserved miRNA families, with the exception of miR5021, are similar in size. Sixteen of these conserved miRNA families are deeply rooted in land plant phylogeny as they are present in mosses and/or lycophytes. Within the other conserved miRNA families, five families (miR1446, miR473, miR479, miR3629, and miR3627) were reported only in tree species (Populustrichocarpa, Citrus trifolia, and Prunus persica). Expression analysis identified several up-regulated or down-regulated miRNAs in winter buds versus young leaves. A search of the peach proteome allowed the prediction of target genes for most of the conserved miRNAs and a large fraction of non-conserved miRNAs. A fraction of predicted targets in peach have not been previously reported in other species. Several conserved and non-conserved miRNAs and miRNA-regulated genes co-localize with Quantitative Trait Loci (QTLs) for chilling requirement (CR-QTL) and bloom date (BD-QTL).ConclusionsIn this work, we identified a large set of conserved and non-conserved miRNAs and describe their evolutionary footprint in angiosperm lineages. Several of these miRNAs were induced in winter buds and co-localized with QTLs for chilling requirement and bloom date thus making their gene targets potential candidates for mediating plant responses to cold stress. Several peach homologs of genes participating in the regulation of vernalization in Arabidopsis were identified as differentially expressed miRNAs targets, potentially linking these gene activities to cold responses in peach dormant buds. The non-conserved miRNAs may regulate cellular, physiological or developmental processes specific to peach and/or other tree species.

Overexpression of Constans Homologs CO1 and CO2 Fails to Alter Normal Reproductive Onset and Fall Bud Set in Woody Perennial Poplar

CONSTANS (CO) is an important flowering-time gene in the photoperiodic flowering pathway of annual Arabidopsis thaliana in which overexpression of CO induces early flowering, whereas mutations in CO cause delayed flowering. The closest homologs of CO in woody perennial poplar (Populus spp.) are CO1 and CO2. A previous report [1] showed that the CO2/FLOWERING LOCUS T1 (FT1) regulon controls the onset of reproduction in poplar, similar to what is seen with the CO/FLOWERING LOCUS T (FT) regulon in Arabidopsis. The CO2/FT1 regulon was also reported to control fall bud set. Our long-term field observations show that overexpression of CO1 and CO2 individually or together did not alter normal reproductive onset, spring bud break, or fall dormancy in poplar, but did result in smaller trees when compared with controls. Transcripts of CO1 and CO2 were normally most abundant in the growing season and rhythmic within a day, peaking at dawn. Our manipulative experiments did not provide evidence for transcriptional regulation being affected by photoperiod, light intensity, temperature, or water stress when transcripts of CO1 and CO2 were consistently measured in the morning. A genetic network analysis using overexpressing trees, microarrays, and computation demonstrated that a majority of functionally known genes downstream of CO1 and CO2 are associated with metabolic processes, which could explain their effect on tree size. In conclusion, the function of CO1 and CO2 in poplar does not appear to overlap with that of CO from Arabidopsis, nor do our data support the involvement of CO1 and CO2 in spring bud break or fall bud set.

Chestnut resistance to the blight disease: insights from transcriptome analysis

BackgroundA century ago, Chestnut Blight Disease (CBD) devastated the American chestnut. Backcross breeding has been underway to introgress resistance from Chinese chestnut into surviving American chestnut genotypes. Development of genomic resources for the family Fagaceae, has focused in this project on Castanea mollissima Blume (Chinese chestnut) and Castaneadentata (Marsh.) Borkh (American chestnut) to aid in the backcross breeding effort and in the eventual identification of blight resistance genes through genomic sequencing and map based cloning. A previous study reported partial characterization of the transcriptomes from these two species. Here, further analyses of a larger dataset and assemblies including both 454 and capillary sequences were performed and defense related genes with differential transcript abundance (GDTA) in canker versus healthy stem tissues were identified.ResultsOver one and a half million cDNA reads were assembled into 34,800 transcript contigs from American chestnut and 48,335 transcript contigs from Chinese chestnut. Chestnut cDNA showed higher coding sequence similarity to genes in other woody plants than in herbaceous species. The number of genes tagged, the length of coding sequences, and the numbers of tagged members within gene families showed that the cDNA dataset provides a good resource for studying the American and Chinese chestnut transcriptomes. In silico analysis of transcript abundance identified hundreds of GDTA in canker versus healthy stem tissues. A significant number of additional DTA genes involved in the defense-response not reported in a previous study were identified here. These DTA genes belong to various pathways involving cell wall biosynthesis, reactive oxygen species (ROS), salicylic acid (SA), ethylene, jasmonic acid (JA), abscissic acid (ABA), and hormone signalling. DTA genes were also identified in the hypersensitive response and programmed cell death (PCD) pathways. These DTA genes are candidates for host resistance to the chestnut blight fungus, Cryphonectria parasitica.ConclusionsOur data allowed the identification of many genes and gene network candidates for host resistance to the chestnut blight fungus, Cryphonectria parasitica. The similar set of GDTAs in American chestnut and Chinese chestnut suggests that the variation in sensitivity to this pathogen between these species may be the result of different timing and amplitude of the response of the two to the pathogen infection. Resources developed in this study are useful for functional genomics, comparative genomics, resistance breeding and phylogenetics in the Fagaceae.

Phylogeny and expression profiling of CAD and CAD-like genes in hybrid Populus ( P. deltoides × P. nigra ): evidence from herbivore damage for subfunctionalization and functional divergence

BackgroundCinnamyl Alcohol Dehydrogenase (CAD) proteins function in lignin biosynthesis and play a critical role in wood development and plant defense against stresses. Previous phylogenetic studies did not include genes from seedless plants and did not reflect the deep evolutionary history of this gene family. We reanalyzed the phylogeny of CAD and CAD-like genes using a representative dataset including lycophyte and bryophyte sequences. Many CAD/CAD-like genes do not seem to be associated with wood development under normal growth conditions. To gain insight into the functional evolution of CAD/CAD-like genes, we analyzed their expression in Populus plant tissues in response to feeding damage by gypsy moth larvae (Lymantria dispar L.). Expression of CAD/CAD-like genes in Populus tissues (xylem, leaves, and barks) was analyzed in herbivore-treated and non-treated plants by real time quantitative RT-PCR.ResultsCAD family genes were distributed in three classes based on sequence conservation. All the three classes are represented by seedless as well as seed plants, including the class of bona fide lignin pathway genes. The expression of some CAD/CAD-like genes that are not associated with xylem development were induced following herbivore damage in leaves, while other genes were induced in only bark or xylem tissues. Five of the CAD/CAD-like genes, however, showed a shift in expression from one tissue to another between non-treated and herbivore-treated plants. Systemic expression of the CAD/CAD-like genes was generally suppressed.ConclusionsOur results indicated a correlation between the evolution of the CAD gene family and lignin and that the three classes of genes may have evolved in the ancestor of land plants. Our results also suggest that the CAD/CAD-like genes have evolved a diversity of expression profiles and potentially different functions, but that they are nonetheless co-regulated under stress conditions.

Comparative and phylogenomic analyses of cinnamoyl-CoA reductase and cinnamoyl-CoA-reductase-like gene family in land plants

The biosynthesis of monolignols, the main components of lignin, involves many intermediates and enzymes. The cinnamoyl-CoA reductase (CCR) enzyme catalyzes the conversion of cinnamoyl-CoAs to cinnamaldehydes, i.e. the first specific step in lignin synthesis. The CCR and CCR-like gene family was studied partially in several plant species. This is a comprehensive study of the CCR and CCR-like gene family including genome organization, gene structure, phylogeny across land plant species, and, expression profiling in Populus. Analysis of amino acid motifs enabled the identification of sequence variations in the CCR catalytic site and annotates CCR and CCR-like genes. CCR and CCR-like genes were distributed in three major phylogenetic classes of which one includes the bona fide CCR genes. The other two classes include CCR and CCR-like, of which several genes present a high similarity to cinnamyl alcohol dehydrogenase, or dihydroflavonol reductase (DFR) genes. All CCR, CCR-like, and DFR classes were deeply rooted in the phylogeny of land plants suggesting that their evolution preceded the evolution of lycophytes. Over two thirds of CCR and CCR-like Populus genes were physically distributed on duplicated regions. This suggests that these duplication/retention processes contributed significantly to the size of the CCR and CCR-like gene family. The Populus CCR and CCR-like genes showed six expression patterns in the tissues studied with a preferential expression of PoptrCCR12 in xylem. The other genes present divergent expression profiles with some preferentially expressed in leaves, bark, or both. Several CCR and CCR-like genes were induced or repressed under various abiotic stresses suggesting that their duplication was followed by the evolution of divergent expression profiles and divergence of functions.