Megabase-scale inversion polymorphism in the wild ancestor of maize

Population genomics highlights structural variations in local adaptation to saline coastal environments in woolly grape

Structural variations (SVs) are a feature of plant genomes that has been largely unexplored despite their significant impact on plant phenotypic traits and local adaptation to abiotic and biotic stress. In this study, we employed woolly grape (Vitis retordii), a species native to the tropical and subtropical regions of East Asia with both coastal and inland habitats, as a valuable model for examining the impact of SVs on local adaptation. We assembled a haplotype-resolved chromosomal reference genome for woolly grape, and conducted population genetic analyses based on whole-genome sequencing (WGS) data from coastal and inland populations. The demographic analyses revealed recent bottlenecks in all populations and asymmetric gene flow from the inland to the coastal population. In total, 1,035 genes associated with plant adaptive regulation for salt stress, radiation, and environmental adaptation were detected underlying local selection by SVs and SNPs in the coastal population, of which 37.29% and 65.26% were detected by SVs and SNPs, respectively. Candidate genes such as FSD2, RGA1, and AAP8 associated with salt tolerance were found to be highly differentiated and selected during the process of local adaptation to coastal habitats in SV regions. Our study highlights the importance of SVs in local adaptation; candidate genes related to salt stress and climatic adaptation to tropical and subtropical environments are important genomic resources for future breeding programs of grapevine and its rootstocks.

Genomic determinants, architecture, and constraints in drought-related traits in Corymbia calophylla

BACKGROUND

Drought adaptation is critical to many tree species persisting under climate change, however our knowledge of the genetic basis for trees to adapt to drought is limited. This knowledge gap impedes our fundamental understanding of drought response and application to forest production and conservation. To improve our understanding of the genomic determinants, architecture, and trait constraints, we assembled a reference genome and detected ~ 6.5 M variants in 432 phenotyped individuals for the foundational tree Corymbia calophylla.

RESULTS

We found 273 genomic variants determining traits with moderate heritability (h2SNP = 0.26-0.64). Significant variants were predominantly in gene regulatory elements distributed among several haplotype blocks across all chromosomes. Furthermore, traits were constrained by frequent epistatic and pleiotropic interactions.

CONCLUSIONS

Our results on the genetic basis for drought traits in Corymbia calophylla have several implications for the ability to adapt to climate change: (1) drought related traits are controlled by complex genomic architectures with large haplotypes, epistatic, and pleiotropic interactions; (2) the most significant variants determining drought related traits occurred in regulatory regions; and (3) models incorporating epistatic interactions increase trait predictions. Our findings indicate that despite moderate heritability drought traits are likely constrained by complex genomic architecture potentially limiting trees response to climate change.

Chromosomal Inversions and the Demography of Speciation in Drosophila montana and Drosophila flavomontana

Chromosomal inversions may play a central role in speciation given their ability to locally reduce recombination and therefore genetic exchange between diverging populations. We analyzed long- and short-read whole-genome data from sympatric and allopatric populations of 2 Drosophila virilis group species, Drosophila montana and Drosophila flavomontana, to understand if inversions have contributed to their divergence. We identified 3 large alternatively fixed inversions on the X chromosome and one on each of the autosomes 4 and 5. A comparison of demographic models estimated for inverted and noninverted (colinear) chromosomal regions suggests that these inversions arose before the time of the species split. We detected a low rate of interspecific gene flow (introgression) from D. montana to D. flavomontana, which was further reduced inside inversions and was lower in allopatric than in sympatric populations. Together, these results suggest that the inversions were already present in the common ancestral population and that gene exchange between the sister taxa was reduced within inversions both before and after the onset of species divergence. Such ancestrally polymorphic inversions may foster speciation by allowing the accumulation of genetic divergence in loci involved in adaptation and reproductive isolation inside inversions early in the speciation process, while gene exchange at colinear regions continues until the evolving reproductive barriers complete speciation. The overlapping X inversions are particularly good candidates for driving the speciation process of D. montana and D. flavomontana, since they harbor strong genetic incompatibilities that were detected in a recent study of experimental introgression.

Two teosintes made modern maize

The origins of maize were the topic of vigorous debate for nearly a century, but neither the current genetic model nor earlier archaeological models account for the totality of available data, and recent work has highlighted the potential contribution of a wild relative, Zea mays ssp. mexicana. Our population genetic analysis reveals that the origin of modern maize can be traced to an admixture between ancient maize and Zea mays ssp. mexicana in the highlands of Mexico some 4000 years after domestication began. We show that variation in admixture is a key component of maize diversity, both at individual loci and for additive genetic variation underlying agronomic traits. Our results clarify the origin of modern maize and raise new questions about the anthropogenic mechanisms underlying dispersal throughout the Americas.

Genomics of plant speciation

Studies of plants have been instrumental for revealing how new species originate. For several decades, botanical research has complemented and, in some cases, challenged concepts on speciation developed via the study of other organisms while also revealing additional ways in which species can form. Now, the ability to sequence genomes at an unprecedented pace and scale has allowed biologists to settle decades-long debates and tackle other emerging challenges in speciation research. Here, we review these recent genome-enabled developments in plant speciation. We discuss complications related to identification of reproductive isolation (RI) loci using analyses of the landscape of genomic divergence and highlight the important role that structural variants have in speciation, as increasingly revealed by new sequencing technologies. Further, we review how genomics has advanced what we know of some routes to new species formation, like hybridization or whole-genome duplication, while casting doubt on others, like population bottlenecks and genetic drift. While genomics can fast-track identification of genes and mutations that confer RI, we emphasize that follow-up molecular and field experiments remain critical. Nonetheless, genomics has clarified the outsized role of ancient variants rather than new mutations, particularly early during speciation. We conclude by highlighting promising avenues of future study. These include expanding what we know so far about the role of epigenetic and structural changes during speciation, broadening the scope and taxonomic breadth of plant speciation genomics studies, and synthesizing information from extensive genomic data that have already been generated by the plant speciation community.

Chromosomal inversion polymorphisms shape the genomic landscape of deer mice

Chromosomal inversions are an important form of structural variation that can affect recombination, chromosome structure and fitness. However, because inversions can be challenging to detect, the prevalence and hence the significance of inversions segregating within species remains largely unknown, especially in natural populations of mammals. Here, by combining population-genomic and long-read sequencing analyses in a single, widespread species of deer mouse (Peromyscus maniculatus), we identified 21 polymorphic inversions that are large (1.5-43.8 Mb) and cause near-complete suppression of recombination when heterozygous (0-0.03 cM Mb-1). We found that inversion breakpoints frequently occur in centromeric and telomeric regions and are often flanked by long inverted repeats (0.5-50 kb), suggesting that they probably arose via ectopic recombination. By genotyping inversions in populations across the species' range, we found that the inversions are often widespread and do not harbour deleterious mutational loads, and many are likely to be maintained as polymorphisms by divergent selection. Comparisons of forest and prairie ecotypes of deer mice revealed 13 inversions that contribute to differentiation between populations, of which five exhibit significant associations with traits implicated in local adaptation. Taken together, these results show that inversion polymorphisms have a significant impact on recombination, genome structure and genetic diversity in deer mice and likely facilitate local adaptation across the widespread range of this species.

