Design and Validation of a High-Throughput Reductive Catalytic Fractionation Method.

Reductive catalytic fractionation (RCF) is a promising method to extract and depolymerize lignin from biomass, and bench-scale studies have enabled considerable progress in the past decade. RCF experiments are typically conducted in pressurized batch reactors with volumes ranging between 50 and 1000 mL, limiting the throughput of these experiments to one to six reactions per day for an individual researcher. Here, we report a high-throughput RCF (HTP-RCF) method in which batch RCF reactions are conducted in 1 mL wells machined directly into Hastelloy reactor plates. The plate reactors can seal high pressures produced by organic solvents by vertically stacking multiple reactor plates, leading to a compact and modular system capable of performing 240 reactions per experiment. Using this setup, we screened solvent mixtures and catalyst loadings for hydrogen-free RCF using 50 mg poplar and 0.5 mL reaction solvent. The system of 1:1 isopropanol/methanol showed optimal monomer yields and selectivity to 4-propyl substituted monomers, and validation reactions using 75 mL batch reactors produced identical monomer yields. To accommodate the low material loadings, we then developed a workup procedure for parallel filtration, washing, and drying of samples and a 1H nuclear magnetic resonance spectroscopy method to measure the RCF oil yield without performing liquid-liquid extraction. As a demonstration of this experimental pipeline, 50 unique switchgrass samples were screened in RCF reactions in the HTP-RCF system, revealing a wide range of monomer yields (21-36%), S/G ratios (0.41-0.93), and oil yields (40-75%). These results were successfully validated by repeating RCF reactions in 75 mL batch reactors for a subset of samples. We anticipate that this approach can be used to rapidly screen substrates, catalysts, and reaction conditions in high-pressure batch reactions with higher throughput than standard batch reactors.

Finger millet: a hero in the making to combat food insecurity.

Climate change and population growth pose challenges to food security. Major crops such as maize, wheat, and rice are expected to face yield reductions due to warming in the coming years, highlighting the need for incorporating climate-resilient crops in agricultural production systems. Finger millet (Eleusine coracana (L.) Gaertn) is a nutritious cereal crop adapted to arid regions that could serve as an alternative crop for sustaining the food supply in low rainfall environments where other crops routinely fail. Despite finger millet's nutritional qualities and climate resilience, it is deemed an "orphan crop," neglected by researchers compared to major crops, which has hampered breeding efforts. However, in recent years, finger millet has entered the genomics era. Next-generation sequencing resources, including a chromosome-scale genome assembly, have been developed to support trait characterization. This review discusses the current genetic and genomic resources available for finger millet while addressing the gaps in knowledge and tools that are still needed to aid breeders in bringing finger millet to its full production potential.

Economic and Sustainability Impacts of Yield and Composition Variation in Bioenergy Crops: Switchgrass (Panicum virgatum L.).

Economically viable production of biobased products and fuels requires high-yielding, high-quality, sustainable process-advantaged crops, developed using bioengineering or advanced breeding approaches. Understanding which crop phenotypic traits have the largest impact on biofuel economics and sustainability outcomes is important for the targeted feedstock crop development. Here, we evaluated biomass yield and cell-wall composition traits across a large natural variant population of switchgrass (Panicum virgatum L.) grown across three common garden sites. Samples from 331 switchgrass genotypes were collected and analyzed for carbohydrate and lignin components. Considering plant survival and biomass after multiple years of growth, we found that 84 of the genotypes analyzed may be suited for commercial production in the southeastern U.S. These genotypes show a range of growth and compositional traits across the population that are apparently independent of each other. We used these data to conduct techno-economic analyses and life cycle assessments evaluating the performance of each switchgrass genotype under a standard cellulosic ethanol process model with pretreatment, added enzymes, and fermentation. We find that switchgrass yield per area is the largest economic driver of the minimum fuel selling price (MSFP), ethanol yield per hectare, global warming potential (GWP), and cumulative energy demand (CED). At any yield, the carbohydrate content is significant but of secondary importance. Water use follows similar trends but has more variability due to an increased dependence on the biorefinery model. Analyses presented here highlight the primary importance of plant yield and the secondary importance of carbohydrate content when selecting a feedstock that is both economical and sustainable.

A single amino acid change led to structural and functional differentiation of PvHd1 to control flowering in switchgrass.

