Epigenomic Diversity in a Global Collection of Arabidopsis thaliana Accessions.

The epigenome orchestrates genome accessibility, functionality, and three-dimensional structure. Because epigenetic variation can impact transcription and thus phenotypes, it may contribute to adaptation. Here, we report 1,107 high-quality single-base resolution methylomes and 1,203 transcriptomes from the 1001 Genomes collection of Arabidopsis thaliana. Although the genetic basis of methylation variation is highly complex, geographic origin is a major predictor of genome-wide DNA methylation levels and of altered gene expression caused by epialleles. Comparison to cistrome and epicistrome datasets identifies associations between transcription factor binding sites, methylation, nucleotide variation, and co-expression modules. Physical maps for nine of the most diverse genomes reveal how transposons and other structural variants shape the epigenome, with dramatic effects on immunity genes. The 1001 Epigenomes Project provides a comprehensive resource for understanding how variation in DNA methylation contributes to molecular and non-molecular phenotypes in natural populations of the most studied model plant.

Conditional GWAS of non-CG transposon methylation in Arabidopsis thaliana reveals major polymorphisms in five genes.

Genome-wide association studies (GWAS) have revealed that the striking natural variation for DNA CHH-methylation (mCHH; H is A, T, or C) of transposons has oligogenic architecture involving major alleles at a handful of known methylation regulators. Here we use a conditional GWAS approach to show that CHG-methylation (mCHG) has a similar genetic architecture-once mCHH is statistically controlled for. We identify five key trans-regulators that appear to modulate mCHG levels, and show that they interact with a previously identified modifier of mCHH in regulating natural transposon mobilization.

Beyond the standard GWAS - a guide for plant biologists.

Classic genome-wide association studies (GWAS) look for associations between individual SNPs and phenotypes of interest. With the rapid progress of high-throughput genotyping and phenotyping technologies, GWAS have become increasingly powerful for detecting genetic determinants and their molecular mechanisms underpinning natural phenotypic variation. However, GWAS frequently yield results with neither expected nor promising loci, nor any significant associations. This is often because associations between SNPs and a single phenotype are confounded, for example with the environment, other traits, or complex genetic structures. Such confounding can mask true genotype-phenotype associations, or inflate spurious associations. To address these problems, numerous methods have been developed that go beyond the standard model. Such advanced GWAS models are flexible and can offer improved statistical power for understanding the genetics underlying complex traits. Despite this advantage, these models have not been widely adopted and implemented compared to the standard GWAS approach, partly because this literature is diverse and often technical. In this review, our aim is to provide an overview of the application and the benefits of various advanced GWAS models for handling complex traits and genetic structures, targeting plant biologists who wish to carry out GWAS more effectively.

A cross-scale approach to unravel the molecular basis of plant phenology in temperate and tropical climates.

Plants have evolved to time their leafing, flowering and fruiting in appropriate seasons for growth, reproduction and resting. As a consequence of their adaptation to geographically different environments, there is a rich diversity in plant phenology from temperate and tropical climates. Recent progress in genetic and molecular studies will provide numerous opportunities to study the genetic basis of phenological traits and the history of adaptation of phenological traits to seasonal and aseasonal environments. Integrating molecular data with long-term phenology and climate data into predictive models will be a powerful tool to forecast future phenological changes in the face of global environmental change. Here, we review the cross-scale approach from genes to plant communities from three aspects: the latitudinal gradient of plant phenology at the community level, the environmental and genetic factors underlying the diversity of plant phenology, and an integrated approach to forecast future plant phenology based on genetically informed knowledge. Synthesizing the latest knowledge about plant phenology from molecular, ecological and mathematical perspectives will help us understand how natural selection can lead to the further evolution of the gene regulatory mechanisms in phenological traits in future forest ecosystems.

Common alleles of CMT2 and NRPE1 are major determinants of CHH methylation variation in Arabidopsis thaliana.

DNA cytosine methylation is an epigenetic mark associated with silencing of transposable elements (TEs) and heterochromatin formation. In plants, it occurs in three sequence contexts: CG, CHG, and CHH (where H is A, T, or C). The latter does not allow direct inheritance of methylation during DNA replication due to lack of symmetry, and methylation must therefore be re-established every cell generation. Genome-wide association studies (GWAS) have previously shown that CMT2 and NRPE1 are major determinants of genome-wide patterns of TE CHH methylation. Here we instead focus on CHH methylation of individual TEs and TE-families, allowing us to identify the pathways involved in CHH methylation simply from natural variation and confirm the associations by comparing them with mutant phenotypes. Methylation at TEs targeted by the RNA-directed DNA methylation (RdDM) pathway is unaffected by CMT2 variation, but is strongly affected by variation at NRPE1, which is largely responsible for the longitudinal cline in this phenotype. In contrast, CMT2-targeted TEs are affected by both loci, which jointly explain 7.3% of the phenotypic variation (13.2% of total genetic effects). There is no longitudinal pattern for this phenotype, however, because the geographic patterns appear to compensate for each other in a pattern suggestive of stabilizing selection.

