Unique genes in plants: specificities and conserved features throughout evolution

Functional analysis of fasciclin-like arabinogalactan in carotenoid synthesis during tomato fruit ripening

Carotenoids are important pigmented nutrients synthesized by tomato fruits during ripening. To reveal the molecular mechanism underlying carotenoid synthesis during tomato fruit ripening, we analyzed carotenoid metabolites and transcriptomes in six development stages of tomato fruits. A total of thirty different carotenoids were detected and quantified in tomato fruits from 10 to 60 DPA. Based on differential gene expression profiles and WGCNA, we explored several genes that were highly significant and negatively correlated with lycopene, all of which encode fasciclin-like arabinogalactan proteins (FLAs). The FLAs are involved in plant signal transduction, however the functional role of these proteins has not been studied in tomato. Genome-wide analysis revealed that cultivated and wild tomato species contained 18 to 22 FLA family members, clustered into four groups, and mainly evolved by means of segmental duplication. The functional characterization of FLAs showed that silencing of SlFLA1, 5, and 13 were found to contribute to the early coloration of tomato fruits, and the expression of carotenoid synthesis-related genes was up-regulated in fruits that changed phenotypically, especially in SlFLA13-silenced plants. Furthermore, the content of multiple carotenoids (including (E/Z)-phytoene, lycopene, γ-carotene, and α-carotene) was significantly increased in SlFLA13-silenced fruits, suggesting that SlFLA13 has a potential inhibitory function in regulating carotenoid synthesis in tomato fruits. The results of the present study broaden the idea of analyzing the biological functions of tomato FLAs and preliminary evidence for the inhibitory role of SlFLA13 in carotenoid synthesis in fruit, providing the theoretical basis and a candidate for improving tomato fruit quality.

Fitness difference between two synonymous mutations of Phytophthora infestans ATP6 gene

BACKGROUND

Sequence variation produced by mutation provides the ultimate source of natural selection for species adaptation. Unlike nonsynonymous mutation, synonymous mutations are generally considered to be selectively neutral but accumulating evidence suggests they also contribute to species adaptation by regulating the flow of genetic information and the development of functional traits. In this study, we analysed sequence characteristics of ATP6, a housekeeping gene from 139 Phytophthora infestans isolates, and compared the fitness components including metabolic rate, temperature sensitivity, aggressiveness, and fungicide tolerance among synonymous mutations.

RESULTS

We found that the housekeeping gene exhibited low genetic variation and was represented by two major synonymous mutants at similar frequency (0.496 and 0.468, respectively). The two synonymous mutants were generated by a single nucleotide substitution but differed significantly in fitness as well as temperature-mediated spatial distribution and expression. The synonymous mutant ending in AT was more common in cold regions and was more expressed at lower experimental temperature than the synonymous mutant ending in GC and vice versa.

CONCLUSION

Our results are consistent with the argument that synonymous mutations can modulate the adaptive evolution of species including pathogens and have important implications for sustainable disease management, especially under climate change.

GERDH: an interactive multi-omics database for cross-species data mining in horticultural crops

Horticultural plants contribute immensely to the quality of human's life. The rapid development of omics studies on horticultural plants has resulted in large volumes of valuable growth- and development-related data. Genes that are essential for growth and development are highly conserved in evolution. Cross-species data mining reduces the impact of species heterogeneity and has been extensively used for conserved gene identification. Owing to the lack of a comprehensive database for cross-species data mining using multi-omics data from all horticultural plant species, the current resources in this field are far from satisfactory. Here, we introduce GERDH (https://dphdatabase.com), a database platform for cross-species data mining among horticultural plants, based on 12 961 uniformly processed publicly available omics libraries from more than 150 horticultural plant accessions, including fruits, vegetables and ornamental plants. Important and conserved genes that are essential for a specific biological process can be obtained by cross-species analysis module with interactive web-based data analysis and visualization. Moreover, GERDH is equipped with seven online analysis tools, including gene expression, in-species analysis, epigenetic regulation, gene co-expression, enrichment/pathway and phylogenetic analysis. By interactive cross-species analysis, we identified key genes contributing to postharvest storage. By gene expression analysis, we explored new functions of CmEIN3 in flower development, which was validated by transgenic chrysanthemum analysis. We believe that GERDH will be a useful resource for key gene identification and will allow for omics big data to be more available and accessible to horticultural plant community members.

