Sequence and analysis of chromosome 2 of the plant Arabidopsis thaliana

Herbicide injury induces DNA methylome alterations in Arabidopsis

The emergence of herbicide-resistant weeds is a major threat facing modern agriculture. Over 470 weedy-plant populations have developed resistance to herbicides. Traditional evolutionary mechanisms are not always sufficient to explain the rapidity with which certain weed populations adapt in response to herbicide exposure. Stress-induced epigenetic changes, such as alterations in DNA methylation, are potential additional adaptive mechanisms for herbicide resistance. We performed methylC sequencing of Arabidopsis thaliana leaves that developed after either mock treatment or two different sub-lethal doses of the herbicide glyphosate, the most-used herbicide in the history of agriculture. The herbicide injury resulted in 9,205 differentially methylated regions (DMRs) across the genome. In total, 5,914 of these DMRs were induced in a dose-dependent manner, wherein the methylation levels were positively correlated to the severity of the herbicide injury, suggesting that plants can modulate the magnitude of methylation changes based on the severity of the stress. Of the 3,680 genes associated with glyphosate-induced DMRs, only 7% were also implicated in methylation changes following biotic or salinity stress. These results demonstrate that plants respond to herbicide stress through changes in methylation patterns that are, in general, dose-sensitive and, at least partially, stress-specific.

Ectopic Expression of Aeluropus littoralis Plasma Membrane Protein Gene AlTMP1 Confers Abiotic Stress Tolerance in Transgenic Tobacco by Improving Water Status and Cation Homeostasis

We report here the isolation and functional analysis of AlTMP1 gene encoding a member of the PMP3 protein family. In Aeluropus littoralis, AlTMP1 is highly induced by abscisic acid (ABA), cold, salt, and osmotic stresses. Transgenic tobacco expressing AlTMP1 exhibited enhanced tolerance to salt, osmotic, H₂O₂, heat and freezing stresses at the seedling stage. Under greenhouse conditions, the transgenic plants showed a higher level of tolerance to drought than to salinity. Noteworthy, AlTMP1 plants yielded two- and five-fold more seeds than non-transgenic plants (NT) under salt and drought stresses, respectively. The leaves of AlTMP1 plants accumulated lower Na⁺ but higher K⁺ and Ca2+ than those of NT plants. Tolerance to osmotic and salt stresses was associated with higher membrane stability, low electrolyte leakage, and improved water status. Finally, accumulation of AlTMP1 in tobacco altered the regulation of some stress-related genes in either a positive (NHX1, CAT1, APX1, and DREB1A) or negative (HKT1 and KT1) manner that could be related to the observed tolerance. These results suggest that AlTMP1 confers stress tolerance in tobacco through maintenance of ion homeostasis, increased membrane integrity, and water status. The observed tolerance may be due to a direct or indirect effect of AlTMP1 on the expression of stress-related genes which could stimulate an adaptive potential not present in NT plants.

Analysis of nuclear mitochondrial DNAs and factors affecting patterns of integration in plant species

Sequences homologous to organellar DNA that have been integrated into nuclear genomes are referred to as nuclear mitochondrial DNAs (NUMTs) and nuclear plastid DNAs (NUPTs). NUMTs in nine plant species were analyzed to reveal the integration patterns and possible factors involved. The cumulative lengths of NUMTs in two-thirds of species analyzed were greater than those of NUPTs observed in a previous study. The age distribution of NUMTs was similar to that of NUPTs, suggesting similar mechanisms for integration and degradation of both NUPTs and NUMTs. Nuclear genome size and the cumulative length of NUMTs showed a significant positive correlation for older but not younger NUMTs. The same correlation was also found between nuclear genome size and older NUPTs in 17 species. These results suggested that genome size is a key factor to determine the cumulative length of relatively older NUPTs/NUMTs. Although the factor(s) determining the cumulative length of younger NUPTs/NUMTs is unclear, these sequences may be more deleterious, which could explain the different manner of determining the cumulative length of younger NUPTs/NUMTs in nuclear genomes. In addition, a relationship between the cumulative length of integrated NUMTs and complexity of mitochondrial genomes (i.e., the number of repeats) was found. The results indicate that the structural complexity of both NUMTs and their original mitochondrial sequences affects integration and degradation processes.

Centromere location in Arabidopsis is unaltered by extreme divergence in CENH3 protein sequence

During cell division, spindle fibers attach to chromosomes at centromeres. The DNA sequence at regional centromeres is fast evolving with no conserved genetic signature for centromere identity. Instead CENH3, a centromere-specific histone H3 variant, is the epigenetic signature that specifies centromere location across both plant and animal kingdoms. Paradoxically, CENH3 is also adaptively evolving. An ongoing question is whether CENH3 evolution is driven by a functional relationship with the underlying DNA sequence. Here, we demonstrate that despite extensive protein sequence divergence, CENH3 histones from distant species assemble centromeres on the same underlying DNA sequence. We first characterized the organization and diversity of centromere repeats in wild-type Arabidopsis thaliana. We show that A. thaliana CENH3-containing nucleosomes exhibit a strong preference for a unique subset of centromeric repeats. These sequences are largely missing from the genome assemblies and represent the youngest and most homogeneous class of repeats. Next, we tested the evolutionary specificity of this interaction in a background in which the native A. thaliana CENH3 is replaced with CENH3s from distant species. Strikingly, we find that CENH3 from Lepidium oleraceum and Zea mays, although specifying epigenetically weaker centromeres that result in genome elimination upon outcrossing, show a binding pattern on A. thaliana centromere repeats that is indistinguishable from the native CENH3. Our results demonstrate positional stability of a highly diverged CENH3 on independently evolved repeats, suggesting that the sequence specificity of centromeres is determined by a mechanism independent of CENH3.

Phosphoribosyl Diphosphate (PRPP): Biosynthesis, Enzymology, Utilization, and Metabolic Significance

Phosphoribosyl diphosphate (PRPP) is an important intermediate in cellular metabolism. PRPP is synthesized by PRPP synthase, as follows: ribose 5-phosphate + ATP → PRPP + AMP. PRPP is ubiquitously found in living organisms and is used in substitution reactions with the formation of glycosidic bonds. PRPP is utilized in the biosynthesis of purine and pyrimidine nucleotides, the amino acids histidine and tryptophan, the cofactors NAD and tetrahydromethanopterin, arabinosyl monophosphodecaprenol, and certain aminoglycoside antibiotics. The participation of PRPP in each of these metabolic pathways is reviewed. Central to the metabolism of PRPP is PRPP synthase, which has been studied from all kingdoms of life by classical mechanistic procedures. The results of these analyses are unified with recent progress in molecular enzymology and the elucidation of the three-dimensional structures of PRPP synthases from eubacteria, archaea, and humans. The structures and mechanisms of catalysis of the five diphosphoryltransferases are compared, as are those of selected enzymes of diphosphoryl transfer, phosphoryl transfer, and nucleotidyl transfer reactions. PRPP is used as a substrate by a large number phosphoribosyltransferases. The protein structures and reaction mechanisms of these phosphoribosyltransferases vary and demonstrate the versatility of PRPP as an intermediate in cellular physiology. PRPP synthases appear to have originated from a phosphoribosyltransferase during evolution, as demonstrated by phylogenetic analysis. PRPP, furthermore, is an effector molecule of purine and pyrimidine nucleotide biosynthesis, either by binding to PurR or PyrR regulatory proteins or as an allosteric activator of carbamoylphosphate synthetase. Genetic analyses have disclosed a number of mutants altered in the PRPP synthase-specifying genes in humans as well as bacterial species.