Revisiting a GWAS peak in Arabidopsis thaliana reveals possible confounding by genetic heterogeneity.

Genome-wide association studies (GWAS) have become a standard approach for exploring the genetic basis of phenotypic variation. However, correlation is not causation, and only a tiny fraction of all associations have been experimentally confirmed. One practical problem is that a peak of association does not always pinpoint a causal gene, but may instead be tagging multiple causal variants. In this study, we reanalyze a previously reported peak associated with flowering time traits in Swedish Arabidopsis thaliana population. The peak appeared to pinpoint the AOP2/AOP3 cluster of glucosinolate biosynthesis genes, which is known to be responsible for natural variation in herbivore resistance. Here we propose an alternative hypothesis, by demonstrating that the AOP2/AOP3 flowering association can be wholly accounted for by allelic variation in two flanking genes with clear roles in regulating flowering: NDX1, a regulator of the main flowering time controller FLC, and GA1, which plays a central role in gibberellin synthesis and is required for flowering under some conditions. In other words, we propose that the AOP2/AOP3 flowering-time association may be yet another example of a spurious, "synthetic" association, arising from trying to fit a single-locus model in the presence of two statistically associated causative loci. We conclude that caution is needed when using GWAS for fine-mapping.

"Missing" G x E Variation Controls Flowering Time in Arabidopsis thaliana.

Understanding how genetic variation interacts with the environment is essential for understanding adaptation. In particular, the life cycle of plants is tightly coordinated with local environmental signals through complex interactions with the genetic variation (G x E). The mechanistic basis for G x E is almost completely unknown. We collected flowering time data for 173 natural inbred lines of Arabidopsis thaliana from Sweden under two growth temperatures (10°C and 16°C), and observed massive G x E variation. To identify the genetic polymorphisms underlying this variation, we conducted genome-wide scans using both SNPs and local variance components. The SNP-based scan identified several variants that had common effects in both environments, but found no trace of G x E effects, whereas the scan using local variance components found both. Furthermore, the G x E effects appears to be concentrated in a small fraction of the genome (0.5%). Our conclusion is that G x E effects in this study are mostly due to large numbers of allele or haplotypes at a small number of loci, many of which correspond to previously identified flowering time genes.

The genetic architecture of freezing tolerance varies across the range of Arabidopsis thaliana.

The capacity to tolerate freezing temperatures limits the geographical distribution of many plants, including several species of agricultural importance. However, the genes involved in freezing tolerance remain largely unknown. Here, we describe the variation in constitutive freezing tolerance that occurs among worldwide accessions of Arabidopsis thaliana. We found that although plants from high latitudes tend to be more freezing tolerant than plants from low latitudes, the environmental factors that shape cold adaptation differ across the species range. Consistent with this, we found that the genetic architecture of freezing tolerance also differs across its range. Conventional genome-wide association studies helped identify a priori and other promising candidate genes. However, simultaneously modelling climate variables and freezing tolerance together pinpointed other excellent a priori candidate genes. This suggests that if the selective factor underlying phenotypic variation is known, multi-trait mixed models may aid in identifying the genes that underlie adaptation.

Whole genotype constellation of prototype feline rotavirus strains FRV-1 and FRV64 and their phylogenetic relationships with feline-like human rotavirus strains.

