Search

New graduate students in biology programs may lack the quantitative skills necessary for their research and professional careers. The acquisition of these skills may be impeded by teaching and mentoring experiences that decrease rather than increase students’ beliefs in their ability to learn and apply quantitative approaches. In this opinion piece, we argue that revising instructional experiences to ensure that both student confidence and quantitative skills are enhanced may improve both educational outcomes and professional success. A few studies suggest that explicitly addressing productive failure in an instructional setting and ensuring effective mentoring may be the most effective routes to simultaneously increasing both quantitative self-efficacy and quantitative skills. However, there is little work that specifically addresses graduate student needs, and more research is required to reach evidence-backed conclusions.

The broad host range of Fusarium virguliforme represents a unique comparative system to identify and define differentially induced responses between an asymptomatic monocot host, maize (Zea mays), and a symptomatic eudicot host, soybean (Glycine max). Using a temporal, comparative transcriptome-based approach, we observed that early gene expression profiles of root tissue from infected maize suggest that pathogen tolerance coincides with the rapid induction of senescence dampening transcriptional regulators, including ANACs (Arabidopsis thaliana NAM/ATAF/CUC protein) and Ethylene-Responsive Factors. In contrast, the expression of senescence-associated processes in soybean was coincident with the appearance of disease symptom development, suggesting pathogen-induced senescence as a key pathway driving pathogen susceptibility in soybean. Based on the analyses described herein, we posit that root senescence is a primary contributing factor underlying colonization and disease progression in symptomatic versus asymptomatic host–fungal interactions. This process also supports the lifestyle and virulence of F. virguliforme during biotrophy to necrotrophy transitions. Further support for this hypothesis lies in comprehensive co-expression and comparative transcriptome analyses, and in total, supports the emerging concept of necrotrophy-activated senescence. We propose that F. virguliforme conditions an environment within symptomatic hosts, which favors susceptibility through transcriptomic reprogramming, and as described herein, the induction of pathways associated with senescence during the necrotrophic stage of fungal development.

Iron (Fe) is a key cofactor in many cellular redox processes, including respiration and photosynthesis. Plant Fe deficiency (-Fe) activates a complex regulatory network which coordinates root Fe uptake and distribution to sink tissues, while avoiding over-accumulation of Fe and other metals to toxic levels. In Arabidopsis (Arabidopsis thaliana), FIT (FER-LIKE FE DEFICIENCY-INDUCED TRANSCRIPTION FACTOR), a bHLH transcription factor (TF), is required for up-regulation of root Fe acquisition genes. However, other root and shoot -Fe-induced genes involved in Fe allocation and signaling are FIT-independent. The cis-regulatory code, i.e. the cis-regulatory elements (CREs) and their combinations that regulate plant -Fe-responses, remains largely elusive. Using Arabidopsis genome and transcriptome data, we identified over 100 putative CREs (pCREs) that were predictive of -Fe-induced up-regulation of genes in root tissue. We used large-scale in vitro TF binding data, association with FIT-dependent or FIT-independent co-expression clusters, positional bias, and evolutionary conservation to assess pCRE properties and possible functions. In addition to bHLH and MYB TFs, also B3, NAC, bZIP, and TCP TFs might be important regulators for -Fe responses. Our approach uncovered IDE1 (Iron Deficiency-responsive Element 1), a -Fe response CRE in grass species, to be conserved in regulating genes for biosynthesis of Fe-chelating compounds also in Arabidopsis. Our findings provide a comprehensive source of cis-regulatory information for -Fe-responsive genes, that advances our mechanistic understanding and informs future efforts in engineering plants with more efficient Fe uptake or transport systems.One sentence summary>100 putative cis-regulatory elements robustly predict Arabidopsis root Fe deficiency-responses in computational models, and shed light on the mechanisms of transcriptional regulation.

The usefulness of genomic prediction in crop and livestock breeding programs has prompted efforts to develop new and improved genomic prediction algorithms, such as artificial neural networks and gradient tree boosting. However, the performance of these algorithms has not been compared in a systematic manner using a wide range of datasets and models. Using data of 18 traits across six plant species with different marker densities and training population sizes, we compared the performance of six linear and six non-linear algorithms. First, we found that hyperparameter selection was necessary for all non-linear algorithms and that feature selection prior to model training was critical for artificial neural networks when the markers greatly outnumbered the number of training lines. Across all species and trait combinations, no one algorithm performed best, however predictions based on a combination of results from multiple algorithms (i.e., ensemble predictions) performed consistently well. While linear and non-linear algorithms performed best for a similar number of traits, the performance of non-linear algorithms vary more between traits. Although artificial neural networks did not perform best for any trait, we identified strategies (i.e., feature selection, seeded starting weights) that boosted their performance to near the level of other algorithms. Our results highlight the importance of algorithm selection for the prediction of trait values.

