Beth A. Rowan

1,135 Genomes Reveal the Global Pattern of Polymorphism in Arabidopsis thaliana

Summary Arabidopsis thaliana serves as a model organism for the study of fundamental physiological, cellular, and molecular processes. It has also greatly advanced our understanding of intraspecific genome variation. We present a detailed map of variation in 1,135 high-quality re-sequenced natural inbred lines representing the native Eurasian and North African range and recently colonized North America. We identify relict populations that continue to inhabit ancestral habitats, primarily in the Iberian Peninsula. They have mixed with a lineage that has spread to northern latitudes from an unknown glacial refugium and is now found in a much broader spectrum of habitats. Insights into the history of the species and the fine-scale distribution of genetic diversity provide the basis for full exploitation of A. thaliana natural variation through integration of genomes and epigenomes with molecular and non-molecular phenotypes.

Species-wide Genetic Incompatibility Analysis Identifies Immune Genes as Hot Spots of Deleterious Epistasis

Intraspecific genetic incompatibilities prevent the assembly of specific alleles into single genotypes and influence genome- and species-wide patterns of sequence variation. A common incompatibility in plants is hybrid necrosis, characterized by autoimmune responses due to epistatic interactions between natural genetic variants. By systematically testing thousands of F1 hybrids of Arabidopsis thaliana strains, we identified a small number of incompatibility hot spots in the genome, often in regions densely populated by nucleotide-binding domain and leucine-rich repeat (NLR) immune receptor genes. In several cases, these immune receptor loci interact with each other, suggestive of conflict within the immune system. A particularly dangerous locus is a highly variable cluster of NLR genes, DM2, which causes multiple independent incompatibilities with genes that encode a range of biochemical functions, including NLRs. Our findings suggest that deleterious interactions of immune receptors limit the combinations of favorable disease resistance alleles accessible to plant genomes.

The genomic landscape of meiotic crossovers and gene conversions in Arabidopsis thaliana

Knowledge of the exact distribution of meiotic crossovers (COs) and gene conversions (GCs) is essential for understanding many aspects of population genetics and evolution, from haplotype structure and long-distance genetic linkage to the generation of new allelic variants of genes. To this end, we resequenced the four products of 13 meiotic tetrads along with 10 doubled haploids derived from Arabidopsis thaliana hybrids. GC detection through short reads has previously been confounded by genomic rearrangements. Rigid filtering for misaligned reads allowed GC identification at high accuracy and revealed an ∼80-kb transposition, which undergoes copy-number changes mediated by meiotic recombination. Non-crossover associated GCs were extremely rare most likely due to their short average length of ∼25–50 bp, which is significantly shorter than the length of CO-associated GCs. Overall, recombination preferentially targeted non-methylated nucleosome-free regions at gene promoters, which showed significant enrichment of two sequence motifs. DOI: http://dx.doi.org/10.7554/eLife.01426.001

Rapid and Inexpensive Whole-Genome Genotyping-by-Sequencing for Crossover Localization and Fine-Scale Genetic Mapping

The reshuffling of existing genetic variation during meiosis is important both during evolution and in breeding. The reassortment of genetic variants relies on the formation of crossovers (COs) between homologous chromosomes. The pattern of genome-wide CO distributions can be rapidly and precisely established by the short-read sequencing of individuals from F2 populations, which in turn are useful for quantitative trait locus (QTL) mapping. Although sequencing costs have decreased precipitously in recent years, the costs of library preparation for hundreds of individuals have remained high. To enable rapid and inexpensive CO detection and QTL mapping using low-coverage whole-genome sequencing of large mapping populations, we have developed a new method for library preparation along with Trained Individual GenomE Reconstruction, a probabilistic method for genotype and CO predictions for recombinant individuals. In an example case with hundreds of F2 individuals from two Arabidopsis thaliana accessions, we resolved most CO breakpoints to within 2 kb and reduced a major flowering time QTL to a 9-kb interval. In addition, an extended region of unusually low recombination revealed a 1.8-Mb inversion polymorphism on the long arm of chromosome 4. We observed no significant differences in the frequency and distribution of COs between F2 individuals with and without a functional copy of the DNA helicase gene RECQ4A. In summary, we present a new, cost-efficient method for large-scale, high-precision genotyping-by-sequencing.

Developmental Genetics and New Sequencing Technologies: The Rise of Nonmodel Organisms

Much of developmental biology in the past decades has been driven by forward genetic studies in a few model organisms. We review recent work with relatives of these species, motivated by a desire to understand the evolutionary and ecological context for morphological innovation. Unfortunately, despite a number of shining examples, progress in nonmodel systems has often been slow. The current revolution in DNA sequencing has, however, enormous potential in extending the reach of genetics. We discuss how developmental biology will benefit from these advances, particularly by increasing the universe of study species.

Genome-wide detection of intervals of genetic heterogeneity associated with complex traits