Selection on Position of Nonsense Codons in Introns

Introns occasionally remain in mature messenger RNAs (mRNAs) due to splicing errors and the translated, aberrant proteins that result represent a metabolic cost and may have other deleterious consequences. The nonsense-mediated decay (NMD) pathway degrades aberrant mRNAs, which it recognizes by the presence of an in-frame premature termination codon (PTC). We investigated whether selection has shaped the location of PTCs in introns to reduce waste and facilitate NMD. We found across seven model organisms, that in both first and last introns, PTCs occur earlier in introns than expected by chance, suggesting that selection favors earlier position. This pattern is more pronounced in species with larger effective population sizes. The pattern does not hold for last introns in the two mammal species, however, perhaps because in these species NMD is not initiated from 3'-terminal introns. We conclude that there is compelling evidence that the location of PTCs is shaped by selection for reduced waste and efficient degradation of aberrant mRNAs.

Genome-wide characterization of non-reference transposons in crops suggests non-random insertion

BACKGROUND

Transposons (transposable elements or TEs) are DNA sequences that can change their position within the genome. A large number of TEs have been identified in reference genome of each crop(named accumulated TEs), which are the important part of genome. However, whether there existed TEs with different insert positions in resequenced crop accession genomes from those of reference genome (named non-reference transposable elements, non-ref TEs), and what the characteristics (such as the number, type and distribution) are. To identify and characterize crop non-ref TEs, we analyzed non-ref TEs in more than 125 accessions from rice (Oryza sativa), maize (Zea mays) and sorghum (Sorghum bicolor) using resequenced data with paired-end mapping methods.

RESULTS

We identified 13,066, 23,866 and 35,679 non-ref TEs in rice, maize and sorghum, respectively. Genome-wide characterization analysis shows that most of non-ref TEs were unique and non-ref TE classes shows different among rice, maize and sorghum. We found that non-ref TEs have a strong positive correlation with gene number and have a bias toward insertion near genes, but with a preference for avoiding coding regions in maize and sorghum. The genes affected by non-ref TE insertion were functionally enriched for stress response mechanisms in all three crops.

CONCLUSIONS

These observations suggest that transposon insertion is not a random event and it makes genomic diversity, which may affect the intraspecific adaption and evolution of crops.

Computational Study of the Binding Mechanism of Actin-Depolymerizing Factor 1 with Actin in Arabidopsis thaliana

Actin is a highly conserved protein. It plays important roles in cellular function and exists either in the monomeric (G-actin) or polymeric form (F-actin). Members of the actin-depolymerizing factor (ADF)/cofilin protein family bind to both G-actin and F-actin and play vital roles in actin dynamics by manipulating the rates of filament polymerization and depolymerization. It has been reported that the S6D and R98A/K100A mutants of actin-depolymerizing factor 1 (ADF1) in Arabidopsis thaliana decreased the binding affinity of ADF for the actin monomer. To investigate the binding mechanism and dynamic behavior of the ADF1–actin complex, we constructed a homology model of the AtADF1–actin complex based on the crystal structure of AtADF1 and the twinfilin C-terminal ADF-H domain in a complex with a mouse actin monomer. The model was then refined for subsequent molecular dynamics simulations. Increased binding energy of the mutated system was observed using the Molecular Mechanics Generalized Born Surface Area and Poisson–Boltzmann Surface Area (MM-GB/PBSA) methods. To determine the residues that make decisive contributions to the ADF1 actin-binding affinity, per-residue decomposition and computational alanine scanning analyses were performed, which provided more detailed information on the binding mechanism. Root-mean-square fluctuation and principal component analyses confirmed that the S6D and R98A/K100A mutants induced an increased conformational flexibility. The comprehensive molecular insight gained from this study is of great importance for understanding the binding mechanism of ADF1 and G-actin.

A computational framework for mapping the timing of vegetative phase change

Phase change plays a prominent role in determining the form of growth and development. Although considerable attention has been focused on identifying the regulatory control mechanisms of phase change, a detailed understanding of the genetic architecture of this phenomenon is still lacking. We address this issue by deriving a computational model. The model is founded on the framework of functional mapping aimed at characterizing the interplay between quantitative trait loci (QTLs) and development through biologically meaningful mathematical equations. A multiphasic growth equation was implemented into functional mapping, which, via a series of hypothesis tests, allows the quantification of how QTLs regulate the timing and pattern of vegetative phase transition between independently regulated, temporally coordinated processes. The model was applied to analyze stem radial growth data of an interspecific hybrid family derived from two Populus species during the first 24 yr of ontogeny. Several key QTLs related to phase change have been characterized, most of which were observed to be in the adjacent regions of candidate genes. The identification of phase transition QTLs, whose expression is regulated by endogenous and environmental signals, may enhance our understanding of the evolution of development in changing environments.

Systematic analysis of plant mitochondrial and chloroplast small RNAs suggests organelle-specific mRNA stabilization mechanisms

Land plant organellar genomes encode a small number of genes, many of which are essential for respiration and photosynthesis. Organellar gene expression is characterized by a multitude of RNA processing events that lead to stable, translatable transcripts. RNA binding proteins (RBPs), have been shown to generate and protect transcript termini and eventually induce the accumulation of short RNA footprints. We applied knowledge of such RBP-derived footprints to develop software (sRNA miner) that enables identification of RBP footprints, or other clusters of small RNAs, in organelles. We used this tool to determine mitochondrial and chloroplast cosRNAs (clustered organellar sRNAs) in Arabidopsis. We found that in mitochondria, cosRNAs coincide with transcript 3'-ends, but are largely absent from 5'-ends. In chloroplasts this bias is absent, suggesting a different mode of 5' processing, possibly owing to different sets of RNases. Furthermore, we identified a large number of cosRNAs that represent silenced insertions of mitochondrial DNA in the nuclear genome of Arabidopsis. Steady-state RNA analyses demonstrate that cosRNAs display differential accumulation during development. Finally, we demonstrate that the chloroplast RBP PPR10 associates in vivo with its cognate cosRNA. A hypothetical role of cosRNAs as competitors of mRNAs for PPR proteins is discussed.

Horizontal Gene Transfer of Pectinases from Bacteria Preceded the Diversification of Stick and Leaf Insects

Genes acquired by horizontal transfer are increasingly being found in animal genomes. Understanding their origin and evolution requires knowledge about the phylogenetic relationships from both source and recipient organisms. We used RNASeq data and respective assembled transcript libraries to trace the evolutionary history of polygalacturonase (pectinase) genes in stick insects (Phasmatodea). By mapping the distribution of pectinase genes on a Polyneoptera phylogeny, we identified the transfer of pectinase genes from known phasmatodean gut microbes into the genome of an early euphasmatodean ancestor that took place between 60 and 100 million years ago. This transfer preceded the rapid diversification of the suborder, enabling symbiont-free pectinase production that would increase the insects’ digestive efficiency and reduce dependence on microbes. Bacteria-to-insect gene transfer was thought to be uncommon, however the increasing availability of large-scale genomic data may change this prevailing notion.