Facing adverse conditions such as nitrogen (N) deprivation, microalgae enter cellular quiescence, a reversible cell cycle arrest with drastic changes in metabolism allowing cells to remain viable. Recovering from N deprivation and quiescence is an active and orderly process as we are showing here for Chlamydomonas reinhardtii . We conducted comparative transcriptomics on this alga to discern processes relevant to quiescence in the context of N deprivation and recovery following refeeding. A mutant with slow recovery from N deprivation, compromised hydrolysis of triacylglycerols7 ( cht7 ), was included to better define the regulatory processes governing the respective transitions. We identified an ordered set of biological processes with expression patterns that showed sequential reversal following N resupply and uncovered acclimation responses specific to the recovery phase. Biochemical assays and microscopy validated selected inferences made based on the transcriptional analyses. These comprise (1) the restoration of N source preference and cellular bioenergetics during the early stage of recovery; (2) flagellum-based motility in the mid to late stage of recovery; and (3) recovery phase-specific gene groups cooperating in the rapid replenishment of chloroplast proteins. In the cht7 mutant, a large number of programmed responses failed to readjust in a timely manner. Finally, evidence is provided for the involvement of the cAMP-protein kinase A pathway in gating the recovery. We conclude that the recovery from N deprivation represents not simply a reversal of processes directly following N deprivation, but a distinct cellular state.

Plants are exposed to a variety of environmental conditions, and their ability to respond to environmental variation depends on the proper regulation of gene expression in an organ-, tissue-, and cell type-specific manner. Although our knowledge of how stress responses are regulated is accumulating, a genome-wide model of how plant transcription factors (TFs) and cis-regulatory elements control spatially specific stress response has yet to emerge. Using Arabidopsis (Arabidopsis thaliana) as a model, we identified a set of 1,894 putative cis-regulatory elements (pCREs) that are associated with high-salinity (salt) up-regulated genes in the root or the shoot. We used these pCREs to develop computational models that can better predict salt up-regulated genes in the root and shoot compared with models based on known TF binding motifs. In addition, we incorporated TF binding sites identified via large-scale in vitro assays, chromatin accessibility, evolutionary conservation, and pCRE combinatorial relationships in machine learning models and found that only consideration of pCRE combinations led to better performance in salt up-regulation prediction in the root and shoot. Our results suggest that the plant organ transcriptional response to high salinity is regulated by a core set of pCREs and provide a genome-wide view of the cis-regulatory code of plant spatial transcriptional responses to environmental stress.

Ancient duplication events and a high rate of retention of extant pairs of duplicate genes have contributed to an abundance of duplicate genes in plant genomes. These duplicates have contributed to the evolution of novel functions, such as the production of floral structures, induction of disease resistance, and adaptation to stress. Additionally, recent whole-genome duplications that have occurred in the lineages of several domesticated crop species, including wheat (Triticum aestivum), cotton (Gossypium hirsutum), and soybean (Glycine max), have contributed to important agronomic traits, such as grain quality, fruit shape, and flowering time. Therefore, understanding the mechanisms and impacts of gene duplication will be important to future studies of plants in general and of agronomically important crops in particular. In this review, we survey the current knowledge about gene duplication, including gene duplication mechanisms, the potential fates of duplicate genes, models explaining duplicate gene retention, the properties that distinguish duplicate from singleton genes, and the evolutionary impact of gene duplication.

Identification of polymorphic transposable elements (TEs) is important because TE polymorphism creates genetic diversity and influences the function of genes in the host genome. However, de novo scanning of polymorphic TEs remains a challenge. Here, we report a novel computational method, called PTEMD (polymorphic TEs and their movement detection), for de novo discovery of genome-wide polymorphic TEs. PTEMD searches highly identical sequences using reads supported breakpoint evidences. Using PTEMD, we identified 14 polymorphic TE families (905 sequences) in rice blast fungus Magnaporthe oryzae, and 68 (10,618 sequences) in maize. We validated one polymorphic TE family experimentally, MoTE-1; all MoTE-1 family members are located in different genomic loci in the three tested isolates. We found that 57.1% (8 of 14) of the PTEMD-detected polymorphic TE families in M. oryzae are active. Furthermore, our data indicate that there are more polymorphic DNA transposons in maize than their counterparts of retrotransposons despite the fact that retrotransposons occupy largest fraction of genomic mass. We demonstrated that PTEMD is an effective tool for identifying polymorphic TEs in M. oryzae and maize genomes. PTEMD and the genome-wide polymorphic TEs in M. oryzae and maize are publically available at http://www.kanglab.cn/blast/PTEMD_V1.02.htm.

