Homoeologous evolution of the allotetraploid genome of Poa annua L.

BACKGROUND: Poa annua (annual bluegrass) is an allotetraploid turfgrass, an agronomically significant weed, and one of the most widely dispersed plant species on earth. Here, we report the chromosome-scale genome assemblies of P. annua's diploid progenitors, P. infirma and P. supina, and use multi-omic analyses spanning all three species to better understand P. annua's evolutionary novelty. RESULTS: We find that the diploids diverged from their common ancestor 5.5 - 6.3 million years ago and hybridized to form P. annua ≤ 50,000 years ago. The diploid genomes are similar in chromosome structure and most notably distinguished by the divergent evolutionary histories of their transposable elements, leading to a 1.7 × difference in genome size. In allotetraploid P. annua, we find biased movement of retrotransposons from the larger (A) subgenome to the smaller (B) subgenome. We show that P. annua's B subgenome is preferentially accumulating genes and that its genes are more highly expressed. Whole-genome resequencing of several additional P. annua accessions revealed large-scale chromosomal rearrangements characterized by extensive TE-downsizing and evidence to support the Genome Balance Hypothesis. CONCLUSIONS: The divergent evolutions of the diploid progenitors played a central role in conferring onto P. annua its remarkable phenotypic plasticity. We find that plant genes (guided by selection and drift) and transposable elements (mostly guided by host immunity) each respond to polyploidy in unique ways and that P. annua uses whole-genome duplication to purge highly parasitized heterochromatic sequences. The findings and genomic resources presented here will enable the development of homoeolog-specific markers for accelerated weed science and turfgrass breeding.

DNA sequence-based mapping and comparative genomics of the St genome of Pseudoroegneria spicata (Pursh) Á. Löve versus wheat (Triticum aestivum L.) and barley (Hordeum vulgare L.).

Bluebunch wheatgrass (referred to as BBWG) [Pseudoroegneria spicata (Pursh) Á. Löve] is an important rangeland Triticeae grass used for forage, conservation, and restoration. This diploid has the basic St genome that occurs also in many polyploid Triticeae species, which serve as a gene reservoir for wheat improvement. Until now, the St genome in diploid species of Pseudoroegneria has not been mapped. Using a double-cross mapping populations, we mapped 230 expressed sequence tag derived simple sequence repeat (EST-SSR) and 3468 genotyping-by-sequencing (GBS) markers to 14 linkage groups (LGs), two each for the seven homologous groups of the St genome. The 227 GBS markers of BBWG that matched those in a previous study helped identify the unclassified seven LGs of the St sub-genome among 21 LGs of Thinopyrum intermedium (Host) Barkworth & D.R. Dewey. Comparisons of GBS sequences in BBWG to whole-genome sequences in bread wheat (Triticum aestivum L.) and barley (Hordeum vulgare L.) revealed that the St genome shared a homology of 35% and 24%, a synteny of 86% and 84%, and a collinearity of 0.85 and 0.86, with ABD and H, respectively. This first-draft molecular map of the St genome will be useful in breeding cereal and forage crops.

Phylogenetic relationships among low-ploidy species of Poa using chloroplast sequences.

Species of the genus Poa are taxonomically and genetically difficult to delineate owing to high and variable polyploidy, aneuploidy, and challenging breeding systems. Approximately 5% of the proposed species in Poa are considered to include or comprise diploids, but very few of those diploids are represented in seed collections. Recent phylogenetic studies of Poa have included some diploid species to elucidate Poa genome relationships. In this study, we build upon that foundation of diploid Poa relationships with additional confirmed diploid species and accessions, and with additional chloroplast sequences. We also include samples of P. pratensis and P. arachnifera to hone in on possible ancestral genomes in these two agronomic and highly polyploidy species. Relative to most species of Poa, Poa section Dioicopoa (P. ligularis, P. iridifolia, and P. arachnifera) contained relatively large chromosomes. Phylogenies were constructed using the TLF gene region and five additional chloroplast genes, and the placement of new species and accessions fit within chloroplast lineages previously reported better than by taxonomic subgenera and sections. Low-ploidy species in the P chloroplast lineage, such as P. iberica and P. remota, grouped closest to P. pratensis.

