Regulatory mechanisms behind the phenotypic plasticity associated with Setaria italica water deficit tolerance.

Drought is one of the main environmental stresses that negatively impacts vegetative and reproductive yield. Water deficit responses are determined by the duration and intensity of the stress, which, together with plant genotype, will define the chances of plant survival. The metabolic adjustments in response to water deficit are complex and involve gene expression modulation regulated by DNA-binding proteins and epigenetic modifications. This last mechanism may also regulate the activity of transposable elements, which in turn impact the expression of nearby loci. Setaria italica plants submitted to five water deficit regimes were analyzed through a phenotypical approach, including growth, physiological, RNA-seq and sRNA-seq analyses. The results showed a progressive reduction in yield as a function of water deficit intensity associated with signaling pathway modulation and metabolic adjustments. We identified a group of loci that were consistently associated with drought responses, some of which were related to water deficit perception, signaling and regulation. Finally, an analysis of the transcriptome and sRNAome allowed us to identify genes putatively regulated by TE- and sRNA-related mechanisms and an intriguing positive correlation between transcript levels and sRNA accumulation in gene body regions. These findings shed light on the processes that allow S. italica to overcome drought and survive under water restrictive conditions.

Genome-wide characterization of the laccase gene family in Setaria viridis reveals members potentially involved in lignification.

MAIN CONCLUSION: Five laccase genes are potentially involved in developmental lignification in the model C4 grass Setaria viridis and their different tissue specificities suggest subfunctionalization events. Plant laccases are copper-containing glycoproteins involved in monolignol oxidation and, therefore, their activity is essential for lignin polymerization. Although these enzymes belong to large multigene families with highly redundant members, not all of them are thought to be involved in lignin metabolism. Here, we report on the genome-wide characterization of the laccase gene family in the model C4 grass Setaria viridis and further identification of the members potentially involved in monolignol oxidation. A total of 52 genes encoding laccases (SvLAC1 to SvLAC52) were found in the genome of S. viridis, and phylogenetic analyses showed that these genes were heterogeneously distributed among the characteristic six subclades of the family and are under relaxed selective constraints. The observed expansion in the total number of genes in this species was mainly caused by tandem duplications within subclade V, which accounts for 68% of the whole family. Comparative phylogenetic analyses showed that the expansion of subclade V is specifically observed for the Paniceae tribe within the Panicoideae subfamily in grasses. Five SvLAC genes (SvLAC9, SvLAC13, SvLAC15, SvLAC50, and SvLAC52) fulfilled the criteria established to identify lignin-related candidates: (1) phylogenetic proximity to previously characterized lignin-related laccases from other species, (2) similar expression pattern to that observed for lignin biosynthetic genes in the S. viridis elongating internode, and (3) high expression in S. viridis tissues undergoing active lignification. In addition, in situ hybridization experiments not only confirmed that these selected SvLAC genes were expressed in lignifying cells, but also that their expression showed different tissue specificities, suggesting subfunctionalization events within the family. These five laccase genes are strong candidates to be involved in lignin polymerization in S. viridis and might be good targets for lignin bioengineering strategies.

The control of endopolygalacturonase expression by the sugarcane RAV transcription factor during aerenchyma formation.

The development of lysigenous aerenchyma starts with cell expansion and degradation of pectin from the middle lamella, leading to cell wall modification, and culminating with cell separation. Here we report that nutritional starvation of sugarcane induced gene expression along sections of the first 5 cm of the root and between treatments. We selected two candidate genes: a RAV transcription factor, from the ethylene response factors superfamily, and an endopolygalacturonase (EPG), a glycosyl hydrolase related to homogalacturonan hydrolysis from the middle lamella. epg1 and rav1 transcriptional patterns suggest they are essential genes at the initial steps of pectin degradation during aerenchyma development in sugarcane. Due to the high complexity of the sugarcane genome, rav1 and epg1 were sequenced from 17 bacterial artificial chromosome clones containing hom(e)ologous genomic regions, and the sequences were compared with those of Sorghum bicolor. We used one hom(e)olog sequence from each gene for transactivation assays in tobacco. rav1 was shown to bind to the epg1 promoter, repressing ß-glucuronidase activity. RAV repression upon epg1 transcription is the first reported link between ethylene regulation and pectin hydrolysis during aerenchyma formation. Our findings may help to elucidate cell wall degradation in sugarcane and therefore contribute to second-generation bioethanol production.