Characterization of the watercress (Nasturtium officinale R. Br.; Brassicaceae) transcriptome using RNASeq and identification of candidate genes for important phytonutrient traits linked to human health

BACKGROUND

Consuming watercress is thought to provide health benefits as a consequence of its phytonutrient composition. However, for watercress there are currently limited genetic resources underpinning breeding efforts for either yield or phytonutritional traits. In this paper, we use RNASeq data from twelve watercress accessions to characterize the transcriptome, perform candidate gene mining and conduct differential expression analysis for two key phytonutritional traits: antioxidant (AO) capacity and glucosinolate (GLS) content.

RESULTS

The watercress transcriptome was assembled to 80,800 transcripts (48,732 unigenes); 71 % of which were annotated based on orthology to Arabidopsis. Differential expression analysis comparing watercress accessions with 'high' and 'low' AO and GLS resulted in 145 and 94 differentially expressed loci for AO capacity and GLS respectively. Differentially expressed loci between high and low AO watercress were significantly enriched for genes involved in plant defence and response to stimuli, in line with the observation that AO are involved in plant stress-response. Differential expression between the high and low GLS watercress identified links to GLS regulation and also novel transcripts warranting further investigation. Additionally, we successfully identified watercress orthologs for Arabidopsis phenylpropanoid, GLS and shikimate biosynthesis pathway genes, and compiled a catalogue of polymorphic markers for future applications.

CONCLUSIONS

Our work describes the first transcriptome of watercress and establishes the foundation for further molecular study by providing valuable resources, including sequence data, annotated transcripts, candidate genes and markers.

Docking, QM/MM, and molecular dynamics simulations of the hexose transporter from Plasmodium falciparum (PfHT)

Malaria is the most prevalent parasitic disease in the world. Currently, an effective vaccine for malaria does not exist, and chemotherapy must be used to treat the disease. Because of increasing resistance to current antimalarial drugs, new treatments must be developed. Among the many potential molecular targets, the hexose transporter of Plasmodium falciparum (PfHT) is particularly promising because it plays a vital role in glucose transport for the parasite. Thus, this study aims to determine the three-dimensional structure of PfHT and to describe the intermolecular interactions between active glycoside derivatives and PfHT. Such information should aid in the development of new antimalarial drugs. The receptor PfHT was constructed from primary sequences deposited in the SWISS MODEL database. Next, molecular docking simulations between O-(undec-10-en)-l-D-glucose and the constructed active site models were performed using Autodock Vina. The glycoside derivative-PfHT complexes were then refined using the hybrid QM/MM (PM3/ff03) method within the AMBER package. The models were then evaluated using Ramachandran plots, which indicated that 93.2% of the residues in the refined PfHT models (P5) were present in favorable regions. Furthermore, graphical plots using ANOLEA showed that the potential energies of interaction for atoms unbonded to P5 were negative. Finally, the O-(undec-10-en)-l-D-glucose-PfHT complex was evaluated using 20-ns Molecular Dynamics simulations with an ff03 force field. Docking and QM/MM studies revealed the amino acids essential for molecular recognition of and activity on glycosides. Inhibition of glucose transporters may prevent the development and metabolism of P. falciparum, so a description of the receptor's structure is a critical step towards rational drug design.

Functional investigation of two 1-aminocyclopropane-1-carboxylate (ACC) synthase-like genes in the moss Physcomitrella patens

Two ACC synthase-like (ACL) proteins in the moss Physcomitrella patens have no ACS activity, and PpACL1 functions as an L -cystine/ L -cysteine C-S lyase. The ethylene biosynthetic pathway has been well characterized in higher plants, and homologs of a key enzyme in this pathway, ACS, have been reported in several algae and mosses, including Physcomitrella patens. However, the function of the ACS homologs in P. patens has not been investigated. In this research, we cloned two putative ACS genes from the P. patens genome, namely PpACS-Like 1 and 2, and investigated whether their encoded proteins had in vitro and in vivo ACS activity. In vitro biochemical assays using purified PpACL1 and PpACL2 showed that neither protein had ACS activity. Subsequently, we generated transgenic Arabidopsis lines expressing 35S:PpACL1 and 35S:PpACL2, and found that the transgenic etiolated seedlings that overexpressed either of these proteins lacked the constitutive triple response phenotype and did not emit excess levels of ethylene, indicating that neither of the PpACS-Like proteins had in vivo ACS activity. Furthermore, we found that PpACL1 functions as a C-S lyase that uses L-cystine and L-cysteine as substrates, rather than as an aminotransferase. Together, these results indicated that PpACL1 and PpACL2 are not true ACS genes as those found in higher plants.

Identification and Comparative Analysis of H2O2-Scavenging Enzymes (Ascorbate Peroxidase and Glutathione Peroxidase) in Selected Plants Employing Bioinformatics Approaches

Among major reactive oxygen species (ROS), hydrogen peroxide (H2O2) exhibits dual roles in plant metabolism. Low levels of H2O2 modulate many biological/physiological processes in plants; whereas, its high level can cause damage to cell structures, having severe consequences. Thus, steady-state level of cellular H2O2 must be tightly regulated. Glutathione peroxidases (GPX) and ascorbate peroxidase (APX) are two major ROS-scavenging enzymes which catalyze the reduction of H2O2 in order to prevent potential H2O2-derived cellular damage. Employing bioinformatics approaches, this study presents a comparative evaluation of both GPX and APX in 18 different plant species, and provides valuable insights into the nature and complex regulation of these enzymes. Herein, (a) potential GPX and APX genes/proteins from 18 different plant species were identified, (b) their exon/intron organization were analyzed, (c) detailed information about their physicochemical properties were provided, (d) conserved motif signatures of GPX and APX were identified, (e) their phylogenetic trees and 3D models were constructed, (f) protein-protein interaction networks were generated, and finally (g) GPX and APX gene expression profiles were analyzed. Study outcomes enlightened GPX and APX as major H2O2-scavenging enzymes at their structural and functional levels, which could be used in future studies in the current direction.

Mobile Bacterial Group II Introns at the Crux of Eukaryotic Evolution

This review focuses on recent developments in our understanding of group II intron function, the relationships of these introns to retrotransposons and spliceosomes, and how their common features have informed thinking about bacterial group II introns as key elements in eukaryotic evolution. Reverse transcriptase-mediated and host factor-aided intron retrohoming pathways are considered along with retrotransposition mechanisms to novel sites in bacteria, where group II introns are thought to have originated. DNA target recognition and movement by target-primed reverse transcription infer an evolutionary relationship among group II introns, non-LTR retrotransposons, such as LINE elements, and telomerase. Additionally, group II introns are almost certainly the progenitors of spliceosomal introns. Their profound similarities include splicing chemistry extending to RNA catalysis, reaction stereochemistry, and the position of two divalent metals that perform catalysis at the RNA active site. There are also sequence and structural similarities between group II introns and the spliceosome's small nuclear RNAs (snRNAs) and between a highly conserved core spliceosomal protein Prp8 and a group II intron-like reverse transcriptase. It has been proposed that group II introns entered eukaryotes during bacterial endosymbiosis or bacterial-archaeal fusion, proliferated within the nuclear genome, necessitating evolution of the nuclear envelope, and fragmented giving rise to spliceosomal introns. Thus, these bacterial self-splicing mobile elements have fundamentally impacted the composition of extant eukaryotic genomes, including the human genome, most of which is derived from close relatives of mobile group II introns.