Genome-wide association study reveals structural chromosome variations with phenotypic effects in wheat (Triticum aestivum L.)

Structural chromosome variations (SCVs) are large-scale genomic variations that can be detected by fluorescence in situ hybridization (FISH). SCVs have played important roles in the genome evolution of wheat (Triticum aestivum L.), but little is known about their genetic effects. In this study, a total of 543 wheat accessions from the Chinese wheat mini-core collection and the Shanxi Province wheat collection were used for chromosome analysis using oligonucleotide probe multiplex FISH. A total of 139 SCVs including translocations, pericentric inversions, presence/absence variations (PAVs), and copy number variations (CNVs) in heterochromatin were identified at 230 loci. The distribution frequency of SCVs varied between ecological regions and between landraces and modern cultivars. Structural analysis using SCVs as markers clearly divided the landraces and modern cultivars into different groups. There are very clear instances illustrating alien introgression and wide application of foreign germplasms improved the chromosome diversity of Chinese modern wheat cultivars. A genome-wide association study (GWAS) identified 29 SCVs associated with 12 phenotypic traits, and five (RT4AS•4AL-1DS/1DL•1DS-4AL, Mg2A-3, Mr3B-10, Mr7B-13, and Mr4A-7) of them were further validated using a doubled haploid population and advanced sib-lines, implying the potential value of these SCVs. Importantly, the number of favored SCVs that were associated with agronomic trait improvement was significantly higher in modern cultivars compared to landraces, indicating positive selection in wheat breeding. This study demonstrates the significant effects of SCVs during wheat breeding and provides an efficient method of mining favored SCVs in wheat and other crops.

Genome sequencing reveals evidence of adaptive variation in the genus Zea

Maize is a globally valuable commodity and one of the most extensively studied genetic model organisms. However, we know surprisingly little about the extent and potential utility of the genetic variation found in wild relatives of maize. Here, we characterize a high-density genomic variation map from 744 genomes encompassing maize and all wild taxa of the genus Zea, identifying over 70 million single-nucleotide polymorphisms. The variation map reveals evidence of selection within taxa displaying novel adaptations. We focus on adaptive alleles in highland teosinte and temperate maize, highlighting the key role of flowering-time-related pathways in their adaptation. To show the utility of variants in these data, we generate mutant alleles for two flowering-time candidate genes. This work provides an extensive sampling of the genetic diversity of Zea, resolving questions on evolution and identifying adaptive variants for direct use in modern breeding.

Temperature-associated selection linked to putative chromosomal inversions in king scallop (Pecten maximus)

The genomic landscape of divergence—the distribution of differences among populations or species across the genome—is increasingly characterized to understand the role that microevolutionary forces such as natural selection and recombination play in causing and maintaining genetic divergence. This line of inquiry has also revealed chromosome structure variation to be an important factor shaping the landscape of adaptive genetic variation. Owing to a high prevalence of chromosome structure variation and the strong pressure for local adaptation necessitated by their sessile nature, bivalve molluscs are an ideal taxon for exploring the relationship between chromosome structure variation and local adaptation. Here, we report a population genomic survey of king scallop (Pecten maximus) across its natural range in the northeastern Atlantic Ocean, using a recent chromosome-level genome assembly. We report the presence of at least three large (12–22 Mb), putative chromosomal inversions associated with sea surface temperature and whose frequencies are in contrast to neutral population structure. These results highlight a potentially large role for recombination-suppressing chromosomal inversions in local adaptation and suggest a hypothesis to explain the maintenance of differences in reproductive timing found at relatively small spatial scales across king scallop populations.

The evolutionary patterns of barley pericentromeric chromosome regions, as shaped by linkage disequilibrium and domestication

The distribution of recombination events along large cereal chromosomes is uneven and is generally restricted to gene-rich telomeric ends. To understand how the lack of recombination affects diversity in the large pericentromeric regions, we analysed deep exome capture data from a final panel of 815 Hordeum vulgare (barley) cultivars, landraces and wild barleys, sampled from across their eco-geographical ranges. We defined and compared variant data across the pericentromeric and non-pericentromeric regions, observing a clear partitioning of diversity both within and between chromosomes and germplasm groups. Dramatically reduced diversity was found in the pericentromeres of both cultivars and landraces when compared with wild barley. We observed a mixture of completely and partially differentiated single-nucleotide polymorphisms (SNPs) between domesticated and wild gene pools, suggesting that domesticated gene pools were derived from multiple wild ancestors. Patterns of genome-wide linkage disequilibrium, haplotype block size and number, and variant frequency within blocks showed clear contrasts among individual chromosomes and between cultivars and wild barleys. Although most cultivar chromosomes shared a single major pericentromeric haplotype, chromosome 7H clearly differentiated the two-row and six-row types associated with different geographical origins. Within the pericentromeric regions we identified 22 387 non-synonymous SNPs, 92 of which were fixed for alternative alleles in cultivar versus wild accessions. Surprisingly, only 29 SNPs found exclusively in the cultivars were predicted to be 'highly deleterious'. Overall, our data reveal an unconventional pericentromeric genetic landscape among distinct barley gene pools, with different evolutionary processes driving domestication and diversification.

A chromosomal inversion contributes to divergence in multiple traits between deer mouse ecotypes

How locally adapted ecotypes are established and maintained within a species is a long-standing question in evolutionary biology. Using forest and prairie ecotypes of deer mice (Peromyscus maniculatus), we characterized the genetic basis of variation in two defining traits-tail length and coat color-and discovered a 41-megabase chromosomal inversion linked to both. The inversion frequency is 90% in the dark, long-tailed forest ecotype; decreases across a habitat transition; and is absent from the light, short-tailed prairie ecotype. We implicate divergent selection in maintaining the inversion at frequencies observed in the wild, despite high levels of gene flow, and explore fitness benefits that arise from suppressed recombination within the inversion. We uncover a key role for a large, previously uncharacterized inversion in the evolution and maintenance of classic mammalian ecotypes.

An accurate assignment test for extremely low-coverage whole-genome sequence data

Genomic assignment tests can provide important diagnostic biological characteristics, such as population of origin or ecotype. Yet, assignment tests often rely on moderate- to high-coverage sequence data that can be difficult to obtain for fields such as molecular ecology and ancient DNA. We have developed a novel approach that efficiently assigns biologically relevant information (i.e., population identity or structural variants such as inversions) in extremely low-coverage sequence data. First, we generate databases from existing reference data using a subset of diagnostic single nucleotide polymorphisms (SNPs) associated with a biological characteristic. Low-coverage alignment files are subsequently compared to these databases to ascertain allelic state, yielding a joint probability for each association. To assess the efficacy of this approach, we assigned haplotypes and population identity in Heliconius butterflies, Atlantic herring, and Atlantic cod using chromosomal inversion sites and whole-genome data. We scored both modern and ancient specimens, including the first whole-genome sequence data recovered from ancient Atlantic herring bones. The method accurately assigns biological characteristics, including population membership, using extremely low-coverage data (as low as 0.0001x) based on genome-wide SNPs. This approach will therefore increase the number of samples in evolutionary, ecological and archaeological research for which relevant biological information can be obtained.

