Thellungilla halophila is more adaptive to salinity than Arabidopsis thaliana at stages of seed germination and seedling establishment

Thellungiella halophila is a salt tolerant relative of Arabidopsis thaliana with high genetic and morphological similarity. In the present study, effects of salinity on germination and seedling growth of T. halophila and A. thaliana were compared. The present results showed that the salinity inhibited seed germination in both species. Unexpectedly, percentages of seed germination in A. thaliana were higher than T. halophila in a range of 0?C200?mM NaCl. Seeds of both species could not germinate when the concentration of NaCl was over 200?mM. However, when compared with A. thaliana, seeds of T. halophila did not suffer ion toxicity, as evidenced by the higher final germination rate after ungerminated seeds pretreated with NaCl were transferred to distilled water. Seedlings of T. halophila were more salt tolerant than those of A. thaliana, e.g., seedlings of T. halophila had better plant growth (root length, fresh and dry mass), higher chlorophyll content, less MDA content and higher proline content and K⁺/Na⁺ ratio under salinity. These results indicate that T. halophila is more salt tolerant than A. thaliana during both seed germination and seedling stages and explain why A. thaliana is excluded from saline locations and T. halophila can survive in saline soils.     

Thellungiella halophila, a salt-tolerant relative of Arabidopsis thaliana, possesses effective mechanisms to discriminate between potassium and sodium

Thellungiella halophila is a salt‐tolerant close relative of Arabidopsis thaliana. Significant mRNA similarity was confirmed by hybridization of T. halophila mRNA with the A. thaliana GeneChip ATH1. To establish a platform for future molecular comparison of the two species several physiological mechanisms, which may confer high salt tolerance to T. halophila, were investigated. Determination of ion content in shoots and roots of A. thaliana and T. halophila indicated different strategies of ion uptake and translocation from root to shoot in the two species. During salt stress T. halophila accumulated less sodium than A. thaliana. Tissue concentrations of sodium and potassium showed negative correlation in A. thaliana but not in T. halophila. Electrophysiological experiments proved high potassium/sodium selectivity of root plasma membrane channels in T. halophila. In particular, voltage‐independent currents were more selective for potassium in T. halophila than in A. thaliana. Single cell sampling of T. halophila leaves during salt exposure revealed increased concentrations of sodium and decreased concentrations of potassium in epidermal cells suggesting that this cell type could function to ensure storage of sodium and exchange of potassium with the rest of leaf. Application of salt resulted in a sharp drop of transpiration in A. thaliana. By contrast, transpiration in T. halophila responded more slowly and was only slightly inhibited by salt treatment, thus maintaining high water uptake and ion transport.     

The ARC1 E3 Ligase Promotes Two Different Self-Pollen Avoidance Traits in Arabidopsis

Flowering plants have evolved various strategies for avoiding self-pollen to drive genetic diversity. These strategies include spatially separated sexual organs (herkogamy), timing differences between male pollen release and female pistil receptivity (dichogamy), and self-pollen rejection. Within the Brassicaceae, these outcrossing systems are the evolutionary default state, and many species display these traits, including Arabidopsis lyrata. In contrast to A. lyrata, closely related Arabidopsis thaliana has lost these self-pollen traits and thus represents an excellent system to test genes for reconstructing these evolutionary traits. We previously demonstrated that the ARC1 E3 ligase is required for self-incompatibility in two diverse Brassicaceae species, Brassica napus and A. lyrata, and is frequently deleted in self-compatible species, including A. thaliana. In this study, we examined ARC1’s requirement for reconstituting self-incompatibility in A. thaliana and uncovered an important role for ARC1 in promoting a strong and stable pollen rejection response when expressed with two other A. lyrata self-incompatibility factors. Furthermore, we discovered that ARC1 promoted an approach herkogamous phenotype in A. thaliana flowers. Thus, ARC1’s expression resulted in two different A. lyrata traits for self-pollen avoidance and highlights the key role that ARC1 plays in the evolution and retention of outcrossing systems.     

Divergent Evolutionary and Expression Patterns between Lineage Specific New Duplicate Genes and Their Parental Paralogs in Arabidopsis thaliana