Welcome to the big leaves: Best practices for improving genome annotation in non-model plant genomes

Premise

Robust standards to evaluate quality and completeness are lacking in eukaryotic structural genome annotation, as genome annotation software is developed using model organisms and typically lacks benchmarking to comprehensively evaluate the quality and accuracy of the final predictions. The annotation of plant genomes is particularly challenging due to their large sizes, abundant transposable elements, and variable ploidies. This study investigates the impact of genome quality, complexity, sequence read input, and method on protein-coding gene predictions.

Methods

The impact of repeat masking, long-read and short-read inputs, and de novo and genome-guided protein evidence was examined in the context of the popular BRAKER and MAKER workflows for five plant genomes. The annotations were benchmarked for structural traits and sequence similarity.

Results

Benchmarks that reflect gene structures, reciprocal similarity search alignments, and mono-exonic/multi-exonic gene counts provide a more complete view of annotation accuracy. Transcripts derived from RNA-read alignments alone are not sufficient for genome annotation. Gene prediction workflows that combine evidence-based and ab initio approaches are recommended, and a combination of short and long reads can improve genome annotation. Adding protein evidence from de novo assemblies, genome-guided transcriptome assemblies, or full-length proteins from OrthoDB generates more putative false positives as implemented in the current workflows. Post-processing with functional and structural filters is highly recommended.

Discussion

While the annotation of non-model plant genomes remains complex, this study provides recommendations for inputs and methodological approaches. We discuss a set of best practices to generate an optimal plant genome annotation and present a more robust set of metrics to evaluate the resulting predictions.

Heterologous expression of a lycophyte protein enhances angiosperm seedling vigor

Seedling vigor is a key agronomic trait that determines juvenile plant performance. Angiosperm seeds develop inside fruits and are connected to the mother plant through vascular tissues. Their formation requires plant-specific genes, such as BREVIS RADIX (BRX) in Arabidopsis thaliana roots. BRX family proteins are found throughout the euphyllophytes but also occur in non-vascular bryophytes and non-seed lycophytes. They consist of four conserved domains, including the tandem BRX domains. We found that bryophyte or lycophyte BRX homologs can only partially substitute for Arabidopsis BRX (AtBRX) because they miss key features in the linker between the BRX domains. Intriguingly, however, expression of a BRX homolog from the lycophyte Selaginella moellendorffii (SmBRX) in an A. thaliana wild-type background confers robustly enhanced root growth vigor that persists throughout the life cycle. This effect can be traced to a substantial increase in seed and embryo size, is associated with enhanced vascular tissue proliferation, and can be reproduced with a modified, SmBRX-like variant of AtBRX. Our results thus suggest that BRX variants can boost seedling vigor and shed light on the activity of ancient, non-angiosperm BRX family proteins.

Bird nests as botanical time capsules: DNA barcoding identifies the contents of contemporary and historical nests

Bird nests in natural history collections are an abundant yet vastly underutilized source of genetic information. We sequenced the nuclear ribosomal internal transcribed spacer to identify plant species used as nest material in two contemporary (2003 and 2018) and two historical (both 1915) nest specimens constructed by Song Sparrows (Melospiza melodia) and Savannah Sparrows (Passerculus sandwichensis). A total of 13 (22%) samples yielded single, strong bands that could be identified using GenBank resources: six plants (Angiospermae), six green algae (Chlorophyta), and one ciliate (Ciliophora). Two native plant species identified in the nests included Festuca microstachys, which was introduced to the nest collection site by restoration practitioners, and Rosa californica, identified in a nest collected from a lost habitat that existed about 100 years ago. Successful sequencing was correlated with higher sample mass and DNA quality, suggesting future studies should select larger pieces of contiguous material from nests and materials that appear to have been fresh when incorporated into the nest. This molecular approach was used to distinguish plant species that were not visually identifiable, and did not require disassembling the nest specimens as is a traditional practice with nest material studies. The many thousands of nest specimens in natural history collections hold great promise as sources of genetic information to address myriad ecological questions.

Highly efficient multiplex editing: one-shot generation of 8× Nicotiana benthamiana and 12× Arabidopsis mutants