Switchgrass, a forage and bioenergy crop, occurs as two main ecotypes with different but overlapping ranges of adaptation. The two ecotypes differ in a range of characteristics, including flowering time. Flowering time determines the duration of vegetative development and therefore biomass accumulation, a key trait in bioenergy crops. No causal variants for flowering time differences between switchgrass ecotypes have, as yet, been identified. In this study, we mapped a robust flowering time quantitative trait locus (QTL) on chromosome 4K in a biparental F2 population and characterized the flowering-associated transcription factor gene PvHd1, an ortholog of CONSTANS in Arabidopsis and Heading date 1 in rice, as the underlying causal gene. Protein modeling predicted that a serine to glycine substitution at position 35 (p.S35G) in B-Box domain 1 greatly altered the global structure of the PvHd1 protein. The predicted variation in protein compactness was supported in vitro by a 4 °C shift in denaturation temperature. Overexpressing the PvHd1-p.35S allele in a late-flowering CONSTANS-null Arabidopsis mutant rescued earlier flowering, whereas PvHd1-p.35G had a reduced ability to promote flowering, demonstrating that the structural variation led to functional divergence. Our findings provide us with a tool to manipulate the timing of floral transition in switchgrass cultivars and, potentially, expand their cultivation range.

Genome analyses reveal population structure and a purple stigma color gene candidate in finger millet.

Finger millet is a key food security crop widely grown in eastern Africa, India and Nepal. Long considered a 'poor man's crop', finger millet has regained attention over the past decade for its climate resilience and the nutritional qualities of its grain. To bring finger millet breeding into the 21st century, here we present the assembly and annotation of a chromosome-scale reference genome. We show that this ~1.3 million years old allotetraploid has a high level of homoeologous gene retention and lacks subgenome dominance. Population structure is mainly driven by the differential presence of large wild segments in the pericentromeric regions of several chromosomes. Trait mapping, followed by variant analysis of gene candidates, reveals that loss of purple coloration of anthers and stigma is associated with loss-of-function mutations in the finger millet orthologs of the maize R1/B1 and Arabidopsis GL3/EGL3 anthocyanin regulatory genes. Proanthocyanidin production in seed is not affected by these gene knockouts.

Host Specificity Controlled by PWL1 and PWL2 Effector Genes in the Finger Millet Blast Pathogen Magnaporthe oryzae in Eastern Africa.

Magnaporthe oryzae, a devastating pathogen of finger millet (Eleusine coracana), secretes effector molecules during infection to manipulate host immunity. This study determined the presence of avirulence effector genes PWL1 and PWL2 in 221 Eleusine blast isolates from eastern Africa. Most Ethiopian isolates carried both PWL1 and PWL2. Kenyan and Ugandan isolates largely lacked both genes, and Tanzanian isolates carried either PWL1 or lacked both. The roles of PWL1 and PWL2 towards pathogenicity on alternative chloridoid hosts, including weeping lovegrass (Eragrostis curvula), were also investigated. PWL1 and PWL2 were cloned from Ethiopian isolate E22 and were transformed separately into Ugandan isolate U34, which lacked both genes. Resulting transformants harboring either gene gained varying degrees of avirulence on Eragrostis curvula but remained virulent on finger millet. Strains carrying one or both PWL1 and PWL2 infected the chloridoid species Sporobolus phyllotrichus and Eleusine tristachya, indicating the absence of cognate resistance (R) genes for PWL1 and PWL2 in these species. Other chloridoid grasses, however, were fully resistant, regardless of the presence of one or both PWL1 and PWL2, suggesting the presence of effective R genes against PWL and other effectors. Partial resistance in some Eragrostis curvula accessions to some blast isolates lacking PWL1 and PWL2 also indicated the presence of other interactions between fungal avirulence (AVR) genes and host resistance (R) genes. Related chloridoid species thus harbor resistance genes that could be useful to improve finger millet for blast resistance. Conversely, loss of AVR genes in the fungus could expand its host range, as demonstrated by the susceptibility of Eragrostis curvula to finger millet blast isolates that had lost PWL1 and PWL2. [Formula: see text] Copyright © 2023 The Author(s). This is an open access article distributed under the CC BY 4.0 International license.

Genome of Paspalum vaginatum and the role of trehalose mediated autophagy in increasing maize biomass.