The genetic architecture of photosynthesis and plant growth-related traits in tomato.

To identify genomic regions involved in the regulation of fundamental physiological processes such as photosynthesis and respiration, a population of Solanum pennellii introgression lines was analyzed. We determined phenotypes for physiological, metabolic, and growth related traits, including gas exchange and chlorophyll fluorescence parameters. Data analysis allowed the identification of 208 physiological and metabolic quantitative trait loci with 33 of these being associated to smaller intervals of the genomic regions, termed BINs. Eight BINs were identified that were associated with higher assimilation rates than the recurrent parent M82. Two and 10 genomic regions were related to shoot and root dry matter accumulation, respectively. Nine genomic regions were associated with starch levels, whereas 12 BINs were associated with the levels of other metabolites. Additionally, a comprehensive and detailed annotation of the genomic regions spanning these quantitative trait loci allowed us to identify 87 candidate genes that putatively control the investigated traits. We confirmed 8 of these at the level of variance in gene expression. Taken together, our results allowed the identification of candidate genes that most likely regulate photosynthesis, primary metabolism, and plant growth and as such provide new avenues for crop improvement.

Building the sugarcane genome for biotechnology and identifying evolutionary trends.

BACKGROUND: Sugarcane is the source of sugar in all tropical and subtropical countries and is becoming increasingly important for bio-based fuels. However, its large (10 Gb), polyploid, complex genome has hindered genome based breeding efforts. Here we release the largest and most diverse set of sugarcane genome sequences to date, as part of an on-going initiative to provide a sugarcane genomic information resource, with the ultimate goal of producing a gold standard genome. RESULTS: Three hundred and seventeen chiefly euchromatic BACs were sequenced. A reference set of one thousand four hundred manually-annotated protein-coding genes was generated. A small RNA collection and a RNA-seq library were used to explore expression patterns and the sRNA landscape. In the sucrose and starch metabolism pathway, 16 non-redundant enzyme-encoding genes were identified. One of the sucrose pathway genes, sucrose-6-phosphate phosphohydrolase, is duplicated in sugarcane and sorghum, but not in rice and maize. A diversity analysis of the s6pp duplication region revealed haplotype-structured sequence composition. Examination of hom(e)ologous loci indicate both sequence structural and sRNA landscape variation. A synteny analysis shows that the sugarcane genome has expanded relative to the sorghum genome, largely due to the presence of transposable elements and uncharacterized intergenic and intronic sequences. CONCLUSION: This release of sugarcane genomic sequences will advance our understanding of sugarcane genetics and contribute to the development of molecular tools for breeding purposes and gene discovery.

Drought Responses in Poaceae: Exploring the Core Components of the ABA Signaling Pathway in Setaria italica and Setaria viridis.

Drought severely impacts plant development and reproduction, reducing biomass and seed number, and altering flowering patterns. Drought-tolerant Setaria italica and Setaria viridis species have emerged as prominent model species for investigating water deficit responses in the Poaceae family, the most important source of food and biofuel biomass worldwide. In higher plants, abscisic acid (ABA) regulates environmental stress responses, and its signaling entails interactions between PYR/PYL/RCAR receptors and clade A PP2C phosphatases, which in turn modulate SnRK2 kinases via reversible phosphorylation to activate ABA-responsive genes. To compare the diversity of PYR/PYL/RCAR, PP2C, and SnRK2 between S. italica and S. viridis, and their involvement in water deficit responses, we examined gene and regulatory region structures, investigated orthology relationships, and analyzed their gene expression patterns under water stress via a meta-analysis approach. Results showed that coding and regulatory sequences of PYR/PYL/RCARs, PP2Cs, and SnRK2s are highly conserved between Setaria spp., allowing us to propose pairs of orthologous genes for all the loci identified. Phylogenetic relationships indicate which clades of Setaria spp. sequences are homologous to the functionally well-characterized Arabidopsis thaliana PYR/PYL/RCAR, PP2C, and SnRK2 genes. Gene expression analysis showed a general downregulation of PYL genes, contrasting with upregulation of PP2C genes, and variable expression modulation of SnRK2 genes under drought stress. This complex network implies that ABA core signaling is a diverse and multifaceted process. Through our analysis, we identified promising candidate genes for further functional characterization, with great potential as targets for drought resistance studies, ultimately leading to advances in Poaceae biology and crop-breeding strategies.