Genome editing by RNA-guided nucleases, such as SpCas9, has been used in numerous different plant species. However, to what extent multiple independent loci can be targeted simultaneously by multiplexing has not been well documented. Here, we developed a toolkit, based on a highly intron-optimized zCas9i gene, which allows assembly of nuclease constructs expressing up to 32 single guide RNAs (sgRNAs). We used this toolkit to explore the limits of multiplexing in two major model species, and report on the isolation of transgene-free octuple (8×) Nicotiana benthamiana and duodecuple (12×) Arabidopsis thaliana mutant lines in a single generation (T1 and T2 , respectively). We developed novel counter-selection markers for N. benthamiana, most importantly Sl-FAST2, comparable to the well-established Arabidopsis seed fluorescence marker, and FCY-UPP, based on the production of toxic 5-fluorouracil in the presence of a precursor. Targeting eight genes with an array of nine different sgRNAs and relying on FCY-UPP for selection of non-transgenic T1 , we identified N. benthamiana mutant lines with astonishingly high efficiencies: All analyzed plants carried mutations in all genes (approximately 112/116 target sites edited). Furthermore, we targeted 12 genes by an array of 24 sgRNAs in A. thaliana. Efficiency was significantly lower in A. thaliana, and our results indicate Cas9 availability is the limiting factor in such higher-order multiplexing applications. We identified a duodecuple mutant line by a combination of phenotypic screening and amplicon sequencing. The resources and results presented provide new perspectives for how multiplexing can be used to generate complex genotypes or to functionally interrogate groups of candidate genes.

Involvement of ethylene biosynthesis and perception during germination of dormant Avena fatua L. caryopses induced by KAR1 or GA3

Germination of primary dormant wild oat caused by KAR₁ or GA₃ is associated with ACC accumulation and increased ethylene production shortly before radicle protrusion as a result of the non-transcriptional and transcriptional activation of ACS and ACO enzymes, respectively. Response to both compounds involves the modulation of ethylene sensitivity through ethylene receptor genes. Harvested Avena fatua caryopses are primary dormant and, therefore, germinated poorly at 20 °C. Karrikin 1 (KAR₁), which action probably requires endogenous gibberellins (GAs), and gibberellin A₃ (GA₃) was found to induce dormant caryopses to germinate. The stimulatory effects were accompanied by the activation of the ethylene biosynthesis pathway and depended on undisturbed ethylene perception. KAR₁ and GA₃ promoted 1-aminocyclopropane-1-carboxylic acid (ACC) accumulation during coleorhizae emergence and ethylene production shortly prior to the radicle protrusion, which resulted from the enhanced activity of two ethylene biosynthesis enzymes, ACC synthase (ACS) and ACC oxidase (ACO). The inhibitor of ACS adversely affected beneficial impacts of both KAR₁ and GA₃ on A. fatua caryopses germination, while the inhibitor of ACO more efficiently impeded the GA₃ effect. The inhibitors of ethylene action markedly lowered germination in response to KAR₁ and GA₃. Gene expression studies preceded by the identification of several genes related to ethylene biosynthesis (AfACS6, AfACO1, and AfACO5) and perception (AfERS1b, AfERS1c, AfERS2, AfETR2, AfETR3, and AfETR4) provided further evidence for the engagement of ethylene in KAR₁ and GA₃ induced germination of A. fatua caryopses. Both AfACO1 and AfACO5 were upregulated, whereas AfACS6 remained unaffected by the treatment. This suggests the existence of different regulatory mechanisms of enzymatic activity, transcriptional for ACO and non-transcriptional for ACS. During imbibition in water, AfERS1b was stronger expressed than other receptor genes. In the presence of KAR₁ or GA₃, the expression of AfETR3 was substantially induced. Differential expression of ethylene receptor genes implies the modulation of caryopses sensitivity adjusted to ethylene availability and suggests the functional diversification of individual receptors.

Comparative transcriptomics uncovers differences in photoautotrophic versus photoheterotrophic modes of nutrition in relation to secondary metabolites biosynthesis in Swertia chirayita

