Positional information resolves structural variations and uncovers an evolutionarily divergent genetic locus in accessions of Arabidopsis thaliana

Genomic Rearrangements in Arabidopsis Considered as Quantitative Traits

Structural Rearrangements can have unexpected effects on quantitative phenotypes. Surprisingly, these rearrangements can also be considered as...

To understand the population genetics of structural variants and their effects on phenotypes, we developed an approach to mapping structural variants that segregate in a population sequenced at low coverage. We avoid calling structural variants directly. Instead, the evidence for a potential structural variant at a locus is indicated by variation in the counts of short-reads that map anomalously to that locus. These structural variant traits are treated as quantitative traits and mapped genetically, analogously to a gene expression study. Association between a structural variant trait at one locus, and genotypes at a distant locus indicate the origin and target of a transposition. Using ultra-low-coverage (0.3×) population sequence data from 488 recombinant inbred Arabidopsis thaliana genomes, we identified 6502 segregating structural variants. Remarkably, 25% of these were transpositions. While many structural variants cannot be delineated precisely, we validated 83% of 44 predicted transposition breakpoints by polymerase chain reaction. We show that specific structural variants may be causative for quantitative trait loci for germination and resistance to infection by the fungus Albugo laibachii, isolate Nc14. Further we show that the phenotypic heritability attributable to read-mapping anomalies differs from, and, in the case of time to germination and bolting, exceeds that due to standard genetic variation. Genes within structural variants are also more likely to be silenced or dysregulated. This approach complements the prevalent strategy of structural variant discovery in fewer individuals sequenced at high coverage. It is generally applicable to large populations sequenced at low-coverage, and is particularly suited to mapping transpositions.

Chromosome fragile sites in Arabidopsis harbor matrix attachment regions that may be associated with ancestral chromosome rearrangement events

Mutations in the BREVIPEDICELLUS (BP) gene of Arabidopsis thaliana condition a pleiotropic phenotype featuring defects in internode elongation, the homeotic conversion of internode to node tissue, and downward pointing flowers and pedicels. We have characterized five mutant alleles of BP, generated by EMS, fast neutrons, x-rays, and aberrant T–DNA insertion events. Curiously, all of these mutagens resulted in large deletions that range from 140 kbp to over 900 kbp just south of the centromere of chromosome 4. The breakpoints of these mutants were identified by employing inverse PCR and DNA sequencing. The south breakpoints of all alleles cluster in BAC T12G13, while the north breakpoint locations are scattered. With the exception of a microhomology at the bp-5 breakpoint, there is no homology in the junction regions, suggesting that double-stranded breaks are repaired via non-homologous end joining. Southwestern blotting demonstrated the presence of nuclear matrix binding sites in the south breakpoint cluster (SBC), which is A/T rich and possesses a variety of repeat sequences. In situ hybridization on pachytene chromosome spreads complemented the molecular analyses and revealed heretofore unrecognized structural variation between the Columbia and Landsberg erecta genomes. Data mining was employed to localize other large deletions around the HY4 locus to the SBC region and to show that chromatin modifications in the region shift from a heterochromatic to euchromatic profile. Comparisons between the BP/HY4 regions of A. lyrata and A. thaliana revealed that several chromosome rearrangement events have occurred during the evolution of these two genomes. Collectively, the features of the region are strikingly similar to the features of characterized metazoan chromosome fragile sites, some of which are associated with karyotype evolution.

Chromosome evolution involves both small-scale (e.g. single nucleotide) changes, as well as large-scale rearrangements such as inversions, translocations, and fusion events. We investigated mutations of the BREVIPEDICELLUS gene of Arabidopsis, which is a master regulator of inflorescence architecture. These mutations are not due to single nucleotide changes, but rather to large deletions, some spanning nearly one million base pairs. Molecular and biochemical analyses reveal that the chromosome breakpoints cluster in an area that is rich in repetitive elements and harbor multiple binding sites for nuclear matrix proteins. Data mining revealed intriguing correlations between the breakpoint cluster and hotspots of genetic recombination, regions of the chromosome that have undergone several rearrangement events during evolution, and changes in histone protein modifications. We propose that these unstable regions are chromosome fragile sites that assist in marking a boundary between gene-poor, transcriptionally repressed centromeric chromatin and a more relaxed chromatin domain that is gene-rich.

Natural variation in GL1 and its effects on trichome density in Arabidopsis thaliana

The ultimate understanding of how biological diversity arises, is maintained, and lost depends on identifying the genes responsible. Although a good deal has been discovered about gene function over the past few decades, far less is understood about gene effects, that is, how natural variation in a gene contributes to natural variation in phenotypes. Trichome density in Arabidopsis thaliana is an ideal trait for studies of natural molecular and phenotypic variation, as trichome initiation is genetically well-characterized and trichome density is highly variable in and among natural populations. Here, we show that variation at GLABRA1 (GL1), an R2R3 MYB transcription factor gene, which has a role in trichome initiation, has qualitative and likely quantitative effects on trichome density in natural accessions of A. thaliana. Specifically, we characterize four independent loss-of-function alleles for GL1, each of which yields a glabrous phenotype. Further, we find that a pattern of common polymorphisms confined to the GL1 locus is associated with quantitative variation for trichome density. While mutations resulting in a glabrous phenotype are primarily coding changes, the pattern resulting in quantitative variation spans both coding and regulatory regions. These data show that GL1 is an important source of trichome density variation within A. thaliana and, along with recent reports, suggest that the TTG1 epidermal cell fate pathway generally may be the key molecular genetic source of natural trichome density variation and an important model for the study of molecular evolution.

Multiple reference genomes and transcriptomes for Arabidopsis thaliana

Genetic differences between Arabidopsis thaliana accessions underlie the plant's extensive phenotypic variation, and until now these have been interpreted largely in the context of the annotated reference accession Col-0. Here we report the sequencing, assembly and annotation of the genomes of 18 natural A. thaliana accessions, and their transcriptomes. When assessed on the basis of the reference annotation, one-third of protein-coding genes are predicted to be disrupted in at least one accession. However, re-annotation of each genome revealed that alternative gene models often restore coding potential. Gene expression in seedlings differed for nearly half of expressed genes and was frequently associated with cis variants within 5 kilobases, as were intron retention alternative splicing events. Sequence and expression variation is most pronounced in genes that respond to the biotic environment. Our data further promote evolutionary and functional studies in A. thaliana, especially the MAGIC genetic reference population descended from these accessions.