Genome sequencing reveals the genetic architecture of heterostyly and domestication history of common buckwheat.

Common buckwheat, Fagopyrum esculentum, is an orphan crop domesticated in southwest China that exhibits heterostylous self-incompatibility. Here we present chromosome-scale assemblies of a self-compatible F. esculentum accession and a self-compatible wild relative, Fagopyrum homotropicum, together with the resequencing of 104 wild and cultivated F. esculentum accessions. Using these genomic data, we report the roles of transposable elements and whole-genome duplications in the evolution of Fagopyrum. In addition, we show that (1) the breakdown of heterostyly occurs through the disruption of a hemizygous gene jointly regulating the style length and female compatibility and (2) southeast Tibet was involved in common buckwheat domestication. Moreover, we obtained mutants conferring the waxy phenotype for the first time in buckwheat. These findings demonstrate the utility of our F. esculentum assembly as a reference genome and promise to accelerate buckwheat research and breeding.

ABCA1 is an extracellular phospholipid translocase.

Production of high density lipoprotein (HDL) requires ATP-binding cassette transporter A1 (ABCA1) to drive phospholipid (PL) from the plasma membrane into extracellular apolipoprotein A-I. Here, we use simulations to show that domains of ABCA1 within the plasma membrane remove PL from the membrane's outer leaflet. In our simulations, after the lipid diffuses into the interior of ABCA1's outward-open cavity, PL extracted by the gateway passes through a ring-shaped domain, the annulus orifice, which forms the base of an elongated hydrophobic tunnel in the transporter's extracellular domain. Engineered mutations in the gateway and annulus strongly inhibit lipid export by ABCA1 without affecting cell-surface expression levels. Our finding that ABCA1 extracts lipid from the outer face of the plasma membrane and forces it through its gateway and annulus into an elongated hydrophobic tunnel contrasts with the alternating access model, which proposes that ABCA1 flops PL substrate from the inner leaflet to the outer leaflet of the membrane. Consistent with our model, ABCA1 lacks the charged amino acid residues in the transmembrane domain found in the floppase members of the ABC transporter family.

Population genomic structure of Eurasian and African foxtail millet landrace accessions inferred from genotyping-by-sequencing.

Foxtail millet [Setaria italica (L.) P. Beauv.] is the second most important millet species globally and is adapted to cultivation in diverse environments. Like its wild progenitor, green foxtail [S. viridis (L.) P. Beauv.], it is a model species for C4 photosynthetic pathways and stress tolerance genes in related bioenergy crops. We addressed questions regarding the evolution and spread of foxtail millet through a population genomic study of landraces from across its cultivated range in Europe, Asia, and Africa. We sought to determine population genomic structure and the relationship of domesticated lineages relative to green foxtail. Further, we aimed to identify genes involved in environmental stress tolerance that have undergone differential selection between geographical and genetic groups. Foxtail millet landrace accessions (n = 328) and green foxtail accessions (n = 12) were sequenced by genotyping-by-sequencing (GBS). After filtering, 5,677 single nucleotide polymorphisms (SNPs) were retained for the combined foxtail millet-green foxtail dataset and 5,020 for the foxtail millet dataset. We extended geographic coverage of green foxtail by including previously published GBS sequence tags, yielding a 4,515-SNP dataset for phylogenetic reconstruction. All foxtail millet samples were monophyletic relative to green foxtail, suggesting a single origin of foxtail millet, although no group of foxtail millet was clearly the most ancestral. Four genetic clusters were found within foxtail millet, each with a distinctive geographical distribution. These results, together with archaeobotanical evidence, suggest plausible routes of spread of foxtail millet. Selection scans identified nine candidate genes potentially involved in environmental adaptations, particularly to novel climates encountered, as domesticated foxtail millet spread to new altitudes and latitudes.

Carbon and nitrogen isotopic variability in foxtail millet (Setaria italica) with watering regime.

