Search

Background: Plants have complex and dynamic immune systems that have evolved to resist pathogens. Humans have worked to enhance these defenses in crops through breeding. However, many crops harbor only a fraction of the genetic diversity present in wild relatives. Increased utilization of diverse germplasm to search for desirable traits, such as disease resistance, is therefore a valuable step towards breeding crops that are adapted to both current and emerging threats. Here, we examine diversity of defense responses across four populations of the long-generation tree crop Theobroma cacao L., as well as four non-cacao Theobroma species, with the goal of identifying genetic elements essential for protection against the oomycete pathogen Phytophthora palmivora., Results: We began by creating a new, highly contiguous genome assembly for the P. palmivora-resistant genotype SCA 6 (Additional file 1: Tables S1-S5), deposited in GenBank under accessions CP139290-CP139299. We then used this high-quality assembly to combine RNA and whole-genome sequencing data to discover several genes and pathways associated with resistance. Many of these are unique, i.e., differentially regulated in only one of the four populations (diverged 40 k-900 k generations). Among the pathways shared across all populations is phenylpropanoid biosynthesis, a metabolic pathway with well-documented roles in plant defense. One gene in this pathway, caffeoyl shikimate esterase (CSE), was upregulated across all four populations following pathogen treatment, indicating its broad importance for cacao's defense response. Further experimental evidence suggests this gene hydrolyzes caffeoyl shikimate to create caffeic acid, an antimicrobial compound and known inhibitor of Phytophthora spp., Conclusions: Our results indicate most expression variation associated with resistance is unique to populations. Moreover, our findings demonstrate the value of using a broad sample of evolutionarily diverged populations for revealing the genetic bases of cacao resistance to P. palmivora. This approach has promise for further revealing and harnessing valuable genetic resources in this and other long-generation plants., (© 2024. The Author(s).)

The ability to detect several types of cancer using a non-invasive, blood-based test holds the potential to revolutionize oncology screening. We mined tumor methylation array data from the Cancer Genome Atlas (TCGA) covering 14 cancer types and identified two novel, broadly-occurring methylation markers at TLX1 and GALR1 . To evaluate their performance as a generalized blood-based screening approach, along with our previously reported methylation biomarker, ZNF154 , we rigorously assessed each marker individually or combined. Utilizing TCGA methylation data and applying logistic regression models within each individual cancer type, we found that the three-marker combination significantly increased the average area under the ROC curve (AUC) across the 14 tumor types compared to single markers ( p = 1.158 × 10 -10 ; Friedman test). Furthermore, we simulated dilutions of tumor DNA into healthy blood cell DNA and demonstrated increased AUC of combined markers across all dilution levels. Finally, we evaluated assay performance in bisulfite sequenced DNA from patient tumors and plasma, including early-stage samples. When combining all three markers, the assay correctly identified nine out of nine lung cancer plasma samples. In patient plasma from hepatocellular carcinoma, ZNF154 alone yielded the highest combined sensitivity and specificity values averaging 68% and 72%, whereas multiple markers could achieve higher sensitivity or specificity, but not both. Altogether, this study presents a comprehensive pipeline for the identification, testing, and validation of multi-cancer methylation biomarkers with a considerable potential for detecting a broad range of cancer types in patient blood samples.

Apples grown in high heat, high light, and low humidity environments are at risk for sun injury disorders like sunburn and associated crop losses. Understanding the physiological and molecular mechanisms underlying sunburn will support improvement of mitigation strategies and breeding for more resilient varieties. Numerous studies have highlighted key biochemical processes involved in sun injury, such as the phenylpropanoid and reactive oxygen species (ROS) pathways, demonstrating both enzyme activities and expression of related genes in response to sunburn conditions. Most previous studies have focused on at-harvest activity of a small number of genes in response to heat stress. Thus, it remains unclear how stress events earlier in the season affect physiology and gene expression. Here, we applied heat stress to mid-season apples in the field and collected tissue along a time course-24, 48, and 72 h following a heat stimulus-to investigate dynamic gene expression changes using a transcriptomic lens. We found a relatively small number of differentially expressed genes (DEGs) and enriched functional terms in response to heat treatments. Only a few of these belonged to pathways previously described to be involved in sunburn, such as the AsA-GSH pathway, while most DEGs had not yet been implicated in sunburn or heat stress in pome fruit., Competing Interests: Conflicts of interest The authors declare no conflict of interest., (Published by Oxford University Press on behalf of The Genetics Society of America 2023.)