Feline rotaviruses, members of the species Rotavirus A, are an infrequent source of zoonotic infections, and were previously shown by RNA-RNA hybridization assays to possess two distinct genomic RNA constellations, represented by strains FRV-1 and FRV64. Due to the lack of whole genome sequence information for FRV-1, human rotavirus strain AU-1 has been used as a surrogate for the genotype constellation of feline rotaviruses. The aim of this study was to determine the whole genome sequence of FRV-1 and FRV64 to help understand the genetic relationships among existing feline rotaviruses from the evolutionary perspective. The genotype constellations of FRV-1 and FRV64 were G3-P[9]-I3-R3-C3-M3-A3-N3-T3-E3-H3 and G3-P[3]-I3-R3-C2-M3-A9-N2-T3-E3-H6, respectively. FRV-1 has a genotype constellation identical to that of the AU-1 strain. Although for individual genes they shared lineages, with the exception of genes encoding VP2, VP6 and VP7, the sequence identity between FRV-1 and AU-1 was considered to be sufficiently high for the AU-1 to be regarded as an example of the direct transmission of a feline rotavirus to a child. On the other hand, the FRV64 strain was not only similar in all the 11 genome segments to another feline rotavirus strain, Cat97, but also to canine rotavirus strains (K9 and CU-1) and feline/canine-like human rotavirus strains (Ro1845 and HCR3A). In conclusion, this study revealed intermingled sharing of genotypes and lineages among feline rotaviruses, suggesting the occurrence of frequent reassortment events over the course of evolution to emerge in four genotype constellations represented by FRV-1, FRV64/Cat97, Cat2 and BA222 strains.

Molecular identification of a G2 rotavirus that provided a G1P[4] mono-reassortant with a DS-1-like genotype constellation.

The segmented nature of the rotavirus genome provides an opportunity for the virus to reassort upon co-infection of more than one rotavirus strain. Previously, two G1P[4] strains isolated from children with diarrhoea, AU64, and AU67, were shown by RNA-RNA hybridization to be a VP7 mono-reassortant possessing a DS-1 genogroup background. However, the origin of the parental G2 strain was not sought for at that time. The aim of this study, therefore, was to identify the G2 strain that provided AU64/AU67 with the DS-1-like genotype constellation. By carefully comparing the genomic RNA migration patterns on polyacrylamide gels of strains circulating at the same and nearby geographic locations around the time of the mono-reassortant detection (1989), a G2 strain designated 88H449 was shown to possess 10 genome segments co-migrating with those of AU64. Sequencing of all genome segments of AU64 and 88H449 showed that those two strains were 99.6-100% identical at the nucleotide level in all except the VP7 gene, indicating that the parental G2 strain that provided AU64 with the 10 genome segments was a sibling of 88H449. Sequencing of the VP7 genes of 36 G1P[8] strains circulating locally between 1981 and 1990 revealed the presence of strains bearing diverse VP7 sequences among which the closest nucleotide identity to AU64 was 98.6% and the closest amino acid identity in the major antigenic regions was 100%. This unusual mono-reassortant was concluded to be the result of a contemporary reassosrtment event occurring between locally co-circulating strains.

Caregivers' understanding of pediatric medication in central Malawi.

OBJECTIVE: To investigate caregivers' understanding of prescriptions of antibiotics and antimalarials for children in Malawi. METHOD: A cross-sectional study using interviews with 513 caregivers was undertaken in three different medical settings in two districts in Malawi for their understanding of medication defined as comprehension of dosage, frequency and duration immediately after receiving medicine, and its association with their sociodemographic characteristics. RESULTS: Nearly half (49%) of the caregivers understood medication on the day of prescription. Understanding varied by type of drug. Sixty-five percent understood medication for antimalarials, LA (lumefantrine and artemether), while only 38% understood medication for amoxicillin syrup. Multivariate analysis showed that caregivers with higher education and older children had greater understanding. CONCLUSIONS: Caregivers' understanding of pediatric medication was differentially related to kinds and types of medicine as well as their education level. The instruction to caregivers on how to administer medication should be evaluated especially when they are less educated.

The linkage method: a novel approach for SNP detection and haplotype reconstruction from a single diploid individual using next-generation sequence data.

When we sequence a diploid individual, the output actually comprises two genomes: one from the paternal parent and the other from the maternal parent. In this study, we introduce a novel heuristic algorithm for distinguishing single-nucleotide polymorphisms (SNPs) from the two parents and phasing them into haplotypes. The algorithm is unique because it simultaneously performs SNP calling and haplotype phasing. This approach can exploit the linkage information of nearby SNPs, which facilitates the efficient removal of haplotypes that originate from incorrectly mapped short reads. Using simulated data we demonstrated that our approach increased the accuracy of SNP calls. The haplotype reconstruction performance depended largely on the density of SNPs. Using current next-generation sequence technology with a relatively short read length, reasonable performance is expected when this approach is applied to species with an average of five heterozygous sites per 1 kb. The algorithm was implemented as the program "linkSNPs."

An NLR paralog Pit2 generated from tandem duplication of Pit1 fine-tunes Pit1 localization and function.