Genomic structural variants constrain and facilitate adaptation in natural populations of Theobroma cacao, the chocolate tree

Genomic structural variants (SVs) can play important roles in adaptation and speciation. Yet the overall fitness effects of SVs are poorly understood, partly because accurate population-level identification of SVs requires multiple high-quality genome assemblies. Here, we use 31 chromosome-scale, haplotype-resolved genome assemblies of Theobroma cacao-an outcrossing, long-lived tree species that is the source of chocolate-to investigate the fitness consequences of SVs in natural populations. Among the 31 accessions, we find over 160,000 SVs, which together cover eight times more of the genome than single-nucleotide polymorphisms and short indels (125 versus 15 Mb). Our results indicate that a vast majority of these SVs are deleterious: they segregate at low frequencies and are depleted from functional regions of the genome. We show that SVs influence gene expression, which likely impairs gene function and contributes to the detrimental effects of SVs. We also provide empirical support for a theoretical prediction that SVs, particularly inversions, increase genetic load through the accumulation of deleterious nucleotide variants as a result of suppressed recombination. Despite the overall detrimental effects, we identify individual SVs bearing signatures of local adaptation, several of which are associated with genes differentially expressed between populations. Genes involved in pathogen resistance are strongly enriched among these candidates, highlighting the contribution of SVs to this important local adaptation trait. Beyond revealing empirical evidence for the evolutionary importance of SVs, these 31 de novo assemblies provide a valuable resource for genetic and breeding studies in T cacao.

A Draft Genome of the Ginger Species Alpinia nigra and New Insights into the Genetic Basis of Flexistyly

Angiosperms possess various strategies to ensure reproductive success, such as stylar polymorphisms that encourage outcrossing. Here, we investigate the genetic basis of one such dimorphism that combines both temporal and spatial separation of sexual function, termed flexistyly. It is a floral strategy characterised by the presence of two morphs that differ in the timing of stylar movement. We performed a de novo assembly of the genome of Alpinia nigra using high-depth genomic sequencing. We then used Pool-seq to identify candidate regions for flexistyly based on allele frequency or coverage differences between pools of anaflexistylous and cataflexistylous morphs. The final genome assembly size was 2 Gb, and showed no evidence of recent polyploidy. The Pool-seq did not reveal large regions with high FST values, suggesting large structural chromosomal polymorphisms are unlikely to underlie differences between morphs. Similarly, no region had a 1:2 mapping depth ratio which would be indicative of hemizygosity. We propose that flexistyly is governed by a small genomic region that might be difficult to detect with Pool-seq, or a complex genomic region that proved difficult to assemble. Our genome will be a valuable resource for future studies of gingers, and provides the first steps towards characterising this complex floral phenotype.

Locally Adaptive Inversions Modulate Genetic Variation at Different Geographic Scales in a Seaweed Fly

Across a species range, multiple sources of environmental heterogeneity, at both small and large scales, create complex landscapes of selection, which may challenge adaptation, particularly when gene flow is high. One key to multidimensional adaptation may reside in the heterogeneity of recombination along the genome. Structural variants, like chromosomal inversions, reduce recombination, increasing linkage disequilibrium among loci at a potentially massive scale. In this study, we examined how chromosomal inversions shape genetic variation across a species range and ask how their contribution to adaptation in the face of gene flow varies across geographic scales. We sampled the seaweed fly Coelopa frigida along a bioclimatic gradient stretching across 10° of latitude, a salinity gradient, and a range of heterogeneous, patchy habitats. We generated a chromosome-level genome assembly to analyze 1,446 low-coverage whole genomes collected along those gradients. We found several large nonrecombining genomic regions, including putative inversions. In contrast to the collinear regions, inversions and low-recombining regions differentiated populations more strongly, either along an ecogeographic cline or at a fine-grained scale. These genomic regions were associated with environmental factors and adaptive phenotypes, albeit with contrasting patterns. Altogether, our results highlight the importance of recombination in shaping adaptation to environmental heterogeneity at local and large scales.

The genomic basis of geographic differentiation and fiber improvement in cultivated cotton

Large-scale genomic surveys of crop germplasm are important for understanding the genetic architecture of favorable traits. The genomic basis of geographic differentiation and fiber improvement in cultivated cotton is poorly understood. Here, we analyzed 3,248 tetraploid cotton genomes and confirmed that the extensive chromosome inversions on chromosomes A06 and A08 underlies the geographic differentiation in cultivated Gossypium hirsutum. We further revealed that the haplotypic diversity originated from landraces, which might be essential for understanding adaptative evolution in cultivated cotton. Introgression and association analyses identified new fiber quality-related loci and demonstrated that the introgressed alleles from two diploid cottons had a large effect on fiber quality improvement. These loci provided the potential power to overcome the bottleneck in fiber quality improvement. Our study uncovered several critical genomic signatures generated by historical breeding effects in cotton and a wealth of data that enrich genomic resources for the research community.

Detection of ploidy and chromosomal aberrations in commercial oil palm using high-throughput SNP markers

Karyotyping using high-density genome-wide SNP markers identified various chromosomal aberrations in oil palm (Elaeis guineensis Jacq.) with supporting evidence from the 2C DNA content measurements (determined using FCM) and chromosome counts. Oil palm produces a quarter of the world's total vegetable oil. In line with its global importance, an initiative to sequence the oil palm genome was carried out successfully, producing huge amounts of sequence information, allowing SNP discovery. High-capacity SNP genotyping platforms have been widely used for marker-trait association studies in oil palm. Besides genotyping, a SNP array is also an attractive tool for understanding aberrations in chromosome inheritance. Exploiting this, the present study utilized chromosome-wide SNP allelic distributions to determine the ploidy composition of over 1,000 oil palms from a commercial F₁ family, including 197 derived from twin-embryo seeds. Our method consisted of an inspection of the allelic intensity ratio using SNP markers. For palms with a shifted or abnormal distribution ratio, the SNP allelic frequencies were plotted along the pseudo-chromosomes. This method proved to be efficient in identifying whole genome duplication (triploids) and aneuploidy. We also detected several loss of heterozygosity regions which may indicate small chromosomal deletions and/or inheritance of identical by descent regions from both parents. The SNP analysis was validated by flow cytometry and chromosome counts. The triploids were all derived from twin-embryo seeds. This is the first report on the efficiency and reliability of SNP array data for karyotyping oil palm chromosomes, as an alternative to the conventional cytogenetic technique. Information on the ploidy composition and chromosomal structural variation can help to better understand the genetic makeup of samples and lead to a more robust interpretation of the genomic data in marker-trait association analyses.