Variation in translation-elongation kinetics along a transcript's coding sequence plays an important role in the maintenance of cellular protein homeostasis by regulating co-translational protein folding, localization, and maturation. Translation-elongation speed is influenced by molecular factors within mRNA and protein sequences. For example, the presence of proline in the ribosome's P- or A-site slows down translation, but the effect of other pairs of amino acids, in the context of all 400 possible pairs, has not been characterized. Here, we study Saccharomyces cerevisiae using a combination of bioinformatics, mutational experiments, and evolutionary analyses, and show that many different pairs of amino acids and their associated tRNA molecules predictably and causally encode translation rate information when these pairs are present in the A- and P-sites of the ribosome independent of other factors known to influence translation speed including mRNA structure, wobble base pairing, tripeptide motifs, positively charged upstream nascent chain residues, and cognate tRNA concentration. The fast-translating pairs of amino acids that we identify are enriched four-fold relative to the slow-translating pairs across Saccharomyces cerevisiae's proteome, while the slow-translating pairs are enriched downstream of domain boundaries. Thus, the chemical identity of amino acid pairs contributes to variability in translation rates, elongation kinetics are causally encoded in the primary structure of proteins, and signatures of evolutionary selection indicate their potential role in co-translational processes., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2020 Elsevier Ltd. All rights reserved.)

Context-dependent biological variation presents a unique challenge to the reproducibility of results in experimental animal research, because organisms' responses to experimental treatments can vary with both genotype and environmental conditions. In March 2019, experts in animal biology, experimental design and statistics convened in Blonay, Switzerland, to discuss strategies addressing this challenge. In contrast to the current gold standard of rigorous standardization in experimental animal research, we recommend the use of systematic heterogenization of study samples and conditions by actively incorporating biological variation into study design through diversifying study samples and conditions. Here we provide the scientific rationale for this approach in the hope that researchers, regulators, funders and editors can embrace this paradigm shift. We also present a road map towards better practices in view of improving the reproducibility of animal research.

Horizontal gene transfer (HGT), the movement and genomic integration of DNA across species boundaries, is commonly associated with bacteria and other microorganisms, but functional HGT (fHGT) is increasingly being recognized in heterotrophic parasitic plants that obtain their nutrients and water from their host plants through direct haustorial feeding. Here, in the holoparasitic stem parasite Cuscuta, we identify 108 transcribed and probably functional HGT events in Cuscuta campestris and related species, plus 42 additional regions with host-derived transposon, pseudogene and non-coding sequences. Surprisingly, 18 Cuscuta fHGTs were acquired from the same gene families by independent HGT events in Orobanchaceae parasites, and the majority are highly expressed in the haustorial feeding structures in both lineages. Convergent retention and expression of HGT sequences suggests an adaptive role for specific additional genes in parasite biology. Between 16 and 20 of the transcribed HGT events are inferred as ancestral in Cuscuta based on transcriptome sequences from species across the phylogenetic range of the genus, implicating fHGT in the successful radiation of Cuscuta parasites. Genome sequencing of C. campestris supports transfer of genomic DNA-rather than retroprocessed RNA-as the mechanism of fHGT. Many of the C. campestris genes horizontally acquired are also frequent sources of 24-nucleotide small RNAs that are typically associated with RNA-directed DNA methylation. One HGT encoding a leucine-rich repeat protein kinase overlaps with a microRNA that has been shown to regulate host gene expression, suggesting that HGT-derived parasite small RNAs may function in the parasite-host interaction. This study enriches our understanding of HGT by describing a parasite-host system with unprecedented gene exchange that points to convergent evolution of HGT events and the functional importance of horizontally transferred coding and non-coding sequences.

Background: Parasitic plants engage in a complex molecular dialog with potential host plants to identify a host and overcome host defenses to initiate development of the parasitic feeding organ, the haustorium, invade host tissues, and withdraw water and nutrients. While one of two critical signaling events in the parasitic plant life cycle (germination via stimulant chemicals) has been relatively well-studied, the signaling event that triggers haustorium formation remains elusive. Elucidation of this poorly understood molecular dialogue will shed light on plant-plant communication, parasitic plant physiology, and the evolution of parasitism in plants., Results: Here we present an experimental framework that develops easily quantifiable contrasts for the facultative generalist parasitic plant, Triphysaria, as it feeds across a broad range of diverse flowering plants. The contrasts, including variable parasite growth form and mortality when grown with different hosts, suggest a dynamic and host-dependent molecular dialogue between the parasite and host. Finally, by comparing transcriptome datasets from attached versus unattached parasites we gain insight into some of the physiological processes that are altered during parasitic behavior including shifts in photosynthesis-related and stress response genes., Conclusions: This work sheds light on Triphysaria's parasitic life habit and is an important step towards understanding the mechanisms of haustorium initiation factor perception, a unique form of plant-plant communication.