Association Analysis for Bacterial Spot Resistance in a Directionally Selected Complex Breeding Population of Tomato.

Bacterial spot of tomato is caused by at least four species of Xanthomonas with multiple physiological races. We developed a complex breeding population for simultaneous discovery of marker-trait linkage, validation of existing quantitative trait loci (QTL), and pyramiding of resistance. Six advanced accessions with resistance from distinct sources were crossed in all combinations and their F1 hybrids were intercrossed. Over 1,100 segregating progeny were evaluated in the field following inoculation with X. euvesicatoria race T1 strains. We selected 5% of the most resistant and 5% of the most susceptible progeny for evaluation as plots in two subsequent replicated field trials inoculated with T1 and T3 (X. perforans) strains. The estimated heritability of T1 resistance was 0.32. In order to detect previously reported resistance genes, as well as novel QTL, we explored methods to correct for population structure and analysis based on single markers or haplotypes. Both single-point and haplotype analyses identified strong associations in the genomic regions known to carry Rx-3 (chromosome 5) and Rx-4/Xv3 (chromosome 11). Accounting for kinship and structure generally improved the fit of statistical models. Detection of known loci was improved by adding kinship or a combination of kinship and structure using a Q matrix from model-based clustering. Additional QTL were detected on chromosomes 1, 4, 6, and 7 for T1 resistance and chromosomes 2, 4, and 6 for T3 resistance (P < 0.01). Haplotype analysis improved our ability to trace the origin of positive alleles. These results demonstrate that both known and novel associations can be identified using complex breeding populations that have experienced directional selection.

Chromosome-Scale Genome Assembly and Annotation of Allotetraploid Annual Bluegrass (Poa annua L.).

Poa annua L. is a globally distributed grass with economic and horticultural significance as a weed and as a turfgrass. This dual significance, and its phenotypic plasticity and ecological adaptation, have made P. annua an intriguing plant for genetic and evolutionary studies. Because of the lack of genomic resources and its allotetraploid (2n = 4x = 28) nature, a reference genome sequence would be a valuable asset to better understand the significance and polyploid origin of P. annua. Here we report a genome assembly with scaffolds representing the 14 haploid chromosomes that are 1.78 Gb in length with an N50 of 112 Mb and 96.7% of BUSCO orthologs. Seventy percent of the genome was identified as repetitive elements, 91.0% of which were Copia- or Gypsy-like long-terminal repeats. The genome was annotated with 76,420 genes spanning 13.3% of the 14 chromosomes. The two subgenomes originating from Poa infirma (Knuth) and Poa supina (Schrad) were sufficiently divergent to be distinguishable but syntenic in sequence and annotation with repetitive elements contributing to the expansion of the P. infirma subgenome.

Mapping and linkage disequilibrium analysis with a genome-wide collection of SNPs that detect polymorphism in cultivated tomato.

The history of tomato (Solanum lycopersicum L.) improvement includes genetic bottlenecks, wild species introgressions, and divergence into distinct market classes. This history makes tomato an excellent model to investigate the effects of selection on genome variation. A combination of linkage mapping in two F(2) populations and physical mapping with emerging genome sequence data was used to position 434 PCR-based markers including SNPs. Three-hundred-and-forty markers were used to genotype 102 tomato lines representing wild species, landraces, vintage cultivars, and contemporary (fresh market and processing) varieties. Principal component analysis confirmed genetic divergence between market classes of cultivated tomato (P <0.0001). A genome-wide survey indicated that linkage disequilibrium (LD) decays over 6-8 cM when all cultivated tomatoes, including vintage and contemporary, were considered together. Within contemporary processing varieties, LD decayed over 6-14 cM, and decay was over 3-16 cM within fresh market varieties. Significant inter-chromosomal (gametic phase) LD was detected in both fresh market and processing varieties between chromosomes 2 and 3, and 2 and 4, but in distinct chromosomal locations for each market class. Additional LD was detected between chromosomes 3 and 4, 3 and 11, and 4 and 6 in fresh market varieties and chromosomes 3 and 12 in processing varieties. These results suggest that breeding practices for market specialization in tomato have led to a genetic divergence between fresh market and processing types.