Gene duplication is an important mechanism for the origination of functional novelties in organisms. We performed a comparative genome analysis to systematically estimate recent lineage specific gene duplication events in Arabidopsis thaliana and further investigate whether and how these new duplicate genes (NDGs) play a functional role in the evolution and adaption of A. thaliana. We accomplished this using syntenic relationship among four closely related species, A. thaliana, A. lyrata, Capsella rubella and Brassica rapa. We identified 100 NDGs, showing clear origination patterns, whose parental genes are located in syntenic regions and/or have clear orthologs in at least one of three outgroup species. All 100 NDGs were transcribed and under functional constraints, while 24% of the NDGs have differential expression patterns compared to their parental genes. We explored the underlying evolutionary forces of these paralogous pairs through conducting neutrality tests with sequence divergence and polymorphism data. Evolution of about 15% of NDGs appeared to be driven by natural selection. Moreover, we found that 3 NDGs not only altered their expression patterns when compared with parental genes, but also evolved under positive selection. We investigated the underlying mechanisms driving the differential expression of NDGs and their parents, and found a number of NDGs had different cis-elements and methylation patterns from their parental genes. Overall, we demonstrated that NDGs acquired divergent cis-elements and methylation patterns and may experience sub-functionalization or neo-functionalization influencing the evolution and adaption of A. thaliana.     

A Population Genomics Study of the Arabidopsis Core Cell Cycle Genes Shows the Signature of Natural Selection

Large-scale comparison of sequence polymorphism and divergence at numerous genomic loci within and between closely related species can reveal signatures of natural selection. Here, we present a population genomics study based on direct sequencing of 61 mitotic cell cycle genes from 30 Arabidopsis thaliana accessions and comparison of the resulting data to the close relative Arabidopsis lyrata. We found that the Arabidopsis core cell cycle (CCC) machinery is not highly constrained but is subject to different modes of selection. We found patterns of purifying selection for the cyclin-dependent kinase (CDK), CDK subunit, retinoblastoma, and WEE1 gene families. Other CCC gene families often showed a mix of one or two constrained genes and relaxed purifying selection on the other genes. We found several large effect mutations in CDKB1;2 that segregate in the species. We found a strong signature of adaptive protein evolution in the Kip-related protein KRP6 and departures from equilibrium at CDKD;1 and CYCA3;3 consistent with the operation of selection in these gene regions. Our data suggest that within Arabidopsis, the genetic robustness of cell cycle–related processes is more due to functional redundancy than high selective constraint.     

Comparative salt tolerance analysis between Arabidopsis thaliana and Thellungiella halophila, with special emphasis on K(+)/Na(+) selectivity and proline accumulation

The eco-physiology of salt tolerance, with an emphasis on K⁺ nutrition and proline accumulation, was investigated in the halophyte Thellungiella halophila and in both wild type and eskimo-1 mutant of the glycophyte Arabidopsis thaliana, which differ in their proline accumulation capacity. Plants cultivated in inert sand were challenged for 3 weeks with up to 500 mM NaCl. Low salinity significantly decreased A. thaliana growth, whereas growth restriction was significant only at salt concentrations equal to or exceeding 300 mM NaCl in T. halophila. Na⁺ content generally increased with the amount of salt added in the culture medium in both species, but T. halophila showed an ability to control Na⁺ accumulation in shoots. The analysis of the relationship between water and Na⁺ contents suggested an apoplastic sodium accumulation in both species; this trait was more pronounced in A. thaliana than in T. halophila. The better NaCl tolerance in the latter was associated with a better K⁺ supply, resulting in higher K⁺/Na⁺ ratios. It was also noteworthy that, despite highly accumulating proline, the A. thaliana eskimo-1 mutant was the most salt-sensitive species. Taken together, our findings indicate that salt tolerance may be partly linked to the plants’ ability to control Na⁺ influx and to ensure appropriate K⁺ nutrition, but is not linked to proline accumulation.     

Genetic transformation of Arabidopsis lyrata: specific expression of the green fluorescent protein (GFP) in pistil tissues

Arabidopsis thaliana has become widely used as a model system for plant biology. Recent phylogenetic studies led to a severe revision of the systematic relationships across species of the Brassicaceae family. This provided an opportunity to examine close relatives of A. thaliana and to study the function and molecular evolution of genes that play roles in ecology and speciation. In this context, developing tools to genetically transform “non-model plants” appears as a major issue to ascertain gene function. Here, we report a method to transform A. lyrata, one of the closest relatives of A. thaliana.     