Functional and comparative analysis of THI1 gene in grasses with a focus on sugarcane.

De novo synthesis of thiamine (vitamin B1) in plants depends on the action of thiamine thiazole synthase, which synthesizes the thiazole ring, and is encoded by the THI1 gene. Here, we investigated the evolution and diversity of THI1 in Poaceae, where C4 and C3 photosynthetic plants co-evolved. An ancestral duplication of THI1 is observed in Panicoideae that remains in many modern monocots, including sugarcane. In addition to the two sugarcane copies (ScTHI1-1 and ScTHI1-2), we identified ScTHI1-2 alleles showing differences in their sequence, indicating divergence between ScTHI1-2a and ScTHI1-2b. Such variations are observed only in the Saccharum complex, corroborating the phylogeny. At least five THI1 genomic environments were found in Poaceae, two in sugarcane, M. sinensis, and S. bicolor. The THI1 promoter in Poaceae is highly conserved at 300 bp upstream of the start codon ATG and has cis-regulatory elements that putatively bind to transcription factors associated with development, growth, development and biological rhythms. An experiment set to compare gene expression levels in different tissues across the sugarcane R570 life cycle showed that ScTHI1-1 was expressed mainly in leaves regardless of age. Furthermore, ScTHI1 displayed relatively high expression levels in meristem and culm, which varied with the plant age. Finally, yeast complementation studies with THI4-defective strain demonstrate that only ScTHI1-1 and ScTHI1-2b isoforms can partially restore thiamine auxotrophy, albeit at a low frequency. Taken together, the present work supports the existence of multiple origins of THI1 harboring genomic regions in Poaceae with predicted functional redundancy. In addition, it questions the contribution of the levels of the thiazole ring in C4 photosynthetic plant tissues or potentially the relevance of the THI1 protein activity.

Identification and functional characterization of the sugarcane (Saccharum spp.) AMT2-type ammonium transporter ScAMT3;3 revealed a presumed role in shoot ammonium remobilization.

Sugarcane (Saccharum spp.) is an important crop for sugar and bioethanol production worldwide. To maintain and increase sugarcane yields in marginal areas, the use of nitrogen (N) fertilizers is essential, but N overuse may result in the leaching of reactive N to the natural environment. Despite the importance of N in sugarcane production, little is known about the molecular mechanisms involved in N homeostasis in this crop, particularly regarding ammonium (NH4 +), the sugarcane's preferred source of N. Here, using a sugarcane bacterial artificial chromosome (BAC) library and a series of in silico analyses, we identified an AMMONIUM TRANSPORTER (AMT) from the AMT2 subfamily, sugarcane AMMONIUM TRANSPORTER 3;3 (ScAMT3;3), which is constitutively and highly expressed in young and mature leaves. To characterize its biochemical function, we ectopically expressed ScAMT3;3 in heterologous systems (Saccharomyces cerevisiae and Arabidopsis thaliana). The complementation of triple mep mutant yeast demonstrated that ScAMT3;3 is functional for NH3/H+ cotransport at high availability of NH4 + and under physiological pH conditions. The ectopic expression of ScAMT3;3 in the Arabidopsis quadruple AMT knockout mutant restored the transport capacity of 15N-NH4 + in roots and plant growth under specific N availability conditions, confirming the role of ScAMT3;3 in NH4 + transport in planta. Our results indicate that ScAMT3;3 belongs to the low-affinity transport system (Km 270.9 µM; Vmax 209.3 µmol g-1 root DW h-1). We were able to infer that ScAMT3;3 plays a presumed role in NH4 + source-sink remobilization in the shoots via phloem loading. These findings help to shed light on the functionality of a novel AMT2-type protein and provide bases for future research focusing on the improvement of sugarcane yield and N use efficiency.