PSEUDO-RESPONSE REGULATORs (PRRs) play overlapping and distinct roles in maintaining circadian rhythms and regulating diverse biological processes, including the photoperiodic control of flowering, growth, and abiotic stress responses. PRRs act as transcriptional repressors and associate with chromatin via their conserved C-terminal CCT (CONSTANS, CONSTANS-like, and TIMING OF CAB EXPRESSION 1 [TOC1/PRR1]) domains by a still-poorly understood mechanism. Here, we identified genome-wide targets of PRR9 using chromatin immunoprecipitation followed by high-throughput sequencing (ChIP-seq) and compared them with PRR7, PRR5, and TOC1/PRR1 ChIP-seq data. We found that PRR binding sites are located within genomic regions of low nucleosome occupancy and high DNase I hypersensitivity. Moreover, conserved noncoding regions among Brassicaceae species are enriched around PRR binding sites, indicating that PRRs associate with functionally relevant cis-regulatory regions. The PRRs shared a significant number of binding regions, and our results indicate that they coordinately restrict the expression of target genes to around dawn. A G-box-like motif was overrepresented at PRR binding regions, and we showed that this motif is necessary for mediating transcriptional regulation of CIRCADIAN CLOCK ASSOCIATED 1 and PRR9 by the PRRs. Our results further our understanding of how PRRs target specific promoters and provide an extensive resource for studying circadian regulatory networks in plants.

Essential genes represent critical cellular components whose disruption results in lethality. Characteristics shared among essential genes have been uncovered in fungal and metazoan model systems. However, features associated with plant essential genes are largely unknown and the full set of essential genes remains to be discovered in any plant species. Here, we show that essential genes in Arabidopsis thaliana have distinct features useful for constructing within- and cross-species prediction models. Essential genes in A. thaliana are often single copy or derived from older duplications, highly and broadly expressed, slow evolving, and highly connected within molecular networks compared with genes with nonlethal mutant phenotypes. These gene features allowed the application of machine learning methods that predicted known lethal genes as well as an additional 1970 likely essential genes without documented phenotypes. Prediction models from A. thaliana could also be applied to predict Oryza sativa and Saccharomyces cerevisiae essential genes. Importantly, successful predictions drew upon many features, while any single feature was not sufficient. Our findings show that essential genes can be distinguished from genes with nonlethal phenotypes using features that are similar across kingdoms and indicate the possibility for translational application of our approach to species without extensive functional genomic and phenomic resources.

Functional divergence between duplicate transcription factors (TFs) has been linked to critical events in the evolution of land plants and can result from changes in patterns of expression, binding site divergence, and/or interactions with other proteins. Although plant TFs tend to be retained post polyploidization, many are lost within tens to hundreds of million years. Thus, it can be hypothesized that some TFs in plant genomes are in the process of becoming pseudogenes. Here, we use a pair of salt tolerance-conferring transcription factors, DWARF AND DELAYED FLOWERING1 (DDF1) and DDF2, that duplicated through paleopolyploidy 50 to 65 million years ago, as examples to illustrate potential mechanisms leading to duplicate retention and loss. We found that the expression patterns of Arabidopsis thaliana (At)DDF1 and AtDDF2 have diverged in a highly asymmetric manner, and AtDDF2 has lost most inferred ancestral stress responses. Consistent with promoter disablement, the AtDDF2 promoter has fewer predicted cis-elements and a methylated repetitive element. Through comparisons of AtDDF1, AtDDF2, and their Arabidopsis lyrata orthologs, we identified significant differences in binding affinities and binding site preference. In particular, an AtDDF2-specific substitution within the DNA-binding domain significantly reduces binding affinity. Cross-species analyses indicate that both AtDDF1 and AtDDF2 are under selective constraint, but among A. thaliana accessions, AtDDF2 has a higher level of nonsynonymous nucleotide diversity compared with AtDDF1. This may be the result of selection in different environments or may point toward the possibility of ongoing functional decay despite retention for millions of years after gene duplication.