NLR family proteins act as intracellular receptors. Gene duplication amplifies the number of NLR genes, and subsequent mutations occasionally provide modifications to the second gene that benefits immunity. However, evolutionary processes after gene duplication and functional relationships between duplicated NLRs remain largely unclear. Here, we report that the rice NLR protein Pit1 is associated with its paralogue Pit2. The two are required for the resistance to rice blast fungus but have different functions: Pit1 induces cell death, while Pit2 competitively suppresses Pit1-mediated cell death. During evolution, the suppression of Pit1 by Pit2 was probably generated through positive selection on two fate-determining residues in the NB-ARC domain of Pit2, which account for functional differences between Pit1 and Pit2. Consequently, Pit2 lost its plasma membrane localization but acquired a new function to interfere with Pit1 in the cytosol. These findings illuminate the evolutionary trajectory of tandemly duplicated NLR genes after gene duplication.

The genetic basis of epigenetic variation and its consequences for adaptation.

Recent population genomic studies in plants have shed new light on natural epigenetic variation by identifying key genetic determinants, "trans modifiers," that influence epigenetic states genome-wide and their interplay with environmental factors. Here, we review this progress by focusing on the epigenetic control of transposition and life-cycle transitions to highlight the ecological consequences of this genetic architecture and its evolutionary significance. This knowledge provides new opportunities to address long-standing questions about the establishment of environment-associated epigenetic variation and its relevance in adaptation.

Role of camalexin, indole glucosinolates, and side chain modification of glucosinolate-derived isothiocyanates in defense of Arabidopsis against Sclerotinia sclerotiorum.

Plant secondary metabolites are known to facilitate interactions with a variety of beneficial and detrimental organisms, yet the contribution of specific metabolites to interactions with fungal pathogens is poorly understood. Here we show that, with respect to aliphatic glucosinolate-derived isothiocyanates, toxicity against the pathogenic ascomycete Sclerotinia sclerotiorum depends on side chain structure. Genes associated with the formation of the secondary metabolites camalexin and glucosinolate were induced in Arabidopsis thaliana leaves challenged with the necrotrophic pathogen S. sclerotiorum. Unlike S. sclerotiorum, the closely related ascomycete Botrytis cinerea was not identified to induce genes associated with aliphatic glucosinolate biosynthesis in pathogen-challenged leaves. Mutant plant lines deficient in camalexin, indole, or aliphatic glucosinolate biosynthesis were hypersusceptible to S. sclerotiorum, among them the myb28 mutant, which has a regulatory defect resulting in decreased production of long-chained aliphatic glucosinolates. The antimicrobial activity of aliphatic glucosinolate-derived isothiocyanates was dependent on side chain elongation and modification, with 8-methylsulfinyloctyl isothiocyanate being most toxic to S. sclerotiorum. This information is important for microbial associations with cruciferous host plants and for metabolic engineering of pathogen defenses in cruciferous plants that produce short-chained aliphatic glucosinolates.

AtCAST, a tool for exploring gene expression similarities among DNA microarray experiments using networks.

The comparison of gene expression profiles among DNA microarray experiments enables the identification of unknown relationships among experiments to uncover the underlying biological relationships. Despite the ongoing accumulation of data in public databases, detecting biological correlations among gene expression profiles from multiple laboratories on a large scale remains difficult. Here, we applied a module (sets of genes working in the same biological action)-based correlation analysis in combination with a network analysis to Arabidopsis data and developed a 'module-based correlation network' (MCN) which represents relationships among DNA microarray experiments on a large scale. We developed a Web-based data analysis tool, 'AtCAST' (Arabidopsis thaliana: DNA Microarray Correlation Analysis Tool), which enables browsing of an MCN or mining of users' microarray data by mapping the data into an MCN. AtCAST can help researchers to find novel connections among DNA microarray experiments, which in turn will help to build new hypotheses to uncover physiological mechanisms or gene functions in Arabidopsis.

Jasmonate-dependent and COI1-independent defense responses against Sclerotinia sclerotiorum in Arabidopsis thaliana: auxin is part of COI1-independent defense signaling.