A number of crop wild relatives can tolerate extreme stress to a degree outside the range observed in their domesticated relatives. However, it is unclear whether or how the molecular mechanisms employed by these species can be translated to domesticated crops. Paspalum (Paspalum vaginatum) is a self-incompatible and multiply stress-tolerant wild relative of maize and sorghum. Here, we describe the sequencing and pseudomolecule level assembly of a vegetatively propagated accession of P. vaginatum. Phylogenetic analysis based on 6,151 single-copy syntenic orthologues conserved in 6 related grass species places paspalum as an outgroup of the maize-sorghum clade. In parallel metabolic experiments, paspalum, but neither maize nor sorghum, exhibits a significant increase in trehalose when grown under nutrient-deficit conditions. Inducing trehalose accumulation in maize, imitating the metabolic phenotype of paspalum, results in autophagy dependent increases in biomass accumulation.

Identification of microRNAs responsive to arbuscular mycorrhizal fungi in Panicum virgatum (switchgrass).

BACKGROUND: MicroRNAs (miRNAs) are important post-transcriptional regulators involved in the control of a range of processes, including symbiotic interactions in plants. MiRNA involvement in arbuscular mycorrhizae (AM) symbiosis has been mainly studied in model species, and our study is the first to analyze global miRNA expression in the roots of AM colonized switchgrass (Panicum virgatum), an emerging biofuel feedstock. AM symbiosis helps plants gain mineral nutrition from the soil and may enhance switchgrass biomass production on marginal lands. Our goals were to identify miRNAs and their corresponding target genes that are controlling AM symbiosis in switchgrass. RESULTS: Through genome-wide analysis of next-generation miRNA sequencing reads generated from switchgrass roots, we identified 122 mature miRNAs, including 28 novel miRNAs. By comparing miRNA expression profiles of AM-inoculated and control switchgrass roots, we identified 15 AM-responsive miRNAs across lowland accession "Alamo", upland accession "Dacotah", and two upland/lowland F1 hybrids. We used degradome sequencing to identify target genes of the AM-responsive miRNAs revealing targets of miRNAs residing on both K and N subgenomes. Notably, genes involved in copper ion binding were targeted by downregulated miRNAs, while upregulated miRNAs mainly targeted GRAS family transcription factors. CONCLUSION: Through miRNA analysis and degradome sequencing, we revealed that both upland and lowland switchgrass genotypes as well as upland-lowland hybrids respond to AM by altering miRNA expression. We demonstrated complex GRAS transcription factor regulation by the miR171 family, with some miR171 family members being AM responsive while others remained static. Copper miRNA downregulation was common amongst the genotypes tested and we identified superoxide dismutases and laccases as targets, suggesting that these Cu-miRNAs are likely involved in ROS detoxification and lignin deposition, respectively. Other prominent targets of the Cu miRNAs were blue copper proteins. Overall, the potential effect of AM colonization on lignin deposition pathways in this biofuel crop highlights the importance of considering AM and miRNA in future biofuel crop development strategies.

Genotyping-by-Sequencing and QTL Mapping of Biomass Yield in Two Switchgrass F1 Populations (Lowland x Coastal and Coastal x Upland).

The prevalence of genetic diversity in switchgrass germplasm can be exploited to capture favorable alleles that increase its range of adaptation and biomass yield. The objectives of the study were to analyze the extent of polymorphism and patterns of segregation distortion in two F1 populations and use the linkage maps to locate QTL for biomass yield. We conducted genotyping-by-sequencing on two populations derived from crosses between the allotetraploid lowland genotype AP13 (a selection from "Alamo") and coastal genotype B6 (a selection from PI 422001) with 285 progeny (AB population) and between B6 and the allotetraploid upland VS16 (a selection from "Summer") with 227 progeny (BV population). As predictable from the Euclidean distance between the parents, a higher number of raw variants was discovered in the coastal × upland BV cross (6 M) compared to the lowland × coastal AB cross (2.5 M). The final number of mapped markers was 3,107 on the BV map and 2,410 on the AB map. More segregation distortion of alleles was seen in the AB population, with 75% distorted loci compared to 11% distorted loci in the BV population. The distortion in the AB population was seen across all chromosomes in both the AP13 and B6 maps and likely resulted from zygotic or post-zygotic selection for increased levels of heterozygosity. Our results suggest lower genetic compatibility between the lowland AP13 and the coastal B6 ecotype than between B6 and the upland ecotype VS16. Four biomass QTLs were mapped in the AB population (LG 2N, 6K, 6N, and 8N) and six QTLs in the BV population [LG 1N (2), 8N (2), 9K, and 9N]. The QTL, with the largest and most consistent effect across years, explaining between 8.4 and 11.5% of the variation, was identified on 6N in the AP13 map. The cumulative effect of all the QTLs explained a sizeable portion of the phenotypic variation in both AB and BV populations and the markers associated with them may potentially be used for the marker-assisted improvement of biomass yield. Since switchgrass improvement is based on increasing favorable allele frequencies through recurrent selection, the transmission bias within individuals and loci needs to be considered as this may affect the genetic gain if the favorable alleles are distorted.