Group II introns: versatile ribozymes and retroelements

Group II introns are catalytic RNAs (ribozymes) and retroelements found in the genomes of bacteria, archaebacteria, and organelles of some eukaryotes. The prototypical retroelement form consists of a structurally conserved RNA and a multidomain reverse transcriptase protein, which interact with each other to mediate splicing and mobility reactions. A wealth of biochemical, cross-linking, and X-ray crystal structure studies have helped to reveal how the two components cooperate to carry out the splicing and mobility reactions. In addition to the standard retroelement form, group II introns have evolved into derivative forms by either losing specific splicing or mobility characteristics, or becoming functionally specialized. Of particular interest are the eukaryotic derivatives-the spliceosome, spliceosomal introns, and non-LTR retroelements-which together make up approximately half of the human genome. On a practical level, the properties of group II introns have been exploited to develop group II intron-based biotechnological tools. WIREs RNA 2016, 7:341-355. doi: 10.1002/wrna.1339 For further resources related to this article, please visit the WIREs website.

Transcriptome Profiling of Resistance to Fusarium oxysporum f. sp. conglutinans in Cabbage (Brassica oleracea) Roots

Fusarium wilt caused by Fusarium oxysporum f. sp. conglutinans (FOC) is a destructive disease of Brassica crops, which results in severe yield losses. There is little information available about the mechanism of disease resistance. To obtain an overview of the transcriptome profiles in roots of R4P1, a Brassica oleracea variety that is highly resistant to fusarium wilt, we compared the transcriptomes of samples inoculated with FOC and samples inoculated with distilled water. RNA-seq analysis generated more than 136 million 100-bp clean reads, which were assembled into 62,506 unigenes (mean size = 741 bp). Among them, 49,959 (79.92%) genes were identified based on sequence similarity searches, including SwissProt (29,050, 46.47%), Gene Ontology (GO) (33,767, 54.02%), Clusters of Orthologous Groups (KOG) (14,721, 23.55%) and Kyoto Encyclopedia of Genes and Genomes Pathway database (KEGG) (12,974, 20.76%) searches; digital gene expression analysis revealed 885 differentially expressed genes (DEGs) between infected and control samples at 4, 12, 24 and 48 hours after inoculation. The DEGs were assigned to 31 KEGG pathways. Early defense systems, including the MAPK signaling pathway, calcium signaling and salicylic acid-mediated hypersensitive response (SA-mediated HR) were activated after pathogen infection. SA-dependent systemic acquired resistance (SAR), ethylene (ET)- and jasmonic (JA)-mediated pathways and the lignin biosynthesis pathway play important roles in plant resistance. We also analyzed the expression of defense-related genes, such as genes encoding pathogenesis-related (PR) proteins, UDP-glycosyltransferase (UDPG), pleiotropic drug resistance, ATP-binding cassette transporters (PDR-ABC transporters), myrosinase, transcription factors and kinases, which were differentially expressed. The results of this study may contribute to efforts to identify and clone candidate genes associated with disease resistance and to uncover the molecular mechanism underlying FOC resistance in cabbage.

Cytogenetic and Sequence Analyses of Mitochondrial DNA Insertions in Nuclear Chromosomes of Maize

The transfer of mitochondrial DNA (mtDNA) into nuclear genomes is a regularly occurring process that has been observed in many species. Few studies, however, have focused on the variation of nuclear-mtDNA sequences (NUMTs) within a species. This study examined mtDNA insertions within chromosomes of a diverse set of Zea mays ssp. mays (maize) inbred lines by the use of fluorescence in situ hybridization. A relatively large NUMT on the long arm of chromosome 9 (9L) was identified at approximately the same position in four inbred lines (B73, M825, HP301, and Oh7B). Further examination of the similarly positioned 9L NUMT in two lines, B73 and M825, indicated that the large size of these sites is due to the presence of a majority of the mitochondrial genome; however, only portions of this NUMT (~252 kb total) were found in the publically available B73 nuclear sequence for chromosome 9. Fiber-fluorescence in situ hybridization analysis estimated the size of the B73 9L NUMT to be ~1.8 Mb and revealed that the NUMT is methylated. Two regions of mtDNA (2.4 kb and 3.3 kb) within the 9L NUMT are not present in the B73 mitochondrial NB genome; however, these 2.4-kb and 3.3-kb segments are present in other Zea mitochondrial genomes, including that of Zea mays ssp. parviglumis, a progenitor of domesticated maize.

Genome-Wide Analysis of Oleosin Gene Family in 22 Tree Species: An Accelerator for Metabolic Engineering of BioFuel Crops and Agrigenomics Industrial Applications?

Trees contribute to enormous plant oil reserves because many trees contain 50%-80% of oil (triacylglycerols, TAGs) in the fruits and kernels. TAGs accumulate in subcellular structures called oil bodies/droplets, in which TAGs are covered by low-molecular-mass hydrophobic proteins called oleosins (OLEs). The OLEs/TAGs ratio determines the size and shape of intracellular oil bodies. There is a lack of comprehensive sequence analysis and structural information of OLEs among diverse trees. The objectives of this study were to identify OLEs from 22 tree species (e.g., tung tree, tea-oil tree, castor bean), perform genome-wide analysis of OLEs, classify OLEs, identify conserved sequence motifs and amino acid residues, and predict secondary and three-dimensional structures in tree OLEs and OLE subfamilies. Data mining identified 65 OLEs with perfect conservation of the "proline knot" motif (PX5SPX3P) from 19 trees. These OLEs contained >40% hydrophobic amino acid residues. They displayed similar properties and amino acid composition. Genome-wide phylogenetic analysis and multiple sequence alignment demonstrated that these proteins could be classified into five OLE subfamilies. There were distinct patterns of sequence conservation among the OLE subfamilies and within individual tree species. Computational modeling indicated that OLEs were composed of at least three α-helixes connected with short coils without any β-strand and that they exhibited distinct 3D structures and ligand binding sites. These analyses provide fundamental information in the similarity and specificity of diverse OLE isoforms within the same subfamily and among the different species, which should facilitate studying the structure-function relationship and identify critical amino acid residues in OLEs for metabolic engineering of tree TAGs.

Rhamnogalacturonan I modifying enzymes: an update

Rhamnogalacturonan I (RGI) modifying enzymes catalyse the degradation of the RGI backbone and encompass enzymes specific for either the α1,2-bond linking galacturonic acid to rhamnose or the α1,4-bond linking rhamnose to galacturonic acid in the RGI backbone. The first microbial enzyme found to be able to catalyse the degradation of the RGI backbone, an endo-hydrolase (EC 3.2.1.171) derived from Aspergillus aculeatus, was discovered 25 years ago. Today the group of RGI modifying enzymes encompasses endo- and exo-hydrolases as well as lyases. The RGI hydrolases, EC 3.2.1.171-EC 3.2.1.174, have been described to be produced by Aspergillus spp. and Bacillus subtilis and are categorized in glycosyl hydrolase families 28 and 105. The RGI lyases, EC 4.2.2.23-EC 4.2.2.24, have been isolated from different fungi and bacterial species and are categorized in polysaccharide lyase families 4 and 11. This review brings together the available knowledge of the RGI modifying enzymes and provides a detailed overview of biocatalytic reaction characteristics, classification, structure-function traits, and analyses the protein properties of these enzymes by multiple sequence alignments in neighbour-joining phylogenetic trees. Some recently detected unique structural features and dependence of calcium for activity of some of these enzymes (notably the lyases) are discussed and newly published results regarding improvement of their thermostability by protein engineering are highlighted. Knowledge of these enzymes is important for understanding microbial plant cell wall degradation and for advancing enzymatic processing and biorefining of pectinaceous plant biomass.

