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.

Arabidopsis FLL2 promotes liquid-liquid phase separation of polyadenylation complexes.

An important component of cellular biochemistry is the concentration of proteins and nucleic acids in non-membranous compartments1,2. These biomolecular condensates are formed from processes that include liquid-liquid phase separation. The multivalent interactions necessary for liquid-liquid phase separation have been extensively studied in vitro1,3. However, the regulation of this process in vivo is poorly understood. Here we identify an in vivo regulator of liquid-liquid phase separation through a genetic screen targeting factors required for Arabidopsis RNA-binding protein FCA function. FCA contains prion-like domains that phase-separate in vitro, and exhibits behaviour in vivo that is consistent with phase separation. The mutant screen identified a functional requirement for FLL2, a coiled-coil protein, in the formation of FCA nuclear bodies. FCA reduces transcriptional read-through by promoting proximal polyadenylation at many sites in the Arabidopsis genome3,4. FLL2 was required to promote this proximal polyadenylation, but not the binding of FCA to target RNA. Ectopic expression of FLL2 increased the size and number of FCA nuclear bodies. Crosslinking with formaldehyde captured in vivo interactions between FLL2, FCA and the polymerase and nuclease modules of the RNA 3'-end processing machinery. These 3' RNA-processing components colocalized with FCA in the nuclear bodies in vivo, which indicates that FCA nuclear bodies compartmentalize 3'-end processing factors to enhance polyadenylation at specific sites. Our findings show that coiled-coil proteins can promote liquid-liquid phase separation, which expands our understanding of the principles that govern the in vivo dynamics of liquid-like bodies.

MSH2 shapes the meiotic crossover landscape in relation to interhomolog polymorphism in Arabidopsis.

During meiosis, DNA double-strand breaks undergo interhomolog repair to yield crossovers between homologous chromosomes. To investigate how interhomolog sequence polymorphism affects crossovers, we sequenced multiple recombinant populations of the model plant Arabidopsis thaliana. Crossovers were elevated in the diverse pericentromeric regions, showing a local preference for polymorphic regions. We provide evidence that crossover association with elevated diversity is mediated via the Class I crossover formation pathway, although very high levels of diversity suppress crossovers. Interhomolog polymorphism causes mismatches in recombining molecules, which can be detected by MutS homolog (MSH) mismatch repair protein heterodimers. Therefore, we mapped crossovers in a msh2 mutant, defective in mismatch recognition, using multiple hybrid backgrounds. Although total crossover numbers were unchanged in msh2 mutants, recombination was remodelled from the diverse pericentromeres towards the less-polymorphic sub-telomeric regions. Juxtaposition of megabase heterozygous and homozygous regions causes crossover remodelling towards the heterozygous regions in wild type Arabidopsis, but not in msh2 mutants. Immunostaining showed that MSH2 protein accumulates on meiotic chromosomes during prophase I, consistent with MSH2 regulating meiotic recombination. Our results reveal a pro-crossover role for MSH2 in regions of higher sequence diversity in A. thaliana.

The alternative reality of plant mitochondrial DNA: One ring does not rule them all.

Plant mitochondrial genomes are usually assembled and displayed as circular maps based on the widely-held view across the broad community of life scientists that circular genome-sized molecules are the primary form of plant mitochondrial DNA, despite the understanding by plant mitochondrial researchers that this is an inaccurate and outdated concept. Many plant mitochondrial genomes have one or more pairs of large repeats that can act as sites for inter- or intramolecular recombination, leading to multiple alternative arrangements (isoforms). Most mitochondrial genomes have been assembled using methods unable to capture the complete spectrum of isoforms within a species, leading to an incomplete inference of their structure and recombinational activity. To document and investigate underlying reasons for structural diversity in plant mitochondrial DNA, we used long-read (PacBio) and short-read (Illumina) sequencing data to assemble and compare mitochondrial genomes of domesticated (Lactuca sativa) and wild (L. saligna and L. serriola) lettuce species. We characterized a comprehensive, complex set of isoforms within each species and compared genome structures between species. Physical analysis of L. sativa mtDNA molecules by fluorescence microscopy revealed a variety of linear, branched, and circular structures. The mitochondrial genomes for L. sativa and L. serriola were identical in sequence and arrangement and differed substantially from L. saligna, indicating that the mitochondrial genome structure did not change during domestication. From the isoforms in our data, we infer that recombination occurs at repeats of all sizes at variable frequencies. The differences in genome structure between L. saligna and the two other Lactuca species can be largely explained by rare recombination events that rearranged the structure. Our data demonstrate that representations of plant mitochondrial genomes as simple, circular molecules are not accurate descriptions of their true nature and that in reality plant mitochondrial DNA is a complex, dynamic mixture of forms.