Functional characterization of the sugarcane (Saccharum spp.) ammonium transporter AMT2;1 suggests a role in ammonium root-to-shoot translocation.

AMMONIUM TRANSPORTER/METHYLAMMONIUM PERMEASE/RHESUS (AMT) family members transport ammonium across membranes in all life domains. Plant AMTs can be categorized into AMT1 and AMT2 subfamilies. Functional studies of AMTs, particularly AMT1-type, have been conducted using model plants but little is known about the function of AMTs from crops. Sugarcane (Saccharum spp.) is a major bioenergy crop that requires heavy nitrogen fertilization but depends on a low carbon-footprint for competitive sustainability. Here, we identified and functionally characterized sugarcane ScAMT2;1 by complementing ammonium uptake-defective mutants of Saccharomyces cerevisiae and Arabidopsis thaliana. Reporter gene driven by the ScAMT2;1 promoter in A. thaliana revealed preferential expression in the shoot vasculature and root endodermis/pericycle according to nitrogen availability and source. Arabidopsis quadruple mutant plants expressing ScAMT2;1 driven by the CaMV35S promoter or by a sugarcane endogenous promoter produced significantly more biomass than mutant plants when grown in NH4 + and showed more 15N-ammonium uptake by roots and nitrogen translocation to shoots. In A. thaliana, ScAMT2;1 displayed a Km of 90.17 µM and Vmax of 338.99 µmoles h-1 g-1 root DW. Altogether, our results suggest that ScAMT2;1 is a functional high-affinity ammonium transporter that might contribute to ammonium uptake and presumably to root-to-shoot translocation under high NH4 + conditions.

Copia retrotransposon in the Zaprionus genus: another case of transposable element sharing with the Drosophila melanogaster subgroup.

Copia is a retrotransposon that appears to be distributed widely among the Drosophilidae subfamily. Evolutionary analyses of regulatory regions have indicated that the Copia retrotransposon evolved through both positive and purifying selection, and that horizontal transfer (HT) could also explain its patchy distribution of the among the subfamilies of the melanogaster subgroup. Additionally, Copia elements could also have transferred between melanogaster subgroup and other species of Drosophilidae-D. willistoni and Z. tuberculatus. In this study, we surveyed seven species of the Zaprionus genus by sequencing the LTR-ULR and reverse transcriptase regions, and by using RT-PCR in order to understand the distribution and evolutionary history of Copia in the Zaprionus genus. The Copia element was detected, and was transcriptionally active, in all species investigated. Structural and selection analysis revealed Zaprionus elements to be closely related to the most ancient subfamily of the melanogaster subgroup, and they seem to be evolving mainly under relaxed purifying selection. Taken together, these results allowed us to classify the Zaprionus sequences as a new subfamily-ZapCopia, a member of the Copia retrotransposon family of the melanogaster subgroup. These findings indicate that the Copia retrotransposon is an ancient component of the genomes of the Zaprionus species and broaden our understanding of the diversity of retrotransposons in the Zaprionus genus.

Evidence for sub-functionalization of tandemly duplicated XPB nucleotide excision repair genes in Arabidopsis thaliana.

Plants are continuously exposed to agents that can generate DNA lesions. Nucleotide Excision Repair (NER) is one of the repair pathways employed by plants to protect their genome, including from sunlight. The Xeroderma Pigmentosum type B (XPB) protein is a DNA helicase shown to be involved in NER and is also an essential subunitof the Transcription Factor IIH (TFIIH) complex. XPB was found to be a single copy gene in eukaryotes, but found as a tandem duplication in the plant Arabidopsis thaliana, AtXPB1 and AtXPB2. We aimed to investigate whether the XPB in tandem duplication was common within members of the Brassicaceae. We analyzed genomic DNA of species from different tribes of the family and the results indicate that the tandem duplication occurred in Camelineae tribe ancestor, of which A. thaliana belongs, at approximately 8 million years ago. Further experiments were devised to study possible functional roles for the A. thaliana AtXPB paralogs. A non-coincident expression profile of the paralogs was observed in various plant organs, developmental and cell cycle stages. AtXPB2 expression was observed in proliferating cells and clustered with the transcription of other components of the TFIIH such as p44, p52 and XPD/UVH6 along the cell cycle. AtXPB1 gene transcription, on the other hand, was enhanced specifically after UV-B irradiation in leaf trichomes. Altogether, our results reported herein suggest a functional specialization for the AtXPB paralogs: while the AtXPB2 paralog may have a role in cell proliferation and repair as XPB of other eukaryotes, the AtXPB1 paralog is most likely implicated in repair functions in highly specialized A. thaliana cells.

