Claudio Varotto

CHLOROPHYLL BINDING TO MONOMERIC LIGHT-HARVESTING COMPLEX : A MUTATION ANALYSIS OF CHROMOPHORE-BINDING RESIDUES

The chromophore binding properties of the higher plant light-harvesting complex II have been studied by site-directed mutagenesis of pigment-binding residues. Mutant apoproteins were overexpressed in Escherichia coli and then refoldedin vitro with purified chromophores to yield holoproteins selectively affected in chlorophyll-binding sites. Biochemical and spectroscopic characterization showed a specific loss of pigments and absorption spectral forms for each mutant, thus allowing identification of the chromophores bound to most of the binding sites. On these bases a map for the occupancy of individual sites by chlorophyll a and chlorophyll b is proposed. In some cases a single mutation led to the loss of more than one chromophore indicating that four chlorophylls and one xanthophyll could be bound by pigment-pigment interactions. Differential absorption spectroscopy allowed identification of the Qy transition energy level for each chlorophyll within the complex. It is shown that not only site selectivity is largely conserved between light-harvesting complex II and CP29 but also the distribution of absorption forms among different protein domains, suggesting conservation of energy transfer pathways within the protein and outward to neighbor subunits of the photosystem.

The neoxanthin binding site of the major light harvesting complex (LHCII) from higher plants

The localisation of the xanthophyll neoxanthin within the structure of the major light harvesting complex (LHCII) of higher plants has been investigated by site‐directed mutagenesis and spectroscopic methods. Mutation analysis performed on pigment binding sites in different helix domains leads to selective loss of neoxanthin for mutations on helix C thus localising this pigment between the helix C and helix A/B domains. Recombinant proteins binding two lutein molecules per polypeptide but lacking neoxanthin have been used in order to determine the contribution of neoxanthin to the absorption and linear dichroism spectra. The data were used to derive the orientation of the neoxanthin transition moment, lying in the polyene chain, which was thus determined to form an angle of 57±1.5° with respect to the normal to the membrane plane where the protein is inserted. On the basis of these results we propose a model for the localisation of the carotenoid site in the LHCII structure which is still unresolved.

isomiRex: Web-based identification of microRNAs, isomiR variations and differential expression using next-generation sequencing datasets

We present an open‐access web platform isomiRex, to identify isomiRs and on the fly graphical visualization of the differentially expressed miRNAs in control as well as treated library. The open‐access web‐platform is not restricted only to NGS sequence dataset from animals and potentially analyzes a wider dataset for plants, animals and viral NGS dataset supporting miRBase (version 19 supporting 193 species). The platform can handle the bloated amount of the read counts and reports the annotated microRNAs from plant, animal and viral NGS datasets. isomiRex also provides an estimation of the the isomiRs, of miRNAs with higher copy number relative to their mature reference sequences indexed in miRBase (version 19 supporting 193 species). Visually enhanced graphs potentially display differentially expressed isomiRs, which will help the user to demonstrate and correlate the abundance of the isomiR as a signature event to the specific condition. An additional module for estimating the differential expression has been implemented allowing the users to postulate the differential expression across the user input samples. The developed web‐platform can be accessed at http://bioinfo1.uni‐plovdiv.bg/isomiRex/.

Spatiotemporal reconstruction of the Aquilegia rapid radiation through next‐generation sequencing of rapidly evolving cpDNA regions

Aquilegia is a well-known model system in the field of evolutionary biology, but obtaining a resolved and well-supported phylogenetic reconstruction for the genus has been hindered by its recent and rapid diversification. Here, we applied 454 next-generation sequencing to PCR amplicons of 21 of the most rapidly evolving regions of the plastome to generate c. 24 kb of sequences from each of 84 individuals from throughout the genus. The resulting phylogeny has well-supported resolution of the main lineages of the genus, although recent diversification such as in the European taxa remains unresolved. By producing a chronogram of the whole Ranunculaceae family based on published data, we inferred calibration points for dating the Aquilegia radiation. The genus originated in the upper Miocene c. 6.9 million yr ago (Ma) in Eastern Asia, and diversification occurred c. 4.8 Ma with the split of two main clades, one colonizing North America, and the other Western Eurasia through the mountains of Central Asia. This was followed by a back-to-Asia migration, originating from the European stock using a North Asian route. These results provide the first backbone phylogeny and spatiotemporal reconstruction of the Aquilegia radiation, and constitute a robust framework to address the adaptative nature of speciation within the group.

Plastome organization and evolution of chloroplast genes in Cardamine species adapted to contrasting habitats