The large size and relative complexity of many plant genomes make creation, quality control, and dissemination of high-quality gene structure annotations challenging. In response, we have developed MAKER-P, a fast and easy-to-use genome annotation engine for plants. Here, we report the use of MAKER-P to update and revise the maize (Zea mays) B73 RefGen_v3 annotation build (5b+) in less than 3 h using the iPlant Cyberinfrastructure. MAKER-P identified and annotated 4,466 additional, well-supported protein-coding genes not present in the 5b+ annotation build, added additional untranslated regions to 1,393 5b+ gene models, identified 2,647 5b+ gene models that lack any supporting evidence (despite the use of large and diverse evidence data sets), identified 104,215 pseudogene fragments, and created an additional 2,522 noncoding gene annotations. We also describe a method for de novo training of MAKER-P for the annotation of newly sequenced grass genomes. Collectively, these results lead to the 6a maize genome annotation and demonstrate the utility of MAKER-P for rapid annotation, management, and quality control of grasses and other difficult-to-annotate plant genomes.

We have optimized and extended the widely used annotation engine MAKER in order to better support plant genome annotation efforts. New features include better parallelization for large repeat-rich plant genomes, noncoding RNA annotation capabilities, and support for pseudogene identification. We have benchmarked the resulting software tool kit, MAKER-P, using the Arabidopsis (Arabidopsis thaliana) and maize (Zea mays) genomes. Here, we demonstrate the ability of the MAKER-P tool kit to automatically update, extend, and revise the Arabidopsis annotations in light of newly available data and to annotate pseudogenes and noncoding RNAs absent from The Arabidopsis Informatics Resource 10 build. Our results demonstrate that MAKER-P can be used to manage and improve the annotations of even Arabidopsis, perhaps the best-annotated plant genome. We have also installed and benchmarked MAKER-P on the Texas Advanced Computing Center. We show that this public resource can de novo annotate the entire Arabidopsis and maize genomes in less than 3 h and produce annotations of comparable quality to those of the current The Arabidopsis Information Resource 10 and maize V2 annotation builds.

The Arabidopsis (Arabidopsis thaliana) genome is the most well-annotated plant genome. However, transcriptome sequencing in Arabidopsis continues to suggest the presence of polyadenylated (polyA) transcripts originating from presumed intergenic regions. It is not clear whether these transcripts represent novel noncoding or protein-coding genes. To understand the nature of intergenic polyA transcription, we first assessed its abundance using multiple messenger RNA sequencing data sets. We found 6,545 intergenic transcribed fragments (ITFs) occupying 3.6% of Arabidopsis intergenic space. In contrast to transcribed fragments that map to protein-coding and RNA genes, most ITFs are significantly shorter, are expressed at significantly lower levels, and tend to be more data set specific. A surprisingly large number of ITFs (32.1%) may be protein coding based on evidence of translation. However, our results indicate that these “translated” ITFs tend to be close to and are likely associated with known genes. To investigate if ITFs are under selection and are functional, we assessed ITF conservation through cross-species as well as within-species comparisons. Our analysis reveals that 237 ITFs, including 49 with translation evidence, are under strong selective constraint and relatively distant from annotated features. These ITFs are likely parts of novel genes. However, the selective pressure imposed on most ITFs is similar to that of randomly selected, untranscribed intergenic sequences. Our findings indicate that despite the prevalence of ITFs, apart from the possibility of genomic contamination, many may be background or noisy transcripts derived from “junk” DNA, whose production may be inherent to the process of transcription and which, on rare occasions, may act as catalysts for the creation of novel genes.

Two-component signaling elements play important roles in plants, including a central role in cytokinin signaling. We characterized two-component elements from the monocot rice (Oryza sativa) using several complementary approaches. Phylogenetic analysis reveals relatively simple orthologous relationships among the histidine kinases in rice and Arabidopsis (Arabidopsis thaliana). In contrast, the histidine-containing phosphotransfer proteins (OsHPs) and response regulators (OsRRs) display a higher degree of lineage-specific expansion. The intracellular localizations of several OsHPs and OsRRs were examined in rice and generally found to correspond to the localizations of their dicot counterparts. The functionality of rice type-B OsRRs was tested in Arabidopsis; one from a clade composed of both monocot and dicot type-B OsRRs complemented an Arabidopsis type-B response regulator mutant, but a type-B OsRR from a monocot-specific subfamily generally did not. The expression of genes encoding two-component elements and proteins involved in cytokinin biosynthesis and degradation was analyzed in rice roots and shoots and in response to phytohormones. Nearly all type-A OsRRs and OsHK4 were up-regulated in response to cytokinin, but other cytokinin signaling elements were not appreciably affected. Furthermore, multiple cytokinin oxidase (OsCKX) genes were up-regulated by cytokinin. Abscisic acid treatment decreased the expression of several genes involved in cytokinin biosynthesis and degradation. Auxin affected the expression of a few genes; brassinosteroid and gibberellin had only modest effects. Our results support a shared role for two-component elements in mediating cytokinin signaling in monocots and dicots and reveal how phytohormones can impact cytokinin function through modulating gene expression.