Dose-dependent interactions between two loci trigger altered shoot growth in BG-5 × Krotzenburg-0 (Kro-0) hybrids of Arabidopsis thaliana.

Hybrids occasionally exhibit genetic interactions resulting in reduced fitness in comparison to their parents. Studies of Arabidopsis thaliana have highlighted the role of immune conflicts, but less is known about the role of other factors in hybrid incompatibility in plants. Here, we present a new hybrid incompatibility phenomenon in this species. We have characterized a new case of F1 hybrid incompatibility from a cross between the A. thaliana accessions Krotzenburg-0 (Kro-0) and BG-5, by conducting transcript, metabolite and hormone analyses, and identified the causal loci through genetic mapping. The F1 hybrids showed arrested growth of the main stem, altered shoot architecture, and altered concentrations of hormones in comparison to parents. The F1 phenotype could be rescued in a developmental-stage-dependent manner by shifting to a higher growth temperature. These F1 phenotypes were linked to two loci, one on chromosome 2 and one on chromosome 3. The F2 generation segregated plants with more severe phenotypes which were linked to the same loci as those in the F1 . This study provides novel insights into how previously unknown mechanisms controlling shoot branching and stem growth can result in hybrid incompatibility.

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.

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.

RecA maintains the integrity of chloroplast DNA molecules in Arabidopsis.

Although our understanding of mechanisms of DNA repair in bacteria and eukaryotic nuclei continues to improve, almost nothing is known about the DNA repair process in plant organelles, especially chloroplasts. Since the RecA protein functions in DNA repair for bacteria, an analogous function may exist for chloroplasts. The effects on chloroplast DNA (cpDNA) structure of two nuclear-encoded, chloroplast-targeted homologues of RecA in Arabidopsis were examined. A homozygous T-DNA insertion mutation in one of these genes (cpRecA) resulted in altered structural forms of cpDNA molecules and a reduced amount of cpDNA, while a similar mutation in the other gene (DRT100) had no effect. Double mutants exhibited a similar phenotype to cprecA single mutants. The cprecA mutants also exhibited an increased amount of single-stranded cpDNA, consistent with impaired RecA function. After four generations, the cprecA mutant plants showed signs of reduced chloroplast function: variegation and necrosis. Double-stranded breaks in cpDNA of wild-type plants caused by ciprofloxacin (an inhibitor of Escherichia coli gyrase, a type II topoisomerase) led to an alteration of cpDNA structure that was similar to that seen in cprecA mutants. It is concluded that the process by which damaged DNA is repaired in bacteria has been retained in their endosymbiotic descendent, the chloroplast.

A multiple-method approach reveals a declining amount of chloroplast DNA during development in Arabidopsis.

BACKGROUND: A decline in chloroplast DNA (cpDNA) during leaf maturity has been reported previously for eight plant species, including Arabidopsis thaliana. Recent studies, however, concluded that the amount of cpDNA during leaf development in Arabidopsis remained constant. RESULTS: To evaluate alternative hypotheses for these two contradictory observations, we examined cpDNA in Arabidopsis shoot tissues at different times during development using several methods: staining leaf sections as well as individual isolated chloroplasts with 4',6-diamidino-2-phenylindole (DAPI), real-time quantitative PCR with DNA prepared from total tissue as well as from isolated chloroplasts, fluorescence microscopy of ethidium-stained DNA molecules prepared in gel from isolated plastids, and blot-hybridization of restriction-digested total tissue DNA. We observed a developmental decline of about two- to three-fold in mean DNA per chloroplast and two- to five-fold in the fraction of cellular DNA represented by chloroplast DNA. CONCLUSION: Since the two- to five-fold reduction in cpDNA content could not be attributed to an artifact of chloroplast isolation, we conclude that DNA within Arabidopsis chloroplasts is degraded in vivo as leaves mature.