BackgroundPlastid genomes, also known as plastomes, are shaped by the selective forces acting on the fundamental cellular functions they code for and thus they are expected to preserve signatures of the adaptive path undertaken by different plant species during evolution. To identify molecular signatures of positive selection associated to adaptation to contrasting ecological niches, we sequenced with Solexa technology the plastomes of two congeneric Brassicaceae species with different habitat preference, Cardamine resedifolia and Cardamine impatiens.ResultsFollowing in-depth characterization of plastome organization, repeat patterns and gene space, the comparison of the newly sequenced plastomes between each other and with 15 fully sequenced Brassicaceae plastomes publically available in GenBank uncovered dynamic variation of the IR boundaries in the Cardamine lineage. We further detected signatures of positive selection in ten of the 75 protein-coding genes of the examined plastomes, identifying a range of chloroplast functions putatively involved in adaptive processes within the family. For instance, the three residues found to be under positive selection in RUBISCO could possibly be involved in the modulation of RUBISCO aggregation/activation and enzymatic specificty in Brassicaceae. In addition, our results points to differential evolutionary rates in Cardamine plastomes.ConclusionsOverall our results support the existence of wider signatures of positive selection in the plastome of C. resedifolia, possibly as a consequence of adaptation to high altitude environments. We further provide a first characterization of the selective patterns shaping the Brassicaceae plastomes, which could help elucidate the driving forces underlying adaptation and evolution in this important plant family.

Fuelling genetic and metabolic exploration of C3 bioenergy crops through the first reference transcriptome of Arundo donax L.