Next Generation Sequencing and Transcriptome Analysis Predicts Biosynthetic Pathway of Sennosides from Senna (Cassia angustifolia Vahl.), a Non-Model Plant with Potent Laxative Properties

Senna (Cassia angustifolia Vahl.) is a world's natural laxative medicinal plant. Laxative properties are due to sennosides (anthraquinone glycosides) natural products. However, little genetic information is available for this species, especially concerning the biosynthetic pathways of sennosides. We present here the transcriptome sequencing of young and mature leaf tissue of Cassia angustifolia using Illumina MiSeq platform that resulted in a total of 6.34 Gb of raw nucleotide sequence. The sequence assembly resulted in 42230 and 37174 transcripts with an average length of 1119 bp and 1467 bp for young and mature leaf, respectively. The transcripts were annotated using NCBI BLAST with 'green plant database (txid 33090)', Swiss Prot, Kyoto Encylcopedia of Genes & Genomes (KEGG), Cluster of Orthologous Gene (COG) and Gene Ontology (GO). Out of the total transcripts, 40138 (95.0%) and 36349 (97.7%) from young and mature leaf, respectively, were annotated by BLASTX against green plant database of NCBI. We used InterProscan to see protein similarity at domain level, a total of 34031 (young leaf) and 32077 (mature leaf) transcripts were annotated against the Pfam domains. All transcripts from young and mature leaf were assigned to 191 KEGG pathways. There were 166 and 159 CDS, respectively, from young and mature leaf involved in metabolism of terpenoids and polyketides. Many CDS encoding enzymes leading to biosynthesis of sennosides were identified. A total of 10,763 CDS differentially expressing in both young and mature leaf libraries of which 2,343 (21.7%) CDS were up-regulated in young compared to mature leaf. Several differentially expressed genes found functionally associated with sennoside biosynthesis. CDS encoding for many CYPs and TF families were identified having probable roles in metabolism of primary as well as secondary metabolites. We developed SSR markers for molecular breeding of senna. We have identified a set of putative genes involved in various secondary metabolite pathways, especially those related to the synthesis of sennosides which will serve as an important platform for public information about gene expression, genomics, and functional genomics in senna.

In silico discovery of a nearly complete mitochondrial genome Numt in the dog (Canis lupus familiaris) nuclear genome

Through bacterial cloning, a non-specific product co-amplified in a previous whole mitochondrial genome study of Canis lupus familiaris was identified as part of a Numt on chromosome 29 of the dog. Even though further analysis confirmed the fidelity of the mitochondrial genome sequencing results, it still highlighted the risk of Numt contamination. A computer-based search of the dog's nuclear genome for segments homologous to the mtDNA sequence revealed the extent of this risk. Over 150 Numts of various sizes were observed throughout all but two chromosomes, covering all positions of the mtDNA. One of the Numts on chromosome 11 even covered over 95 % of the entire dog mtDNA sequence. This comprehensive list of Numts was provided to assist researchers with the evaluation of dog mtDNA sequencing protocols for Numt co-amplification.

Determining mitochondrial transcript termini for the study of transcription start sites and transcript 5' end maturation

Mitochondrial gene expression in plants is considerably more complex than in animals or fungi. In plants, mitochondrial transcripts are generated from transcription initiation at numerous, poorly conserved promoters located throughout the mitochondrial genome. Most genes have more than one transcription start site. Posttranscriptional RNA 5' end maturation contributes to the diversity of transcripts produced from each mitochondrial gene. Understanding transcriptional mechanisms and transcript maturation requires knowledge on transcription start sites and processing sites. This chapter describes two different, complementary experimental approaches for determining these sites in mitochondrial genomes through mapping of transcript 5' ends. In order to distinguish 5' ends deriving from transcription initiation, both strategies exploit the presence of triphosphates at these specific 5' termini.

Evolution of group II introns

Present in the genomes of bacteria and eukaryotic organelles, group II introns are an ancient class of ribozymes and retroelements that are believed to have been the ancestors of nuclear pre-mRNA introns. Despite long-standing speculation, there is limited understanding about the actual pathway by which group II introns evolved into eukaryotic introns. In this review, we focus on the evolution of group II introns themselves. We describe the different forms of group II introns known to exist in nature and then address how these forms may have evolved to give rise to spliceosomal introns and other genetic elements. Finally, we summarize the structural and biochemical parallels between group II introns and the spliceosome, including recent data that strongly support their hypothesized evolutionary relationship.

Enhanced seed oil content by overexpressing genes related to triacylglyceride synthesis

Oilseed rape (Brassica napus) is one of the most important oilseed crops globally. To meet increasing demand for oil-based products, the ability to enhance desirable oil content in the seed is required. This study assessed the capability of five genes in the triacylglyceride (TAG) synthesis pathway to enhance oil content. The genes BnGPDH, BnGPAT, BnDGAT, ScGPDH and ScLPAAT were overexpressed separately in a tobacco (Nicotiana benthamiana) model system, and simultaneously by pyramiding in B. napus, under the control of a seed specific Napin promoter. ScLPAAT transgenic plants showed a significant increase of 6.84% to 8.55% in oil content in tobacco seeds, while a ~4% increase was noted for BnGPDH and BnGPAT transgenic seeds. Seed-specific overexpression of all four genes in B. napus resulted in as high a 12.57% to 14.46% increased in seed oil content when compared to WT, equaling close to the sum of the single-gene overexpression increases in tobacco. Taken together, our study demonstrates that BnGPDH, BnGPAT and ScLPAAT may effectively increase seed oil content, and that simultaneous overexpression of these in transgenic B. napus may further enhance the desirable oil content relative to single-gene overexpressors.

Transcriptional consequence and impaired gametogenesis with high-grade aneuploidy in Arabidopsis thaliana

Aneuploidy features a numerical chromosome variant that the number of chromosomes in the nucleus of a cell is not an exact multiple of the haploid number, which may have an impact on morphology and gene expression. Here we report a tertiary trisomy uncovered by characterizing a T-DNA insertion mutant (aur2-1/+) in the Arabidopsis (Arabidopsis thaliana) AURORA2 locus. Whole-genome analysis with DNA tiling arrays revealed a chromosomal translocation linked to the aur2-1 allele, which collectively accounted for a tertiary trisomy 2. Morphologic, cytogenetic and genetic analyses of aur2-1 progeny showed impaired male and female gametogenesis to various degrees and a tight association of the aur2-1 allele with the tertiary trisomy that was preferentially inherited. Transcriptome analysis showed overlapping and distinct gene expression profiles between primary and tertiary trisomy 2 plants, particularly genes involved in response to stress and various types of external and internal stimuli. Additionally, transcriptome and gene ontology analyses revealed an overrepresentation of nuclear-encoded organelle-related genes functionally involved in plastids, mitochondria and peroxisomes that were differentially expressed in at least three if not all Arabidopsis trisomics. These observations support a previous hypothesis that aneuploid cells have higher energy requirement to overcome the detrimental effects of an unbalanced genome. Moreover, our findings extend the knowledge of the complex nature of the T-DNA insertion event influencing plant genomic integrity by creating high-grade trisomy. Finally, gene expression profiling results provide useful information for future research to compare primary and tertiary trisomics for the effects of aneuploidy on plant cell physiology.