Molecular mapping of hypersensitive resistance from tomato 'Hawaii 7981' to Xanthomonas perforans race T3.

Bacterial spot of tomato (Solanum lycopersicum) is caused by four species of Xanthomonas. The disease causes significant yield losses and a reduction in fruit quality. Physiological races have been described with tomato race 3 (T3) corresponding to strains of Xanthomonas perforans. The breeding line Hawaii 7981 (hereafter H7981) shows a hypersensitive reaction (HR) to race T3 strains conditioned by the interaction of the host resistance locus Xv3 and the bacterial effector avrXv3. The Xv3 gene is required for H7981-derived resistance to be effective under field conditions, though its expression is subject to genetic background. The segregation of HR in F(2) populations derived from H7981 crossed to processing tomato parents OH88119 and OH7870 was studied in 331 progeny, with the two independent crosses providing validation. We screened 453 simple-sequence repeat, insertion/deletion, and single-nucleotide polymorphism markers and identified 44 polymorphic markers each for the OH88119 and OH7870 populations covering 84.6 and 73.3% of the genome, respectively, within 20 centimorgans (cM). Marker-trait analysis using all polymorphic markers demonstrated that Xv3-mediated resistance maps to chromosome 11 in the two independent crosses. Allelism tests were conducted in crosses between lines carrying Xv3 derived from H7981, Rx-4 derived from plant introduction (PI) 128216, and resistance derived from PI 126932. These allelism tests suggested that the loci conditioning HR to race T3 strains are linked within 0.1 cM, are allelic, or are the same gene.

Stopped-flow kinetic analysis using Hadamard transform time-of-flight mass spectrometry.

A home-built stopped-flow apparatus is interfaced to a Hadamard transform time-of-flight mass spectrometer, which permits study of reaction kinetics with a time between reaction initiation and observation as short as about 100 ms and a sampling rate of chemical change that can approach 1 ms. This technique is applied to the trypsin-catalyzed hydrolysis of several peptides and is validated by comparing the results with literature values as well as to optical data obtained with the present stopped-flow apparatus. In addition, we report a kinetic study of the action of trypsin on a peptide having more than one cleavage site.

Oligonucleotide array discovery of polymorphisms in cultivated tomato (Solanum lycopersicum L.) reveals patterns of SNP variation associated with breeding.

BACKGROUND: Cultivated tomato (Solanum lycopersicum L.) has narrow genetic diversity that makes it difficult to identify polymorphisms between elite germplasm. We explored array-based single feature polymorphism (SFP) discovery as a high-throughput approach for marker development in cultivated tomato. RESULTS: Three varieties, FL7600 (fresh-market), OH9242 (processing), and PI114490 (cherry) were used as a source of genomic DNA for hybridization to oligonucleotide arrays. Identification of SFPs was based on outlier detection using regression analysis of normalized hybridization data within a probe set for each gene. A subset of 189 putative SFPs was sequenced for validation. The rate of validation depended on the desired level of significance (alpha) used to define the confidence interval (CI), and ranged from 76% for polymorphisms identified at alpha or= 2 SNPs per locus. We used a subset of validated SNPs for genetic diversity analysis of 92 tomato varieties and accessions. Pairwise estimation of theta (Fst) suggested significant differentiation between collections of fresh-market, processing, vintage, Latin American (landrace), and S. pimpinellifolium accessions. The fresh-market and processing groups displayed high genetic diversity relative to vintage and landrace groups. Furthermore, the patterns of SNP variation indicated that domestication and early breeding practices have led to progressive genetic bottlenecks while modern breeding practices have reintroduced genetic variation into the crop from wild species. Finally, we examined the ratio of non-synonymous (Ka) to synonymous substitutions (Ks) for 20 loci with multiple SNPs (>or= 4 per locus). Six of 20 loci showed ratios of Ka/Ks >or= 0.9. CONCLUSION: Array-based SFP discovery was an efficient method to identify a large number of molecular markers for genetics and breeding in elite tomato germplasm. Patterns of sequence variation across five major tomato groups provided insight into to the effect of human selection on genetic variation.