CRISPR/Cas-mediated chromosome engineering: opening up a new avenue for plant breeding

The advent of powerful site-specific nucleases, particularly the clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein (Cas) system, which enables precise genome manipulation, has revolutionized plant breeding. Until recently, the main focus of researchers has been to simply knock-in or knock-out single genes, or to induce single base changes, but constant improvements of this technology have enabled more ambitious applications that aim to improve plant productivity or other desirable traits. One long-standing aim has been the induction of targeted chromosomal rearrangements (crossovers, inversions, or translocations). The feasibility of this technique has the potential to transform plant breeding, because natural rearrangements, like inversions, for example, typically present obstacles to the breeding process. In this way, genetic linkages between traits could be altered to combine or separate favorable and deleterious genes, respectively. In this review, we discuss recent breakthroughs in the field of chromosome engineering in plants and their potential applications in the field of plant breeding. In the future, these approaches might be applicable in shaping plant chromosomes in a directed manner, based on plant breeding needs.

Extensive haplotypes are associated with population differentiation and environmental adaptability in Upland cotton (Gossypium hirsutum)

Three extensive eco-haplotypes associated with population differentiation and environmental adaptability in Upland cotton were identified, with A06_85658585, A08_43734499 and A06_113104285 considered the eco-loci for environmental adaptability. Population divergence is suggested to be the primary force driving the evolution of environmental adaptability in various species. Chromosome inversion increases reproductive isolation between subspecies and accelerates population divergence to adapt to new environments. Although modern cultivated Upland cotton (Gossypium hirsutum L.) has spread worldwide, the noticeable phenotypic differences still existed among cultivars grown in different areas. In recent years, the long-distance migration of cotton cultivation areas throughout China has demanded that breeders better understand the genetic basis of environmental adaptability in Upland cotton. Here, we integrated the genotypes of 419 diverse accessions, long-term environment-associated variables (EAVs) and environment-associated traits (EATs) to evaluate subgroup differentiation and identify adaptive loci in Upland cotton. Two highly divergent genomic regions were found on chromosomes A06 and A08, which likely caused by extensive chromosome inversions. The subgroups could be geographically classified based on distinct haplotypes in the divergent regions. A genome-wide association study (GWAS) also confirmed that loci located in these regions were significantly associated with environmental adaptability in Upland cotton. Our study first revealed the cause of population divergence in Upland cotton, as well as the consequences of variation in its environmental adaptability. These findings provide new insights into the genetic basis of environmental adaptability in Upland cotton, which could accelerate the development of molecular markers for adaptation to climate change in future cotton breeding.

Changing local recombination patterns in Arabidopsis by CRISPR/Cas mediated chromosome engineering

Chromosomal inversions are recurrent rearrangements that occur between different plant isolates or cultivars. Such inversions may underlie reproductive isolation in evolution and represent a major obstacle for classical breeding as no crossovers can be observed between inverted sequences on homologous chromosomes. The heterochromatic knob (hk4S) on chromosome 4 is the most well-known inversion of Arabidopsis. If a knob carrying accession such as Col-0 is crossed with a knob-less accession such as Ler-1, crossovers cannot be recovered within the inverted region. Our work shows that by egg-cell specific expression of the Cas9 nuclease from Staphylococcus aureus, a targeted reversal of the 1.1 Mb long hk4S-inversion can be achieved. By crossing Col-0 harbouring the rearranged chromosome 4 with Ler-1, meiotic crossovers can be restored into a region with previously no detectable genetic exchange. The strategy of somatic chromosome engineering for breaking genetic linkage has huge potential for application in plant breeding.

Changing local recombination patterns in Arabidopsis by CRISPR/Cas mediated chromosome engineering

Chromosomal inversions are recurrent rearrangements that occur between different plant isolates or cultivars. Such inversions may underlie reproductive isolation in evolution and represent a major obstacle for classical breeding as no crossovers can be observed between inverted sequences on homologous chromosomes. The heterochromatic knob (hk4S) on chromosome 4 is the most well-known inversion of Arabidopsis. If a knob carrying accession such as Col-0 is crossed with a knob-less accession such as Ler-1, crossovers cannot be recovered within the inverted region. Our work shows that by egg-cell specific expression of the Cas9 nuclease from Staphylococcus aureus, a targeted reversal of the 1.1 Mb long hk4S-inversion can be achieved. By crossing Col-0 harbouring the rearranged chromosome 4 with Ler-1, meiotic crossovers can be restored into a region with previously no detectable genetic exchange. The strategy of somatic chromosome engineering for breaking genetic linkage has huge potential for application in plant breeding.

Footprints of local adaptation span hundreds of linked genes in the Atlantic silverside genome

The study of local adaptation in the presence of ongoing gene flow is the study of natural selection in action, revealing the functional genetic diversity most relevant to contemporary pressures. In addition to individual genes, genome‐wide architecture can itself evolve to enable adaptation. Distributed across a steep thermal gradient along the east coast of North America, Atlantic silversides (Menidia menidia) exhibit an extraordinary degree of local adaptation in a suite of traits, and the capacity for rapid adaptation from standing genetic variation, but we know little about the patterns of genomic variation across the species range that enable this remarkable adaptability. Here, we use low‐coverage, whole‐transcriptome sequencing of Atlantic silversides sampled along an environmental cline to show marked signatures of divergent selection across a gradient of neutral differentiation. Atlantic silversides sampled across 1371 km of the southern section of its distribution have very low genome‐wide differentiation (median F_ST = 0.006 across 1.9 million variants), consistent with historical connectivity and observations of recent migrants. Yet almost 14,000 single nucleotide polymorphisms (SNPs) are nearly fixed (F_ST > 0.95) for alternate alleles. Highly differentiated SNPs cluster into four tight linkage disequilibrium (LD) blocks that span hundreds of genes and several megabases. Variants in these LD blocks are disproportionately nonsynonymous and concentrated in genes enriched for multiple functions related to known adaptations in silversides, including variation in lipid storage, metabolic rate, and spawning behavior. Elevated levels of absolute divergence and demographic modeling suggest selection maintaining divergence across these blocks under gene flow. These findings represent an extreme case of heterogeneity in levels of differentiation across the genome, and highlight how gene flow shapes genomic architecture in continuous populations. Locally adapted alleles may be common features of populations distributed along environmental gradients, and will likely be key to conserving variation to enable future responses to environmental change.

QTL for horticulturally important traits associated with pleiotropic andromonoecy and carpel number loci, and a paracentric inversion in cucumber