RATIONALE: Carbonised plant remains are analysed for reconstruction of past climates and agricultural regimes. Several recent studies have used C4 plants to address related questions, and correlations between modern C4 plant δ13 C values and rainfall have been found. The millets were important food crops in prehistoric Eurasia, yet little is known about causes of isotopic variation within millet species. Previous research has shown there to be significant isotopic variation between millet accessions. Here we compare isotope ratios from plants grown under different watering regimes. This allows for a consideration of whether or not Setaria italica is a good proxy for environmental reconstruction. METHODS: We compare stable isotope ratios of Setaria italica plants grown in a controlled environment chamber with different watering regimes. We compare the carbon isotope ratios of leaves and grains, and the nitrogen isotope ratios of grains, from 12 accessions of Setaria italica. RESULTS: We find significant isotopic variability between watering regimes. Carbon isotope ratios are positively correlated with water availability, and on average vary by 1.9‰ and 1.7‰ for leaves and grains, respectively. Grain nitrogen isotope ratios also vary with watering regime; however, the highest isotope ratios are found with the 130-mL watering regime. CONCLUSIONS: The carbon isotope ratios of Setaria italica are strongly correlated with water availability. However, the correlation is the opposite to that seen in studies of C3 plants. The difference in isotopic ratio due to watering regime is comparable with that seen between different accessions; thus distinguishing between changing varieties of Setaria italica and changing climate is problematic. In terms of grain nitrogen isotope ratios, the highest δ15 N values were not associated with the lowest watering regime. Again, δ15 N variation is comparable with that which would be expected from an aridity effect or a manuring effect, and thus distinguishing between these factors is probably problematic.

Cooked starchy food in hearths ca. 120 kya and 65 kya (MIS 5e and MIS 4) from Klasies River Cave, South Africa.

Plant carbohydrates currently constitute 55-80% of the modern human diet (FAO and WHO, 1997) and some of today's key global health issues are associated with excessive carbohydrate consumption. However, starch carbohydrate is still a poorly understood element of modern human diet and our past starch diet may provide insights for future research. Despite an archaeological narrative that links our early hominin ancestors to a diet that is rich in roots and tubers, there is little deep time archaeological evidence of human plant starch consumption. Geneticists hypothesise that the duplication of starch digestion genes in early Homo sapiens (∼300 kya), is an adaptive response to an increased starch diet. Here we offer the earliest evidence of identified fragments of charred starch plant tissue (parenchyma) from cave and rock shelter hearths dated to Marine Isotope Stage (MIS) 5e and MIS 4, from the Middle Stone Age (MSA) site of Klasies River main site, South Africa (34.06°S, 24.24°E).

Genetic evidence for a western Chinese origin of broomcorn millet (Panicum miliaceum).

Broomcorn millet (Panicum miliaceum) is a key domesticated cereal that has been associated with the north China centre of agricultural origins. Early archaeobotanical evidence for this crop has generated two major debates. First, its contested presence in pre-7000 cal. BP sites in eastern Europe has admitted the possibility of a western origin. Second, its occurrence in the 7th and 8th millennia cal. BP in diverse regions of northern China is consistent with several possible origin foci, associated with different Neolithic cultures. We used microsatellite and granule-bound starch synthase I (GBSSI) genotype data from 341 landrace samples across Eurasia, including 195 newly genotyped samples from China, to address these questions. A spatially explicit discriminative modelling approach favours an eastern Eurasian origin for the expansion of broomcorn millet. This is consistent with recent archaeobotanical and chronological re-evaluations, and stable isotopic data. The same approach, together with the distribution of GBSSI alleles, is also suggestive that the origin of broomcorn millet expansion was in western China. This second unexpected finding stimulates new questions regarding the ecology of wild millet and vegetation dynamics in China prior to the mid-Holocene domestication of millet. The chronological relationship between population expansion and domestication is unclear, but our analyses are consistent with the western Loess Plateau being at least one region of primary domestication of broomcorn millet. Patterns of genetic variation indicate that this region was the source of populations to the west in Eurasia, which broomcorn probably reached via the Inner Asia Mountain Corridor from the 3rd millennium BC. A secondary westward expansion along the steppe may have taken place from the 2nd millennium BC.

Correction: Journey to the east: Diverse routes and variable flowering times for wheat and barley en route to prehistoric China.

[This corrects the article DOI: 10.1371/journal.pone.0187405.].

Barley heads east: Genetic analyses reveal routes of spread through diverse Eurasian landscapes.

