Spatio-temporal expression patterns of Arabidopsis thaliana and Medicago truncatula defensin-like genes

Plant genomes contain several hundred defensin-like (DEFL) genes that encode short cysteine-rich proteins resembling defensins, which are well known antimicrobial polypeptides. Little is known about the expression patterns or functions of many DEFLs because most were discovered recently and hence are not well represented on standard microarrays. We designed a custom Affymetrix chip consisting of probe sets for 317 and 684 DEFLs from Arabidopsis thaliana and Medicago truncatula, respectively for cataloging DEFL expression in a variety of plant organs at different developmental stages and during symbiotic and pathogenic associations. The microarray analysis provided evidence for the transcription of 71% and 90% of the DEFLs identified in Arabidopsis and Medicago, respectively, including many of the recently annotated DEFL genes that previously lacked expression information. Both model plants contain a subset of DEFLs specifically expressed in seeds or fruits. A few DEFLs, including some plant defensins, were significantly up-regulated in Arabidopsis leaves inoculated with Alternaria brassicicola or Pseudomonas syringae pathogens. Among these, some were dependent on jasmonic acid signaling or were associated with specific types of immune responses. There were notable differences in DEFL gene expression patterns between Arabidopsis and Medicago, as the majority of Arabidopsis DEFLs were expressed in inflorescences, while only a few exhibited root-enhanced expression. By contrast, Medicago DEFLs were most prominently expressed in nitrogen-fixing root nodules. Thus, our data document salient differences in DEFL temporal and spatial expression between Arabidopsis and Medicago, suggesting distinct signaling routes and distinct roles for these proteins in the two plant species

Using the R Package crlmm for Genotyping and Copy Number Estimation

Genotyping platforms such as Affymetrix can be used to assess genotype-phenotype as well as copy number-phenotype associations at millions of markers. While genotyping algorithms are largely concordant when assessed on HapMap samples, tools to assess copy number changes are more variable and often discordant. One explanation for the discordance is that copy number estimates are susceptible to systematic differences between groups of samples that were processed at different times or by different labs. Analysis algorithms that do not adjust for batch effects are prone to spurious measures of association. The R package crlmm implements a multilevel model that adjusts for batch effects and provides allele-specific estimates of copy number. This paper illustrates a workflow for the estimation of allele-specific copy number and integration of the marker-level estimates with complimentary Bioconductor software for inferring regions of copy number gain or loss. All analyses are performed in the statistical environment R.