Swertia chirayita is a high-value medicinal herb exhibiting antidiabetic, hepatoprotective, anticancer, antiediematogenic and antipyretic properties. Scarcity of its plant material has necessitated in vitro production of therapeutic metabolites; however, their yields were low compared to field grown plants. Possible reasons for this could be differences in physiological and biochemical processes between plants grown in photoautotrophic versus photoheterotrophic modes of nutrition. Comparative transcriptomes of S. chirayita were generated to decipher the crucial molecular components associated with the secondary metabolites biosynthesis. Illumina HiSeq sequencing yielded 57,460 and 43,702 transcripts for green house grown (SCFG) and tissue cultured (SCTC) plants, respectively. Biological role analysis (GO and COG assignments) revealed major differences in SCFG and SCTC transcriptomes. KEGG orthology mapped 351 and 341 transcripts onto secondary metabolites biosynthesis pathways for SCFG and SCTC transcriptomes, respectively. Nineteen out of 30 genes from primary metabolism showed higher in silico expression (FPKM) in SCFG versus SCTC, possibly indicating their involvement in regulating the central carbon pool. In silico data were validated by RT-qPCR using a set of 16 genes, wherein 10 genes showed similar expression pattern across both the methods. Comparative transcriptomes identified differentially expressed transcription factors and ABC-type transporters putatively associated with secondary metabolism in S. chirayita. Additionally, functional classification was performed using NCBI Biosystems database. This study identified the molecular components implicated in differential modes of nutrition (photoautotrophic vs. photoheterotrophic) in relation to secondary metabolites production in S. chirayita.

Genome-Based Selection and Characterization of Fusarium circinatum-Specific Sequences

Fusarium circinatum is an important pathogen of pine trees and its management in the commercial forestry environment relies largely on early detection, particularly in seedling nurseries. The fact that the entire genome of this pathogen is available opens new avenues for the development of diagnostic tools for this fungus. In this study we identified open reading frames (ORFs) unique to F. circinatum and determined that they were specific to the pathogen. The ORF identification process involved bioinformatics-based screening of all the putative F. circinatum ORFs against public databases. This was followed by functional characterization of ORFs found to be unique to F. circinatum. We used PCR- and hybridization-based approaches to confirm the presence of selected unique genes in different strains of F. circinatum and their absence from other Fusarium species for which genome sequence data are not yet available. These included species that are closely related to F. circinatum as well as those that are commonly encountered in the forestry environment. Thirty-six ORFs were identified as potentially unique to F. circinatum. Nineteen of these encode proteins with known domains while the other 17 encode proteins of unknown function. The results of our PCR analyses and hybridization assays showed that three of the selected genes were present in all of the strains of F. circinatum tested and absent from the other Fusarium species screened. These data thus indicate that the selected genes are common and unique to F. circinatum. These genes thus could be good candidates for use in rapid, in-the-field diagnostic assays specific to F. circinatum. Our study further demonstrates how genome sequence information can be mined for the identification of new diagnostic markers for the detection of plant pathogens.

Gene Duplicability of Core Genes Is Highly Consistent across All Angiosperms

Gene duplication is an important mechanism for adding to genomic novelty. Hence, which genes undergo duplication and are preserved following duplication is an important question. It has been observed that gene duplicability, or the ability of genes to be retained following duplication, is a nonrandom process, with certain genes being more amenable to survive duplication events than others. Primarily, gene essentiality and the type of duplication (small-scale versus large-scale) have been shown in different species to influence the (long-term) survival of novel genes. However, an overarching view of "gene duplicability" is lacking, mainly due to the fact that previous studies usually focused on individual species and did not account for the influence of genomic context and the time of duplication. Here, we present a large-scale study in which we investigated duplicate retention for 9178 gene families shared between 37 flowering plant species, referred to as angiosperm core gene families. For most gene families, we observe a strikingly consistent pattern of gene duplicability across species, with gene families being either primarily single-copy or multicopy in all species. An intermediate class contains gene families that are often retained in duplicate for periods extending to tens of millions of years after whole-genome duplication, but ultimately appear to be largely restored to singleton status, suggesting that these genes may be dosage balance sensitive. The distinction between single-copy and multicopy gene families is reflected in their functional annotation, with single-copy genes being mainly involved in the maintenance of genome stability and organelle function and multicopy genes in signaling, transport, and metabolism. The intermediate class was overrepresented in regulatory genes, further suggesting that these represent putative dosage-balance-sensitive genes.

Single-copy gene based 50 K SNP chip for genetic studies and molecular breeding in rice