One of the world's most important crops, barley, was domesticated in the Near East around 11,000 years ago. Barley is a highly resilient crop, able to grown in varied and marginal environments, such as in regions of high altitude and latitude. Archaeobotanical evidence shows that barley had spread throughout Eurasia by 2,000 BC. To further elucidate the routes by which barley cultivation was spread through Eurasia, simple sequence repeat (SSR) analysis was used to determine genetic diversity and population structure in three extant barley taxa: domesticated barley (Hordeum vulgare L. subsp. vulgare), wild barley (H. vulgare subsp. spontaneum) and a six-rowed brittle rachis form (H. vulgare subsp. vulgare f. agriocrithon (Åberg) Bowd.). Analysis of data using the Bayesian clustering algorithm InStruct suggests a model with three ancestral genepools, which captures a major split in the data, with substantial additional resolution provided under a model with eight genepools. Our results indicate that H. vulgare subsp. vulgare f. agriocrithon accessions and Tibetan Plateau H. vulgare subsp. spontaneum are closely related to the H. vulgare subsp. vulgare in their vicinity, and are therefore likely to be feral derivatives of H. vulgare subsp. vulgare. Under the eight genepool model, cultivated barley is split into six ancestral genepools, each of which has a distinct distribution through Eurasia, along with distinct morphological features and flowering time phenotypes. The distribution of these genepools and their phenotypic characteristics is discussed together with archaeological evidence for the spread of barley eastwards across Eurasia.

Buckwheat: a crop from outside the major Chinese domestication centres? A review of the archaeobotanical, palynological and genetic evidence.

The two cultivated species of buckwheat, Fagopyrum esculentum (common buckwheat) and F. tataricum (Tartary buckwheat) are Chinese domesticates whose origins are usually thought to lie in upland southwestern China, outside the major centres of agricultural origins associated with rice and millet. Synthesis of the macro- and microfossil evidence for buckwheat cultivation in China found just 26 records across all time periods, of which the majority were pollen finds. There are few or no identifying criteria distinguishing F. esculentum and F. tataricum for any sample type. The earliest plausibly agricultural Fagopyrum occurs in northern China from the mid 6th millennium cal bp. The archaeobotanical record requires reconciliation with biogeographic and genetic inferences of a southwestern Chinese origin for buckwheat. Scrutiny of the genetic data indicates limitations related to sampling, molecular markers and analytical approaches. Common buckwheat may have been domesticated at the range margins of its wild progenitor before its cultivation expanded in the north, mediated by changing ranges of wild species during the Holocene and/or by cultural exchange or movement of early agriculturalists between southwest China, the Chengdu Plain and the southern Loess Plateau. Buckwheat probably became a pan-Eurasian crop by the 3rd millennium cal bp, with the pattern of finds suggesting a route of westward expansion via the southern Himalaya to the Caucasus and Europe.

Journey to the east: Diverse routes and variable flowering times for wheat and barley en route to prehistoric China.

Today, farmers in many regions of eastern Asia sow their barley grains in the spring and harvest them in the autumn of the same year (spring barley). However, when it was first domesticated in southwest Asia, barley was grown between the autumn and subsequent spring (winter barley), to complete their life cycles before the summer drought. The question of when the eastern barley shifted from the original winter habit to flexible growing schedules is of significance in terms of understanding its spread. This article investigates when barley cultivation dispersed from southwest Asia to regions of eastern Asia and how the eastern spring barley evolved in this context. We report 70 new radiocarbon measurements obtained directly from barley grains recovered from archaeological sites in eastern Eurasia. Our results indicate that the eastern dispersals of wheat and barley were distinct in both space and time. We infer that barley had been cultivated in a range of markedly contrasting environments by the second millennium BC. In this context, we consider the distribution of known haplotypes of a flowering-time gene in barley, Ppd-H1, and infer that the distributions of those haplotypes may reflect the early dispersal of barley. These patterns of dispersal resonate with the second and first millennia BC textual records documenting sowing and harvesting times for barley in central/eastern China.

A consensus model of human apolipoprotein A-I in its monomeric and lipid-free state.