Identification and expression analysis of cold and freezing stress responsive genes of Brassica oleracea

Cold and freezing stress is a major environmental constraint to the production of Brassica crops. Enhancement of tolerance by exploiting cold and freezing tolerance related genes offers the most efficient approach to address this problem. Cold-induced transcriptional profiling is a promising approach to the identification of potential genes related to cold and freezing stress tolerance. In this study, 99 highly expressed genes were identified from a whole genome microarray dataset of Brassica rapa. Blast search analysis of the Brassica oleracea database revealed the corresponding homologous genes. To validate their expression, pre-selected cold tolerant and susceptible cabbage lines were analyzed. Out of 99 BoCRGs, 43 were differentially expressed in response to varying degrees of cold and freezing stress in the contrasting cabbage lines. Among the differentially expressed genes, 18 were highly up-regulated in the tolerant lines, which is consistent with their microarray expression. Additionally, 12 BoCRGs were expressed differentially after cold stress treatment in two contrasting cabbage lines, and BoCRG54, 56, 59, 62, 70, 72 and 99 were predicted to be involved in cold regulatory pathways. Taken together, the cold-responsive genes identified in this study provide additional direction for elucidating the regulatory network of low temperature stress tolerance and developing cold and freezing stress resistant Brassica crops.

Dynamic transposable element accumulation in the nascent sex chromosomes of papaya

From their inception, Y chromosomes in plants and animals are subjected to the powerful effects of Müller's ratchet, a process spurred by suppression of recombination that results in a rapid accumulation of mutations and repetitive elements. These mutations eventually lead to gene loss and degeneration of the Y chromosome. Y chromosomes in mammals are ancient, whereas most sex chromosomes in plants and many in insects and fish evolved recently. Sex type in papaya is controlled by a pair of nascent sex chromosomes that evolved around 7 million years ago. The papaya X and Y^h were recently sequenced, providing valuable insight into the early stages of sex chromosome evolution. Here we discuss the fruits of this work with a focus on the repeat accumulation, gene trafficking and promiscuous DNA sequences found in the slowly degenerating Y^h chromosome of papaya.

A skellam model to identify differential patterns of gene expression induced by environmental signals

Background

RNA-seq, based on deep-sequencing techniques, has been widely employed to precisely measure levels of transcripts and their isoforms expressed under different conditions. However, robust statistical tools used to analyze these complex datasets are lacking. By grouping genes with similar expression profiles across treatments, cluster analysis provides insight into gene functions and networks that have become increasingly important.

Results

We proposed and verified a cluster algorithm based on a skellam model for grouping genes into distinct groups based on the pattern of gene expression in response to changing conditions or in different tissues. This algorithm capitalizes on the skellam distribution to capture the count property of RNA-seq data and clusters genes in different environments. A two-stage hierarchical expectation-maximization (EM) algorithm was implemented to estimate the optimal number of groups and mean expression levels of each group across two environments. A procedure was formulated to test whether and how a given group shows a plastic response to environmental changes. The model was used to analyze an RNA-seq dataset measured from reciprocal crosses of early Arabidopsis thaliana embryos that respond differently based on the extent of maternal and paternal genome contributions, from which genes associated with maternal and paternal contributions were identified. Simulation studies were also performed to validate the statistical behavior of the model.

Conclusions

This model is a useful tool for clustering gene expression data by RNA-seq, thus facilitating our understanding of gene functions and networks.

The coffee genome provides insight into the convergent evolution of caffeine biosynthesis

Coffee is a valuable beverage crop due to its characteristic flavor, aroma, and the stimulating effects of caffeine. We generated a high-quality draft genome of the species Coffea canephora, which displays a conserved chromosomal gene order among asterid angiosperms. Although it shows no sign of the whole-genome triplication identified in Solanaceae species such as tomato, the genome includes several species-specific gene family expansions, among them N-methyltransferases (NMTs) involved in caffeine production, defense-related genes, and alkaloid and flavonoid enzymes involved in secondary compound synthesis. Comparative analyses of caffeine NMTs demonstrate that these genes expanded through sequential tandem duplications independently of genes from cacao and tea, suggesting that caffeine in eudicots is of polyphyletic origin.

Identification and characterization of LIM gene family in Brassica rapa

BACKGROUND

LIM (Lin-11, Isl-1 and Mec-3 domains) genes have been reported to trigger the formation of actin bundles, a major higher-order cytoskeletal assembly, in higher plants; however, the stress resistance related functions of these genes are still not well known. In this study, we collected 22 LIM genes designated as Brassica rapa LIM (BrLIM) from the Brassica database, analyzed the sequences, compared them with LIM genes of other plants and analyzed their expression after applying biotic and abiotic stresses in Chinese cabbage.

RESULTS

Upon sequence analysis these genes were confirmed as LIM genes and found to have a high degree of homology with LIM genes of other species. These genes showed distinct clusters when compared to other recognized LIM proteins upon phylogenetic analysis. Additionally, organ specific expression of these genes was observed in Chinese cabbage plants, with BrPLIM2a, b, c, BrDAR1, BrPLIM2e, f and g only being expressed in flower buds. Furthermore, the expression of these genes (except for BrDAR1 and BrPLIM2e) was high in the early flowering stages. The remaining genes were expressed in almost all organs tested. All BrDAR genes showed higher expression in flower buds compared to other organs. These organ specific expressions were clearly correlated with the phylogenetic grouping. In addition, BrWLIM2c and BrDAR4 responded to Fusarium oxysporum f. sp. conglutinans infection, while commonly two BrDARs and eight BrLIMs responded to cold, ABA and pH (pH5, pH7 and pH9) stress treatments in Chinese cabbage plants.

CONCLUSION

Taken together, the results of this study indicate that BrLIM and BrDAR genes may be involved in resistance against biotic and abiotic stresses in Brassica.

Identification and stress-induced expression of three 3β-hydroxysteroid dehydrogenases from Erysimum crepidifolium Rchb. and their putative role in cardenolide biosynthesis