A high-density linkage map of finger millet provides QTL for blast resistance and other agronomic traits.

Finger millet [Eleusine coracana (L.) Gaertn.] is a critical subsistence crop in eastern Africa and southern Asia but has few genomic resources and modern breeding programs. To aid in the understanding of finger millet genomic organization and genes underlying disease resistance and agronomically important traits, we generated a F2:3 population from a cross between E. coracana (L.) Gaertn. subsp. coracana accession ACC 100007 and E. coracana (L.) Gaertn. subsp. africana , accession GBK 030647. Phenotypic data on morphology, yield, and blast (Magnaporthe oryzae) resistance traits were taken on a subset of the F2:3 population in a Kenyan field trial. The F2:3 population was genotyped via genotyping-by-sequencing (GBS) and the UGbS-Flex pipeline was used for sequence alignment, nucleotide polymorphism calling, and genetic map construction. An 18-linkage-group genetic map consisting of 5,422 markers was generated that enabled comparative genomic analyses with rice (Oryza sativa L.), foxtail millet [Setaria italica (L.) P. Beauv.], and sorghum [Sorghum bicolor (L.) Moench]. Notably, we identified conserved acrocentric homoeologous chromosomes (4A and 4B in finger millet) across all species. Significant quantitative trait loci (QTL) were discovered for flowering date, plant height, panicle number, and blast incidence and severity. Sixteen putative candidate genes that may underlie trait variation were identified. Seven LEUCINE-RICH REPEAT-CONTAINING PROTEIN genes, with homology to nucleotide-binding site leucine-rich repeat (NBS-LRR) disease resistance proteins, were found on three chromosomes under blast resistance QTL. This high-marker-density genetic map provides an important tool for plant breeding programs and identifies genomic regions and genes of critical interest for agronomic traits and blast resistance.

Aegilops tauschii genome assembly Aet v5.0 features greater sequence contiguity and improved annotation.

Aegilops tauschii is the donor of the D subgenome of hexaploid wheat and an important genetic resource. The reference-quality genome sequence Aet v4.0 for Ae. tauschii acc. AL8/78 was therefore an important milestone for wheat biology and breeding. Further advances in sequencing acc. AL8/78 and release of the Aet v5.0 sequence assembly are reported here. Two new optical maps were constructed and used in the revision of pseudomolecules. Gaps were closed with Pacific Biosciences long-read contigs, decreasing the gap number by 38,899. Transposable elements and protein-coding genes were reannotated. The number of annotated high-confidence genes was reduced from 39,635 in Aet v4.0 to 32,885 in Aet v5.0. A total of 2245 biologically important genes, including those affecting plant phenology, grain quality, and tolerance of abiotic stresses in wheat, was manually annotated and disease-resistance genes were annotated by a dedicated pipeline. Disease-resistance genes encoding nucleotide-binding site domains, receptor-like protein kinases, and receptor-like proteins were preferentially located in distal chromosome regions, whereas those encoding transmembrane coiled-coil proteins were dispersed more evenly along the chromosomes. Discovery, annotation, and expression analyses of microRNA (miRNA) precursors, mature miRNAs, and phasiRNAs are reported, including miRNA target genes. Other small RNAs, such as hc-siRNAs and tRFs, were characterized. These advances enhance the utility of the Ae. tauschii genome sequence for wheat genetics, biotechnology, and breeding.

Quantitative trait locus mapping combined with variant and transcriptome analyses identifies a cluster of gene candidates underlying the variation in leaf wax between upland and lowland switchgrass ecotypes.