Reproduction in social insect societies reflects a delicate balance between cooperation and conflict over offspring production, and worker reproduction is widespread even in species showing strong reproductive skew in favor of the queen. To navigate these conflicts, workers are predicted to develop the means to estimate the queen's fecundity - potentially through behavioral and/or chemical cues - and to adjust their reproduction to maximize their fitness. Here, we introduced bumble bee, Bombus impatiens, workers to queens of different mating and reproductive status and examined worker reproduction and expression levels of two genes which were previously shown to be sensitive to the presence of the queen, vitellogenin and Krüppel-homolog 1. We further explored whether the queen's chemical secretion alone is sufficient to regulate worker reproduction, aggression and gene expression. We found that worker ovary activation was inhibited only in the presence of egg-laying queens, regardless of their mating status. Workers reared in the presence of newly-mated queens showed intermediate vitellogenin expression levels relative to workers reared with mated egg-laying and virgin queens. However, none of the whole-body chemical extracts of any of the queen treatment groups affected ovary activation, aggressive behavior, or gene expression in workers. Our findings indicate that only the presence of a freely-behaving, egg-laying queen can fully inhibit worker reproduction. It remains to be determined if workers detect differences in queen mating status and fecundity through differences in the queens' behavior alone or through the queen's behavior in concert with fertility signals.

Whole Exome Sequencing (WES) is a powerful clinical diagnostic tool for discovering the genetic basis of many diseases. A major shortcoming of WES is uneven coverage of sequence reads over the exome targets contributing to many low coverage regions, which hinders accurate variant calling. In this study, we devised two novel metrics, Cohort Coverage Sparseness (CCS) and Unevenness (U E ) Scores for a detailed assessment of the distribution of coverage of sequence reads. Employing these metrics we revealed non-uniformity of coverage and low coverage regions in the WES data generated by three different platforms. This non-uniformity of coverage is both local (coverage of a given exon across different platforms) and global (coverage of all exons across the genome in the given platform). The low coverage regions encompassing functionally important genes were often associated with high GC content, repeat elements and segmental duplications. While a majority of the problems associated with WES are due to the limitations of the capture methods, further refinements in WES technologies have the potential to enhance its clinical applications.

Horizontal gene transfer (HGT) is the transfer of genetic material across species boundaries and has been a driving force in prokaryotic evolution. HGT involving eukaryotes appears to be much less frequent, and the functional implications of HGT in eukaryotes are poorly understood. We test the hypothesis that parasitic plants, because of their intimate feeding contacts with host plant tissues, are especially prone to horizontal gene acquisition. We sought evidence of HGTs in transcriptomes of three parasitic members of Orobanchaceae, a plant family containing species spanning the full spectrum of parasitic capabilities, plus the free-living Lindenbergia Following initial phylogenetic detection and an extensive validation procedure, 52 high-confidence horizontal transfer events were detected, often from lineages of known host plants and with an increasing number of HGT events in species with the greatest parasitic dependence. Analyses of intron sequences in putative donor and recipient lineages provide evidence for integration of genomic fragments far more often than retro-processed RNA sequences. Purifying selection predominates in functionally transferred sequences, with a small fraction of adaptively evolving sites. HGT-acquired genes are preferentially expressed in the haustorium-the organ of parasitic plants-and are strongly biased in predicted gene functions, suggesting that expression products of horizontally acquired genes are contributing to the unique adaptive feeding structure of parasitic plants., Competing Interests: The authors declare no conflict of interest.

Whereas de novo assemblies of RNA-Seq data are being published for a growing number of species across the tree of life, there are currently no broadly accepted methods for evaluating such assemblies. Here we present a detailed comparison of 99 transcriptome assemblies, generated with 6 de novo assemblers including CLC, Trinity, SOAP, Oases, ABySS and NextGENe. Controlled analyses of de novo assemblies for Arabidopsis thaliana and Oryza sativa transcriptomes provide new insights into the strengths and limitations of transcriptome assembly strategies. We find that the leading assemblers generate reassuringly accurate assemblies for the majority of transcripts. At the same time, we find a propensity for assemblers to fail to fully assemble highly expressed genes. Surprisingly, the instance of true chimeric assemblies is very low for all assemblers. Normalized libraries are reduced in highly abundant transcripts, but they also lack 1000s of low abundance transcripts. We conclude that the quality of de novo transcriptome assemblies is best assessed through consideration of a combination of metrics: 1) proportion of reads mapping to an assembly 2) recovery of conserved, widely expressed genes, 3) N50 length statistics, and 4) the total number of unigenes. We provide benchmark Illumina transcriptome data and introduce SCERNA, a broadly applicable modular protocol for de novo assembly improvement. Finally, our de novo assembly of the Arabidopsis leaf transcriptome revealed ~20 putative Arabidopsis genes lacking in the current annotation.