3β-Hydroxysteroid dehydrogenases (3βHSD) are supposed to be involved in cardenolide biosynthesis in plants. Erysimum crepidifolium Rchb., a member of the Brassicaceae accumulating cardenolides, is a close relative to Arabidopsis thaliana. Full length cDNAs encoding for three individual 3βHSDs (EcHSD1, EcHSD2, EcHSD3) were isolated from E. crepidifolium leaves. EcHSD1 and EcHSD2 encode proteins assembled from 257 amino acids whereas EcHSD3 encodes a protein assembled from 260 amino acids. All three proteins qualify as members of the short-chain dehydrogenases/reductases family of proteins (SDRs). EcHSD1 and EcHSD2 shared a high amino acid sequence identity of about 86% and 91% with putative 3βHSDs of A. thaliana (AT2G47140 and AT2G47130). EcHSD3 showed high homology to the A. thaliana SDRs AT2G47150 (74%) and AT2G47120 (81%). All three EcHSD genes were expressed in Escherichia coli and the recombinant enzymes were characterized biochemically. All three recombinant EcHSDs catalyzed the dehydrogenation of pregnenolone and the 3-reduction of 5α/β-pregnane-3,20-dione when NAD and NADH were used as cosubstrates, respectively. After exposure to different stress conditions, no increased transcription was seen for EcHSD1 whereas EcHSD2 was expressed four times higher under osmotic stress than under control conditions. EcHSD3 expression was 10 times and 6 times higher after osmotic stress and MeJA treatment, respectively, than in controls.

Genomic features and regulatory roles of intermediate-sized non-coding RNAs in Arabidopsis

Recent advances in genome-wide techniques allowed the identification of thousands of non-coding RNAs with various sizes in eukaryotes, some of which have further been shown to serve important functions in many biological processes. However, in model plant Arabidopsis, novel intermediate-sized ncRNAs (im-ncRNAs) (50~300 nt) have very limited information. By using a modified isolation strategy combined with deep-sequencing technology, we identified 838 im-ncRNAs in Arabidopsis globally. More than half (58%) are new ncRNA species, mostly evolutionary divergent. Interestingly, annotated protein-coding genes with 5'-UTR-derived novel im-ncRNAs tend to be highly expressed. For intergenic im-ncRNAs, their average abundances were comparable to mRNAs in seedlings, but subsets exhibited significantly lower expression in senescing leaves. Further, intergenic im-ncRNAs were regulated by similar genetic and epigenetic mechanisms to those of protein-coding genes, and some showed developmentally regulated expression patterns. Large-scale reverse genetic screening showed that the down-regulation of a number of im-ncRNAs resulted in either obvious molecular changes or abnormal developmental phenotypes in vivo, indicating the functional importance of im-ncRNAs in plant growth and development. Together, our results demonstrate that novel Arabidopsis im-ncRNAs are developmentally regulated and functional components discovered in the transcriptome.

Identification, classification and differential expression of oleosin genes in tung tree (Vernicia fordii)

Triacylglycerols (TAG) are the major molecules of energy storage in eukaryotes. TAG are packed in subcellular structures called oil bodies or lipid droplets. Oleosins (OLE) are the major proteins in plant oil bodies. Multiple isoforms of OLE are present in plants such as tung tree (Vernicia fordii), whose seeds are rich in novel TAG with a wide range of industrial applications. The objectives of this study were to identify OLE genes, classify OLE proteins and analyze OLE gene expression in tung trees. We identified five tung tree OLE genes coding for small hydrophobic proteins. Genome-wide phylogenetic analysis and multiple sequence alignment demonstrated that the five tung OLE genes represented the five OLE subfamilies and all contained the "proline knot" motif (PX5SPX3P) shared among 65 OLE from 19 tree species, including the sequenced genomes of Prunus persica (peach), Populus trichocarpa (poplar), Ricinus communis (castor bean), Theobroma cacao (cacao) and Vitis vinifera (grapevine). Tung OLE1, OLE2 and OLE3 belong to the S type and OLE4 and OLE5 belong to the SM type of Arabidopsis OLE. TaqMan and SYBR Green qPCR methods were used to study the differential expression of OLE genes in tung tree tissues. Expression results demonstrated that 1) All five OLE genes were expressed in developing tung seeds, leaves and flowers; 2) OLE mRNA levels were much higher in seeds than leaves or flowers; 3) OLE1, OLE2 and OLE3 genes were expressed in tung seeds at much higher levels than OLE4 and OLE5 genes; 4) OLE mRNA levels rapidly increased during seed development; and 5) OLE gene expression was well-coordinated with tung oil accumulation in the seeds. These results suggest that tung OLE genes 1-3 probably play major roles in tung oil accumulation and/or oil body development. Therefore, they might be preferred targets for tung oil engineering in transgenic plants.

Comparative transcriptomics of early meiosis in Arabidopsis and maize

Though sexually reproductive plants share the same principle and most processes in meiosis, there are distinct features detectable. To address the similarities and differences of early meiosis transcriptomes from the dicot model system Arabidopsis and monocot model system maize, we performed comparative analyses of RNA-seq data of isolated meiocytes, anthers and seedlings from both species separately and via orthologous genes. Overall gene expression showed similarities, such as an increased number of reads mapping to unannotated features, and differences, such as the amount of differentially expressed genes. We detected major similarities and differences in functional annotations of genes up-regulated in meiocytes, which point to conserved features as well as unique features. Transcriptional regulation seems to be quite similar in Arabidopsis and maize, and we could reveal known and novel transcription factors and cis-regulatory elements acting in early meiosis. Taken together, meiosis between Arabidopsis and maize is conserved in many ways, but displays key distinctions that lie in the patterns of gene expression.

Characterization and expression analysis of dirigent family genes related to stresses in Brassica

The dirigent (DIR) genes are playing a vital role in enhancing stress resistance in different crop plants. In this study, we collected 29 DIR like genes, two from a Brassica rapa cv. Osome full length cDNA library and 27 from the B. rapa database designated as B. rapa Dirigent (BrDIR) like genes. Sequence analysis and a comparison study of these genes confirmed that seven were dirigent and the remaining 22 were dirigent like genes. Expression analysis revealed an organ specific expression of these genes. BrDIR2 showed differential responses after Fusarium oxysporum f.sp. conglutinans infection in cabbage. Four Brassica oleracea dirigent like genes highly homologous to BrDIR2 also showed similar responses in cabbage plants infected with this fungus. Moreover, several BrDIR like genes showed significant responses after water, ABA and cold stress treatments in Chinese cabbage. Under water stress, most responsive genes showed the highest expression at 24 h, at which time the acid soluble lignin content of samples under the same stress condition were also highest, indicating a possible relationship between BrDIR like genes and lignin content. Taken together, our results indicate a protective role of BrDIR genes against biotic and abiotic stresses in Brassica.

The Arabidopsis stress responsive gene database

Plants in nature may face a wide range of favorable or unfavorable biotic and abiotic factors during their life cycle. Any of these factors may cause stress in plants; therefore, they have to be more adaptable to stressful environments and must acquire greater response to different stresses. The objective of this study is to retrieve and arrange data from the literature in a standardized electronic format for the development of information resources on potential stress responsive genes in Arabidopsis thaliana. This provides a powerful mean for manipulation, comparison, search, and retrieval of records describing the nature of various stress responsive genes in Arabidopsis thaliana. The database is based exclusively on published stress tolerance genes associated with plants.

IFN-inducible GTPases in host cell defense