Single nucleotide polymorphism (SNP) is the most abundant DNA sequence variation present in plant genomes. Here, we report the design and validation of a unique genic-SNP genotyping chip for genetic and evolutionary studies as well as molecular breeding applications in rice. The chip incorporates 50,051 SNPs from 18,980 different genes spanning 12 rice chromosomes, including 3,710 single-copy (SC) genes conserved between wheat and rice, 14,959 SC genes unique to rice, 194 agronomically important cloned rice genes and 117 multi-copy rice genes. Assays with this chip showed high success rate and reproducibility because of the SC gene based array with no sequence redundancy and cross-hybridisation problems. The usefulness of the chip in genetic diversity and phylogenetic studies of cultivated and wild rice germplasm was demonstrated. Furthermore, its efficacy was validated for analysing background recovery in improved mega rice varieties with submergence tolerance developed through marker-assisted backcross breeding.

Explosive tandem and segmental duplications of multigenic families in Eucalyptus grandis

Plant organisms contain a large number of genes belonging to numerous multigenic families whose evolution size reflects some functional constraints. Sequences from eight multigenic families, involved in biotic and abiotic responses, have been analyzed in Eucalyptus grandis and compared with Arabidopsis thaliana. Two transcription factor families APETALA 2 (AP2)/ethylene responsive factor and GRAS, two auxin transporter families PIN-FORMED and AUX/LAX, two oxidoreductase families (ascorbate peroxidases [APx] and Class III peroxidases [CIII Prx]), and two families of protective molecules late embryogenesis abundant (LEA) and DNAj were annotated in expert and exhaustive manner. Many recent tandem duplications leading to the emergence of species-specific gene clusters and the explosion of the gene numbers have been observed for the AP2, GRAS, LEA, PIN, and CIII Prx in E. grandis, while the APx, the AUX/LAX and DNAj are conserved between species. Although no direct evidence has yet demonstrated the roles of these recent duplicated genes observed in E. grandis, this could indicate their putative implications in the morphological and physiological characteristics of E. grandis, and be the key factor for the survival of this nondormant species. Global analysis of key families would be a good criterion to evaluate the capabilities of some organisms to adapt to environmental variations.

Genome-wide analysis reveals diverged patterns of codon bias, gene expression, and rates of sequence evolution in picea gene families

The recent sequencing of several gymnosperm genomes has greatly facilitated studying the evolution of their genes and gene families. In this study, we examine the evidence for expression-mediated selection in the first two fully sequenced representatives of the gymnosperm plant clade (Picea abies and Picea glauca). We use genome-wide estimates of gene expression (>50,000 expressed genes) to study the relationship between gene expression, codon bias, rates of sequence divergence, protein length, and gene duplication. We found that gene expression is correlated with rates of sequence divergence and codon bias, suggesting that natural selection is acting on Picea protein-coding genes for translational efficiency. Gene expression, rates of sequence divergence, and codon bias are correlated with the size of gene families, with large multicopy gene families having, on average, a lower expression level and breadth, lower codon bias, and higher rates of sequence divergence than single-copy gene families. Tissue-specific patterns of gene expression were more common in large gene families with large gene expression divergence than in single-copy families. Recent family expansions combined with large gene expression variation in paralogs and increased rates of sequence evolution suggest that some Picea gene families are rapidly evolving to cope with biotic and abiotic stress. Our study highlights the importance of gene expression and natural selection in shaping the evolution of protein-coding genes in Picea species, and sets the ground for further studies investigating the evolution of individual gene families in gymnosperms.

An easy-to-use primer design tool to address paralogous loci and T-DNA insertion sites in the genome of Arabidopsis thaliana

BACKGROUND

More than 90% of the Arabidopsis thaliana genes are members of multigene families. DNA sequence similarities present in such related genes can cause trouble, e.g. when molecularly analysing mutant alleles of these genes. Also, flanking-sequence-tag (FST) based predictions of T-DNA insertion positions are often located within paralogous regions of the genome. In such cases, the prediction of the correct insertion site must include careful sequence analyses on the one hand and a paralog specific primer design for experimental confirmation of the prediction on the other hand.

RESULTS