Once a biochemical method has been devised to sample RNA or DNA of interest, sequencing can be used to identify the sampled molecules with high fidelity and low bias. High-throughput sequencing has therefore become the primary data acquisition method for many genomics studies and is being used more and more to address molecular biology questions. By applying principles of statistical experimental design, sequencing experiments can be made more sensitive to the effects under study as well as more biologically sound, hence more replicable.

The capacity of coral-dinoflagellate mutualisms to adapt to a changing climate relies in part on standing variation in host and symbiont populations, but rarely have the interactions between symbiotic partners been considered at the level of individuals. Here, we tested the importance of inter-individual variation with respect to the physiology of coral holobionts. We identified six genetically distinct Acropora palmata coral colonies that all shared the same isoclonal Symbiodinium 'fitti' dinoflagellate strain. No other Symbiodinium could be detected in host tissues. We exposed fragments of each colony to extreme cold and found that the stress-induced change in symbiont photochemical efficiency varied up to 3.6-fold depending on host genetic background. The S. 'fitti' strain was least stressed when associating with hosts that significantly altered the expression of 184 genes under cold shock; it was most stressed in hosts that only adjusted 14 genes. Key expression differences among hosts were related to redox signaling and iron availability pathways. Fine-scale interactions among unique host colonies and symbiont strains provide an underappreciated source of raw material for natural selection in coral symbioses.

Motivation: In high-dimensional testing problems π0, the proportion of null hypotheses that are true is an important parameter. For discrete test statistics, the P values come from a discrete distribution with finite support and the null distribution may depend on an ancillary statistic such as a table margin that varies among the test statistics. Methods for estimating π0 developed for continuous test statistics, which depend on a uniform or identical null distribution of P values, may not perform well when applied to discrete testing problems., Results: This article introduces a number of π0 estimators, the regression and 'T' methods that perform well with discrete test statistics and also assesses how well methods developed for or adapted from continuous tests perform with discrete tests. We demonstrate the usefulness of these estimators in the analysis of high-throughput biological RNA-seq and single-nucleotide polymorphism data., Availability and Implementation: implemented in R., (© The Author 2015. Published by Oxford University Press. All rights reserved. For Permissions, please e-mail: journals.permissions@oup.com.)