KEY MESSAGE: Mapping combined with expression and variant analyses in switchgrass, a crop with complex genetics, identified a cluster of candidate genes for leaf wax in a fast-evolving region of chromosome 7K. Switchgrass (Panicum virgatum L.) is a promising warm-season candidate energy crop. It occurs in two ecotypes, upland and lowland, which vary in a number of phenotypic traits, including leaf glaucousness. To initiate trait mapping, two F2 mapping populations were developed by crossing two different F1 sibs derived from a cross between the tetraploid lowland genotype AP13 and the tetraploid upland genotype VS16, and high-density linkage maps were generated. Quantitative trait locus (QTL) analyses of visually scored leaf glaucousness and of hydrophobicity of the abaxial leaf surface measured using a drop shape analyzer identified highly significant colocalizing QTL on chromosome 7K (Chr07K). Using a multipronged approach, we identified a cluster of genes including Pavir.7KG077009, which encodes a Type III polyketide synthase-like protein, and Pavir.7KG013754 and Pavir.7KG030500, two highly similar genes that encode putative acyl-acyl carrier protein (ACP) thioesterases, as strong candidates underlying the QTL. The lack of homoeologs for any of the three genes on Chr07N, the relatively low level of identity with other switchgrass KCS proteins and thioesterases, as well as the organization of the surrounding region suggest that Pavir.7KG077009 and Pavir.7KG013754/Pavir.7KG030500 were duplicated into a fast-evolving chromosome region, which led to their neofunctionalization. Furthermore, sequence analyses showed all three genes to be absent in the two upland compared to the two lowland accessions analyzed. This study provides an example of and practical guide for trait mapping and candidate gene identification in a complex genetic system by combining QTL mapping, transcriptomics and variant analysis.

Integration of high-density genetic mapping with transcriptome analysis uncovers numerous agronomic QTL and reveals candidate genes for the control of tillering in sorghum.

Phenotypes such as branching, photoperiod sensitivity, and height were modified during plant domestication and crop improvement. Here, we perform quantitative trait locus (QTL) mapping of these and other agronomic traits in a recombinant inbred line (RIL) population derived from an interspecific cross between Sorghum propinquum and Sorghum bicolor inbred Tx7000. Using low-coverage Illumina sequencing and a bin-mapping approach, we generated ∼1920 bin markers spanning ∼875 cM. Phenotyping data were collected and analyzed from two field locations and one greenhouse experiment for six agronomic traits, thereby identifying a total of 30 QTL. Many of these QTL were penetrant across environments and co-mapped with major QTL identified in other studies. Other QTL uncovered new genomic regions associated with these traits, and some of these were environment-specific in their action. To further dissect the genetic underpinnings of tillering, we complemented QTL analysis with transcriptomics, identifying 6189 genes that were differentially expressed during tiller bud elongation. We identified genes such as Dormancy Associated Protein 1 (DRM1) in addition to various transcription factors that are differentially expressed in comparisons of dormant to elongating tiller buds and lie within tillering QTL, suggesting that these genes are key regulators of tiller elongation in sorghum. Our study demonstrates the usefulness of this RIL population in detecting domestication and improvement-associated genes in sorghum, thus providing a valuable resource for genetic investigation and improvement to the sorghum community.

A rapid thioacidolysis method for biomass lignin composition and tricin analysis.

BACKGROUND: Biomass composition varies from plant to plant and greatly affects biomass utilization. Lignin is a heterogeneous phenolic polymer derived mainly from p-coumaryl, coniferyl, and sinapyl alcohols and makes up to 10-25% of lignocellulosic biomass. Recently, tricin, an O-methylated flavone, was identified as a lignin monomer in many grass species. Tricin may function as a nucleation site for lignification and is advocated as a novel target for lignin engineering to reduce lignin content and improve biomass digestibility in grasses. Thioacidolysis is an analytical method that can be adapted to analyze both lignin monomeric composition and tricin content in the lignin polymer. However, the original thioacidolysis procedure is complex, laborious, and time consuming, making it difficult to be adopted for large-scale screening in biomass research. In this study, a modified, rapid higher throughput thioacidolysis method was developed. RESULTS: In combination with gas chromatography-mass spectrometry (GC-MS) and liquid chromatography-mass spectrometry (LC-MS), the modified thioacidolysis method can be used to simultaneously characterize the lignin composition and tricin content using 2-5 mg of dry samples. The modified method eliminates the solvent extraction and drastically improves the throughput; 80 samples can be processed in one day per person. Our results indicate that there is no significant difference in the determination of lignin S/G ratio and tricin content between the original and modified methods. CONCLUSIONS: A modified thioacidolysis protocol was established. The results demonstrate that the modified method can be used for rapid, high-throughput, and reliable lignin composition and tricin content analyses for screening transgenic plants for cell wall modifications or in large-scale genome-wide association studies (GWAS).