Apolipoprotein (apo)A-I is an organizing scaffold protein that is critical to high-density lipoprotein (HDL) structure and metabolism, probably mediating many of its cardioprotective properties. However, HDL biogenesis is poorly understood, as lipid-free apoA-I has been notoriously resistant to high-resolution structural study. Published models from low-resolution techniques share certain features but vary considerably in shape and secondary structure. To tackle this central issue in lipoprotein biology, we assembled a team of structural biologists specializing in apolipoproteins and set out to build a consensus model of monomeric lipid-free human apoA-I. Combining novel and published cross-link constraints, small-angle X-ray scattering (SAXS), hydrogen-deuterium exchange (HDX) and crystallography data, we propose a time-averaged model consistent with much of the experimental data published over the last 40 years. The model provides a long-sought platform for understanding and testing details of HDL biogenesis, structure and function.

Selective sweep on human amylase genes postdates the split with Neanderthals.

Humans have more copies of amylase genes than other primates. It is still poorly understood, however, when the copy number expansion occurred and whether its spread was enhanced by selection. Here we assess amylase copy numbers in a global sample of 480 high coverage genomes and find that regions flanking the amylase locus show notable depression of genetic diversity both in African and non-African populations. Analysis of genetic variation in these regions supports the model of an early selective sweep in the human lineage after the split of humans from Neanderthals which led to the fixation of multiple copies of AMY1 in place of a single copy. We find evidence of multiple secondary losses of copy number with the highest frequency (52%) of a deletion of AMY2A and associated low copy number of AMY1 in Northeast Siberian populations whose diet has been low in starch content.

Intraspecific carbon and nitrogen isotopic variability in foxtail millet (Setaria italica).

RATIONALE: Isotopic palaeodietary studies generally focus on bone collagen from human and/or animal remains. While plant remains are rarely analysed, it is known that plant isotope values can vary as a result of numerous factors, including soil conditions, the environment and type of plant. The millets were important food crops in prehistoric Eurasia, yet little is known about the isotopic differences within millet species. METHODS: Here we compare the stable isotope ratios within and between Setaria italica plants grown in a controlled environment chamber. Using homogenised samples, we compare carbon isotope ratios of leaves and grains, and nitrogen isotope ratios of grains, from 29 accessions of Setaria italica. RESULTS: We find significant isotopic variability within single leaves and panicles, and between leaves and panicles within the same plant, which must be considered when undertaking plant isotope studies. We find that the leaves and grains from the different accessions have a ca 2‰ range in δ(13) C values, while the nitrogen isotope values in the grains have a ca 6‰ range. We also find an average offset of 0.9‰ between leaves and grains in their δ(13) C values. CONCLUSIONS: The variation found is large enough to have archaeological implications and within- and between-plant isotope variability should be considered in isotope studies. The range in δ(15) N values is particularly significant as it is larger than the typical values quoted for a trophic level enrichment, and as such may lead to erroneous interpretations of the amount of animal protein in human or animal diets. It is therefore necessary to account for the variability in plant stable isotope values during palaeodietary reconstructions. Copyright © 2016 John Wiley & Sons, Ltd.

An Evaluation of the Crystal Structure of C-terminal Truncated Apolipoprotein A-I in Solution Reveals Structural Dynamics Related to Lipid Binding.

Apolipoprotein (apo) A-I mediates many of the anti-atherogenic functions attributed to high density lipoprotein. Unfortunately, efforts toward a high resolution structure of full-length apoA-I have not been fruitful, although there have been successes with deletion mutants. Recently, a C-terminal truncation (apoA-I(Δ185-243)) was crystallized as a dimer. The structure showed two helical bundles connected by a long, curved pair of swapped helical domains. To compare this structure to that existing under solution conditions, we applied small angle x-ray scattering and isotope-assisted chemical cross-linking to apoA-I(Δ185-243) in its dimeric and monomeric forms. For the dimer, we found evidence for the shared domains and aspects of the N-terminal bundles, but not the molecular curvature seen in the crystal. We also found that the N-terminal bundles equilibrate between open and closed states. Interestingly, this movement is one of the transitions proposed during lipid binding. The monomer was consistent with a model in which the long shared helix doubles back onto the helical bundle. Combined with the crystal structure, these data offer an important starting point to understand the molecular details of high density lipoprotein biogenesis.

Surface Density-Induced Pleating of a Lipid Monolayer Drives Nascent High-Density Lipoprotein Assembly.