Multiple invasions of Gypsy and Micropia retroelements in genus Zaprionus and melanogaster subgroup of the genus Drosophila.

BACKGROUND: The Zaprionus genus shares evolutionary features with the melanogaster subgroup, such as space and time of origin. Although little information about the transposable element content in the Zaprionus genus had been accumulated, some of their elements appear to be more closely related with those of the melanogaster subgroup, indicating that these two groups of species were involved in horizontal transfer events during their evolution. Among these elements, the Gypsy and the Micropia retroelements were chosen for screening in seven species of the two Zaprionus subgenera, Anaprionus and Zaprionus. RESULTS: Screening allowed the identification of diverse Gypsy and Micropia retroelements only in species of the Zaprionus subgenus, showing that they are transcriptionally active in the sampled species. The sequences of each retroelement were closely related to those of the melanogaster species subgroup, and the most parsimonious hypothesis would be that 15 horizontal transfer events shaped their evolution. The Gypsy retroelement of the melanogaster subgroup probably invaded the Zaprionus genomes about 11 MYA. In contrast, the Micropia retroelement may have been introduced into the Zaprionus subgenus and the melanogaster subgroup from an unknown donor more recently (~3 MYA). CONCLUSION: Gypsy and Micropia of Zaprionus and melanogaster species share similar evolutionary patterns. The sharing of evolutionary, ecological and ethological features probably allowed these species to pass through a permissive period of transposable element invasion, explaining the proposed waves of horizontal transfers.

Estimating genomic instability mediated by Alu retroelements in breast cancer.

Alu-PCR is a relatively simple technique that can be used to investigate genomic instability in cancer. This technique allows identification of the loss, gain or amplification of gene sequences based on the analysis of segments between two Alu elements coupled with quantitative and qualitative analyses of the profiles obtained from tumor samples, surgical margins and blood. In this work, we used Alu-PCR to identify gene alterations in ten patients with invasive ductal breast cancer. Several deletions and insertions were identified, indicating genomic instability in the tumor and adjacent normal tissue. Although not associated with specific genes, the alterations, which involved chromosomal bands 1p36.23, 1q41, 11q14.3, 13q14.2, occurred in areas of well-known genomic instability in breast and other types of cancer. These results indicate the potential usefulness of Alu-PCR in identifying altered gene sequences in breast cancer. However, caution is required in its application since the Alu primer can produce non-specific amplification.

The population genetic structure approach adds new insights into the evolution of plant LTR retrotransposon lineages.

Long terminal repeat retrotransposons (LTR-RTs) in plant genomes differ in abundance, structure and genomic distribution, reflecting the large number of evolutionary lineages. Elements within lineages can be considered populations, in which each element is an individual in its genomic environment. In this way, it would be reasonable to apply microevolutionary analyses to understand transposable element (TE) evolution, such as those used to study the genetic structure of natural populations. Here, we applied a Bayesian method to infer genetic structure of populations together with classical phylogenetic and dating tools to analyze LTR-RT evolution using the monocot Setaria italica as a model species. In contrast to a phylogeny, the Bayesian clusterization method identifies populations by assigning individuals to one or more clusters according to the most probabilistic scenario of admixture, based on genetic diversity patterns. In this work, each LTR-RT insertion was considered to be one individual and each LTR-RT lineage was considered to be a single species. Nine evolutionary lineages of LTR-RTs were identified in the S. italica genome that had different genetic structures with variable numbers of clusters and levels of admixture. Comprehensive analysis of the phylogenetic, clusterization and time of insertion data allowed us to hypothesize that admixed elements represent sequences that harbor ancestral polymorphic sequence signatures. In conclusion, application of microevolutionary concepts in genome evolution studies is suitable as a complementary approach to phylogenetic analyses to address the evolutionary history and functional features of TEs.