The legendary cucumber inbred line WI2757 possesses a rare combination of resistances against nine pathogens, which is an important germplasm for cucumber breeding. However, WI2757 flowers late and does not perform well under field conditions. The genetic basis for horticulturally important traits other than disease resistances in WI2757 is largely unknown. In this study, we conducted QTL mapping using F₂ and recombinant inbred line (RIL) populations from the WI2757 × True Lemon cross that were segregating for multiple traits. Phenotypic data were collected in replicated field trials across multiple years for seven traits including fruit carpel number (CN) and sex expression. A high-density SNP-based genetic map was developed with genotyping by sequencing of the RIL population, which revealed a region on chromosome 1 with strong recombination suppression. The reduced recombination in this region was due to a ~ 10-Mbp paracentric inversion in WI2757 that was confirmed with additional segregation and cytological (FISH) analyses. Thirty-six QTL were detected for flowering time, fruit length (FL), fruit diameter (FD), fruit shape (LD), fruit number (FN), CN, and powdery mildew resistance. Five moderate- or major-effect QTL for FL, FD, LD, and FN inside the inversion are likely the pleiotropic effects of the andromonoecy (m), or the cn locus. The major-effect flowering time QTL ft1.1 was also mapped inside the inversion, which seems to be different from the previously assigned delayed flowering in WI2757. Implications of these findings on the use of WI2757 in cucumber breeding are discussed.

Repetitive sequences and structural chromosome alterations promote intraspecific variations in Zea mays L. karyotype

LTR-retrotransposons, knobs and structural chromosome alterations contribute to shape the structure and organization of the Zea mays karyotype. Our initial nuclear DNA content data of Z. mays accessions revealed an intraspecific variation (2 C = 2.00 pg to 2 C = 6.10 pg), suggesting differences in their karyotypes. We aimed to compare the karyotypes of three Z. mays accessions in search of the differences and similarities among them. Karyotype divergences were demonstrated among the accessions, despite their common chromosome number (2n = 20) and ancestral origin. Cytogenomic analyses showed that repetitive sequences and structural chromosome alterations play a significant role in promoting intraspecific nuclear DNA content variation. In addition, heterozygous terminal deletion in chromosome 3 was pointed out as a cause of lower nuclear 2 C value. Besides this, translocation was also observed in the short arm of chromosome 1. Differently, higher 2 C value was associated with the more abundant distribution of LTR-retrotransposons from the family Grande in the karyotype. Moreover, heteromorphism involving the number and position of the 180-bp knob sequence was found among the accessions. Taken together, we provide insights on the pivotal role played by repetitive sequences and structural chromosome alterations in shaping the karyotype of Z. mays.

The genetic architecture of the maize progenitor, teosinte, and how it was altered during maize domestication

The genetics of domestication has been extensively studied ever since the rediscovery of Mendel's law of inheritance and much has been learned about the genetic control of trait differences between crops and their ancestors. Here, we ask how domestication has altered genetic architecture by comparing the genetic architecture of 18 domestication traits in maize and its ancestor teosinte using matched populations. We observed a strongly reduced number of QTL for domestication traits in maize relative to teosinte, which is consistent with the previously reported depletion of additive variance by selection during domestication. We also observed more dominance in maize than teosinte, likely a consequence of selective removal of additive variants. We observed that large effect QTL have low minor allele frequency (MAF) in both maize and teosinte. Regions of the genome that are strongly differentiated between teosinte and maize (high FST) explain less quantitative variation in maize than teosinte, suggesting that, in these regions, allelic variants were brought to (or near) fixation during domestication. We also observed that genomic regions of high recombination explain a disproportionately large proportion of heritable variance both before and after domestication. Finally, we observed that about 75% of the additive variance in both teosinte and maize is "missing" in the sense that it cannot be ascribed to detectable QTL and only 25% of variance maps to specific QTL. This latter result suggests that morphological evolution during domestication is largely attributable to very large numbers of QTL of very small effect.

A universal chromosome identification system for maize and wild Zea species

Maize was one of the first eukaryotic species in which individual chromosomes can be identified cytologically, which made maize one of the oldest models for genetics and cytogenetics research. Nevertheless, consistent identification of all 10 chromosomes from different maize lines as well as from wild Zea species remains a challenge. We developed a new technique for maize chromosome identification based on fluorescence in situ hybridization (FISH). We developed two oligonucleotide-based probes that hybridize to 24 chromosomal regions. Individual maize chromosomes show distinct FISH signal patterns, which allow universal identification of all chromosomes from different Zea species. We developed karyotypes from three Zea mays subspecies and two additional wild Zea species based on individually identified chromosomes. A paracentric inversion was discovered on the long arm of chromosome 4 in Z. nicaraguensis and Z. luxurians based on modifications of the FISH signal patterns. Chromosomes from these two species also showed distinct distribution patterns of terminal knobs compared with other Zea species. These results support that Z. nicaraguensis and Z. luxurians are closely related species.

Frequency, Origins, and Evolutionary Role of Chromosomal Inversions in Plants

Chromosomal inversions have the potential to play an important role in evolution by reducing recombination between favorable combinations of alleles. Until recently, however, most evidence for their likely importance derived from dipteran flies, whose giant larval salivary chromosomes aided early cytogenetic studies. The widespread application of new genomic technologies has revealed that inversions are ubiquitous across much of the plant and animal kingdoms. Here we review the rapidly accumulating literature on inversions in the plant kingdom and discuss what we have learned about their establishment and likely evolutionary role. We show that inversions are prevalent across a wide range of plant groups. We find that inversions are often associated with locally favored traits, as well as with traits that contribute to assortative mating, suggesting that they may be key to adaptation and speciation in the face of gene flow. We also discuss the role of inversions in sex chromosome formation, and explore possible parallels with inversion establishment on autosomes. The identification of inversion origins, as well as the causal variants within them, will advance our understanding of chromosomal evolution in plants.

Frequency, Origins, and Evolutionary Role of Chromosomal Inversions in Plants

Chromosomal inversions have the potential to play an important role in evolution by reducing recombination between favorable combinations of alleles. Until recently, however, most evidence for their likely importance derived from dipteran flies, whose giant larval salivary chromosomes aided early cytogenetic studies. The widespread application of new genomic technologies has revealed that inversions are ubiquitous across much of the plant and animal kingdoms. Here we review the rapidly accumulating literature on inversions in the plant kingdom and discuss what we have learned about their establishment and likely evolutionary role. We show that inversions are prevalent across a wide range of plant groups. We find that inversions are often associated with locally favored traits, as well as with traits that contribute to assortative mating, suggesting that they may be key to adaptation and speciation in the face of gene flow. We also discuss the role of inversions in sex chromosome formation, and explore possible parallels with inversion establishment on autosomes. The identification of inversion origins, as well as the causal variants within them, will advance our understanding of chromosomal evolution in plants.

From gene editing to genome engineering: restructuring plant chromosomes via CRISPR/Cas