The development of inexpensive and highly productive biomass sources of biofuel is a priority in global climate change biology. Arundo donax, also known as the giant reed, is recognized as one of the most promising nonfood bioenergy crops in Europe. Despite its relevance, to date no genomic resources are available to support the characterization of the developmental, adaptive and metabolic traits underlying the high productivity of this nonmodel species. We hereby present the first report on the de novo assembly of bud, culm, leaf and root transcriptomes of A. donax, which can be accessed through a customized BLAST server (http://ecogenomics.fmach.it/arundo/) for mining and exploring the genetic potential of this species. Based on functional annotation and homology comparison to 19 prospective biofuel Poaceae species, we provide the first genomic view of this so far unexplored crop and indicate the model species with highest potential for comparative genomics approaches. The analysis of the transcriptome reveals strong differences in the enrichment of the Gene Ontology categories and the relative expression among different organs, which can guide future efforts for functional genomics or genetic improvement of A. donax. A set of homologs to key genes involved in lignin, cellulose, starch, lipid metabolism and in the domestication of other crops is discussed to provide a platform for possible enhancement of productivity and saccharification efficiency in A. donax.

Evolution of MIR168 paralogs in Brassicaceae

BackgroundIn plants, expression of ARGONAUTE1 (AGO1), the catalytic subunit of the RNA-Induced Silencing Complex responsible for post-transcriptional gene silencing, is controlled through a feedback loop involving the miR168 microRNA. This complex auto-regulatory loop, composed of miR168-guided AGO1-catalyzed cleavage of AGO1 mRNA and AGO1-mediated stabilization of miR168, was shown to ensure the maintenance of AGO1 homeostasis that is pivotal for the correct functioning of the miRNA pathway.ResultsWe applied different approaches to studying the genomic organization and the structural and functional evolution of MIR168 homologs in Brassicaeae. A whole genome comparison of Arabidopsis and poplar, phylogenetic footprinting and phylogenetic reconstruction were used to date the duplication events originating MIR168 homologs in these genomes. While orthology was lacking between Arabidopsis and poplar MIR168 genes, we successfully isolated orthologs of both loci present in Arabidopsis (MIR168a and MIR168b) from all the Brassicaceae species analyzed, including the basal species Aethionema grandiflora, thus indicating that (1) independent duplication events took place in Arabidopsis and poplar lineages and (2) the origin of MIR168 paralogs predates both the Brassicaceae radiation and the Arabidopsis alpha polyploidization. Different phylogenetic footprints, corresponding to known functionally relevant regions (transcription starting site and double-stranded structures responsible for microRNA biogenesis and function) or for which functions could be proposed, were found to be highly conserved among MIR168 homologs. Comparative predictions of the identified microRNAs also indicate extreme conservation of secondary structure and thermodynamic stability.ConclusionWe used a comparative phylogenetic footprinting approach to identify the structural and functional constraints that shaped MIR168 evolution in Brassicaceae. Although their duplication happened at least 40 million years ago, we found evidence that both MIR168 paralogs have been maintained throughout the evolution of Brassicaceae, most likely functionally as indicated by the extremely high conservation of functionally relevant regions, predicted secondary structure and thermodynamic profile. Interestingly, the expression patterns observed in Arabidopsis indicate that MIR168b underwent partial subfunctionalization as determined by the experimental characterization of its expression pattern provided in this study. We found further evolutionary evidence that pre-miR168 lower stem (the RNA-duplex structure adjacent to the miR-miR* stem) is significantly longer than animal lower stems and probably plays a relevant role in multi-step miR168 biogenesis.

Rates of evolution in stress-related genes are associated with habitat preference in two Cardamine lineages

BackgroundElucidating the selective and neutral forces underlying molecular evolution is fundamental to understanding the genetic basis of adaptation. Plants have evolved a suite of adaptive responses to cope with variable environmental conditions, but relatively little is known about which genes are involved in such responses. Here we studied molecular evolution on a genome-wide scale in two species of Cardamine with distinct habitat preferences: C. resedifolia, found at high altitudes, and C. impatiens, found at low altitudes. Our analyses focussed on genes that are involved in stress responses to two factors that differentiate the high- and low-altitude habitats, namely temperature and irradiation.ResultsHigh-throughput sequencing was used to obtain gene sequences from C. resedifolia and C. impatiens. Using the available A. thaliana gene sequences and annotation, we identified nearly 3,000 triplets of putative orthologues, including genes involved in cold response, photosynthesis or in general stress responses. By comparing estimated rates of molecular substitution, codon usage, and gene expression in these species with those of Arabidopsis, we were able to evaluate the role of positive and relaxed selection in driving the evolution of Cardamine genes. Our analyses revealed a statistically significant higher rate of molecular substitution in C. resedifolia than in C. impatiens, compatible with more efficient positive selection in the former. Conversely, the genome-wide level of selective pressure is compatible with more relaxed selection in C. impatiens. Moreover, levels of selective pressure were heterogeneous between functional classes and between species, with cold responsive genes evolving particularly fast in C. resedifolia, but not in C. impatiens.ConclusionsOverall, our comparative genomic analyses revealed that differences in effective population size might contribute to the differences in the rate of protein evolution and in the levels of selective pressure between the C. impatiens and C. resedifolia lineages. The within-species analyses also revealed evolutionary patterns associated with habitat preference of two Cardamine species. We conclude that the selective pressures associated with the habitats typical of C. resedifolia may have caused the rapid evolution of genes involved in cold response.

Dissection of early transcriptional responses to water stress in Arundo donax L. by unigene-based RNA-seq

BackgroundArundo donax L. (Poaceae) is considered one of the most promising energy crops in the Mediterranean region because of its high biomass yield and low input requirements, but to date no information on its transcriptional responses to water stress is available.ResultsWe obtained by Illumina-based RNA-seq the whole root and shoot transcriptomes of young A. donax plants subjected to osmotic/water stress with 10 and 20xa0% polyethylene glycol (PEG; 3 biological replicates/organ/condition corresponding to 18 RNA-Seq libraries), and identified a total of 3034 differentially expressed genes. Blast-based mining of stress-related genes indicated the higher responsivity of roots compared to shoots at the early stages of water stress especially under the milder PEG treatment, with a majority of genes responsive to salt, oxidative, and dehydration stress. Analysis of gene ontology terms underlined the qualitatively different responses between root and shoot tissues. Among the most significantly enriched metabolic pathways identified using a Fisher’s exact test with FDR correction, a crucial role was played in both shoots and roots by genes involved in the signaling cascade of abscisic acid. We further identified relatively large organ-specific differences in the patterns of drought-related transcription factor AP2-EREBP, AUX/IAA, MYB, bZIP, C2H2, and GRAS families, which may underlie the transcriptional reprogramming differences between organs. Through comparative analyses with major Poaceae species based on Blast, we finally identified a set of 53 orthologs that can be considered as a core of evolutionary conserved genes important to mediate water stress responses in the family.ConclusionsThis study provides the first characterization of A. donax transcriptome in response to water stress, thus shedding novel light at the molecular level on the mechanisms of stress response and adaptation in this emerging bioenergy species. The inventory of early-responsive genes to water stress identified could constitute useful markers of the physiological status of A. donax and be a basis for the improvement of its productivity under water limitation. The full water-stressed A. donax transcriptome is available for Blast-based homology searches through a dedicated web server (http://ecogenomics.fmach.it/arundo/).

MULTIPLE LIGHT-HARVESTING II POLYPEPTIDES FROM MAIZE MESOPHYLL CHLOROPLASTS ARE DISTINCT GENE PRODUCTS

The major light-harvesting complex of photosystem II in higher plants is known as LHCII. It is composed of a number of chlorophyll-binding proteins sharing epitopes with each other. The number of apoproteins resolved by fully denaturing sodium dodecylsulfate polyacrylamide gel electrophoresis varies in different species. In order to know if this heterogeneity is caused by the expression of a number of homologous genes or if it is the product of post-translational modifications, we have resolved the six major apoproteins of Zea mays LHCII. Each protein is purified to homogeneity, subjected to direct protein sequencing and the sequences compared with those deduced from lhcb genes in maize and other organisms. All of the six proteins are distinct gene products, since they show differences in their primary structure. Three apoproteins are identified as products of type I lhcb genes and one each as type II and type III gene products. A sixth protein does not fit the requirements for any of the lhcb genes so far cloned and is therefore probably the product of an lhcb gene type not yet described. Our results clearly show that the major source of LHCII protein heterogeneity is the expression of many lhcb genes. Fractionation of maize LHCII by non-denaturing flat-bed isoelectric focusing resolves at least five major isoforms showing distinct differences in their polypeptide composition and also differing in their spectroscopic properties, thus suggesting that individual Lhcb gene products have distinct pigment-binding properties.