Like many microalgae, Chlamydomonas reinhardtii forms lipid droplets rich in triacylglycerols when nutrient deprived. To begin studying the mechanisms underlying this process, nitrogen (N) deprivation was used to induce triacylglycerol accumulation and changes in developmental programs such as gametogenesis. Comparative global analysis of transcripts under induced and noninduced conditions was applied as a first approach to studying molecular changes that promote or accompany triacylglycerol accumulation in cells encountering a new nutrient environment. Towards this goal, high-throughput sequencing technology was employed to generate large numbers of expressed sequence tags of eight biologically independent libraries, four for each condition, N replete and N deprived, allowing a statistically sound comparison of expression levels under the two tested conditions. As expected, N deprivation activated a subset of control genes involved in gametogenesis while down-regulating protein biosynthesis. Genes for components of photosynthesis were also down-regulated, with the exception of the PSBS gene. N deprivation led to a marked redirection of metabolism: the primary carbon source, acetate, was no longer converted to cell building blocks by the glyoxylate cycle and gluconeogenesis but funneled directly into fatty acid biosynthesis. Additional fatty acids may be produced by membrane remodeling, a process that is suggested by the changes observed in transcript abundance of putative lipase genes. Inferences on metabolism based on transcriptional analysis are indirect, but biochemical experiments supported some of these deductions. The data provided here represent a rich source for the exploration of the mechanism of oil accumulation in microalgae.

Pseudogenes (Ψ) are nonfunctional genomic sequences resembling functional genes. Knowledge of Ψs can improve genome annotation and our understanding of genome evolution. However, there has been relatively little systemic study of Ψs in plants. In this study, we characterized the evolution and expression patterns of Ψs in Arabidopsis (Arabidopsis thaliana) and rice (Oryza sativa). In contrast to animal Ψs, many plant Ψs experienced much stronger purifying selection. In addition, plant Ψs experiencing stronger selective constraints tend to be derived from relatively ancient duplicates, suggesting that they were functional for a relatively long time but became Ψs recently. Interestingly, the regions 5′ to the first stops in the Ψs have experienced stronger selective constraints compared with 3′ regions, suggesting that the 5′ regions were functional for a longer period of time after the premature stops appeared. We found that few Ψs have expression evidence, and their expression levels tend to be lower compared with annotated genes. Furthermore, Ψs with expressed sequence tags tend to be derived from relatively recent duplication events, indicating that Ψ expression may be due to insufficient time for complete degeneration of regulatory signals. Finally, larger protein domain families have significantly more Ψs in general. However, while families involved in environmental stress responses have a significant excess of Ψs, transcription factors and receptor-like kinases have lower than expected numbers of Ψs, consistent with their elevated retention rate in plant genomes. Our findings illustrate peculiar properties of plant Ψs, providing additional insight into the evolution of duplicate genes and benefiting future genome annotation.

Gene duplication is an important mechanism for evolution of new genes. In plants, a special group of transposable elements, called Pack-MULEs or transduplicates, is able to duplicate and amplify genes or gene fragments on a large scale. Despite the abundance of Pack-MULEs, the functionality of these duplicates is not clear. Here, we present a comprehensive analysis of expression and purifying selection on 2809 Pack-MULEs in rice (Oryza sativa), which are derived from 1501 parental genes. At least 22% of the Pack-MULEs are transcribed, and 28 Pack-MULEs have direct evidence of translation. Chimeric Pack-MULEs, which contain gene fragments from multiple genes, are much more frequently expressed than those derived only from a single gene. In addition, Pack-MULEs are frequently associated with small RNAs. The presence of these small RNAs is associated with a reduction in expression of both the Pack-MULEs and their parental genes. Furthermore, an assessment of the selection pressure on the Pack-MULEs using the ratio of nonsynonymous (Ka) and synonymous (Ks) substitution rates indicates that a considerable number of Pack-MULEs likely have been under selective constraint. The Ka/Ks values of Pack-MULE and parental gene pairs are lower among Pack-MULEs that are expressed in sense orientations. Taken together, our analysis suggests that a significant number of Pack-MULEs are expressed and subjected to purifying selection, and some are associated with small RNAs. Therefore, at least a subset of Pack-MULEs are likely functional and have great potential in regulating gene expression as well as providing novel coding capacities.