The shrunken genome of Arabidopsis thaliana

This paper examines macro and micro-level patterns of genome size evolution in the Brassicaceae. A phylogeny of 25 relatives of Arabidopsis thaliana was reconstructed using four molecular markers under both parsimony and Bayesian methods. Reconstruction of genome size (C value) evolution as a discrete character and as a continuous character was also performed. In addition, size dynamics in small chromosomal regions were assessed by comparing genomic clones generated for Arabidopsis lyrata and for Boechera stricta to the fully sequenced genome of A. thaliana. The results reveal a sevenfold variation in genome size among the taxa investigated and that the small genome size of A. thaliana is derived. Our results also indicate that the genome is free to increase or decrease in size across these evolutionary lineages without a directional bias. These changes are accomplished by insertions and deletions at both large and small-scales occurring mostly in intergenic regions, with repetitive sequences and transposable elements implicated in genome size increases. The focus upon taxa relatively closely related to the model organism A. thaliana, and the combination of complementary approaches, allows for unique insights into the processes driving genome size changes.     

Codon bias of the gene for chloroplast glycerol-3-phosphate acyltransferase in Camellia sinensis (L.) O. Kuntze

The analysis on codon usage bias of GPAT gene of Camellia sinensis (L.) O. Kuntze may provide a basis for understanding the evolution relationship of C. sinensis and for selecting appropriate host expression systems to improve the expression of target genes. In the present study, the coding sequence of CsGPAT was analyzed with CodonW, CHIPS and CUSP programs, and compared with the genome of C. sinensis and GPAT genes of other 9 plant species. Our results showed that the cluster tree based on CDs could reveal the evolutional relations among the 10 plant species, whereas the cluster tree based on relative synonymous codon usage (RSCU) could not. There were 31 codons showing distinct usage differences between CsGPAT and genome of Escherichia coli, 21 between CsGPAT and yeast, but 13 between CsGPAT and Arabidopsis thaliana. But there were slightly fewer differences in codon usage between CsGPAT and A. thaliana. Therefore, the A. thaliana expression system may be more suitable for the expression of CsGPAT. These results may improve our understanding of the codon usage bias and functional studies of CsGPAT.     

Arabidopsis lyrata Small RNAs: Transient MIRNA and Small Interfering RNA Loci within the Arabidopsis Genus

Twenty-one-nucleotide microRNAs (miRNAs) and 24-nucleotide Pol IV-dependent small interfering RNAs (p4-siRNAs) are the most abundant types of small RNAs in angiosperms. Some miRNAs are well conserved among different plant lineages, whereas others are less conserved, and it is not clear whether less-conserved miRNAs have the same functionality as the well conserved ones. p4-siRNAs are broadly produced in the Arabidopsis genome, sometimes from active hot spot loci, but it is unknown whether individual p4-siRNA hot spots are retained as hot spots between plant species. In this study, we compare small RNAs in two closely related species (Arabidopsis thaliana and Arabidopsis lyrata) and find that less-conserved miRNAs have high rates of divergence in MIRNA hairpin structures, mature miRNA sequences, and target-complementary sites in the other species. The fidelity of miRNA biogenesis from many less-conserved MIRNA hairpins frequently deteriorates in the sister species relative to the species of first discovery. We also observe that p4-siRNA occupied loci have a slight tendency to be retained as p4-siRNA loci between species, but the most active A. lyrata p4-siRNA hot spots are generally not syntenic to the most active p4-siRNA hot spots of A. thaliana. Altogether, our findings indicate that many MIRNAs and most p4-siRNA hot spots are rapidly changing and evolutionarily transient within the Arabidopsis genus.     

Genetic evolution and codon usage analysis of NKX-2.5 gene governing heart development in some mammals