The Halophyte Seashore Paspalum Uses Adaxial Leaf Papillae for Sodium Sequestration.

Salinity is a growing issue worldwide, with nearly 30% of arable land predicted to be lost due to soil salinity in the next 30 years. Many grass crops that are vital to sustain the world's caloric intake are salt sensitive. Studying mechanisms of salt tolerance in halophytic grasses, plants that thrive in salt conditions, may be an effective approach to ultimately improve salt-sensitive grass crops. Seashore paspalum (Paspalum vaginatum) is a halophytic Panicoid grass able to grow in salt concentrations near that of seawater. Despite its widespread cultivation as a sustainable turfgrass, the mechanism underlying its ability to retain high Na+ concentrations in photosynthetic tissue while maintaining growth remains unknown. We examined the leaf structure and ion content in P. vaginatum 'HI10', which shows increased growth under saline conditions, and Paspalum distichum 'Spence', which shows reduced growth under salt, to better understand the superior salt tolerance of cv HI10. A striking difference between cv HI10 and cv Spence was the high steady-state level of K+ in cv HI10. Imaging further showed that the adaxial surface of both cv HI10 and cv Spence contained dense costal ridges of papillae. However, these unicellular extensions of the epidermis were significantly larger in cv HI10 than in cv Spence. The cv HI10 papillae were shown to act as Na+ sinks when plants were grown under saline conditions. We provide evidence that leaf papillae function as specialized structures for Na+ sequestration in P. vaginatum, illustrating a possible path for biotechnological improvement of salt-sensitive Panicoid crops with analogous leaf structures.

Italian weedy rice-A case of de-domestication?

Weedy rice is a representative of the extensive group of feral weeds that derive from crops, but has returned to the lifestyle of a wild species. These weeds develop either from a hybridization of crops with wild relatives (exoferality), or by mutation of crops to weedy forms (endoferality). Due to the close relation of weed and crop, the methods for weed-targeted containment are limited to date. A deeper understanding of the development of such weeds might help to design more efficient and sustainable approaches for weed management. Weedy rice poses a serious threat to rice yields worldwide. It is widely accepted that weedy rice has originated independently in different regions all over the world. However, details of its evolution have remained elusive. In the current study, we investigated the history of weedy rice in northern Italy, the most important rice-growing area in Europe. Our approach was to analyze genes related to weedy traits (SD1, sh4, Rc) in weedy rice accessions compared to cultivars, and to integrate these results with phenotypic and physiological data, as well as historical information about rice farming in Italy. We arrive at a working model for the timeline of evolution of weedy rice in Italy indicating that both exoferality and endoferality acted as forces driving the development of the diverse weedy rice populations found in the region today. Models of weed evolution can help to predict the direction which weed development might take and to develop new, sustainable methods to control feral weeds.

Strategies for eQTL mapping in allopolyploid organisms.

KEY MESSAGE: This study uses simulations to explore statistical power and false-positive rates for eQTL mapping in allopolyploid organisms and provides guidelines to apply eQTL mapping in these organisms. In recent years, RNA-seq has become the dominant technology for eQTL studies. However, most work has been in diploid organisms. Many species of economic and environmental importance are polyploid, and approaches for eQTL mapping in polyploids are not well developed. High similarity between duplicated genes in polyploids will cause misassignment of sequence reads and may cause false-positive results and/or lack of power to detect eQTL. In this paper, we first explore the similarity of homoeologous transcripts in polyploid organisms. We find that 5-20% of genes (varying with organism) in important agricultural plants such as wheat, soybean, and switchgrass are not sufficiently diverged between duplicated genomes to allow unambiguous assignment of reads. Second, we examine the impact of misassigned reads on eQTL mapping and show that both false-positive and false-negative rates can be greatly inflated. Third, we compare four strategies for dealing with ambiguous reads: (1) dividing ambiguous reads evenly between homoeologous transcripts, (2) assigning them proportionally, (3) using all reads for all genes, and (4) discarding ambiguous reads. We find that the strategy of discarding ambiguous reads gives the best balance of false-positive and false-negative rates for most genes. However, for genes that are very similar between genomes, using all reads is the only choice. This leads to reduced power, but false-positive rates will be maintained. We also discuss QTL mapping in polyploids using allele-specific expression (ASE) and show how the proportion of ASE-informative reads varies according to the divergence between homoeologous genes.