Helena and BS: Two Travellers between the Genera Drosophila and Zaprionus.

The frequency of horizontal transfers of transposable elements (HTTs) varies among the types of elements according to the transposition mode and the geographical and temporal overlap of the species involved in the transfer. The drosophilid species of the genus Zaprionus and those of the melanogaster, obscura, repleta, and virilis groups of the genus Drosophila investigated in this study shared space and time at some point in their evolutionary history. This is particularly true of the subgenus Zaprionus and the melanogaster subgroup, which overlapped both geographically and temporally in Tropical Africa during their period of origin and diversification. Here, we tested the hypothesis that this overlap may have facilitated the transfer of retrotransposons without long terminal repeats (non-LTRs) between these species. We estimated the HTT frequency of the non-LTRs BS and Helena at the genome-wide scale by using a phylogenetic framework and a vertical and horizontal inheritance consistence analysis (VHICA). An excessively low synonymous divergence among distantly related species and incongruities between the transposable element and species phylogenies allowed us to propose at least four relatively recent HTT events of Helena and BS involving ancestors of the subgroup melanogaster and ancestors of the subgenus Zaprionus during their concomitant diversification in Tropical Africa, along with older possible events between species of the subgenera Drosophila and Sophophora. This study provides the first evidence for HTT of non-LTRs retrotransposons between Drosophila and Zaprionus, including an in-depth reconstruction of the time frame and geography of these events.

Transposon display supports transpositional activity of P elements in species of the saltans group of Drosophila.

Mobilization of two P element subfamilies (canonical and O-type) from Drosophila sturtevanti and D. saltans was evaluated for copy number and transposition activity using the transposon display (TD) technique. Pairwise distances between strains regarding the insertion polymorphism profile were estimated. Amplification of the P element based on copy number estimates was highly variable among the strains (D. sturtevanti, canonical 20.11, O-type 9.00; D. saltans, canonical 16.4, O-type 12.60 insertions, on average). The larger values obtained by TD compared to our previous data by Southern blotting support the higher sensitivity of TD over Southern analysis for estimating transposable element copy numbers. The higher numbers of the canonical P element and the greater divergence in its distribution within the genome of D. sturtevanti (24.8%) compared to the O-type (16.7%), as well as the greater divergence in the distribution of the canonical P element, between the D. sturtevanti (24.8%) and the D. saltans (18.3%) strains, suggest that the canonical element occupies more sites within the D. sturtevanti genome, most probably due to recent transposition activity. These data corroborate the hypothesis that the O-type is the oldest subfamily of P elements in the saltans group and suggest that the canonical P element is or has been transpositionally active until more recently in D. sturtevanti.

Analysis of Three Sugarcane Homo/Homeologous Regions Suggests Independent Polyploidization Events of Saccharum officinarum and Saccharum spontaneum.

Whole genome duplication has played an important role in plant evolution and diversification. Sugarcane is an important crop with a complex hybrid polyploid genome, for which the process of adaptation to polyploidy is still poorly understood. In order to improve our knowledge about sugarcane genome evolution and the homo/homeologous gene expression balance, we sequenced and analyzed 27 BACs (Bacterial Artificial Chromosome) of sugarcane R570 cultivar, containing the putative single-copy genes LFY (seven haplotypes), PHYC (four haplotypes), and TOR (seven haplotypes). Comparative genomic approaches showed that these sugarcane loci presented a high degree of conservation of gene content and collinearity (synteny) with sorghum and rice orthologous regions, but were invaded by transposable elements (TE). All the homo/homeologous haplotypes of LFY, PHYC, and TOR are likely to be functional, because they are all under purifying selection (dN/dS ≪ 1). However, they were found to participate in a nonequivalently manner to the overall expression of the corresponding gene. SNPs, indels, and amino acid substitutions allowed inferring the S. officinarum or S. spontaneum origin of the TOR haplotypes, which further led to the estimation that these two sugarcane ancestral species diverged between 2.5 and 3.5 Ma. In addition, analysis of shared TE insertions in TOR haplotypes suggested that two autopolyploidization may have occurred in the lineage that gave rise to S. officinarum, after its divergence from S. spontaneum.