NKX-2.5 gene is responsible for cardiac development and its targeted disruption apprehends cardiac development at the linear heart tube stage. Bioinformatic analysis was employed to investigate the codon usage pattern and dN/dS of mammalian NKX-2.5 gene. The relative synonymous codon usage analysis revealed variation in codon usage and two synonymous codons namely ATA (Ile) and GTA (Val) were absent in NKX-2.5 gene across selected mammalian species suggesting that these two codons were possibly selected against during evolution. Parity rule 2 analysis of two and four fold amino acids showed CT bias whereas six-fold amino acids revealed GA bias. Neutrality analysis suggests that selection played a prominent role while mutation had a minor role. The dN/dS analysis suggests synonymous substitution played a significant role and it negatively correlated with p-distance of the gene. Purifying natural selection played a dominant role in the genetic evolution of NKX-2.5 gene in mammals.     

Analysis of synonymous codon usage patterns of HPRT1 gene across twelve mammalian species

Genetic changes in Hypoxanthine guanine phosphoribosyltransferace (HPRT1) gene can alter the expression of the dopamine neurotransmitter leads to abnormal neuron function, a disease called Lesch-Nyhan syndrome (LNS). Although different studies were conducted on LNS, information on codon usage bias (CUB) of HPRT1 gene is limited. The present study examines the genetic determinants of CUB in HPRT1 gene using twelve mammalian species. In the coding sequence of HPRT1 genes, A/T ending codons was most frequently used. A higher ENC value was observed indicating lower HPRT1 gene expression in the selected mammalian species. Correlation analysis indicates that compositional constraints under mutation pressure can involve in CUB of HPRT1 genes among the selected mammalian species. Relative synonymous codon usage (RSCU) value revealed that the codons such as ACT, AGG, ATT and AGC were over-represented in each of the mammalian species. Result from the analysis of the RSCU indicates that compositional constraint is a key driver for the variation in codon usage. Ratio of nonsynonymous (dN) and synonymous (dS) substitution further suggested that purifying selection occurs among the HPRT1 gene of studied mammals to maintain its protein function under the process of evolution. Our findings report an insight into the codon usage patterns of HPRT1 gene and will be useful for LNS management.     

Expression of GFP-mTalin reveals an actin-related role for the Arabidopsis Class II formin AtFH12

Formins (FH2 proteins) are implicated in F-actin nucleation and other aspects of cytoskeletal organization. Plants possess two formin clades, relatively well-described Class I formins and so far poorly characterized Class II formins. Comparison of Class II formin genes of two Arabidopsis species, A. thaliana and A. lyrata, indicates dynamic evolution within the Class II formin clade. Disruption of an outlier A. thaliana Class II formin gene, AtFH12 (At1g42980), whose expression is induced by NaCl, produced only negligible phenotypic effects under a variety of conditions, including salt stress, suggesting functional redundancy among Class II formins. However, the same mutation massively aggravated toxic effects of the expression of a fluorescent actin marker, GFP-tagged mouse talin (GFP-mTalin), known to interfere with normal actin dynamics. Abnormal actin structures were observed in atfh12 mutants expressing GFP-mTalin as compared to wild type. This not only demonstrates an actin-associated function for AtFH12, but also documents the feasibility of using the heterologous actin marker to “stress-test” the actin cytoskeleton in phenotyping “weak“ actin related mutant alleles.     

Natural variation in Arabidopsis lyrata vernalization requirement conferred by a FRIGIDA indel polymorphism

Species share homologous genes to a large extent, but it isnot yet known to what degree the same loci have been targetsfor natural selection in different species. Natural variationin flowering time is determined to a large degree by 2 genes,FLOWERING LOCUS C and FRIGIDA, in Arabidopsis thaliana. Here,we examine whether FRIGIDA has a role in differences in floweringtime between and within natural populations of Arabidopsis lyrata,a close outcrossing perennial relative of A. thaliana. We found2 FRIGIDA sequence variants producing potentially functionalproteins but with a length difference of 14 amino acids. Thesevariants conferred a 15-day difference in flowering time inan association experiment in 2 Scandinavian populations. Thedifference in flowering time between alleles was confirmed withtransformation to A. thaliana. Because the north European late-floweringpopulations harbor both late- and early sequence variants atintermediate frequencies and the late-flowering variant is mostfrequent in the southern early flowering European population,other genetic factors must be responsible for the floweringtime differences between the populations. The length polymorphismoccurs at high frequencies also in several North American populations.The occurrence of functional variants at intermediate frequenciesin several populations suggests that the variation may be maintainedby balancing selection. This is in contrast to A. thaliana,where independent loss-of-function mutations at the FRIGIDAgene are responsible for differences between populations andlocal adaptation.     