Plants have substantially higher gene duplication rates compared with most other eukaryotes. These plant gene duplicates are mostly derived from whole genome and/or tandem duplications. Earlier studies have shown that a large number of duplicate genes are retained over a long evolutionary time, and there is a clear functional bias in retention. However, the influence of duplication mechanism, particularly tandem duplication, on duplicate retention has not been thoroughly investigated. We have defined orthologous groups (OGs) between Arabidopsis (Arabidopsis thaliana) and three other land plants to examine the functional bias of retained duplicate genes during vascular plant evolution. Based on analysis of Gene Ontology categories, it is clear that genes in OGs that expanded via tandem duplication tend to be involved in responses to environmental stimuli, while those that expanded via nontandem mechanisms tend to have intracellular regulatory roles. Using Arabidopsis stress expression data, we further demonstrated that tandem duplicates in expanded OGs are significantly enriched in genes that are up-regulated by biotic stress conditions. In addition, tandem duplication of genes in an OG tends to be highly asymmetric. That is, expansion of OGs with tandem genes in one organismal lineage tends to be coupled with losses in the other. This is consistent with the notion that these tandem genes have experienced lineage-specific selection. In contrast, OGs with genes duplicated via nontandem mechanisms tend to experience convergent expansion, in which similar numbers of genes are gained in parallel. Our study demonstrates that the expansion of gene families and the retention of duplicates in plants exhibit substantial functional biases that are strongly influenced by the mechanism of duplication. In particular, genes involved in stress responses have an elevated probability of retention in a single-lineage fashion following tandem duplication, suggesting that these tandem duplicates are likely important for adaptive evolution to rapidly changing environments.

Selective ubiquitination of proteins is directed by diverse families of ubiquitin-protein ligases (or E3s) in plants. One important type uses Cullin-3 as a scaffold to assemble multisubunit E3 complexes containing one of a multitude of bric-a-brac/tramtrack/broad complex (BTB) proteins that function as substrate recognition factors. We previously described the 80-member BTB gene superfamily in Arabidopsis thaliana. Here, we describe the complete BTB superfamily in rice (Oryza sativa spp japonica cv Nipponbare) that contains 149 BTB domain–encoding genes and 43 putative pseudogenes. Amino acid sequence comparisons of the rice and Arabidopsis superfamilies revealed a near equal repertoire of putative substrate recognition module types. However, phylogenetic comparisons detected numerous gene duplication and/or loss events since the rice and Arabidopsis BTB lineages split, suggesting possible functional specialization within individual BTB families. In particular, a major expansion and diversification of a subset of BTB proteins containing Meprin and TRAF homology (MATH) substrate recognition sites was evident in rice and other monocots that likely occurred following the monocot/dicot split. The MATH domain of a subset appears to have evolved significantly faster than those in a smaller core subset that predates flowering plants, suggesting that the substrate recognition module in many monocot MATH-BTB E3s are diversifying to ubiquitinate a set of substrates that are themselves rapidly changing. Intriguing possibilities include pathogen proteins attempting to avoid inactivation by the monocot host.

There are 2 ways to infer selection pressures in the evolution of protein-coding genes, the nonsynonymous and synonymous substitution rate ratio (K(A)/K(S)) and the radical and conservative amino acid replacement rate ratio (K(R)/K(C)). Because the K(R)/K(C) ratio depends on the definition of radical and conservative changes in the classification of amino acids, we develop an amino acid classification that maximizes the correlation between K(A)/K(S) and K(R)/K(C). An analysis of 3,375 orthologous gene groups among 5 mammalian species shows that our classification gives a significantly higher correlation coefficient between the 2 ratios than those of existing classifications. However, there are many orthologous gene groups with a low K(A)/K(S) but a high K(R)/K(C) ratio. Examining the functions of these genes, we found an overrepresentation of functional categories related to development. To determine if the overrepresentation is stage specific, we examined the expression patterns of these genes at different developmental stages of the mouse. Interestingly, these genes are highly expressed in the early middle stage of development (blastocyst to amnion). It is commonly thought that developmental genes tend to be conservative in evolution, but some molecular changes in developmental stages should have contributed to morphological divergence in adult mammals. Therefore, we propose that the relaxed pressures indicated by the K(R)/K(C) ratio but not by K(A)/K(S) in the early middle stage of development may be important for the morphological divergence of mammals at the adult stage, whereas purifying selection detected by K(A)/K(S) occurs in the early middle developmental stage.