Intra- and interspecific diversity analyses in the genus Eremurus in Iran using genotyping-by-sequencing reveal geographic population structure.

Eremurus species, better known as 'Foxtail Lily' or 'Desert Candle', are important worldwide in landscaping and the cut-flower industry. One of the centers of highest diversity of the genus Eremurus is Iran, which has seven species. However, little is known about the genetic diversity within the genus Eremurus. With the advent of genotyping-by-sequencing (GBS), it is possible to develop and employ single nucleotide polymorphism (SNP) markers in a cost-efficient manner in any species, regardless of its ploidy level, genome size or availability of a reference genome. Population structure and phylogeographic analyses of the genus Eremurus in Iran using a minimum of 3002 SNP markers identified either at the genus level or at the species level from GBS data showed longitudinal geographic structuring at the country scale for the genus and for the species E. spectabilis and E. luteus, and at the regional scale for E. olgae. Our analyses furthermore showed a close genetic relatedness between E. olgae and E. stenophyllus to the extent that they should be considered subspecies within an E. olgae/stenophyllus species complex. Their close genetic relatedness may explain why crosses between these two (sub)species have been found in the wild and are exploited extensively as ornamentals. Last, current species identification, while robust, relies on flower morphology. A subset of seven SNPs with species-specific (private) alleles were selected that differentiate the seven Eremurus species. The markers will be especially useful for cultivar protection and in hybrid production, where true hybrids could be identified at the seedling stage.

High Density Genetic Maps of Seashore Paspalum Using Genotyping-By-Sequencing and Their Relationship to The Sorghum Bicolor Genome.

As a step towards trait mapping in the halophyte seashore paspalum (Paspalum vaginatum Sw.), we developed an F1 mapping population from a cross between two genetically diverse and heterozygous accessions, 509022 and HI33. Progeny were genotyped using a genotyping-by-sequencing (GBS) approach and sequence reads were analyzed for single nucleotide polymorphisms (SNPs) using the UGbS-Flex pipeline. More markers were identified that segregated in the maternal parent (HA maps) compared to the paternal parent (AH maps), suggesting that 509022 had overall higher levels of heterozygosity than HI33. We also generated maps that consisted of markers that were heterozygous in both parents (HH maps). The AH, HA and HH maps each comprised more than 1000 markers. Markers formed 10 linkage groups, corresponding to the ten seashore paspalum chromosomes. Comparative analyses showed that each seashore paspalum chromosome was syntenic to and highly colinear with a single sorghum chromosome. Four inversions were identified, two of which were sorghum-specific while the other two were likely specific to seashore paspalum. These high-density maps are the first available genetic maps for seashore paspalum. The maps will provide a valuable tool for plant breeders and others in the Paspalum community to identify traits of interest, including salt tolerance.

Structural Characterization of ABCB1, the Gene Underlying the d2 Dwarf Phenotype in Pearl Millet, Cenchrus Americanus (L.) Morrone.

Pearl millet is an important food crop in arid and semi-arid regions of South Asia and sub-Saharan Africa and is grown in Australia and the United States as a summer fodder crop. The d2 dwarf germplasm has been widely used in the last half-century to develop high-performing pearl millet hybrids. We previously mapped the d2 phenotype to a 1.6 cM region in linkage group (LG) 4 and identified the ABCB1 gene as a candidate underlying the trait. Here, we report the sequence, structure and expression of ABCB1 in tall (D2D2) and d2 dwarf (d2d2) germplasm. The ABCB1 allele in d2 dwarfs differs from that in tall inbreds by the presence of two different high copy transposable elements, one in the coding region and the second located 664 bp upstream of the ATG start codon. These transposons were present in all d2 dwarfs tested that were reported to be of independent origin and absent in the analyzed wild-type tall germplasm. We also compared the expression profile of this gene in different organs of multiple tall and d2 dwarf inbreds, including the near-isogenic inbreds at the d2 locus, Tift 23B (D2D2) and Tift 23DB (d2d2). Heterologous transformation of the tall (Ca_ABCB1) and the d2 dwarf (Ca_abcb1) pearl millet alleles in the Arabidopsis double mutant abcb1abcb19 showed that the pearl millet D2 but not the d2 allele complements the Arabidopsis abcb1 mutation. Our studies also show the importance of the COOH-terminal 22 amino acids of the ABCB1 protein in either protein function or stability.