In the last years, tremendous progress has been achieved in the field of gene editing in plants. By the induction of single site-specific double-strand breaks (DSBs), the knockout of genes by non-homologous end joining has become routine in many plant species. Recently, the efficiency of inducing pre-planned mutations by homologous recombination has also been improved considerably. However, very little effort has been undertaken until now to achieve more complex changes in plant genomes by the simultaneous induction of several DSBs. Several reports have been published on the efficient induction of deletions. However, the induction of intrachromosomal inversions and interchromosomal recombination by the use of CRISPR/Cas has only recently been reported. In this review, we want to sum up these results and put them into context with regards to what is known about natural chromosome rearrangements in plants. Moreover, we review the recent progress in CRISPR/Cas-based mammalian chromosomal rearrangements, which might be inspiring for plant biologists. In the long run, the controlled restructuring of plant genomes should enable us to link or break linkage of traits at will, thus defining a new area of plant breeding.

Genomics of Long- and Short-Term Adaptation in Maize and Teosintes

Maize is an excellent model for the study of plant adaptation. Indeed, post domestication maize quickly adapted to a host of new environments across the globe. And work over the last decade has begun to highlight the role of the wild relatives of maize-the teosintes Zea mays ssp. parviglumis and ssp. mexicana-as excellent models for dissecting long-term local adaptation.Although human-driven selection associated with maize domestication has been extensively studied, the genetic basis of natural variation is still poorly understood. Here we review studies on the genetic basis of adaptation and plasticity in maize and its wild relatives. We highlight a range of different processes that contribute to adaptation and discuss evidence from natural, cultivated, and experimental populations. From an applied perspective, understanding the genetic bases of adaptation and the contribution of plasticity will provide us with new tools to both better understand and mitigate the effect of climate changes on natural and cultivated populations.

Common gardens in teosintes reveal the establishment of a syndrome of adaptation to altitude

In plants, local adaptation across species range is frequent. Yet, much has to be discovered on its environmental drivers, the underlying functional traits and their molecular determinants. Genome scans are popular to uncover outlier loci potentially involved in the genetic architecture of local adaptation, however links between outliers and phenotypic variation are rarely addressed. Here we focused on adaptation of teosinte populations along two elevation gradients in Mexico that display continuous environmental changes at a short geographical scale. We used two common gardens, and phenotyped 18 traits in 1664 plants from 11 populations of annual teosintes. In parallel, we genotyped these plants for 38 microsatellite markers as well as for 171 outlier single nucleotide polymorphisms (SNPs) that displayed excess of allele differentiation between pairs of lowland and highland populations and/or correlation with environmental variables. Our results revealed that phenotypic differentiation at 10 out of the 18 traits was driven by local selection. Trait covariation along the elevation gradient indicated that adaptation to altitude results from the assembly of multiple co-adapted traits into a complex syndrome: as elevation increases, plants flower earlier, produce less tillers, display lower stomata density and carry larger, longer and heavier grains. The proportion of outlier SNPs associating with phenotypic variation, however, largely depended on whether we considered a neutral structure with 5 genetic groups (73.7%) or 11 populations (13.5%), indicating that population stratification greatly affected our results. Finally, chromosomal inversions were enriched for both SNPs whose allele frequencies shifted along elevation as well as phenotypically-associated SNPs. Altogether, our results are consistent with the establishment of an altitudinal syndrome promoted by local selective forces in teosinte populations in spite of detectable gene flow. Because elevation mimics climate change through space, SNPs that we found underlying phenotypic variation at adaptive traits may be relevant for future maize breeding.

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.

Environmental Association Identifies Candidates for Tolerance to Low Temperature and Drought

Barley (Hordeum vulgare ssp. vulgare) is cultivated from the equator to the Arctic Circle. The wild progenitor species, Hordeum vulgare ssp. spontaneum, occupies a relatively narrow latitudinal range (∼30 - 40° N) primarily at low elevation (< 1,500 m). Adaptation to the range of cultivation has occurred over ∼8,000 years. The genetic basis of adaptation is amenable to study through environmental association. An advantage of environmental association in a well-characterized crop is that many loci that contribute to climatic adaptation and abiotic stress tolerance have already been identified. This provides the opportunity to determine if environmental association approaches effectively identify these loci of large effect. Using published genotyping from 7,864 SNPs in 803 barley landraces, we examined allele frequency differentiation across multiple partitions of the data and mixed model associations relative to bioclimatic variables. Using newly generated resequencing data from a subset of these landraces, we tested for linkage disequilibrium (LD) between SNPs queried in genotyping and SNPs in neighboring loci. Six loci previously reported to contribute to adaptive differences in flowering time and abiotic stress in barley and six loci previously identified in other plant species were identified in our analyses. In many cases, patterns of LD are consistent with the causative variant occurring in the immediate vicinity of the queried SNP. The identification of barley orthologs to well-characterized genes may provide a new understanding of the nature of adaptive variation and could permit a more targeted use of potentially adaptive variants in barley breeding and germplasm improvement.

Population Genetics of the Highly Polymorphic RPP8 Gene Family

Plant nucleotide-binding domain and leucine-rich repeat containing (NLR) genes provide some of the most extreme examples of polymorphism in eukaryotic genomes, rivalling even the vertebrate major histocompatibility complex. Surprisingly, this is also true in Arabidopsis thaliana, a predominantly selfing species with low heterozygosity. Here, we investigate how gene duplication and intergenic exchange contribute to this extraordinary variation. RPP8 is a three-locus system that is configured chromosomally as either a direct-repeat tandem duplication or as a single copy locus, plus a locus 2 Mb distant. We sequenced 48 RPP8 alleles from 37 accessions of A. thaliana and 12 RPP8 alleles from Arabidopsis lyrata to investigate the patterns of interlocus shared variation. The tandem duplicates display fixed differences and share less variation with each other than either shares with the distant paralog. A high level of shared polymorphism among alleles at one of the tandem duplicates, the single-copy locus and the distal locus, must involve both classical crossing over and intergenic gene conversion. Despite these polymorphism-enhancing mechanisms, the observed nucleotide diversity could not be replicated under neutral forward-in-time simulations. Only by adding balancing selection to the simulations do they approach the level of polymorphism observed at RPP8. In this NLR gene triad, genetic architecture, gene function and selection all combine to generate diversity.

Hybrid Decay: A Transgenerational Epigenetic Decline in Vigor and Viability Triggered in Backcross Populations of Teosinte with Maize

In the course of generating populations of maize with teosinte chromosomal introgressions, an unusual sickly plant phenotype was noted in individuals from crosses with two teosinte accessions collected near Valle de Bravo, Mexico. The plants of these Bravo teosinte accessions appear phenotypically normal themselves and the F1 plants appear similar to typical maize × teosinte F1s. However, upon backcrossing to maize, the BC1 and subsequent generations display a number of detrimental characteristics including shorter stature, reduced seed set, and abnormal floral structures. This phenomenon is observed in all BC individuals and there is no chromosomal segment linked to the sickly plant phenotype in advanced backcross generations. Once the sickly phenotype appears in a lineage, normal plants are never again recovered by continued backcrossing to the normal maize parent. Whole-genome shotgun sequencing reveals a small number of genomic sequences, some with homology to transposable elements, that have increased in copy number in the backcross populations. Transcriptome analysis of seedlings, which do not have striking phenotypic abnormalities, identified segments of 18 maize genes that exhibit increased expression in sickly plants. A de novo assembly of transcripts present in plants exhibiting the sickly phenotype identified a set of 59 upregulated novel transcripts. These transcripts include some examples with sequence similarity to transposable elements and other sequences present in the recurrent maize parent (W22) genome as well as novel sequences not present in the W22 genome. Genome-wide profiles of gene expression, DNA methylation, and small RNAs are similar between sickly plants and normal controls, although a few upregulated transcripts and transposable elements are associated with altered small RNA or methylation profiles. This study documents hybrid incompatibility and genome instability triggered by the backcrossing of Bravo teosinte with maize. We name this phenomenon "hybrid decay" and present ideas on the mechanism that may underlie it.