The origin of novel traits is recognized as an important process underlying many major evolutionary radiations. We studied the genetic basis for the evolution of haustoria, the novel feeding organs of parasitic flowering plants, using comparative transcriptome sequencing in three species of Orobanchaceae. Around 180 genes are upregulated during haustorial development following host attachment in at least two species, and these are enriched in proteases, cell wall modifying enzymes, and extracellular secretion proteins. Additionally, about 100 shared genes are upregulated in response to haustorium inducing factors prior to host attachment. Collectively, we refer to these newly identified genes as putative "parasitism genes." Most of these parasitism genes are derived from gene duplications in a common ancestor of Orobanchaceae and Mimulus guttatus, a related nonparasitic plant. Additionally, the signature of relaxed purifying selection and/or adaptive evolution at specific sites was detected in many haustorial genes, and may play an important role in parasite evolution. Comparative analysis of gene expression patterns in parasitic and nonparasitic angiosperms suggests that parasitism genes are derived primarily from root and floral tissues, but with some genes co-opted from other tissues. Gene duplication, often taking place in a nonparasitic ancestor of Orobanchaceae, followed by regulatory neofunctionalization, was an important process in the origin of parasitic haustoria., (© The Author 2014. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Background: Two-channel (or two-color) microarrays are cost-effective platforms for comparative analysis of gene expression. They are traditionally analysed in terms of the log-ratios (M-values) of the two channel intensities at each spot, but this analysis does not use all the information available in the separate channel observations. Mixed models have been proposed to analyse intensities from the two channels as separate observations, but such models can be complex to use and the gain in efficiency over the log-ratio analysis is difficult to quantify. Mixed models yield test statistics for the null distributions can be specified only approximately, and some approaches do not borrow strength between genes., Results: This article reformulates the mixed model to clarify the relationship with the traditional log-ratio analysis, to facilitate information borrowing between genes, and to obtain an exact distributional theory for the resulting test statistics. The mixed model is transformed to operate on the M-values and A-values (average log-expression for each spot) instead of on the log-expression values. The log-ratio analysis is shown to ignore information contained in the A-values. The relative efficiency of the log-ratio analysis is shown to depend on the size of the intraspot correlation. A new separate channel analysis method is proposed that assumes a constant intra-spot correlation coefficient across all genes. This approach permits the mixed model to be transformed into an ordinary linear model, allowing the data analysis to use a well-understood empirical Bayes analysis pipeline for linear modeling of microarray data. This yields statistically powerful test statistics that have an exact distributional theory. The log-ratio, mixed model and common correlation methods are compared using three case studies. The results show that separate channel analyses that borrow strength between genes are more powerful than log-ratio analyses. The common correlation analysis is the most powerful of all., Conclusions: The common correlation method proposed in this article for separate-channel analysis of two-channel microarray data is no more difficult to apply in practice than the traditional log-ratio analysis. It provides an intuitive and powerful means to conduct analyses and make comparisons that might otherwise not be possible.

Coral populations have declined worldwide largely due to increased sea surface temperatures. Recovery of coral populations depends in part upon larval recruitment. Many corals reproduce during the warmest time of year when further increases in temperature can lead to low fertilization rates of eggs and high larval mortality. Microarray experiments were designed to capture and assess variability in the thermal stress responses of Acropora palmata larvae from Puerto Rico. Transcription profiles showed a striking acceleration of normal developmental gene expression patterns with increased temperature. The transcriptional response to heat suggested rapid depletion of larval energy stores via peroxisomal lipid oxidation and included key enzymes that indicated the activation of the glyoxylate cycle. High temperature also resulted in expression differences in key developmental signalling genes including the conserved WNT pathway that is critical for pattern formation and tissue differentiation in developing embryos. Expression of these and other important developmental and thermal stress genes such as ferritin, heat shock proteins, cytoskeletal components, cell adhesion and autophagy proteins also varied among larvae derived from different parent colonies. Disruption of normal developmental and metabolic processes will have negative impacts on larval survival and dispersal as temperatures rise. However, it appears that variation in larval response to high temperature remains despite the dramatic population declines. Further research is needed to determine whether this variation is heritable or attributable to maternal effects., (© 2013 Blackwell Publishing Ltd.)

Background: Orobanchaceae is the only plant family with members representing the full range of parasitic lifestyles plus a free-living lineage sister to all parasitic lineages, Lindenbergia. A generalist member of this family, and an important parasitic plant model, Triphysaria versicolor regularly feeds upon a wide range of host plants. Here, we compare de novo assembled transcriptomes generated from laser micro-dissected tissues at the host-parasite interface to uncover details of the largely uncharacterized interaction between parasitic plants and their hosts., Results: The interaction of Triphysaria with the distantly related hosts Zea mays and Medicago truncatula reveals dramatic host-specific gene expression patterns. Relative to above ground tissues, gene families are disproportionally represented at the interface including enrichment for transcription factors and genes of unknown function. Quantitative Real-Time PCR of a T. versicolor β-expansin shows strong differential (120x) upregulation in response to the monocot host Z. mays; a result that is concordant with our read count estimates. Pathogenesis-related proteins, other cell wall modifying enzymes, and orthologs of genes with unknown function (annotated as such in sequenced plant genomes) are among the parasite genes highly expressed by T. versicolor at the parasite-host interface., Conclusions: Laser capture microdissection makes it possible to sample the small region of cells at the epicenter of parasite host interactions. The results of our analysis suggest that T. versicolor's generalist strategy involves a reliance on overlapping but distinct gene sets, depending upon the host plant it is parasitizing. The massive upregulation of a T. versicolor β-expansin is suggestive of a mechanism for parasite success on grass hosts. In this preliminary study of the interface transcriptomes, we have shown that T. versicolor, and the Orobanchaceae in general, provide excellent opportunities for the characterization of plant genes with unknown functions.

The origin and rapid diversification of the angiosperms (Darwin's "Abominable Mystery") has engaged generations of researchers. Here, we examine the floral genetic programs of phylogenetically pivotal angiosperms (water lily, avocado, California poppy, and Arabidopsis) and a nonflowering seed plant (a cycad) to obtain insight into the origin and subsequent evolution of the flower. Transcriptional cascades with broadly overlapping spatial domains, resembling the hypothesized ancestral gymnosperm program, are deployed across morphologically intergrading organs in water lily and avocado flowers. In contrast, spatially discrete transcriptional programs in distinct floral organs characterize the more recently derived angiosperm lineages represented by California poppy and Arabidopsis. Deep evolutionary conservation in the genetic programs of putatively homologous floral organs traces to those operating in gymnosperm reproductive cones. Female gymnosperm cones and angiosperm carpels share conserved genetic features, which may be associated with the ovule developmental program common to both organs. However, male gymnosperm cones share genetic features with both perianth (sterile attractive and protective) organs and stamens, supporting the evolutionary origin of the floral perianth from the male genetic program of seed plants.

Current understanding of floral developmental genetics comes primarily from the core eudicot model Arabidopsis thaliana. Here, we explore the floral transcriptome of the basal angiosperm, Nuphar advena (water lily), for insights into the ancestral developmental program of flowers. We identify several thousand Nuphar genes with significantly upregulated floral expression, including homologs of the well-known ABCE floral regulators, deployed in broadly overlapping transcriptional programs across floral organ categories. Strong similarities in the expression profiles of different organ categories in Nuphar flowers are shared with the magnoliid Persea americana (avocado), in contrast to the largely organ-specific transcriptional cascades evident in Arabidopsis, supporting the inference that this is the ancestral condition in angiosperms. In contrast to most eudicots, floral organs are weakly differentiated in Nuphar and Persea, with staminodial intermediates between stamens and perianth in Nuphar, and between stamens and carpels in Persea. Consequently, the predominantly organ-specific transcriptional programs that characterize Arabidopsis flowers (and perhaps other eudicots) are derived, and correlate with a shift towards morphologically distinct floral organs, including differentiated sepals and petals, and a perianth distinct from stamens and carpels. Our findings suggest that the genetic regulation of more spatially discrete transcriptional programs underlies the evolution of floral morphology., (© 2010 The Authors. The Plant Journal © 2010 Blackwell Publishing Ltd.)

Background: Molecular genetic studies of floral development have concentrated on several core eudicots and grasses (monocots), which have canalized floral forms. Basal eudicots possess a wider range of floral morphologies than the core eudicots and grasses and can serve as an evolutionary link between core eudicots and monocots, and provide a reference for studies of other basal angiosperms. Recent advances in genomics have enabled researchers to profile gene activities during floral development, primarily in the eudicot Arabidopsis thaliana and the monocots rice and maize. However, our understanding of floral developmental processes among the basal eudicots remains limited., Results: Using a recently generated expressed sequence tag (EST) set, we have designed an oligonucleotide microarray for the basal eudicot Eschscholzia californica (California poppy). We performed microarray experiments with an interwoven-loop design in order to characterize the E. californica floral transcriptome and to identify differentially expressed genes in flower buds with pre-meiotic and meiotic cells, four floral organs at preanthesis stages (sepals, petals, stamens and carpels), developing fruits, and leaves., Conclusions: Our results provide a foundation for comparative gene expression studies between eudicots and basal angiosperms. We identified whorl-specific gene expression patterns in E. californica and examined the floral expression of several gene families. Interestingly, most E. californica homologs of Arabidopsis genes important for flower development, except for genes encoding MADS-box transcription factors, show different expression patterns between the two species. Our comparative transcriptomics study highlights the unique evolutionary position of E. californica compared with basal angiosperms and core eudicots.

Brain structure and function experience dramatic changes from embryonic to postnatal development. Microarray analyses have detected differential gene expression at different stages and in disease models, but gene expression information during early brain development is limited. We have generated >27 million reads to identify mRNAs from the mouse cortex for >16,000 genes at either embryonic day 18 (E18) or postnatal day 7 (P7), a period of significant synaptogenesis for neural circuit formation. In addition, we devised strategies to detect alternative splice forms and uncovered more splice variants. We observed differential expression of 3,758 genes between the 2 stages, many with known functions or predicted to be important for neural development. Neurogenesis-related genes, such as those encoding Sox4, Sox11, and zinc-finger proteins, were more highly expressed at E18 than at P7. In contrast, the genes encoding synaptic proteins such as synaptotagmin, complexin 2, and syntaxin were up-regulated from E18 to P7. We also found that several neurological disorder-related genes were highly expressed at E18. Our transcriptome analysis may serve as a blueprint for gene expression pattern and provide functional clues of previously unknown genes and disease-related genes during early brain development.

The debate on the origin and evolution of flowers has recently entered the field of developmental genetics, with focus on the design of the ancestral floral regulatory program. Flowers can differ dramatically among angiosperm lineages, but in general, male and female reproductive organs surrounded by a sterile perianth of sepals and petals constitute the basic floral structure. However, the basal angiosperm lineages exhibit spectacular diversity in the number, arrangement, and structure of floral organs, whereas the evolutionarily derived monocot and eudicot lineages share a far more uniform floral ground plan. Here we show that broadly overlapping transcriptional programs characterize the floral transcriptome of the basal angiosperm Persea americana (avocado), whereas floral gene expression domains are considerably more organ specific in the model eudicot Arabidopsis thaliana. Our findings therefore support the "fading borders" model for organ identity determination in basal angiosperm flowers and extend it from the action of regulatory genes to downstream transcriptional programs. Furthermore, the declining expression of components of the staminal transcriptome in central and peripheral regions of Persea flowers concurs with elements of a previous hypothesis for developmental regulation in a gymnosperm "floral progenitor." Accordingly, in contrast to the canalized organ-specific regulatory apparatus of Arabidopsis, floral development may have been originally regulated by overlapping transcriptional cascades with fading gradients of influence from focal to bordering organs.

We develop methods to compare the positions of quantitative trait loci (QTL) with a set of genes selected by other methods, such as microarray experiments, from a sequenced genome. We apply our methods to QTL for addictive behavior in mouse, and a set of genes upregulated in a region of the brain associated with addictive behavior, the nucleus accumbens (NA). The association between the QTL and NA genes is not significantly stronger than expected by chance. However, chromosomes 2 and 16 do show strong associations suggesting that genes on these chromosomes might be associated with addictive behavior. The statistical methodology developed for this study can be applied to similar studies to assess the mutual information in microarray and QTL analyses.

The PlantTribes database (http://fgp.huck.psu.edu/tribe.html) is a plant gene family database based on the inferred proteomes of five sequenced plant species: Arabidopsis thaliana, Carica papaya, Medicago truncatula, Oryza sativa and Populus trichocarpa. We used the graph-based clustering algorithm MCL [Van Dongen (Technical Report INS-R0010 2000) and Enright et al. (Nucleic Acids Res. 2002; 30: 1575-1584)] to classify all of these species' protein-coding genes into putative gene families, called tribes, using three clustering stringencies (low, medium and high). For all tribes, we have generated protein and DNA alignments and maximum-likelihood phylogenetic trees. A parallel database of microarray experimental results is linked to the genes, which lets researchers identify groups of related genes and their expression patterns. Unified nomenclatures were developed, and tribes can be related to traditional gene families and conserved domain identifiers. SuperTribes, constructed through a second iteration of MCL clustering, connect distant, but potentially related gene clusters. The global classification of nearly 200 000 plant proteins was used as a scaffold for sorting approximately 4 million additional cDNA sequences from over 200 plant species. All data and analyses are accessible through a flexible interface allowing users to explore the classification, to place query sequences within the classification, and to download results for further study.

Through multifaceted genome-scale research involving phylogenomics, targeted gene surveys, and gene expression analyses in diverse basal lineages of angiosperms, our studies provide insights into the most recent common ancestor of all extant flowering plants. MADS-box gene duplications have played an important role in the origin and diversification of angiosperms. Furthermore, early angiosperms possessed a diverse tool kit of floral genes and exhibited developmental 'flexibility', with broader patterns of expression of key floral organ identity genes than are found in eudicots. In particular, homologs of B-function MADS-box genes are more broadly expressed across the floral meristem in basal lineages. These results prompted formulation of the 'fading borders' model, which states that the gradual transitions in floral organ morphology observed in some basal angiosperms (e.g. Amborella) result from a gradient in the level of expression of floral organ identity genes across the developing floral meristem.

The loop design of Kerr and Churchill is a clever application of incomplete blocks of size 2 to two-channel microarray experiments. In this paper, we extend the loop design to include more replicates, biological and technical replication, multi-factor experiments, and blocking. Loop and extended loop designs are shown to be more efficient than the reference design for any given number of arrays. We also show that adding new treatments to a loop design requires the same number of additional arrays as adding treatments to a reference design, with a greater gain in power. Given the flexibility of extended loop designs and their power, we propose that these should be the designs of choice for most experiments using two-channel microarrays.

Combined blood pool and delayed images produced by use of 99mTc-methylene diphosphonate (99mTcMDP) were evaluated as an objective measurement of the response of equine joints with osteochondral defects to postoperative exercise and intra-articularly administered polysulfated glycosaminoglycan (PSGAG). Osteochondral defects (approx 2.4 x 0.9 cm) were induced arthroscopically in the dorsodistal radial carpal bones of 18 ponies. These ponies were randomized (while balancing for age [range 2 to 15; median, 5.0; mean, 5.1 years]) to 2 treatment groups. Nine ponies were assigned to be exercised, and 9 were stall-rested. Six ponies in each group were administered PSGAG (250 mg) in 1 joint (medicated) and lactated Ringer's solution (LRS) in the contralateral joint. The 3 remaining ponies in each group were administered LRS in both joints (nonmedicated). Medication was given at surgery, then weekly for 4 weeks. The exercise protocol (begun at postoperative day 6 and conducted twice daily) started with 30 minutes walking (approx 0.7 m/s), and, by postoperative month 3, the ponies were being walked for 15 minutes and trotted (approx 1.6 m/s) for 25 minutes. Simultaneous dorsal images of both carpi were made 2 to 3 minutes after IV administration of 99mTcMDP (blood pool image) and 90 to 120 minutes later (delayed image). Scintimetry, in counts per minute per pixel per millicurie, was done before, and at 1, 2, 4, 8, 10, 13, and 17 weeks after surgery, prior to euthanasia. Radionuclide uptake on blood pool images decreased faster than that on delayed images, in which uptake remained high for 17 weeks. This indicated that bone was metabolically active for at least 17 weeks after surgery. Exercise significantly (P < 0.05) decreased uptake on the blood pool images of medicated joints up to 1 month after surgery. Thus, exercise (in the presence of PSGAG) probably had a transient, beneficial effect on soft tissues of the joint. Exercise, without PSGAG, promoted increased bone remodeling, because the highest uptake on delayed images was observed in exercised, nonmedicated ponies up to 3 months after surgery. This was consistent with development of osteoarthritis in these ponies. Medication alone stimulated bone remodeling, and data indicated that an identical effect may take place in contralateral LRS-injected joints, because of systemic circulation of the drug. However, the combination of exercise and medication appeared to moderate the independent effects of each. The combination of exercise and medication in individual joints resulted in notably (P < 0.05) decreased bone remodeling.(ABSTRACT TRUNCATED AT 400 WORDS)

The objective of this work was to compare between and within sample differences in aerobic plate counts (APC) and Salmonella in chicken held at abusive temperatures under air and modified atmosphere (MA). Three independent samples and three subsamples of ground raw chicken meat were inoculated with a nalidixic acid-resistant Salmonella enteritidis var. Rochester, packaged in barrier films containing air or 75% CO2, 20% N2, and 5% O2, and stored at 13 and 27 degrees C for 120 h. Samples were quantitatively enumerated for APC and S. enteritidis colony forming units (cfu/g). Gompertz function parameters for microbial growth in each sample and subsample were statistically compared in order to determine if within sample variation was different to between sample variation. At 13 degrees C, APC organisms averaged 1.27 log units lower in the MA compared to air while Salmonella counts averaged 1.44 log units higher. At 27 degrees C APC averaged 0.46 log units lower in the MA compared to air while Salmonella averaged 0.30 log units lower. The data fit the Gompertz function such that no systematic deviations from the curves were detected and the confidence bands were relatively narrow. Except for initial microbial load, there were no significant differences in variance of Gompertz parameters between the subsamples and the independent samples. This indicates that samples from different chickens with the same initial APC can be expected to be as similar in growth patterns and final counts as subsamples from the same chicken.

The development of Listeria monocytogenes Scott A, serotype 4 and aerobic plate counts on minced raw chicken were determined independently at 4, 10 and 27 degrees C. Samples were packaged in flexible film under two modified atmospheres (one containing oxygen and one containing no oxygen) or air. The anaerobic modified atmosphere (75:25, CO2:N2) resulted in the failure of both the aerobic plate counts and L. monocytogenes to grow at all temperatures. Both the L. monocytogenes and aerobic plate counts grew in air at all temperatures. The aerobic modified atmosphere (72.5:22.5:5, CO2:N2:O2), which more closely duplicates commercial practice, inhibited the increase in aerobic plate counts by more than 4 log10 cfu/g compared to air at 4 degrees C. However, the L. monocytogenes was not affected by this atmosphere and increased in numbers by nearly 6 log10 cfu/g at 4 degrees C in 21 days. Regression analysis of the log10 growth and 95% confidence intervals showed that the differences between aerobic plate counts and L. monocytogenes in modified atmosphere were large. The ability of L. monocytogenes to grow in the aerobic modified atmosphere was not affected by level of the L. monocytogenes inoculum nor by the initial level of aerobic plate counts. These data show that modified atmosphere packaging of raw chicken (and probably other meats) can substantially inhibit the aerobic spoilage flora while allowing pathogenic L. monocytogenes to increase.