The jasmonate receptor COI1 is known to facilitate plant defense responses against necrotrophic pathogens, including the ascomycete Sclerotinia sclerotiorum. However, it is not known to what extent jasmonates contribute to defense nor have COI1-independent defense pathways been sufficiently characterized. Here we show that the susceptibility to S. sclerotiorum of the aos mutant, deficient in biosynthesis of jasmonic acid (JA) and its precursor 12-oxophytadienoic acid, was elevated to a level reminiscent of that of hypersusceptible coi1 mutants. In contrast, susceptibility of the JA-deficient opr3 mutant was comparable with that of the wild type. A set of 99 genes responded similarly to infection with S. sclerotiorum in wild-type and coi1 mutant leaves. Expression of this COI1-independent gene set correlated with known differences in gene expression between wild-type plants and a mutant in the transcriptional repressor auxin response factor 2 (arf2). Susceptibility to S. sclerotiorum was reduced in two arf2 mutants early during infection, implicating ARF2 as a negative regulator of defense responses against this pathogen. Hypersusceptibility of an axr1 mutant to S. sclerotiorum confirmed the contribution of auxin action to defense responses against this fungal pathogen.

AtMetExpress development: a phytochemical atlas of Arabidopsis development.

Plants possess many metabolic genes for the production of a wide variety of phytochemicals in a tissue-specific manner. However, the metabolic systems behind the diversity and tissue-dependent regulation still remain unknown due to incomplete characterization of phytochemicals produced in a single plant species. Thus, having a metabolome dataset in addition to the genome and transcriptome information resources would enrich our knowledge of plant secondary metabolism. Here we analyzed phytochemical accumulation during development of the model plant Arabidopsis (Arabidopsis thaliana) using liquid chromatography-mass spectrometry in samples covering many growth stages and organs. We also obtained tandem mass spectrometry spectral tags of many metabolites as a resource for elucidation of metabolite structure. These are part of the AtMetExpress metabolite accumulation atlas. Based on the dataset, we detected 1,589 metabolite signals from which the structures of 167 metabolites were elucidated. The integrated analyses with transcriptome data demonstrated that Arabidopsis produces various phytochemicals in a highly tissue-specific manner, which often accompanies the expression of key biosynthesis-related genes. We also found that a set of biosynthesis-related genes is coordinately expressed among the tissues. These data suggested that the simple mode of regulation, transcript to metabolite, is an origin of the dynamics and diversity of plant secondary metabolism.

Improved social interaction and increased anterior cingulate metabolism after group reminiscence with reality orientation approach for vascular dementia.

A group reminiscence approach (GRA) with reality orientation (RO) is widely used as a psychosocial intervention for dementia. Since clinical effectiveness was reported for the intervention, interest has been directed toward areas of the neuronal network that might be being stimulated. We hypothesized that the frontal lobe associated with social interaction was being stimulated. To test this hypothesis, we studied 24 patients with vascular dementia. In addition to conventional care, a 1-h session of GRA with RO was provided once a week for 3 months in the GRA-RO arm (n=12). Only supportive care was provided in the control arm (n=12). Before and after the interventions, cognitive function, depressive state, and social activities were assessed. Since glucose metabolism is associated with brain function, cerebral glucose metabolism was measured by positron emission tomography (PET). Regarding behavioral improvement, 10 patients in the GRA-RO arm showed improvement compared with only two patients in the control arm, a significant difference. PET demonstrated that metabolism in the anterior cingulate was increased in the GRA-RO arm, whereas no significant changes were observed in the control arm. These results suggest that GRA-RO stimulates the anterior cingulate and has a positive effect on social interaction.

Auxin biosynthesis inhibitors, identified by a genomics-based approach, provide insights into auxin biosynthesis.

Despite its importance in plant growth and development, the auxin biosynthetic pathway has remained elusive. In this study, we analyzed hormone series transcriptome data from AtGenExpress in Arabidopsis and found that aminoethoxyvinylglycine (AVG) had the strongest anti-auxin activity. We also identified other effective compounds such as L-amino-oxyphenylpropionic acid (AOPP) through additional screening. These inhibitors shared characteristics in that they inhibited pyridoxal enzymes and/or aminotransferases. They reduced endogenous IAA levels in both monocots and dicots. L-AOPP inhibited root development of Arabidopsis in main root elongation, gravitropism, root skewing and root hair formation. This inhibition was generally recovered after exogenous IAA treatment, and the recovery was almost completely to the level of non-inhibited seedlings. The compounds inhibited conversion from tryptophan to indole-3-pyruvic acid in enzyme extracts from Arabidopsis and wheat. Our data collectively suggest that the inhibitors directly blocked auxin biosynthesis, and that the major target site was tryptophan aminotransferase. This enzyme probably makes up one of the major biosynthesis pathways conserved among higher plants. Each inhibitor, however, demonstrated a different action spectrum in shoot and root of rice and tomato, indicating diversity in biosynthesis pathways between organs and species. Our results provide novel insights into auxin biosynthesis and action, and uncover structural characteristics of auxin biosynthesis inhibitors.