The ethylene binding domain (EBD) of the Arabidopsis thaliana ETR1 receptor is modeled as three membrane-spanning helices. We surveyed ethylene binding activity in different kingdoms and performed a bioinformatic analysis of the EBD. Ethylene binding is confined to land plants, Chara, and a group of cyanobacteria but is largely absent in other organisms, consistent with our finding that EBD-like sequences are overrepresented among plant and cyanobacterial species. We made amino acid substitutions in 37 partially or completely conserved residues of the EBD and assayed their effects on ethylene binding and signaling. Mutations primarily in residues in Helices I and II midregions eliminated ethylene binding and conferred constitutive signaling, consistent with the inverse-agonist model of ethylene receptor signaling and indicating that these residues define the ethylene binding pocket. The largest class of mutations, clustered near the cytoplasmic ends of Helices I and III, gave normal ethylene binding activity yet still conferred constitutive signaling. Therefore, these residues may play a role in turning off the signal transmitter domain of the receptor. By contrast, only two mutations were loss of function with respect to signaling. These findings yield insight into the structure and function of the EBD and suggest a conserved role of the EBD as a negative regulator of the signal transmitter domain.

Transcription factors (TFs), which are central to the regulation of gene expression, are usually members of multigene families. In plants, they are involved in diverse processes such as developmental control and elicitation of defense and stress responses. To investigate if differences exist in the expansion patterns of TF gene families between plants and other eukaryotes, we first used Arabidopsis (Arabidopsis thaliana) TFs to identify TF DNA-binding domains. These DNA-binding domains were then used to identify related sequences in 25 other eukaryotic genomes. Interestingly, among 19 families that are shared between animals and plants, more than 14 are larger in plants than in animals. After examining the lineage-specific expansion of TF families in two plants, eight animals, and two fungi, we found that TF families shared among these organisms have undergone much more dramatic expansion in plants than in other eukaryotes. Moreover, this elevated expansion rate of plant TF is not simply due to higher duplication rates of plant genomes but also to a higher degree of expansion compared to other plant genes. Further, in many Arabidopsis-rice (Oryza sativa) TF orthologous groups, the degree of lineage-specific expansion in Arabidopsis is correlated with that in rice. This pattern of parallel expansion is much more pronounced than the whole-genome trend in rice and Arabidopsis. The high rate of expansion among plant TF genes and their propensity for parallel expansion suggest frequent adaptive responses to selection pressure common among higher plants.

The Arabidopsis proline-rich extensin-like receptor kinase (PERK) family consists of 15 predicted receptor kinases. A comprehensive expression analysis was undertaken to identify overlapping and unique expression patterns within this family relative to their phylogeny. Three different approaches were used to study AtPERK gene family expression, and included analyses of the EST, MPSS and NASCArrays databases as well as experimental RNA blot analyses. Some of the AtPERK members were identified as tissue-specific genes while others were more broadly expressed. While in some cases there was a good association between these different expression patterns and the position of the AtPERK members in the kinase phylogeny, in other cases divergence of expression patterns was seen. The PERK expression data identified by the bioinformatics and experimental approaches were found generally to show similar trends and supported the use of data from large-scale expression studies for obtaining preliminary expression data. Thus, the bioinformatics survey for ESTs and microarrays is a powerful comprehensive approach for obtaining a genome-wide view of genes in a multigene family.