Large-scale metabolite quantitative trait locus analysis provides new insights for high-quality maize improvement

It is generally recognized that many favorable genes which were lost during domestication, including those related to both nutritional value and stress resistance, remain hidden in wild relatives. To uncover such genes in teosinte, an ancestor of maize, we conducted metabolite profiling in a BC₂ F₇ population generated from a cross between the maize wild relative (Zea mays ssp. mexicana) and maize inbred line Mo17. In total, 65 primary metabolites were quantified in four tissues (seedling-stage leaf, grouting-stage leaf, young kernel and mature kernel) with clear tissue-specific patterns emerging. Three hundred and fifty quantitative trait loci (QTLs) for these metabolites were obtained, which were distributed unevenly across the genome and included two QTL hotspots. Metabolite concentrations frequently increased in the presence of alleles from the teosinte genome while the opposite was observed for grain yield and shape trait QTLs. Combination of the multi-tissue transcriptome and metabolome data provided considerable insight into the metabolic variations between maize and its wild relatives. This study thus identifies favorable genes hidden in the wild relative which should allow us to balance high yield and quality in future modern crop breeding programs.

Divergence with gene flow is driven by local adaptation to temperature and soil phosphorus concentration in teosinte subspecies (Zea mays parviglumis and Zea mays mexicana)

Patterns of genomic divergence between hybridizing taxa can be heterogeneous along the genome. Both differential introgression and local adaptation may contribute to this pattern. Here, we analysed two teosinte subspecies, Zea mays ssp. parviglumis and ssp. mexicana, to test whether their divergence has occurred in the face of gene flow and to infer which environmental variables have been important drivers of their ecological differentiation. We generated 9,780 DArTseqTM SNPs for 47 populations, and used an additional data set containing 33,454 MaizeSNP50 SNPs for 49 populations. With these data, we inferred features of demographic history and performed genome wide scans to determine the number of outlier SNPs associated with climate and soil variables. The two data sets indicate that divergence has occurred or been maintained despite continuous gene flow and/or secondary contact. Most of the significant SNP associations were to temperature and to phosphorus concentration in the soil. A large proportion of these candidate SNPs were located in regions of high differentiation that had been identified previously as putative inversions. We therefore propose that genomic differentiation in teosintes has occurred by a process of adaptive divergence, with putative inversions contributing to reduced gene flow between locally adapted populations.

Efficient induction of heritable inversions in plant genomes using the CRISPR/Cas system

During the evolution of plant genomes, sequence inversions occurred repeatedly making the respective regions inaccessible for meiotic recombination and thus for breeding. Therefore, it is important to develop technologies that allow the induction of inversions within chromosomes in a directed and efficient manner. Using the Cas9 nuclease from Staphylococcus aureus (SaCas9), we were able to obtain scarless heritable inversions with high efficiency in the model plant Arabidopsis thaliana. Via deep sequencing, we defined the patterns of junction formation in wild-type and in the non-homologous end-joining (NHEJ) mutant ku70-1. Surprisingly, in plants deficient of KU70, inversion induction is enhanced, indicating that KU70 is required for tethering the local broken ends together during repair. However, in contrast to wild-type, most junctions are formed by microhomology-mediated NHEJ and thus are imperfect with mainly deletions, making this approach unsuitable for practical applications. Using egg-cell-specific expression of Cas9, we were able to induce heritable inversions at different genomic loci and at intervals between 3 and 18 kb, in the percentage range, in the T1 generation. By screening individual lines, inversion frequencies of up to the 10% range were found in T2. Most of these inversions had scarless junctions and were without any sequence change within the inverted region, making the technology attractive for use in crop plants. Applying our approach, it should be possible to reverse natural inversions and induce artificial ones to break or fix linkages between traits at will.

GOOGA: A platform to synthesize mapping experiments and identify genomic structural diversity

Understanding genomic structural variation such as inversions and translocations is a key challenge in evolutionary genetics. We develop a novel statistical approach to comparative genetic mapping to detect large-scale structural mutations from low-level sequencing data. The procedure, called Genome Order Optimization by Genetic Algorithm (GOOGA), couples a Hidden Markov Model with a Genetic Algorithm to analyze data from genetic mapping populations. We demonstrate the method using both simulated data (calibrated from experiments on Drosophila melanogaster) and real data from five distinct crosses within the flowering plant genus Mimulus. Application of GOOGA to the Mimulus data corrects numerous errors (misplaced sequences) in the M. guttatus reference genome and confirms or detects eight large inversions polymorphic within the species complex. Finally, we show how this method can be applied in genomic scans to improve the accuracy and resolution of Quantitative Trait Locus (QTL) mapping.

Interhomolog polymorphism shapes meiotic crossover within the Arabidopsis RAC1 and RPP13 disease resistance genes

During meiosis, chromosomes undergo DNA double-strand breaks (DSBs), which can be repaired using a homologous chromosome to produce crossovers. Meiotic recombination frequency is variable along chromosomes and tends to concentrate in narrow hotspots. We mapped crossover hotspots located in the Arabidopsis thaliana RAC1 and RPP13 disease resistance genes, using varying haplotypic combinations. We observed a negative non-linear relationship between interhomolog divergence and crossover frequency within the hotspots, consistent with polymorphism locally suppressing crossover repair of DSBs. The fancm, recq4a recq4b, figl1 and msh2 mutants, or lines with increased HEI10 dosage, are known to show increased crossovers throughout the genome. Surprisingly, RAC1 crossovers were either unchanged or decreased in these genetic backgrounds, showing that chromosome location and local chromatin environment are important for regulation of crossover activity. We employed deep sequencing of crossovers to examine recombination topology within RAC1, in wild type, fancm, recq4a recq4b and fancm recq4a recq4b backgrounds. The RAC1 recombination landscape was broadly conserved in the anti-crossover mutants and showed a negative relationship with interhomolog divergence. However, crossovers at the RAC1 5'-end were relatively suppressed in recq4a recq4b backgrounds, further indicating that local context may influence recombination outcomes. Our results demonstrate the importance of interhomolog divergence in shaping recombination within plant disease resistance genes and crossover hotspots.

Modification of meiotic recombination by natural variation in plants