Biogenesis of high-density lipoproteins (HDL) is coupled to the transmembrane protein, ATP-binding cassette transporter A1 (ABCA1), which transports phospholipid (PL) from the inner to the outer membrane monolayer. Using a combination of computational and experimental approaches, we show that increased outer lipid monolayer surface density, driven by excess PL or membrane insertion of amphipathic helices, results in pleating of the outer monolayer to form membrane-attached discoidal bilayers. Apolipoprotein (apo)A-I accelerates and stabilizes the pleats. In the absence of apoA-I, pleats collapse to form vesicles. These results mimic cells overexpressing ABCA1 that, in the absence of apoA-I, form and release vesicles. We conclude that the basic driving force for nascent discoidal HDL assembly is a PL pump-induced surface density increase that produces lipid monolayer pleating. We then argue that ABCA1 forms an extracellular reservoir containing an isolated pressurized lipid monolayer decoupled from the transbilayer density buffering of cholesterol.

A robust all-atom model for LCAT generated by homology modeling.

LCAT is activated by apoA-I to form cholesteryl ester. We combined two structures, phospholipase A2 (PLA2) that hydrolyzes the ester bond at the sn-2 position of oxidized (short) acyl chains of phospholipid, and bacteriophage tubulin PhuZ, as C- and N-terminal templates, respectively, to create a novel homology model for human LCAT. The juxtaposition of multiple structural motifs matching experimental data is compelling evidence for the general correctness of many features of the model: i) The N-terminal 10 residues of the model, required for LCAT activity, extend the hydrophobic binding trough for the sn-2 chain 15-20 Å relative to PLA2. ii) The topography of the trough places the ester bond of the sn-2 chain less than 5 Å from the hydroxyl of the catalytic nucleophile, S181. iii) A ß-hairpin resembling a lipase lid separates S181 from solvent. iv) S181 interacts with three functionally critical residues: E149, that regulates sn-2 chain specificity, and K128 and R147, whose mutations cause LCAT deficiency. Because the model provides a novel explanation for the complicated thermodynamic problem of the transfer of hydrophobic substrates from HDL to the catalytic triad of LCAT, it is an important step toward understanding the antiatherogenic role of HDL in reverse cholesterol transport.

An experimentally robust model of monomeric apolipoprotein A-I created from a chimera of two X-ray structures and molecular dynamics simulations.

High-density lipoprotein (HDL) retards atherosclerosis by accepting cholesterol from the artery wall. However, the structure of the proposed acceptor, monomeric apolipoprotein A-I (apoA-I), the major protein of HDL, is poorly understood. Two published models for monomeric apoA-I used cross-linking distance constraints to derive best fit conformations. This approach has limitations. (i) Cross-linked peptides provide no information about secondary structure. (ii) A protein chain can be folded in multiple ways to create a best fit. (iii) Ad hoc folding of a secondary structure is unlikely to produce a stable orientation of hydrophobic and hydrophilic residues. To address these limitations, we used a different approach. We first noted that the dimeric apoA-I crystal structure, (Δ185-243)apoA-I, is topologically identical to a monomer in which helix 5 forms a helical hairpin, a monomer with a hydrophobic cleft running the length of the molecule. We then realized that a second crystal structure, (Δ1-43)apoA-I, contains a C-terminal structure that fits snuggly via aromatic and hydrophobic interactions into the hydrophobic cleft. Consequently, we combined these crystal structures into an initial model that was subjected to molecular dynamics simulations. We tested the initial and simulated models and the two previously published models in three ways: against two published data sets (domains predicted to be helical by H/D exchange and six spin-coupled residues) and against our own experimentally determined cross-linking distance constraints. We note that the best fit simulation model, superior by all tests to previously published models, has dynamic features of a molten globule with interesting implications for the functions of apoA-I.

Wheat in the Mediterranean revisited--tetraploid wheat landraces assessed with elite bread wheat Single Nucleotide Polymorphism markers.