From plants to humans, the ability to control infection at the level of an individual cell-a process termed cell-autonomous immunity-equates firmly with survival of the species. Recent work has begun to unravel this programmed cell-intrinsic response and the central roles played by IFN-inducible GTPases in defending the mammalian cell's interior against a diverse group of invading pathogens. These immune GTPases regulate vesicular traffic and protein complex assembly to stimulate oxidative, autophagic, membranolytic, and inflammasome-related antimicrobial activities within the cytosol, as well as on pathogen-containing vacuoles. Moreover, human genome-wide association studies and disease-related transcriptional profiling have linked mutations in the Immunity-Related GTPase M (IRGM) locus and altered expression of guanylate binding proteins (GBPs) with tuberculosis susceptibility and Crohn's colitis.

A novel C-S lyase from the latex-producing plant Taraxacum brevicorniculatum displays alanine aminotransferase and l-cystine lyase activity

We isolated a novel pyridoxal-5-phosphate-dependent l-cystine lyase from the dandelion Taraxacum brevicorniculatum. Real time qPCR analysis showed that C-S lyase from Taraxacum brevicorniculatum (TbCSL) mRNA is expressed in all plant tissues, although at relatively low levels in the latex and pedicel. The 1251 bp TbCSL cDNA encodes a protein with a calculated molecular mass of 46,127 kDa. It is homologous to tyrosine and alanine aminotransferases (AlaATs) as well as to an Arabidopsis thaliana carbon-sulfur lyase (C-S lyase) (SUR1), which has a role in glucosinolate metabolism. TbCSL displayed in vitrol-cystine lyase and AlaAT activities of 4 and 19nkatmg(-1) protein, respectively. However, we detected no in vitro tyrosine aminotransferase (TyrAT) activity and RNAi knockdown of the enzyme had no effect on phenotype, showing that TbCSL substrates might be channeled into redundant pathways. TbCSL is in vivo localized in the cytosol and functions as a C-S lyase or an aminotransferase in planta, but the purified enzyme converts at least two substrates specifically, and can thus be utilized for further in vitro applications.

Molecular characterization of stress resistance-related chitinase genes of Brassica rapa

Brassica is an important vegetable group worldwide that is impacted by biotic and abiotic stresses. Molecular biology techniques offer the most efficient approach to address these concerns. Inducible plant defense responses include the production of pathogenesis-related (PR) proteins, and chitinases are very important PR proteins. We collected 30 chitinase like genes, three from our full-length cDNA library of Brassica rapa cv. Osome and 27 from Brassica databases. Sequence analysis and comparison study confirmed that they were all class I-V and VII chitinase genes. These genes also showed a high degree of homology with other biotic stress resistance-related plant chitinases. An organ-specific expression of these genes was observed and among these, seven genes showed significant responses after infection with Fusarium oxysporum f.sp. conglutinans in cabbage and sixteen genes showed responsive expression after abiotic stress treatments in Chinese cabbage. BrCLP1, 8, 10, 17 and 18 responded commonly after biotic and abiotic stress treatments indicating their higher potentials. Taken together, the results presented herein suggest that these chitinase genes may be useful resources in the development of stress resistant Brassica.

Advances in plant genome sequencing

The study of plant biology in the 21st century is, and will continue to be, vastly different from that in the 20th century. One driver for this has been the use of genomics methods to reveal the genetic blueprints for not one but dozens of plant species, as well as resolving genome differences in thousands of individuals at the population level. Genomics technology has advanced substantially since publication of the first plant genome sequence, that of Arabidopsis thaliana, in 2000. Plant genomics researchers have readily embraced new algorithms, technologies and approaches to generate genome, transcriptome and epigenome datasets for model and crop species that have permitted deep inferences into plant biology. Challenges in sequencing any genome include ploidy, heterozygosity and paralogy, all which are amplified in plant genomes compared to animal genomes due to the large genome sizes, high repetitive sequence content, and rampant whole- or segmental genome duplication. The ability to generate de novo transcriptome assemblies provides an alternative approach to bypass these complex genomes and access the gene space of these recalcitrant species. The field of genomics is driven by technological improvements in sequencing platforms; however, software and algorithm development has lagged behind reductions in sequencing costs, improved throughput, and quality improvements. It is anticipated that sequencing platforms will continue to improve the length and quality of output, and that the complementary algorithms and bioinformatic software needed to handle large, repetitive genomes will improve. The future is bright for an exponential improvement in our understanding of plant biology.

A reevaluation of rice mitochondrial evolution based on the complete sequence of male-fertile and male-sterile mitochondrial genomes

Plant mitochondrial genomes have features that distinguish them radically from their animal counterparts: a high rate of rearrangement, of uptake and loss of DNA sequences, and an extremely low point mutation rate. Perhaps the most unique structural feature of plant mitochondrial DNAs is the presence of large repeated sequences involved in intramolecular and intermolecular recombination. In addition, rare recombination events can occur across shorter repeats, creating rearrangements that result in aberrant phenotypes, including pollen abortion, which is known as cytoplasmic male sterility (CMS). Using next-generation sequencing, we pyrosequenced two rice (Oryza sativa) mitochondrial genomes that belong to the indica subspecies. One genome is normal, while the other carries the wild abortive-CMS. We find that numerous rearrangements in the rice mitochondrial genome occur even between close cytotypes during rice evolution. Unlike maize (Zea mays), a closely related species also belonging to the grass family, integration of plastid sequences did not play a role in the sequence divergence between rice cytotypes. This study also uncovered an excellent candidate for the wild abortive-CMS-encoding gene; like most of the CMS-associated open reading frames that are known in other species, this candidate was created via a rearrangement, is chimeric in structure, possesses predicted transmembrane domains, and coopted the promoter of a genuine mitochondrial gene. Our data give new insights into rice mitochondrial evolution, correcting previous reports.

Unscrambling genetic information at the RNA level

Genomics aims at unraveling the blueprint of life; however, DNA sequence alone does not always reveal the proteins and structural RNAs encoded by the genome. The reason is that genetic information is often encrypted. Recognizing the logic of encryption, and understanding how living cells decode hidden information--at the level of DNA, RNA or protein--is challenging. RNA-level decryption includes topical RNA editing and more 'macroscopic' transcript rearrangements. The latter events involve the four types of introns recognized to date, notably spliceosomal, group I, group II, and archaeal/tRNA splicing. Intricate variants, such as alternative splicing and trans-splicing, have been reported for each intron type, but the biological significance has not always been confirmed. Novel RNA-level unscrambling processes were recently discovered in mitochondria of dinoflagellates and diplonemids, and potentially euglenids. These processes seem not to rely on known introns, and the corresponding molecular mechanisms remain to be elucidated.

Genome-wide patterns of genetic variation in worldwide Arabidopsis thaliana accessions from the RegMap panel

Arabidopsis thaliana is native to Eurasia and naturalized across the world due to human disturbance. Its easy propagation and immense phenotypic variability make it an ideal model system for functional, ecological and evolutionary genetics. To date, analyses of its natural variation have involved small numbers of individuals or genetic markers. Here we genotype 1,307 world-wide accessions, including several regional samples, at 250K SNPs, enabling us to describe the global pattern of genetic variation with high resolution. Three complementary tests applied to these data reveal novel targets of selection. Furthermore, we characterize the pattern of historical recombination and observe an enrichment of hotspots in intergenic regions and repetitive DNA, consistent with the pattern observed for humans but strikingly different from other plant species. We are making seeds for this Regional Mapping (RegMap) panel publicly available; they comprise the largest genomic mapping resource available for a naturally occurring, non-human, species.