Phytochrome Interacting Factors (PIFs) in Solanum lycopersicum: Diversity, Evolutionary History and Expression Profiling during Different Developmental Processes.

Although the importance of light for tomato plant yield and edible fruit quality is well known, the PHYTOCHROME INTERACTING FACTORS (PIFs), main components of phytochrome-mediated light signal transduction, have been studied almost exclusively in Arabidopsis thaliana. Here, the diversity, evolution and expression profile of PIF gene subfamily in Solanum lycopersicum was characterized. Eight tomato PIF loci were identified, named SlPIF1a, SlPIF1b, SlPIF3, SlPIF4, SlPIF7a, SlPIF7b, SlPIF8a and SlPIF8b. The duplication of SlPIF1, SlPIF7 and SlPIF8 genes were dated and temporally coincided with the whole-genome triplication event that preceded tomato and potato divergence. Different patterns of mRNA accumulation in response to light treatments were observed during seedling deetiolation, dark-induced senescence, diel cycle and fruit ripening. SlPIF4 showed similar expression profile as that reported for A. thaliana homologs, indicating an evolutionary conserved function of PIF4 clade. A comprehensive analysis of the evolutionary and transcriptional data allowed proposing that duplicated SlPIFs have undergone sub- and neofunctionalization at mRNA level, pinpointing the importance of transcriptional regulation for the maintenance of duplicated genes. Altogether, the results indicate that genome polyploidization and functional divergence have played a major role in diversification of the Solanum PIF gene subfamily.

Plant degreening: evolution and expression of tomato (Solanum lycopersicum) dephytylation enzymes.

Chlorophyll is the most abundant pigment on earth and even though it is known that its high photo-excitability necessitates a tight regulation of its degradation pathway, to date there are still several steps in chlorophyll breakdown that remain obscure. In order to better understand the 'degreening' processes that accompany leaf senescence and fruit ripening, we characterized the enzyme-encoding genes involved in dephytylation from tomato (Solanum lycopersicum). A single pheophytinase (PPH) gene and four chlorophyllase (CLH) genes were identified in the tomato genome. A phenetic analysis revealed two groups of CLHs in eudicot species and further evolutionary analysis indicated that these enzymes are under diverse selection pressures. A comprehensive expression profile analysis also suggested functional specificity for these dephytylating enzymes. The integrated analysis allows us to propose three general roles for chlorophyll dephytylation: i) PPH, which is under high selective constraint, is responsible for chlorophyll degradation during developmentally programed physiological processes; ii) Group I CLHs, which are under relaxed selection constraint, respond to environmental and hormonal stimuli and play a role in plant adaptation plasticity; and iii) Group II CLHs, which are also under high selective constraint, are mostly involved in chlorophyll recycling.

Virus-like attachment sites and plastic CpG islands:landmarks of diversity in plant Del retrotransposons.

Full-length Del elements from ten angiosperm genomes, 5 monocot and 5 dicot, were retrieved and putative attachment (att) sites were identified. In the 2432 Del elements, two types of U5 att sites and a single conserved type of U3 att site were identified. Retroviral att sites confer specificity to the integration process, different att sites types therefore implies lineage specificity. While some features are common to all Del elements, CpG island patterns within the LTRs were particular to lineage specific clusters. All eudicot copies grouped into one single clade while the monocots harbour a more diverse collection of elements. Furthermore, full-length Del elements and truncated copies were unevenly distributed amongst chromosomes. Elements of Del lineage are organized in plants into three clusters and each cluster is composed of elements with distinct LTR features. Our results suggest that the Del lineage efficiently amplified in the monocots and that one branch is probably a newly emerging sub-lineage. Finally, sequences in all groups are under purifying selection. These results show the LTR region is dynamic and important in the evolution of LTR-retrotransposons, we speculate that it is a trigger for retrotransposon diversification.