GABI-Kat is a large A. thaliana insertion line resource, which uses in-house confirmation to provide highly reliable access to T-DNA insertion alleles. To offer trustworthy mutant alleles of paralogous loci, we considered multiple insertion site predictions for single FSTs and implemented this 1-to-N relation in our database. The resulting paralogous predictions were addressed experimentally and the correct insertion locus was identified in most cases, including cases in which there were multiple predictions with identical prediction scores. A newly developed primer design tool that takes paralogous regions into account was developed to streamline the confirmation process for paralogs. The tool is suitable for all parts of the genome and is freely available at the GABI-Kat website. Although the tool was initially designed for the analysis of T-DNA insertion mutants, it can be used for any experiment that requires locus-specific primers for the A. thaliana genome. It is easy to use and also able to design amplimers with two genome-specific primers as required for genotyping segregating families of insertion mutants when looking for homozygous offspring.

CONCLUSIONS

The paralog-aware confirmation process significantly improved the reliability of the insertion site assignment when paralogous regions of the genome were affected. An automatic online primer design tool that incorporates experience from the in-house confirmation of T-DNA insertion lines has been made available. It provides easy access to primers for the analysis of T-DNA insertion alleles, but it is also beneficial for other applications as well.

Arabidopsis thaliana organellar DNA polymerase IB mutants exhibit reduced mtDNA levels with a decrease in mitochondrial area density

Plant organelle genomes are complex and the mechanisms for their replication and maintenance remain unclear. Arabidopsis thaliana has two DNA polymerase genes, DNA polymerase IA (polIA) and polIB, that are dual targeted to mitochondria and chloroplasts and are differentially expressed in primary plant tissues. PolIB gene expression occurs at higher levels in tissues not primary for photosynthesis. Arabidopsis T-DNA polIB mutants have a 30% reduction in relative mitochondrial DNA (mtDNA) levels, but also exhibit a 70% increase in polIA gene expression. The polIB mutant shows an increase in mitochondrial numbers but a significant decrease in mitochondrial area density within the hypocotyl epidermis, shoot apex and root tips. Chloroplast numbers are not significantly different in mesophyll protoplasts. These mutants do not have a significant difference in total dark mitorespiration levels but exhibit a difference in light respiration levels and photosynthesis capacity. Organelle-encoded genes for components of respiration and photosynthesis are upregulated in polIB mutants. The mutants exhibited slow growth in conjunction with a decreased rate of cell expansion and other secondary phenotypic effects. Evidence suggests that early plastid development and DNA levels are directly affected by a polIB mutation but are resolved to wild-type levels over time. However, mitochondria numbers and DNA levels never reach wild-type levels in the polIB mutant. We propose that both polIA and polIB are required for mtDNA replication. The results suggest that polIB mutants undergo an adjustment in cell homeostasis, enabling them to maintain functional mitochondria at the cost of normal cell expansion and plant growth.

GABI-DUPLO: a collection of double mutants to overcome genetic redundancy in Arabidopsis thaliana

Owing to duplication events in its progenitor, more than 90% of the genes in the Arabidopsis thaliana genome are members of multigene families. A set of 2108 gene families, each consisting of precisely two unlinked paralogous genes, was identified in the nuclear genome of A. thaliana on the basis of sequence similarity. A systematic method for the creation of double knock-out lines for such gene pairs, designated as DUPLO lines, was established and 200 lines are now publicly available. Their initial phenotypic characterisation led to the identification of seven lines with defects that emerge only in the adult stage. A further six lines display seedling lethality and 23 lines were lethal before germination. Another 14 lines are known to show phenotypes under non-standard conditions or at the molecular level. Knock-out of gene pairs with very similar coding sequences or expression profiles is more likely to produce a mutant phenotype than inactivation of gene pairs with dissimilar profiles or sequences. High coding sequence similarity and highly similar expression profiles are only weakly correlated, implying that promoter and coding regions of these gene pairs display different degrees of diversification.

Computational predictions provide insights into the biology of TAL effector target sites