The loss of DNA from chloroplasts as leaves mature: fact or artefact?

In this review, the controversy regarding the preservation or degradation of chloroplast DNA (cpDNA) as chloroplasts develop their photosynthetic capacity and leaves reach maturity is addressed. A constant amount of cpDNA during maturity might be expected in order to support photosynthesis over the lifespan of the leaf. Nevertheless, a decline in cpDNA during leaf development was found for all seven plant species investigated. Initial measurements showed that Arabidopsis was similar to the other seven. The controversy arose with two recent studies concluding that the amount of cpDNA remains constant as Arabidopsis leaves mature. These authors proposed that the observation of Arabidopsis chloroplasts with undetectable levels of DNA was an artefact, although the most recent data support the original findings. If the amount of cpDNA remains constant, then Arabidopsis is atypical and would not serve as a good model for chloroplast development. It is shown that the apparently contradictory data may be attributed to methodology and the choice of leaves to be compared. Thus, it is concluded that the controversy can be resolved, Arabidopsis can serve as a representative model, and cpDNA degradation is a common event in chloroplast development.

An Ultra High-Density Arabidopsis thaliana Crossover Map That Refines the Influences of Structural Variation and Epigenetic Features.

Many environmental, genetic, and epigenetic factors are known to affect the frequency and positioning of meiotic crossovers (COs). Suppression of COs by large, cytologically visible inversions and translocations has long been recognized, but relatively little is known about how smaller structural variants (SVs) affect COs. To examine fine-scale determinants of the CO landscape, including SVs, we used a rapid, cost-effective method for high-throughput sequencing to generate a precise map of >17,000 COs between the Col-0 and Ler-0 accessions of Arabidopsis thaliana COs were generally suppressed in regions with SVs, but this effect did not depend on the size of the variant region, and was only marginally affected by the variant type. CO suppression did not extend far beyond the SV borders and CO rates were slightly elevated in the flanking regions. Disease resistance gene clusters, which often exist as SVs, exhibited high CO rates at some loci, but there was a tendency toward depressed CO rates at loci where large structural differences exist between the two parents. Our high-density map also revealed in fine detail how CO positioning relates to genetic (DNA motifs) and epigenetic (chromatin structure) features of the genome. We conclude that suppression of COs occurs over a narrow region spanning large- and small-scale SVs, representing an influence on the CO landscape in addition to sequence and epigenetic variation along chromosomes.

Linked-read sequencing of gametes allows efficient genome-wide analysis of meiotic recombination.

Meiotic crossovers (COs) ensure proper chromosome segregation and redistribute the genetic variation that is transmitted to the next generation. Large populations and the demand for genome-wide, fine-scale resolution challenge existing methods for CO identification. Taking advantage of linked-read sequencing, we develop a highly efficient method for genome-wide identification of COs at kilobase resolution in pooled recombinants. We first test this method using a pool of Arabidopsis F2 recombinants, and recapitulate results obtained from the same plants using individual whole-genome sequencing. By applying this method to a pool of pollen DNA from an F1 plant, we establish a highly accurate CO landscape without generating or sequencing a single recombinant plant. The simplicity of this approach enables the simultaneous generation and analysis of multiple CO landscapes, accelerating the pace at which mechanisms for the regulation of recombination can be elucidated through efficient comparisons of genotypic and environmental effects on recombination.

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.

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.

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.

Isolation, quantification, and analysis of chloroplast DNA.

Many areas of chloroplast research require methods that can assess the quality and quantity of chloroplast DNA (cpDNA). The study of chloroplast functions that depend on the proper maintenance and expression of the chloroplast genome, understanding cpDNA replication and repair, and the development of technologies for chloroplast transformation are just some of the disciplines that require the isolation of high-quality cpDNA. Arabidopsis thaliana offers several advantages for studying these processes because of the sizeable collection of mutants and natural varieties (accessions) available from stock centers and a broad community of researchers that has developed many other genetic resources. Several approaches for the isolation and quantification of cpDNA have been developed, but little consideration has been given to the strengths and weaknesses and the type of information obtained by each method, especially with respect to A. thaliana. Here, we provide protocols for obtaining high-quality cpDNA for PCR and other applications, and we evaluate several different isolation and analytical methods in order to build a robust framework for the study of cpDNA with this model organism.

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.