BACKGROUND: Single Nucleotide Polymorphism (SNP) panels recently developed for the assessment of genetic diversity in wheat are primarily based on elite varieties, mostly those of bread wheat. The usefulness of such SNP panels for studying wheat evolution and domestication has not yet been fully explored and ascertainment bias issues can potentially affect their applicability when studying landraces and tetraploid ancestors of bread wheat. We here evaluate whether population structure and evolutionary history can be assessed in tetraploid landrace wheats using SNP markers previously developed for the analysis of elite cultivars of hexaploid wheat. RESULTS: We genotyped more than 100 tetraploid wheat landraces and wild emmer wheat accessions, some of which had previously been screened with SSR markers, for an existing SNP panel and obtained publically available genotypes for the same SNPs for hexaploid wheat varieties and landraces. Results showed that quantification of genetic diversity can be affected by ascertainment bias but that the effects of ascertainment bias can at least partly be alleviated by merging SNPs to haplotypes. Analyses of population structure and genetic differentiation show strong subdivision between the tetraploid wheat subspecies, except for durum and rivet that are not separable. A more detailed population structure of durum landraces could be obtained than with SSR markers. The results also suggest an emmer, rather than durum, ancestry of bread wheat and with gene flow from wild emmer. CONCLUSIONS: SNP markers developed for elite cultivars show great potential for inferring population structure and can address evolutionary questions in landrace wheat. Issues of marker genome specificity and mapping need, however, to be addressed. Ascertainment bias does not seem to interfere with the ability of a SNP marker system developed for elite bread wheat accessions to detect population structure in other types of wheat.

Reticulate evolution in Panicum (Poaceae): the origin of tetraploid broomcorn millet, P. miliaceum.

Panicum miliaceum (broomcorn millet) is a tetraploid cereal, which was among the first domesticated crops, but is now a minor crop despite its high water use efficiency. The ancestors of this species have not been determined; we aimed to identify likely candidates within the genus, where phylogenies are poorly resolved. Nuclear and chloroplast DNA sequences from P. miliaceum and a range of diploid and tetraploid relatives were used to develop phylogenies of the diploid and tetraploid species. Chromosomal in situ hybridization with genomic DNA as a probe was used to characterize the genomes in the tetraploid P. miliaceum and a tetraploid accession of P. repens. In situ hybridization showed that half the chromosomes of P. miliaceum hybridized more strongly with labelled genomic DNA from P. capillare, and half with labelled DNA from P. repens. Genomic DNA probes differentiated two sets of 18 chromosomes in the tetraploid P. repens. Our phylogenetic data support the allotetraploid origin of P. miliaceum, with the maternal ancestor being P. capillare (or a close relative) and the other genome being shared with P. repens. Our P. repens accession was also an allotetraploid with two dissimilar but closely related genomes, the maternal genome being similar to P. sumatrense. Further collection of Panicum species, particularly from the Old World, is required. It is important to identify why the water-efficient P. miliaceum is now of minimal importance in agriculture, and it may be valuable to exploit the diversity in this species and its ancestors.

The structure of human apolipoprotein A-IV as revealed by stable isotope-assisted cross-linking, molecular dynamics, and small angle x-ray scattering.

Apolipoprotein (apo)A-IV plays important roles in dietary lipid and glucose metabolism, and knowledge of its structure is required to fully understand the molecular basis of these functions. However, typical of the entire class of exchangeable apolipoproteins, its dynamic nature and affinity for lipid has posed challenges to traditional high resolution structural approaches. We previously reported an x-ray crystal structure of a dimeric truncation mutant of apoA-IV, which showed a unique helix-swapping molecular interface. Unfortunately, the structures of the N and C termini that are important for lipid binding were not visualized. To build a more complete model, we used chemical cross-linking to derive distance constraints across the full-length protein. The approach was enhanced with stable isotope labeling to overcome ambiguities in determining molecular span of the cross-links given the remarkable similarities in the monomeric and dimeric apoA-IV structures. Using 51 distance constraints, we created a starting model for full-length monomeric apoA-IV and then subjected it to two modeling approaches: (i) molecular dynamics simulations and (ii) fitting to small angle x-ray scattering data. This resulted in the most detailed models yet for lipid-free monomeric or dimeric apoA-IV. Importantly, these models were of sufficient detail to direct the experimental identification of new functional residues that participate in a "clasp" mechanism to modulate apoA-IV lipid affinity. The isotope-assisted cross-linking approach should prove useful for further study of this family of apolipoproteins in both the lipid-free and -bound states.