Transcription activator-like (TAL) effectors are injected into host plant cells by Xanthomonas bacteria to function as transcriptional activators for the benefit of the pathogen. The DNA binding domain of TAL effectors is composed of conserved amino acid repeat structures containing repeat-variable diresidues (RVDs) that determine DNA binding specificity. In this paper, we present TALgetter, a new approach for predicting TAL effector target sites based on a statistical model. In contrast to previous approaches, the parameters of TALgetter are estimated from training data computationally. We demonstrate that TALgetter successfully predicts known TAL effector target sites and often yields a greater number of predictions that are consistent with up-regulation in gene expression microarrays than an existing approach, Target Finder of the TALE-NT suite. We study the binding specificities estimated by TALgetter and approve that different RVDs are differently important for transcriptional activation. In subsequent studies, the predictions of TALgetter indicate a previously unreported positional preference of TAL effector target sites relative to the transcription start site. In addition, several TAL effectors are predicted to bind to the TATA-box, which might constitute one general mode of transcriptional activation by TAL effectors. Scrutinizing the predicted target sites of TALgetter, we propose several novel TAL effector virulence targets in rice and sweet orange. TAL-mediated induction of the candidates is supported by gene expression microarrays. Validity of these targets is also supported by functional analogy to known TAL effector targets, by an over-representation of TAL effector targets with similar function, or by a biological function related to pathogen infection. Hence, these predicted TAL effector virulence targets are promising candidates for studying the virulence function of TAL effectors. TALgetter is implemented as part of the open-source Java library Jstacs, and is freely available as a web-application and a command line program.

Convergent gene loss following gene and genome duplications creates single-copy families in flowering plants

The importance of gene gain through duplication has long been appreciated. In contrast, the importance of gene loss has only recently attracted attention. Indeed, studies in organisms ranging from plants to worms and humans suggest that duplication of some genes might be better tolerated than that of others. Here we have undertaken a large-scale study to investigate the existence of duplication-resistant genes in the sequenced genomes of 20 flowering plants. We demonstrate that there is a large set of genes that is convergently restored to single-copy status following multiple genome-wide and smaller scale duplication events. We rule out the possibility that such a pattern could be explained by random gene loss only and therefore propose that there is selection pressure to preserve such genes as singletons. This is further substantiated by the observation that angiosperm single-copy genes do not comprise a random fraction of the genome, but instead are often involved in essential housekeeping functions that are highly conserved across all eukaryotes. Furthermore, single-copy genes are generally expressed more highly and in more tissues than non-single-copy genes, and they exhibit higher sequence conservation. Finally, we propose different hypotheses to explain their resistance against duplication.

Ascorbate peroxidase-related (APx-R) is not a duplicable gene

Phylogenetic, genomic and functional analyses have allowed the identification of a new class of putative heme peroxidases, so called APx-R (APx-Related). These new class, mainly present in the green lineage (including green algae and land plants), can also be detected in other unicellular chloroplastic organisms. Except for recent polyploid organisms, only single-copy of APx-R gene was detected in each genome, suggesting that the majority of the APx-R extra-copies were lost after chromosomal or segmental duplications. In a similar way, most APx-R co-expressed genes in Arabidopsis genome do not have conserved extra-copies after chromosomal duplications and are predicted to be localized in organelles, as are the APx-R. The member of this gene network can be considered as unique gene, well conserved through the evolution due to a strong negative selection pressure and a low evolution rate.

Genome-wide investigation reveals high evolutionary rates in annual model plants

BACKGROUND

Rates of molecular evolution vary widely among species. While significant deviations from molecular clock have been found in many taxa, effects of life histories on molecular evolution are not fully understood. In plants, annual/perennial life history traits have long been suspected to influence the evolutionary rates at the molecular level. To date, however, the number of genes investigated on this subject is limited and the conclusions are mixed. To evaluate the possible heterogeneity in evolutionary rates between annual and perennial plants at the genomic level, we investigated 85 nuclear housekeeping genes, 10 non-housekeeping families, and 34 chloroplast genes using the genomic data from model plants including Arabidopsis thaliana and Medicago truncatula for annuals and grape (Vitis vinifera) and popular (Populus trichocarpa) for perennials.

RESULTS

According to the cross-comparisons among the four species, 74-82% of the nuclear genes and 71-97% of the chloroplast genes suggested higher rates of molecular evolution in the two annuals than those in the two perennials. The significant heterogeneity in evolutionary rate between annuals and perennials was consistently found both in nonsynonymous sites and synonymous sites. While a linear correlation of evolutionary rates in orthologous genes between species was observed in nonsynonymous sites, the correlation was weak or invisible in synonymous sites. This tendency was clearer in nuclear genes than in chloroplast genes, in which the overall evolutionary rate was small. The slope of the regression line was consistently lower than unity, further confirming the higher evolutionary rate in annuals at the genomic level.

CONCLUSIONS