Desorption electrospray ionization: achieving rapid sampling rates.

The sampling rate and imaging capabilities of desorption electrospray ionization (DESI) are examined using a rotating sample platform combined with Hadamard transform time-of-flight mass spectrometry (HTTOFMS), a multiplexed time-of-flight technique that allows for millisecond acquisition of full mass-to-charge ratio scans. DESI-compatible dyes are used to produce spatially defined sample patterns on poly(methyl methacrylate) discs. Control of disk rotation rate sets the residence time of the sample spots in the DESI plume, and thus the sampling rate. Surface patterns of alternating analytes are spectrally resolved up to 80 samples/s and single-analyte spots up to 50 samples/s. The rapid movement of the surface under the DESI plume allows for high DESI solution flow rates without blurring the chemical information on the surface. Data from multiple rotations can be additively combined, generating a chemical image of the surface with improved signal-to-noise characteristics. This multipass data enables analysis of the rising and falling edges of the analyte signal, placing a lower limit on both the temporal resolution of DESI and the maximum achievable sampling rate. Multipass analysis is proposed as a method for DESI surface imaging.

Characterization of hypersensitive resistance to bacterial spot race T3 (Xanthomonas perforans) from tomato accession PI 128216.

Bacterial spot of tomato is caused by four species of Xanthomonas. The accession PI 128216 (Solanum pimpinellifolium) displays a hypersensitive reaction (HR) to race T3 strains (predominantely Xanthomonas perforans). We developed an inbred backcross (IBC) population (BC(2)S(5), 178 families) derived from PI 128216 and OH88119 (S. lycopersicum) as the susceptible recurrent parent for simultaneous introgression and genetic analysis of the HR response. These IBC families were evaluated in the greenhouse for HR to race T3 strain Xcv761. The IBC population was genotyped with molecular markers distributed throughout the genome in order to identify candidate loci conferring resistance. We treated the IBC population as a hypothesis forming generation to guide validation in subsequent crosses. Nonparametric analysis identified an association between HR and markers clustered on chromosome 11 (P < 0.05 to 0.0001) and chromosome 6 (0.04 > P > 0.002). Further analysis of the IBC population suggested that markers on chromosome 6 and 11 failed to assort independently, a phenomenon known as gametic phase disequilibrium. Therefore, to validate marker-trait linkages, resistant IBC plants were crossed with OH88119 and BC(3)F(2) progeny were evaluated for HR in the greenhouse. In these subsequent populations, the HR response was associated with the chromosome 11 markers (P < 0.0002) but not with the markers on chromosome 6 (P > 0.25). Independent F(2) families were developed by crossing resistant IBC lines to OH8245, OH88119, and OH7530. These populations were genotyped, organized into classes based on chromosome 11 markers, and evaluated for resistance in the field. The PI 128216 locus on chromosome 11 provided resistance that was dependent on gene dosage and genetic background. These results define a single locus, Rx-4, from PI 128216, which provides resistance to bacterial spot race T3, has additive gene action, and is located on chromosome 11.

Continuous time-of-flight ion imaging: application to fragmentation.

We have designed and constructed a continuous imaging reflectron time-of-flight mass spectrometer (TOFMS) that provides a mass spectrum at every pixel of a two-dimensional image with a 100% duty cycle. The technique is based on pseudorandom ion beam modulation and three-dimensional ( x, y, t) ion imaging. We use a multichannel plate detector with a delay-line anode that provides x, y positions and flight times t of every ion arrival event. The precision of the peak heights in the 100% duty cycle mass spectra is shown to be enhanced even at short (10 ms) acquisition times, which should prove useful for the study of solution kinetics or fast chromatographic separations. As a demonstration of the system's capability, we have imaged the fragmented ions that underwent surface-induced dissociation inside the reflectron and the ions that fragmented spontaneously through postsource decay.