Receptor-like kinases (RLKs) belong to the large RLK/Pelle gene family, and it is known that the Arabidopsis thaliana genome contains >600 such members, which play important roles in plant growth, development, and defense responses. Surprisingly, we found that rice (Oryza sativa) has nearly twice as many RLK/Pelle members as Arabidopsis does, and it is not simply a consequence of a larger predicted gene number in rice. From the inferred phylogeny of all Arabidopsis and rice RLK/Pelle members, we estimated that the common ancestor of Arabidopsis and rice had >440 RLK/Pelles and that large-scale expansions of certain RLK/Pelle members and fusions of novel domains have occurred in both the Arabidopsis and rice lineages since their divergence. In addition, the extracellular domains have higher nonsynonymous substitution rates than the intracellular domains, consistent with the role of extracellular domains in sensing diverse signals. The lineage-specific expansions in Arabidopsis can be attributed to both tandem and large-scale duplications, whereas tandem duplication seems to be the major mechanism for recent expansions in rice. Interestingly, although the RLKs that are involved in development seem to have rarely been duplicated after the Arabidopsis–rice split, those that are involved in defense/disease resistance apparently have undergone many duplication events. These findings led us to hypothesize that most of the recent expansions of the RLK/Pelle family have involved defense/resistance-related genes.

Tyrosine kinases are important components of metazoan signaling pathways, and their mutant forms are implicated in various malignancies. Searching the sequences from the genomes of 28 eukaryotes and the GenBank, we found tyrosine kinases not only in metazoans but also in the green algae Chlamydomonas reinhardtii, the potato late blight pathogen Phytophthora infestans, and the protozoan pathogen Entamoeba histolytica, contrary to the current view that tyrosine kinases are animal-specific. Based on a phylogenetic analysis, we divided this gene family into 43 subfamilies and found that at least 19 tyrosine kinases were likely present in the common ancestor of chordates, arthropods, and nematodes. Interestingly, most of the subfamilies have conserved domain organizations among subfamily members but have undergone different degrees of expansion during the evolution of metazoans. In particular, a large number of duplications occurred in the lineage leading to the common ancestor of Tagifugu and mammals after its split from the Ciona lineage about 450 to 550 MYA. The timing of expansion coincides with proposed large-scale duplication event in the chordate lineage. Furthermore, gene losses have occurred in most subfamilies. Interestingly, different subfamilies have similar net gain rates in the chordates studied. However, the tyrosine kinases in mouse and human or in fruit fly and mosquito mostly have a one-to-one relationship between species, indicating that static periods of 90 Myr or longer in tyrosine kinase evolution have followed large expansion events.

The Arabidopsis genome was searched to identify predicted proteins containing armadillo (ARM) repeats, a motif known to mediate protein-protein interactions in a number of different animal proteins. Using domain database predictions and models generated in this study, 108 Arabidopsis proteins were identified that contained a minimum of two ARM repeats with the majority of proteins containing four to eight ARM repeats. Clustering analysis showed that the 108 predicted Arabidopsis ARM repeat proteins could be divided into multiple groups with wide differences in their domain compositions and organizations. Interestingly, 41 of the 108 Arabidopsis ARM repeat proteins contained a U-box, a motif present in a family of E3 ligases, and these proteins represented the largest class of Arabidopsis ARM repeat proteins. In 14 of these U-box/ARM repeat proteins, there was also a novel conserved domain identified in the N-terminal region. Based on the phylogenetic tree, representative U-box/ARM repeat proteins were selected for further study. RNA-blot analyses revealed that these U-box/ARM proteins are expressed in a variety of tissues in Arabidopsis. In addition, the selected U-box/ARM proteins were found to be functional E3 ubiquitin ligases. Thus, these U-box/ARM proteins represent a new family of E3 ligases in Arabidopsis.

The branched-chain amino acids (BCAAs) Leu, Ile, and Val are among nine essential amino acids that must be obtained from the diet of humans and other animals, and can be nutritionally limiting in plant foods. Despite genetic evidence of its importance in regulating seed amino acid levels, the full BCAA catabolic network is not completely understood in plants, and limited information is available regarding its regulation. In this study, transcript coexpression analyses revealed positive correlations among BCAA catabolism genes in stress, development, diurnal/circadian, and light data sets. A core subset of BCAA catabolism genes, including those encoding putative branched-chain ketoacid dehydrogenase subunits, is highly expressed during the night in plants on a diel cycle and in prolonged darkness. Mutants defective in these subunits accumulate higher levels of BCAAs in mature seeds, providing genetic evidence for their function in BCAA catabolism. In addition, prolonged dark treatment caused the mutants to undergo senescence early and overaccumulate leaf BCAAs. These results extend the previous evidence that BCAAs can be catabolized and serve as respiratory substrates at multiple steps. Moreover, comparison of amino acid profiles between mature seeds and dark-treated leaves revealed differences in amino acid accumulation when BCAA catabolism is perturbed. Together, these results demonstrate the consequences of blocking BCAA catabolism during both normal growth conditions and under energy-limited conditions.