The higher evolutionary rate in annuals than in perennials appears to be a universal phenomenon both in nuclear and chloroplast genomes in the four dicot model plants we investigated. Therefore, such heterogeneity in evolutionary rate should result from factors that have genome-wide influence, most likely those associated with annual/perennial life history. Although we acknowledge current limitations of this kind of study, mainly due to a small sample size available and a distant taxonomic relationship of the model organisms, our results indicate that the genome-wide survey is a promising approach toward further understanding of the mechanism determining the molecular evolutionary rate at the genomic level.

Natural diversity in flowering responses of Arabidopsis thaliana caused by variation in a tandem gene array

Tandemly arrayed genes that belong to gene families characterize genomes of many organisms. Gene duplication and subsequent relaxation of selection can lead to the establishment of paralogous cluster members that may evolve along different trajectories. Here, we report on the structural variation in MADS AFFECTING FLOWERING 2 (MAF2) gene, one member of the tandemly duplicated cluster of MADS-box-containing transcription factors in Arabidopsis thaliana. The altered gene structure at the MAF2 locus is present as a moderate-frequency polymorphism in Arabidopsis and leads to the extensive diversity in transcript patterns due to alternative splicing. Rearrangements at the MAF2 locus are associated with an early flowering phenotype in BC(5) lines. The lack of suppression of flowering time in a MAF2-insertion line expressing the MAF2-specific artificial miRNA suggests that these MAF2 variants are behaving as loss-of-function alleles. The variation in gene architecture is also associated with segregation distortion, which may have facilitated the spread and the establishment of the corresponding alleles throughout the Eurasian range of the A. thaliana population.

TC-motifs at the TATA-box expected position in plant genes: a novel class of motifs involved in the transcription regulation

BACKGROUND

The TATA-box and TATA-variants are regulatory elements involved in the formation of a transcription initiation complex. Both have been conserved throughout evolution in a restricted region close to the Transcription Start Site (TSS). However, less than half of the genes in model organisms studied so far have been found to contain either one of these elements. Indeed different core-promoter elements are involved in the recruitment of the TATA-box-binding protein. Here we assessed the possibility of identifying novel functional motifs in plant genes, sharing the TATA-box topological constraints.

RESULTS

We developed an ab-initio approach considering the preferential location of motifs relative to the TSS. We identified motifs observed at the TATA-box expected location and conserved in both Arabidopsis thaliana and Oryza sativa promoters. We identified TC-elements within non-TA-rich promoters 30 bases upstream of the TSS. As with the TATA-box and TATA-variant sequences, it was possible to construct a unique distance graph with the TC-element sequences. The structural and functional features of TC-element-containing genes were distinct from those of TATA-box- or TATA-variant-containing genes. Arabidopsis thaliana transcriptome analysis revealed that TATA-box-containing genes were generally those showing relatively high levels of expression and that TC-element-containing genes were generally those expressed in specific conditions.

CONCLUSIONS

Our observations suggest that the TC-elements might constitute a class of novel regulatory elements participating towards the complex modulation of gene expression in plants.

FIDEL-a retrovirus-like retrotransposon and its distinct evolutionary histories in the A- and B-genome components of cultivated peanut

In this paper, we describe a Ty3-gypsy retrotransposon from allotetraploid peanut (Arachis hypogaea) and its putative diploid ancestors Arachis duranensis (A-genome) and Arachis ipaënsis (B-genome). The consensus sequence is 11,223 bp. The element, named FIDEL (Fairly long Inter-Dispersed Euchromatic LTR retrotransposon), is more frequent in the A- than in the B-genome, with copy numbers of about 3,000 (±950, A. duranensis), 820 (±480, A. ipaënsis), and 3,900 (±1,500, A. hypogaea) per haploid genome. Phylogenetic analysis of reverse transcriptase sequences showed distinct evolution of FIDEL in the ancestor species. Fluorescent in situ hybridization revealed disperse distribution in euchromatin and absence from centromeres, telomeric regions, and the nucleolar organizer region. Using paired sequences from bacterial artificial chromosomes, we showed that elements appear less likely to insert near conserved ancestral genes than near the fast evolving disease resistance gene homologs. Within the Ty3-gypsy elements, FIDEL is most closely related with the Athila/Calypso group of retrovirus-like retrotransposons. Putative transmembrane domains were identified, supporting the presence of a vestigial envelope gene. The results emphasize the importance of FIDEL in the evolution and divergence of different Arachis genomes and also may serve as an example of the role of retrotransposons in the evolution of legume genomes in general.