Computer-controlled, variable-frequency power supply for driving multipole ion guides.

A high voltage, variable-frequency driver circuit for powering resonant multipole ion guides is presented. Two key features of this design are (1) the use of integrated circuits in the driver stage and (2) the use a stepper motor for tuning a large variable capacitor in the resonant stage. In the present configuration the available frequency range spans a factor of 2. The actual values of the minimum and maximum frequencies depend on the chosen inductor and the capacitance of the ion guide. Feedback allows for stabilized, computer-adjustable rf amplitudes over the range of 5-500 V. The rf power supply was characterized over the range of 350-750 kHz and evaluated by driving a quadrupole ion guide in an electrospray time-of-flight mass spectrometer.

Simple template-based method to produce Bradbury-Nielsen gates.

A Bradbury-Nielsen gate (BNG) consists of two interleaved and electrically isolated sets of wires and can transmit or deflect charged particles by applying a varying voltage difference across the two wire sets. We present a simple template-based method to fabricate BNGs with wire spacings as small as 50 microm with minimal use of a microscope. The small wire spacing allows modulation rates at tens of megahertz. Using this method, we have fabricated four BNGs with wire spacings of 500, 200, 100, and 50 microm using 10 microm gold-coated tungsten wires. The performance of the four BNGs is characterized using an imaging detector and compared with theoretical predictions.

Association of candidate genes with heading date in a diverse Dactylis glomerata population.

Flowering occurs in response to cues from both temperature and photoperiod elicitors in cool-season, long-day forage grasses, and genes involved in sensing the elicitors and inducing downstream flowering responses have been associated with heading date and flowering time in perennial forage grasses as well as cereal grasses. In this study we test for association between orchardgrass (Dactylis glomerata L.) heading date and polymorphisms in the CONSTANS (DgCO1), FLOWERING TIME (DgFT1), a VRN1 like MADS-box (DgMADS), and PHOTOPERIOD (DgPPD1-like) containing genes. A diverse population of 150 genotypes was measured for heading date across three years, genotyped, and candidate genes sequenced. Although pairwise population kinship values were generally low, the genotypes fit into a two-group structure model. Linkage disequilibrium decayed rapidly, reaching r2 levels below 0.2 within the 500bp of each gene. SNPs significantly associated with heading date were detected in equal-dose and tetraploid dosage models. The DgCO1 gene had the most significant polymorphisms and those with the largest effects, while DgMADS had several significant polymorphisms in its first intron with smaller effects. These polymorphisms can be used for further validation, selection, and development of breeding lines of orchardgrass.

Molecular and morphological evidence for Penstemon luculentus (Plantaginaceae): a replacement name for Penstemon fremontii var. glabrescens.

Penstemon luculentus R.L.Johnson & M.R.Stevens, nom. nov. replaces Penstemon fremontii var. glabrescens Dorn & Lichvar. The varietal name glabrescens was not elevated because it was already occupied by Penstemon glabrescens Pennell, a different species. This new arrangement is supported by molecular and morphological evidence. An analysis of genetic diversity in populations of both varieties of Penstemon fremontii Torr. & A. Gray (glabrescens and fremontii) from the Piceance Basin, Colorado, using SSR (simple sequences repeats) or microsatellites markers, revealed significant genetic differentiation between the two. Penstemon fremontii var. glabrescens was also genetically different from Penstemon gibbensii Dorn and Penstemon scariosus var. garrettii (Pennell) N.H. Holmgren. The combination of hirtellous stems, glabrous leaves, non-glandular inflorescence, and long anther hairs distinguish Penstemon luculentus from other morphologically similar species.

Duty cycle and modulation efficiency of two-channel Hadamard transform time-of-flight mass spectrometry.