Gaining Insights into the Codon Usage Patterns of TP53 Gene across Eight Mammalian Species

TP53 gene is known as the “guardian of the genome” as it plays a vital role in regulating cell cycle, cell proliferation, DNA damage repair, initiation of programmed cell death and suppressing tumor growth. Non uniform usage of synonymous codons for a specific amino acid during translation of protein known as codon usage bias (CUB) is a unique property of the genome and shows species specific deviation. Analysis of codon usage bias with compositional dynamics of coding sequences has contributed to the better understanding of the molecular mechanism and the evolution of a particular gene. In this study, the complete nucleotide coding sequences of TP53 gene from eight different mammalian species were used for CUB analysis. Our results showed that the codon usage patterns in TP53 gene across different mammalian species has been influenced by GC bias particularly GC₃ and a moderate bias exists in the codon usage of TP53 gene. Moreover, we observed that nature has highly favored the most over represented codon CTG for leucine amino acid but selected against the ATA codon for isoleucine in TP53 gene across all mammalian species during the course of evolution.     

Non-uniqueness of factors constraint on the codon usage in Bombyx mori

Background

The analysis of codon usage is a good way to understand the genetic and evolutionary characteristics of an organism. However, there are only a few reports related with the codon usage of the domesticated silkworm, Bombyx mori (B. mori). Hence, the codon usage of B. mori was analyzed here to reveal the constraint factors and it could be helpful to improve the bioreactor based on B. mori.

Results

A total of 1,097 annotated mRNA sequences from B. mori were analyzed, revealing there is only a weak codon bias. It also shows that the gene expression level is related to the GC content, and the amino acids with higher general average hydropathicity (GRAVY) and aromaticity (Aromo). And the genes on the primary axis are strongly positively correlated with the GC content, and GC3s. Meanwhile, the effective number of codons (ENc) is strongly correlated with codon adaptation index (CAI), gene length, and Aromo values. However, the ENc values are correlated with the second axis, which indicates that the codon usage in B. mori is affected by not only mutation pressure and natural selection, but also nucleotide composition and the gene expression level. It is also associated with Aromo values, and gene length. Additionally, B. mori has a greater relative discrepancy in codon preferences with Drosophila melanogaster (D. melanogaster) or Saccharomyces cerevisiae (S. cerevisiae) than with Arabidopsis thaliana (A. thaliana), Escherichia coli (E. coli), or Caenorhabditis elegans (C. elegans).

Conclusions

The codon usage bias in B. mori is relatively weak, and many influence factors are found here, such as nucleotide composition, mutation pressure, natural selection, and expression level. Additionally, it is also associated with Aromo values, and gene length. Among them, natural selection might play a major role. Moreover, the “optimal codons” of B. mori are all encoded by G and C, which provides useful information for enhancing the gene expression in B. mori through codon optimization.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1596-z) contains supplementary material, which is available to authorized users.     

Codon usage <Emphasis Type="BoldItalic">vis-a-vis</Emphasis> start and stop codon context analysis of three dicot species

To understand the variation in genomic composition and its effect on codon usage, we performed the comparative analysis of codon usage and nucleotide usage in the genes of three dicots, Glycine max, Arabidopsis thaliana and Medicago truncatula. The dicot genes were found to be A/T rich and have predominantly A-ending and/or T-ending codons. GC3s directly mimic the usage pattern of global GC content. Relative synonymous codon usage analysis suggests that the high usage frequency of A/T over G/C mononucleotide containing codons in AT-rich dicot genome is due to compositional constraint as a factor of codon usage bias. Odds ratio analysis identified the dinucleotides TpG, TpC, GpA, CpA and CpT as over-represented, where, CpG and TpA as under-represented dinucleotides. The results of (NcExp?NcObs)/NcExp plot suggests that selection pressure other than mutation played a significant role in influencing the pattern of codon usage in these dicots. PR2 analysis revealed the significant role of selection pressure on codon usage. Analysis of varience on codon usage at start and stop site showed variation in codon selection in these sites. This study provides evidence that the dicot genes were subjected to compositional selection pressure.