Motivation: Genetic heterogeneity, the fact that several sequence variants give rise to the same phenotype, is a phenomenon that is of the utmost interest in the analysis of complex phenotypes. Current approaches for finding regions in the genome that exhibit genetic heterogeneity suffer from at least one of two shortcomings: (i) they require the definition of an exact interval in the genome that is to be tested for genetic heterogeneity, potentially missing intervals of high relevance, or (ii) they suffer from an enormous multiple hypothesis testing problem due to the large number of potential candidate intervals being tested, which results in either many false positives or a lack of power to detect true intervals. Results: Here, we present an approach that overcomes both problems: it allows one to automatically find all contiguous sequences of single nucleotide polymorphisms in the genome that are jointly associated with the phenotype. It also solves both the inherent computational efficiency problem and the statistical problem of multiple hypothesis testing, which are both caused by the huge number of candidate intervals. We demonstrate on Arabidopsis thaliana genome-wide association study data that our approach can discover regions that exhibit genetic heterogeneity and would be missed by single-locus mapping. Conclusions: Our novel approach can contribute to the genome-wide discovery of intervals that are involved in the genetic heterogeneity underlying complex phenotypes. Availability and implementation: The code can be obtained at: http://www.bsse.ethz.ch/mlcb/research/bioinformatics-and-computational-biology/sis.html. Contact: felipe.llinares@bsse.ethz.ch Supplementary information: Supplementary data are available at Bioinformatics online.

Methods for Genotyping-by-Sequencing

A major goal for biologists is to understand the connection between genes and phenotypic traits, and genetic mapping in experimental populations remains a powerful approach for discovering the causal genes underlying phenotypes. For genetic mapping, the process of genotyping was previously a major rate-limiting step. Modern sequencing technology has greatly improved the resolution and speed of genetic mapping by reducing the time, labor, and cost per genotyping marker. In addition, the ability to perform genotyping-by-sequencing (GBS) has facilitated large-scale population genetic analyses by providing a simpler way to survey segregating genetic variation in natural populations. Here we present two protocols for GBS, using the Illumina platform, that can be applied to a wide range of genotyping projects in different species. The first protocol is for genotyping a subset of marker positions genome-wide using restriction digestion, and the second is for preparing inexpensive paired-end whole-genome libraries. We discuss the suitability of each approach for different genotyping applications and provide notes for adapting these protocols for use with a liquid-handling robot.

The Plant Immunity Regulating F-Box Protein CPR1 Supports Plastid Function in Absence of Pathogens

The redox imbalanced 6 mutant (rimb6) of Arabidopsis thaliana was isolated in a genetic screening approach for mutants with defects in chloroplast-to-nucleus redox signaling. It has an atypically low activation status of the 2-Cys peroxiredoxin-A promoter in the seedling stage. rimb6 shows wildtype-like germination, seedling development and greening, but slower growth and reduced biomass in the rosette stage. Mapping of the casual mutation revealed that rimb6 carries a single nucleotide polymorphism in the gene encoding CONSTITUTIVE EXPRESSER OF PATHOGENESIS RELATED (PR) GENES 1, CPR1 (At4g12560), leading to a premature stop codon. CPR1 is known as a repressor of pathogen signaling and regulator of microtubule organization. Allelism of rimb6 and cpr1 revealed a function of CPR1 in chloroplast stress protection. Expression studies in pathogen signaling mutants demonstrated that CPR1-mediated activation of genes for photosynthesis and chloroplast antioxidant protection is, in contrast to activation of pathogen responses, regulated independently from PAD4-controlled salicylic acid (SA) accumulation. We conclude that the support of plastid function is a basic, SA-independent function of CPR1.

A rainfall-manipulation experiment with 517 Arabidopsis thaliana accessions

The gold standard for studying natural selection and adaptation in the wild is to quantify lifetime fitness of individuals from natural populations that have been grown together in a common garden, or that have been reciprocally transplanted. By combining fitness values with species traits and genome sequences, one can infer selection coefficients at the genetic level. Here we present a rainfall-manipulation experiment with 517 whole-genome sequenced natural accessions of the plant Arabidopsis thaliana spanning the global distribution of the species. The experiments were conducted in two field stations in contrasting climates, in the Mediterranean and in Central Europe, where we built rainout shelters and simulated high and low rainfall. Using custom image analysis we quantified fitness- and phenology-related traits for 23,154 pots, which contained about 14,500 plants growing independently, and over 310,000 plants growing in small populations (max. 30 plants). This large field experiment dataset, which associates fitness and ecologically-relevant traits with genomes, will provide an important resource to test eco-evolutionary genetic theories and to understand the potential evolutionary impacts of future climates on an important plant model species.

A single haplotype hyposensitive to light and requiring strong vernalization dominates Arabidopsis thaliana populations in Patagonia, Argentina

The growing collection of sequenced or genotyped Arabidopsis thaliana accessions includes mostly individuals from the native Eurasian and N. African range and introduced North American populations. Here, we describe the genetic and phenotypic diversity, along with habitats and life history, of A. thaliana plants collected at the southernmost end of its worldwide distribution. Seed samples were harvested from plants growing in four sites within a ~3500‐km2‐area in Patagonia, Argentina, and represent the first germplasm to be collected in South America for this species. Whole‐genome resequencing revealed that plants from the four sites and a Patagonia herbarium specimen collected in 1967 formed a single haplogroup (Pat), indicating that the phenotypic variation observed in the field reflected plastic responses to the environment. admixture and principal components analyses suggest that the ancestor of the Pat haplogroup either came from Italy or the Balkan/Caucasus regions of Eurasia. In the laboratory, plants from the Pat haplogroup were hyposensitive to continuous red (Rc) and shade light, with corresponding changes in the expression of phytochrome signalling genes. Pat had higher PIF3 and PIF5 and lower HY5 expression under Rc light; and lower expression of PIL1, ATHB2 and HFR1 under shade compared to Col‐0. In addition, Pat plants had a strong vernalization requirement associated with high levels of FLC expression. We conclude that including Pat in studies of natural variation and in comparison with other introduced populations will provide additional information for association studies and allow for a more detailed assessment of the demographic events following colonization.