Hadamard transform time-of-flight mass spectrometry (HT-TOFMS) is based on the pseudorandom gating of ion packets into a time-of-flight mass-to-charge analyzer. In its typical implementation, the technique is able to monitor continuous ion sources with a 50% duty cycle, independent of all other figures of merit. Recently, we have demonstrated that the duty cycle can be extended to 100% using patterned, two-channel detection. Two-channel HT-TOFMS involves the simultaneous optimization of paired one-channel experiments and imposes more stringent conditions to achieve high-quality spectra. An ion modulation device, known as Bradbury-Nielson Gate (BNG), is central to HT-TOFMS. It is an ideal deflection plate, capable of transmitting or deflecting an ion beam according to a known binary sequence without changing the times-of-flight of the ions. Analytical equations are derived that accurately describe the ion modulation process of the BNG as confirmed by good agreement with SimIon simulations and ion beam imaging experiments. From these expressions, the duty cycle and ion modulation efficiency were calculated for various BNG parameters, ion beam characteristics, and detector dimensions, which permit the optimum conditions to be chosen for the two-channel experiment. We conclude that the outer detector should be three times the maximum deflection angle to detect all deflected ions (100% duty cycle) and that the difference between the modulated ion counts in the sequence elements 0 and 1 should be maximized to achieve high modulation efficiency. This condition is best achieved by tight focusing of the ion beam in the center of the inner detector. When both channels are optimized, the two-channel advantage can be exploited to achieve a further improvement over a single-channel experiment.

Tomato Analyzer: a useful software application to collect accurate and detailed morphological and colorimetric data from two-dimensional objects.

Measuring fruit morphology and color traits of vegetable and fruit crops in an objective and reproducible way is important for detailed phenotypic analyses of these traits. Tomato Analyzer (TA) is a software program that measures 37 attributes related to two-dimensional shape in a semi-automatic and reproducible manner. Many of these attributes, such as angles at the distal and proximal ends of the fruit and areas of indentation, are difficult to quantify manually. The attributes are organized in ten categories within the software: Basic Measurement, Fruit Shape Index, Blockiness, Homogeneity, Proximal Fruit End Shape, Distal Fruit End Shape, Asymmetry, Internal Eccentricity, Latitudinal Section and Morphometrics. The last category requires neither prior knowledge nor predetermined notions of the shape attributes, so morphometric analysis offers an unbiased option that may be better adapted to high-throughput analyses than attribute analysis. TA also offers the Color Test application that was designed to collect color measurements from scanned images and allow scanning devices to be calibrated using color standards. TA provides several options to export and analyze shape attribute, morphometric, and color data. The data may be exported to an excel file in batch mode (more than 100 images at one time) or exported as individual images. The user can choose between output that displays the average for each attribute for the objects in each image (including standard deviation), or an output that displays the attribute values for each object on the image. TA has been a valuable and effective tool for indentifying and confirming tomato fruit shape Quantitative Trait Loci (QTL), as well as performing in-depth analyses of the effect of key fruit shape genes on plant morphology. Also, TA can be used to objectively classify fruit into various shape categories. Lastly, fruit shape and color traits in other plant species as well as other plant organs such as leaves and seeds can be evaluated with TA.

Micromachined Bradbury-Nielsen gates.

Bradbury-Nielsen gates (BNGs) are a standard way for gating or steering beams of charged particles in ion mobility spectrometry and time-of-flight mass spectrometry. They consist of a pair of interleaved electrodes that when at the same potential allow ions to pass through the electrodes undeflected and, when a voltage is applied, cause the ions to be deflected from their propagation axis. Previous efforts to construct such devices have relied on mechanical assembly by winding wires across an aperture. We describe a micromachining method for making monolithic BNGs using deep reactive ion etching of silicon-on-insulator wafers. This method enables the creation of electrodes with spacings ranging from 25 to 100 microm with a thickness of 20 microm, covering a 5 mm by 5 mm active area. We characterize the performance of these micromachined BNGs by ion imaging in a pseudorandom time-of-flight mass spectrometer.