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.

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.

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.

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.

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.

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.

GWAS with Heterogeneous Data: Estimating the Fraction of Phenotypic Variation Mediated by Gene Expression Data.

Intermediate phenotypes such as gene expression values can be used to elucidate the mechanisms by which genetic variation causes phenotypic variation, but jointly analyzing such heterogeneous data are far from trivial. Here we extend a so-called mediation model to handle the confounding effects of genetic background, and use it to analyze flowering time variation in Arabidopsis thaliana, focusing in particular on the central role played by the key regulator FLOWERING TIME LOCUS C (FLC). FLC polymorphism and FLC expression are both strongly correlated with flowering time variation, but the effect of the former is only partly mediated through the latter. Furthermore, the latter also reflects genetic background effects. We demonstrate that it is possible to partition these effects, shedding light on the complex regulatory network that underlies flowering time variation.

Multiple alleles at a single locus control seed dormancy in Swedish Arabidopsis.

Seed dormancy is a complex life history trait that determines the timing of germination and is crucial for local adaptation. Genetic studies of dormancy are challenging, because the trait is highly plastic and strongly influenced by the maternal environment. Using a combination of statistical and experimental approaches, we show that multiple alleles at the previously identified dormancy locus DELAY OF GERMINATION1 jointly explain as much as 57% of the variation observed in Swedish Arabidopsis thaliana, but give rise to spurious associations that seriously mislead genome-wide association studies unless modeled correctly. Field experiments confirm that the major alleles affect germination as well as survival under natural conditions, and demonstrate that locally adaptive traits can sometimes be dissected genetically.

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.

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.

Computational analyses of ancient pathogen DNA from herbarium samples: challenges and prospects.

The application of DNA sequencing technology to the study of ancient DNA has enabled the reconstruction of past epidemics from genomes of historically important plant-associated microbes. Recently, the genome sequences of the potato late blight pathogen Phytophthora infestans were analyzed from 19th century herbarium specimens. These herbarium samples originated from infected potatoes collected during and after the Irish potato famine. Herbaria have therefore great potential to help elucidate past epidemics of crops, date the emergence of pathogens, and inform about past pathogen population dynamics. DNA preservation in herbarium samples was unexpectedly good, raising the possibility of a whole new research area in plant and microbial genomics. However, the recovered DNA can be extremely fragmented resulting in specific challenges in reconstructing genome sequences. Here we review some of the challenges in computational analyses of ancient DNA from herbarium samples. We also applied the recently developed linkage method to haplotype reconstruction of diploid or polyploid genomes from fragmented ancient DNA.

"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.

Strigolactone regulates anthocyanin accumulation, acid phosphatases production and plant growth under low phosphate condition in Arabidopsis.

Phosphate is an essential macronutrient in plant growth and development; however, the concentration of inorganic phosphate (Pi) in soil is often suboptimal for crop performance. Accordingly, plants have developed physiological strategies to adapt to low Pi availability. Here, we report that typical Pi starvation responses in Arabidopsis are partially dependent on the strigolactone (SL) signaling pathway. SL treatment induced root hair elongation, anthocyanin accumulation, activation of acid phosphatase, and reduced plant weight, which are characteristic responses to phosphate starvation. Furthermore, the expression profile of SL-response genes correlated with the expression of genes induced by Pi starvation. These results suggest a potential overlap between SL signaling and Pi starvation signaling pathways in plants.

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.

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.

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.

Intracellular diversity of the V4 and V9 regions of the 18S rRNA in marine protists (radiolarians) assessed by high-throughput sequencing.

Metabarcoding is a powerful tool for exploring microbial diversity in the environment, but its accurate interpretation is impeded by diverse technical (e.g. PCR and sequencing errors) and biological biases (e.g. intra-individual polymorphism) that remain poorly understood. To help interpret environmental metabarcoding datasets, we investigated the intracellular diversity of the V4 and V9 regions of the 18S rRNA gene from Acantharia and Nassellaria (radiolarians) using 454 pyrosequencing. Individual cells of radiolarians were isolated, and PCRs were performed with generalist primers to amplify the V4 and V9 regions. Different denoising procedures were employed to filter the pyrosequenced raw amplicons (Acacia, AmpliconNoise, Linkage method). For each of the six isolated cells, an average of 541 V4 and 562 V9 amplicons assigned to radiolarians were obtained, from which one numerically dominant sequence and several minor variants were found. At the 97% identity, a diversity metrics commonly used in environmental surveys, up to 5 distinct OTUs were detected in a single cell. However, most amplicons grouped within a single OTU whereas other OTUs contained very few amplicons. Different analytical methods provided evidence that most minor variants forming different OTUs correspond to PCR and sequencing artifacts. Duplicate PCR and sequencing from the same DNA extract of a single cell had only 9 to 16% of unique amplicons in common, and alignment visualization of V4 and V9 amplicons showed that most minor variants contained substitutions in highly-conserved regions. We conclude that intracellular variability of the 18S rRNA in radiolarians is very limited despite its multi-copy nature and the existence of multiple nuclei in these protists. Our study recommends some technical guidelines to conservatively discard artificial amplicons from metabarcoding datasets, and thus properly assess the diversity and richness of protists in the environment.

QTL map meets population genomics: an application to rice.

Genes involved in the transition from wild to cultivated crop species should be of great agronomic importance. Population genomic approaches utilizing genome resequencing data have been recently applied for this purpose, although it only reports a large list of candidate genes with no biological information. Here, by resequencing more than 30 genomes altogether of wild rice Oryza rufipogon and cultivated rice O. sativa, we identified a number of regions with clear footprints of selection during the domestication process. We then focused on identifying candidate domestication genes in these regions by utilizing the wealth of QTL information in rice. We were able to identify a number of interesting candidates such as transcription factors that should control key domestication traits such as shattering, awn length, and seed dormancy. Other candidates include those that might have been related to the improvement of grain quality and those that might have been involved in the local adaptation to dry conditions and colder environments. Our study shows that population genomic approaches and QTL mapping information can be used together to identify genes that might be of agronomic importance.

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."