Meiosis is a specialized cell division that produces haploid gametes required for sexual reproduction. During the first meiotic division, homologous chromosomes pair and undergo reciprocal crossing over, which recombines linked sequence variation. Meiotic recombination frequency varies extensively both within and between species. In this review, we will examine the molecular basis of meiotic recombination rate variation, with an emphasis on plant genomes. We first consider cis modification caused by polymorphisms at the site of recombination, or elsewhere on the same chromosome. We review cis effects caused by mismatches within recombining joint molecules, the effect of structural hemizygosity, and the role of specific DNA sequence motifs. In contrast, trans modification of recombination is exerted by polymorphic loci encoding diffusible molecules, which are able to modulate recombination on the same and/or other chromosomes. We consider trans modifiers that act to change total recombination levels, hotspot locations, or interactions between homologous and homeologous chromosomes in polyploid species. Finally, we consider the significance of genetic variation that modifies meiotic recombination for adaptation and evolution of plant species.

Using RNA-Seq for Genomic Scaffold Placement, Correcting Assemblies, and Genetic Map Creation in a Common Brassica rapa Mapping Population

Brassica rapa is a model species for agronomic, ecological, evolutionary, and translational studies. Here, we describe high-density SNP discovery and genetic map construction for a B. rapa recombinant inbred line (RIL) population derived from field collected RNA sequencing (RNA-Seq) data. This high-density genotype data enables the detection and correction of putative genome misassemblies and accurate assignment of scaffold sequences to their likely genomic locations. These assembly improvements represent 7.1-8.0% of the annotated B. rapa genome. We demonstrate how using this new resource leads to a significant improvement for QTL analysis over the current low-density genetic map. Improvements are achieved by the increased mapping resolution and by having known genomic coordinates to anchor the markers for candidate gene discovery. These new molecular resources and improvements in the genome annotation will benefit the Brassicaceae genomics community and may help guide other communities in fine-tuning genome annotations.

Connecting genomic patterns of local adaptation and niche suitability in teosintes

The central abundance hypothesis predicts that local adaptation is a function of the distance to the centre of a species' geographic range. To test this hypothesis, we gathered genomic diversity data from 49 populations, 646 individuals and 33,464 SNPs of two wild relatives of maize, the teosintes Zea mays ssp. parviglumis and Zea. mays. ssp. mexicana. We examined the association between the distance to their climatic and geographic centroids and the enrichment of SNPs bearing signals of adaptation. We identified candidate adaptive SNPs in each population by combining neutrality tests and cline analyses. By applying linear regression models, we found that the number of candidate SNPs is positively associated with niche suitability, while genetic diversity is reduced at the limits of the geographic distribution. Our results suggest that overall, populations located at the limit of the species' niches are adapting locally. We argue that local adaptation to this limit could initiate ecological speciation processes and facilitate adaptation to global change.

A novel Zea mays ssp. mexicana L. MYC-type ICE-like transcription factor gene ZmmICE1, enhances freezing tolerance in transgenic Arabidopsis thaliana

The annual Zea mays ssp. mexicana L., a member of the teosinte group, is a close wild relative of maize and thus can be effectively used in maize improvement. In this study, an ICE-like gene, ZmmICE1, was isolated from a cDNA library of RNA-Seq from cold-treated seedling tissues of Zea mays ssp. mexicana L. The deduced protein of ZmmICE1 contains a highly conserved basic helix-loop-helix (bHLH) domain and C-terminal region of ICE-like proteins. The ZmmICE1 protein localizes to the nucleus and shows sumoylation when expressed in an Escherichia coli reconstitution system. In addition, yeast one hybrid assays indicated that ZmmICE1 has transactivation activities. Moreover, ectopic expression of ZmmICE1 in the Arabidopsis ice1-2 mutant increased freezing tolerance. The ZmmICE1 overexpressed plants showed lower electrolyte leakage (EL), reduced contents of malondialdehyde (MDA). The expression of downstream cold related genes of Arabidopsis C-repeat-binding factors (AtCBF1, AtCBF2 and AtCBF3), cold-responsive genes (AtCOR15A and AtCOR47), kinesin-1 member gene (AtKIN1) and responsive to desiccation gene (AtRD29A) was significantly induced when compared with wild type under low temperature treatment. Taken together, these results indicated that ZmmICE1 is the homolog of Arabidopsis inducer of CBF expression genes (AtICE1/2) and plays an important role in the regulation of freezing stress response.

Population structure and local selection yield high genomic variation in Mimulus guttatus

Across western North America, Mimulus guttatus exists as many local populations adapted to site-specific environmental challenges. Gene flow between locally adapted populations will affect genetic diversity both within demes and across the larger metapopulation. Here, we analyse 34 whole-genome sequences from the intensively studied Iron Mountain population (IM) in conjunction with sequences from 22 Mimulus individuals sampled from across western North America. Three striking features of these data address hypotheses about migration and selection in a locally adapted population. First, we find very high levels of intrapopulation polymorphism (synonymous π = 0.033). Variation outside of genes is likely even higher but difficult to estimate because excessive divergence reduces the efficiency of read mapping. Second, IM exhibits a significantly positive genomewide average for Tajima's D. This indicates allele frequencies are typically more intermediate than expected from neutrality, opposite the pattern observed in many other species. Third, IM exhibits a distinctive haplotype structure with a genomewide excess of positive associations between rarer alleles at linked loci. This suggests an important effect of gene flow from other Mimulus populations, although a residual effect of population founding might also contribute. The combination of multiple analyses, including a novel tree-based analytic method, illustrates how the balance of local selection, limited dispersal and metapopulation dynamics manifests across the genome. The overall genomic pattern of sequence diversity suggests successful gene flow of divergent immigrant genotypes into IM. However, many loci show patterns indicative of local adaptation, particularly at SNPs associated with chromosomal inversions.

Molecular, genetic and evolutionary analysis of a paracentric inversion in Arabidopsis thaliana

Chromosomal inversions can provide windows onto the cytogenetic, molecular, evolutionary and demographic histories of a species. Here we investigate a paracentric 1.17-Mb inversion on chromosome 4 of Arabidopsis thaliana with nucleotide precision of its borders. The inversion is created by Vandal transposon activity, splitting an F-box and relocating a pericentric heterochromatin segment in juxtaposition with euchromatin without affecting the epigenetic landscape. Examination of the RegMap panel and the 1001 Arabidopsis genomes revealed more than 170 inversion accessions in Europe and North America. The SNP patterns revealed historical recombinations from which we infer diverse haplotype patterns, ancient introgression events and phylogenetic relationships. We find a robust association between the inversion and fecundity under drought. We also find linkage disequilibrium between the inverted region and the early flowering Col-FRIGIDA allele. Finally, SNP analysis elucidates the origin of the inversion to South-Eastern Europe approximately 5000 years ago and the FRI-Col allele to North-West Europe, and reveals the spreading of a single haplotype to North America during the 17th to 19th century. The 'American haplotype' was identified from several European localities, potentially due to return migration.