Search

The low phytic acid (lpa) trait in soybeans can be conferred by loss-of-function mutations in genes encoding myo-inositol phosphate synthase and two epistatically interacting genes encoding multidrug-resistance protein ATP-binding cassette (ABC) transporters. However, perturbations in phytic acid biosynthesis are associated with poor seed vigor. Since the benefits of the lpa trait, in terms of end-use quality and sustainability, far outweigh the negatives associated with poor seed performance, a fuller understanding of the molecular basis behind the negatives will assist crop breeders and engineers in producing variates with lpa and better germination rate. The gene regulatory network (GRN) for developing low and normal phytic acid soybean seeds was previously constructed, with genes modulating a variety of processes pertinent to phytic acid metabolism and seed viability being identified. In this study, a comparative time series analysis of low and normal phytic acid soybeans was carried out to investigate the transcriptional regulatory elements governing the transitional dynamics from dry seed to germinated seed. GRNs were reverse engineered from time series transcriptomic data of three distinct genotypic subsets composed of lpa soybean lines and their normal phytic acid sibling lines. Using a robust unsupervised network inference scheme, putative regulatory interactions were inferred for each subset of genotypes. These interactions were further validated by published regulatory interactions found in Arabidopsis thaliana and motif sequence analysis. Results indicate that lpa seeds have increased sensitivity to stress, which could be due to changes in phytic acid levels, disrupted inositol phosphate signaling, disrupted phosphate ion (Pi) homeostasis, and altered myo-inositol metabolism. Putative regulatory interactions were identified for the latter two processes. Changes in abscisic acid (ABA) signaling candidate transcription factors (TFs) putatively regulating genes in this process were identified as well. Analysis of the GRNs reveal altered regulation in processes that may be affecting the germination of lpa soybean seeds. Therefore, this work contributes to the ongoing effort to elucidate molecular mechanisms underlying altered seed viability, germination and field emergence of lpa crops, understanding of which is necessary in order to mitigate these problems. Virginia Soybean Board projects [630, 719, 783]; Agricultural Experiment Station Hatch Program at Virginia Tech; Open Access Subvention Fund at Virginia Tech; John Lee Pratt Fellowship Program at Virginia Tech Published version This project was funded in part by the Virginia Soybean Board projects 630, 719, and 783. Additional funding was provided by the Agricultural Experiment Station Hatch Program and Open Access Subvention Fund - both at Virginia Tech. LD was funded in part by the John Lee Pratt Fellowship Program at Virginia Tech.

In this study, four recombinant inbred line (RIL) soybean populations were screened for their response to infection by Pythium sylvaticum, Pythium irregulare, Pythium oopapillum, and Pythium torulosum. The parents, PI 424237A, PI 424237B, PI 408097, and PI 408029, had higher levels of resistance to these species in a preliminary screening and were crossed with “Williams,” a susceptible cultivar. A modified seed rot assay was used to evaluate RIL populations for their response to specific Pythium species selected for a particular population based on preliminary screenings. Over 2500 single-nucleotide polymorphism (SNP) markers were used to construct chromosomal maps to identify regions associated with resistance to Pythium species. Several minor and large effect quantitative disease resistance loci (QDRL) were identified including one large effect QDRL on chromosome 8 in the population of PI 408097 × Williams. It was identified by two different disease reaction traits in P. sylvaticum, P. irregulare, and P. torulosum. Another large effect QDRL was identified on chromosome 6 in the population of PI 408029 × Williams, and conferred resistance to P. sylvaticum and P. irregulare. These large effect QDRL will contribute toward the development of improved soybean cultivars with higher levels of resistance to these common soil-borne pathogens. Published version

Phytophthora sojae causes Phytophthora root and stem rot of soybean and has been primarily managed through deployment of qualitative Resistance to P. sojae genes (Rps genes). The effectiveness of each individual or combination of Rps gene(s) depends on the diversity and pathotypes of the P. sojae populations present. Due to the complex nature of P. sojae populations, identification of more novel Rps genes is needed. In this study, phenotypic data from previous studies of 16 panels of plant introductions (PIs) were analyzed. Panels 1 and 2 consisted of 448 Glycine max and 520 G. soja, which had been evaluated for Rps gene response with a combination of P. sojae isolates. Panels 3 and 4 consisted of 429 and 460 G. max PIs, respectively, which had been evaluated using individual P. sojae isolates with complex virulence pathotypes. Finally, Panels 5–16 (376 G. max PIs) consisted of data deposited in the USDA Soybean Germplasm Collection from evaluations with 12 races of P. sojae. Using these panels, genome-wide association (GWA) analyses were carried out by combining phenotypic and SoySNP50K genotypic data. GWA models identified two, two, six, and seven novel Rps loci with Panels 1, 2, 3, and 4, respectively. A total of 58 novel Rps loci were identified using Panels 5–16. Genetic and phenotypic dissection of these loci may lead to the characterization of novel Rps genes that can be effectively deployed in new soybean cultivars against diverse P. sojae populations. Published (Publication status)

Two low-phytate soybean (Glycine max(L.) Merr.) mutant lines- V99-5089 (mipsmutation on chromosome 11) and CX-1834 (mrp-landmrp-nmutations on chromosomes 19 and 3, respectively) have proven to be valuable resources for breeding of low-phytate, high-sucrose, and low-raffinosaccharide soybeans, traits that are highly desirable from a nutritional and environmental standpoint. A recombinant inbred population derived from the cross CX1834 x V99-5089 provides an opportunity to study the effect of different combinations of these three mutations on soybean phytate and oligosaccharides levels. Of the 173 recombinant inbred lines tested, 163 lines were homozygous for various combinations of MIPS and two MRP loci alleles. These individuals were grouped into eight genotypic classes based on the combination of SNP alleles at the three mutant loci. The two genotypic classes that were homozygousmrp-l/mrp-nand either homozygous wild-type or mutant at themipslocus (MIPS/mrp-l/mrp-normips/mrp-l/mrp-n) displayed relatively similar similar to 55% reductions in seed phytate, 6.94 mg g(-1)and 6.70 mg g(-1)respectively, as compared with 15.2 mg g(-1)in the wild-type MIPS/MRP-L/MRP-N seed. Therefore, in the presence of the double mutantmrp-l/mrp-n, themipsmutation did not cause a substantially greater decrease in seed phytate level. However, the nutritionally-desirable high-sucrose/low-stachyose/low-raffinose seed phenotype originally observed in soybeans homozygous for themipsallele was reversed in the presence ofmrp-l/mrp-nmutations: homozygousmips/mrp-l/mrp-nseed displayed low-sucrose (7.70%), high-stachyose (4.18%), and the highest observed raffinose (0.94%) contents per gram of dry seed. Perhaps the block in phytic acid transport from its cytoplasmic synthesis site to its storage site, conditioned bymrp-l/mrp-n, alters myo-inositol flux inmipsseeds in a way that restores to wild-type levels themipsconditioned reductions in raffinosaccharides. Overall this study determined the combinatorial effects of three low phytic acid causing mutations on regulation of seed phytate and oligosaccharides in soybean. United Soybean Board; Virginia Soybean Board; Agricultural Experiment Station Hatch Program at Virginia Tech; Open Access Subvention Fund at Virginia Tech This project was funded in part by the United Soybean Board and Virginia Soybean Board. Additional funding was provided by the Agricultural Experiment Station Hatch Program and Open Access Subvention Fund - both at Virginia Tech. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

A dominant loss of function mutation in myo-inositol phosphate synthase (MIPS) gene and recessive loss of function mutations in two multidrug resistant protein type-ABC transporter genes not only reduce the seed phytic acid levels in soybean, but also affect the pathways associated with seed development, ultimately resulting in low emergence. To understand the regulatory mechanisms and identify key genes that intervene in the seed development process in low phytic acid crops, we performed computational inference of gene regulatory networks in low and normal phytic acid soybeans using a time course transcriptomic data and multiple network inference algorithms. We identified a set of putative candidate transcription factors and their regulatory interactions with genes that have functions in myo-inositol biosynthesis, auxin-ABA signaling, and seed dormancy. We evaluated the performance of our unsupervised network inference method by comparing the predicted regulatory network with published regulatory interactions in Arabidopsis. Some contrasting regulatory interactions were observed in low phytic acid mutants compared to non-mutant lines. These findings provide important hypotheses on expression regulation of myo-inositol metabolism and phytohormone signaling in developing low phytic acid soybeans. The computational pipeline used for unsupervised network learning in this study is provided as open source software and is freely available at https://lilabatvt.github.io/LPANetwork/. Published version

Meeting the food, feed, and fiber needs of a growing world population represents one of the signature challenges of this century. The UN FAO estimates that food production alone must increase by 70 percent by 2050 to meet the needs of a projected nine billion people. Demand will also increase for livestock feed, biofuel feedstocks, fiber for paper products, and construction materials. Meeting these demands will require implementation of advanced technologies, sustainable management of natural resources, and coordination of political forces. We must think beyond the boundaries of traditional agriculture, to integrate breakthroughs in plant science, engineering, environmental sciences, and other disciplines...

Rsv3 is one of three genetic loci conferring strain-specific resistance to Soybean mosaic virus (SMV). The Rsv3 locus has been mapped to a 154-kb region on chromosome 14, containing a cluster of five nucleotide-binding leucine-rich repeat (NB-LRR) resistance genes. High sequence similarity between the Rsv3 candidate genes challenges fine mapping of the locus. Among the five, Glyma14g38533 showed the highest transcript abundance in 1 to 3 h of SMV-G7 inoculation. Comparative sequence analyses were conducted with the five Rsv3 candidate NB-LRR genes from susceptible (rsv-type) soybean [Glycine max (L.) Merr.] cultivar Williams 82, resistant (Rsv3-type) cultivar Hwangkeum, and resistant lines L29 and RRR. Sequence comparisons revealed that Glyma14g38533 had far more polymorphisms than the other candidate genes. Interestingly, Glyma14g38533 gene from Rsv3-type lines exhibited 150 single-nucleotide polymorphism (SNP and six insertion–deletion (InDel) markers relative to rsv-type line, Furthermore, the polymorphisms identified in three Rsv3-type lines were highly conserved. Several polymorphisms were validated in 18 Rsv3-type resistant and six rsv-type susceptible lines and were found associated with their disease response. The majority of the polymorphisms were located in LRR domain encoding region, which is involved in pathogen recognition via protein–protein interactions. These findings associating Glyma14g38533 with Rsv3-type resistance to SMV suggest it is the most likely candidate gene for Rsv3.

Background Low phytic acid (lpa) crops are potentially eco-friendly alternative to conventional normal phytic acid (PA) crops, improving mineral bioavailability in monogastric animals as well as decreasing phosphate pollution. The lpa crops developed to date carry mutations that are directly or indirectly associated with PA biosynthesis and accumulation during seed development. These lpa crops typically exhibit altered carbohydrate profiles, increased free phosphate, and lower seedling emergence, the latter of which reduces overall crop yield, hence limiting their large-scale cultivation. Improving lpa crop yield requires an understanding of the downstream effects of the lpa genotype on seed development. Towards that end, we present a comprehensive comparison of gene-expression profiles between lpa and normal PA soybean lines (Glycine max) at five stages of seed development using RNA-Seq approaches. The lpa line used in this study carries single point mutations in a myo-inositol phosphate synthase gene along with two multidrug-resistance protein ABC transporter genes. Results RNA sequencing data of lpa and normal PA soybean lines from five seed-developmental stages (total of 30 libraries) were used for differential expression and functional enrichment analyses. A total of 4235 differentially expressed genes, including 512-transcription factor genes were identified. Eighteen biological processes such as apoptosis, glucan metabolism, cellular transport, photosynthesis and 9 transcription factor families including WRKY, CAMTA3 and SNF2 were enriched during seed development. Genes associated with apoptosis, glucan metabolism, and cellular transport showed enhanced expression in early stages of lpa seed development, while those associated with photosynthesis showed decreased expression in late developmental stages. The results suggest that lpa-causing mutations play a role in inducing and suppressing plant defense responses during early and late stages of seed development, respectively. Conclusions This study provides a global perspective of transcriptomal changes during soybean seed development in an lpa mutant. The mutants are characterized by earlier expression of genes associated with cell wall biosynthesis and a decrease in photosynthetic genes in late stages. The biological processes and transcription factors identified in this study are signatures of lpa-causing mutations. Published version

Background Root system architecture is important for water acquisition and nutrient acquisition for all crops. In soybean breeding programs, wild soybean alleles have been used successfully to enhance yield and seed composition traits, but have never been investigated to improve root system architecture. Therefore, in this study, high-density single-feature polymorphic markers and simple sequence repeats were used to map quantitative trait loci (QTLs) governing root system architecture in an inter-specific soybean mapping population developed from a cross between Glycine max and Glycine soja. Results Wild and cultivated soybean both contributed alleles towards significant additive large effect QTLs on chromosome 6 and 7 for a longer total root length and root distribution, respectively. Epistatic effect QTLs were also identified for taproot length, average diameter, and root distribution. These root traits will influence the water and nutrient uptake in soybean. Two cell division-related genes (D type cyclin and auxin efflux carrier protein) with insertion/deletion variations might contribute to the shorter root phenotypes observed in G. soja compared with cultivated soybean. Based on the location of the QTLs and sequence information from a second G. soja accession, three genes (slow anion channel associated 1 like, Auxin responsive NEDD8-activating complex and peroxidase), each with a non-synonymous single nucleotide polymorphism mutation were identified, which may also contribute to changes in root architecture in the cultivated soybean. In addition, Apoptosis inhibitor 5-like on chromosome 7 and slow anion channel associated 1-like on chromosome 15 had epistatic interactions for taproot length QTLs in soybean. Conclusion Rare alleles from a G. soja accession are expected to enhance our understanding of the genetic components involved in root architecture traits, and could be combined to improve root system and drought adaptation in soybean. Published version

Soybean is an integral oilseed crop that is used all over the world. Despite this, soybean has been experiencing a decrease in seed quality seen by the decreasing trend in test weight, and is limited in its uses due to the antinutritional factor, trypsin inhibitor (TI). This project first worked to identify single nucleotide polymorphisms (SNPs) associated with high test weight through a genome-wide association study (GWAS) as well as determine the relationship between important seed composition traits and test weight. Additionally, it focused on identifying and knocking out KTI genes specifically expressed in soybean seed tissue through CRISPR/Cas9. Lastly, this project investigated further use of soybean in aquaculture by evaluating the performance of Rainbow trout fed diets with low-TI, 'VT Barrack' soy meal. Nine SNPs on chromosome 15 were found to be significantly associated with high test weight as and eight potential gene candidates were identified. Test weight was found to be significantly, and negatively related to seed oil content and had inconsistent correlations with protein and sugar content. The KTI genes KTI1, Glyma01g095000, and KTI3, Glyma08g341500, were found to be only expressed in seed tissues and multiple KTI knock-out soybean plants with decreased seed TI content and activity were developed. Lastly, we found that our low-TI soymeal could replace 30% of fish meal in aquafeeds with no negative effects on trout growth and health. Doctor of Philosophy Soybean is an essential and lucrative crop due to its countless uses. Despite its importance and flexibility, it contains antinutritional factors, such as trypsin inhibitors (TI), and has been experiencing a decline in test weight, the amount of soybean seed in kilograms per hectoliter, that can limit the profits that growers receive. This dissertation seeks to contribute to the development of improved soybean varieties in addition to evaluating current low-TI varieties' application in aquaculture as a substitute for fish meal. A genome-wide association study (GWAS) was performed to identify molecular markers associated with high test weight and the relationship between high test weight and important seed composition traits was determined. Gene editing was conducted to knock-out expression of Kunitz trypsin inhibitor (KTI) genes specifically expressed in seed tissue. A fish feeding trial was conducted to evaluate the performance of Rainbow trout-fed diets using our low-TI soybean cultivar, 'VT Barrack'. We were able to successfully identify multiple molecular markers associated with high test weight as well as identify multiple candidate genes. We were able to determine that high test weight has a strong, negative relationship with seed oil content. We identified two KTI genes for gene knockout and were able to successfully develop a low-TI soybean line that has limited TI in just seed tissue. Lastly, we found that our low-TI soymeal could replace 30% of fish meal in aquafeeds with no negative effects on trout growth and health.

The field of functional genomics aims to understand the complex relationship between genotype and phenotype by integrating genome-wide approaches, such as transcriptomics, proteomics, and metabolomics. Large-scale "-omics" research has been made widely possible by the advent of high-throughput techniques, such as next-generation sequencing and mass-spectrometry. The vast data generated from such studies provide a wealth of information on the biological dynamics underlying phenotypes. Though functional genomics approaches are used extensively in human disease research, their use also spans organisms as miniscule as mycoplasmas to as great as sperm whales. In particular, functional genomics is instrumental in agricultural advancements for the improvement of productivity and sustainability in crop and livestock production. Improvement in soybean production is especially imperative, as soybeans are a primary source of oil and protein for human and livestock consumption, respectively. The research presented here employs functional genomics approaches – transcriptomics and metabolomics – to discern the transcriptional regulation and metabolic events underlying two economically important agronomic traits in soybean: seed phytic acid content and Soybean mosaic virus resistance. At normal levels, seed phytic acid content inhibits mineral absorption in humans and livestock, acting as an antinutrient and contributing to phosphorus pollution; however, the development of low phytic acid soybeans has helped mitigate these issues, as their seeds increase nutrient bioavailability and reduce environmental impact. Despite these desirable qualities, low phytic acid soybeans exhibit poor seed performance, which negatively affects germination rates and yield and has prevented their large-scale commercial production. Thus, part of the focus of this research was investigating the effects of mutations conferring the low phytic acid phenotype on seed germination. Comparative studies between low and normal phytic acid soybean seeds were carried out and revealed distinct differences in metabolite profiles and in the transcriptional regulation of biological pathways that may be vital for successful seed germination. The final part of this research concerns Rsv3-mediated extreme resistance, a unique mode of resistance that is effective against the most virulent strains of Soybean mosaic virus. The molecular mechanisms governing this type of resistance are poorly characterized. Therefore, the research presented here attempts to elucidate the regulatory elements responsible for the induction of the Rsv3-mediated extreme resistance response. Utilizing a comparative transcriptomic time series approach on Soybean mosaic virus-inoculated Rsv3 (resistant) and rsv3 (susceptible) soybean lines, this final study provides gene candidates putatively functioning in the regulation of biological pathways demonstrated to be crucial for Rsv3-mediated resistance. Doctor of Philosophy Soybeans are a crop of great economic importance, being a primary source of oil and protein for human and livestock consumption, respectively. Increasing demand for soybean calls for improvement in its production. An emerging field that has had tremendous impact on this endeavor is the field of functional genomics. Functional genomics approaches generate large-scale biological data that can aid in discerning how specific processes are regulated and controlled in an organism. The research presented in this work utilizes functional genomics approaches to elucidate the biological mechanisms underlying two economically important traits in soybean: seed phytic acid content and Soybean mosaic virus resistance. Phytic acid is a compound found in soybean seeds that causes nutrient deficiencies and phosphorus pollution. Soybeans with reduced to phytic acid content have been developed to mitigate these problems; they have poor seed germination and emergence. The studies in this work employ functional genomics approaches to compare unique sets of low and normal phytic acid soybeans to help establish the relationship between seed phytic acid content and seed performance. These studies resulted in new and promising hypotheses for future studies on investigating the low phytic acid trait. The final focus of this work used a functional genomics approach to discern the molecular mechanisms underlying a unique mode of resistance to Soybean mosaic virus. The study identified genes in soybean that are potentially critical to resistance against Soybean mosaic virus.

The PROTEOLYSIS 6 (PRT6) branch of the N-end rule pathway is a well-characterized negative regulator of flooding and low oxygen tolerance in plants. This study investigated the role of this pathway in adaptation to salinity stress in Arabidopsis and maize via physiological and molecular characterization of Arabidopsis prt6-1 and maize prt6 MU insertion mutants, respectively. Our study demonstrated that the loss of function mutation of prt6 in Arabidopsis activated hormonal and transcriptional responses associated with adaptation to salinity stress, enhancing high salt tolerance at seed germination, seedling, and adult plant stages. Our data also indicated that salinity tolerance conferred by the prt6 mutation is attributed to increased mRNA abundance of key transcriptional factors in ABA-dependent (AREB/ABFs) and independent (DREBs) pathways, together with the dominant expression of downstream dehydrins. Furthermore, this study revealed that the prt6 mutation enhances ethylene and brassinosteroid responses, resulting in restricted Na+ accumulation in roots and shoots as well as increased expression of dehydrin genes such as RD29A and RD29B. Maize prt6 mutant plants, contrary to our observation in Arabidopsis, showed lower seed germination, primary root elongation, and shoot biomass growth along with increased malondialdehyde (MDA) accumulation under high salt. Moreover, maize prt6 mutants exhibited reduced grain yield and yield-related components under high salt. These results indicate that PRT6 functions as a negative regulator for salinity tolerance in Arabidopsis, whereas this gene plays a positive role in salinity tolerance in maize. In wheat, we compared two genotypes with contrasting nitrogen-use-efficiency (NUE), VA08MAS-369 and VA07W-415, to dissect physiological and molecular mechanisms underlying NUE regulation. Our agronomic data revealed that line 369 maintained yield and yield-related parameters and exhibited greater NUE indexes relative to line 415 under N deficient conditions. Furthermore, our analyses suggested that the significantly higher nitrogen use efficiency (NUE) in line 369 could be attributed to the greater N uptake efficiency in this genotype. In fact, line 369 was able to maintain the development of root systems under N limitation. Consistently, genes encoding high-affinity nitrate transporters such as TaNRT2.1 and TaNRT2.2 were expressed more abundantly in the roots of line 369 than line 415 at limited N. Overall, the results of this study characterized physiological and molecular phenotypes associated with high N uptake efficiency in line 369. This is useful information for the development of new wheat accessions with improved NUE. Doctor of Philosophy In coastal areas, sea-level rise increases the chances of saltwater intrusion into cultivable lands, making a hostile environment for crop growth and production by imposing flooding and salinity stresses simultaneously. Identification of central regulators that regulate the adaptation to both flooding and salinity is a critical step for the development of new crop genotypes with enhanced tolerance to these stresses. Previous studies have characterized the function of the PROTEOLYSIS 6 (PRT6) gene in adaptation to flooding stress in plants. This study assessed whether this gene is involved in adaptation to salinity stress in Arabidopsis and maize by evaluating the growth and survival of their respective prt6 mutants under high salt. Consistent with the flooding tolerance data, our study showed that the PRT6 gene also functions as a negative regulator of salinity stress tolerance in Arabidopsis. The prt6 mutation in Arabidopsis activated the key transcriptional and hormone response pathways associated with adaptation to both salinity/osmotic stress and sodium toxicity, expressed as enhanced tolerance to excess salt at seed germination, seedling, and adult plant stages. In maize, disruption of the PRT6 gene decreased seed germination, primary root elongation, and shoot biomass growth under high salt, which is opposite to our observations in Arabidopsis. Additionally, the maize mutant plants encountered more oxidative stress, as demonstrated by the higher accumulation of malondialdehyde (MDA) under high salt. Moreover, maize prt6 mutants exhibited reduced grain yield under high salt. Overall, these results indicate that disruption of the PRT6 gene confers increased tolerance to high salt in Arabidopsis, whereas it conversely reduced salinity tolerance in maize. In wheat, we compared two genotypes with distinct nitrogen use efficiency (NUE), VA08MAS-369 and VA07W-415, to determine critical traits involved in NUE regulation. Our study showed that grain yield and yield-related parameters were significantly higher in line 369 than line 415 under low N. Moreover, high NUE in line 369 was attributed to efficient N uptake in this genotype under limited N. Our root architecture analysis demonstrated that line 369 was able to maintain root depth, volume, and thickness even under N limitation. Consistently, line 369 highly induced expression of genes associated with nitrogen transport at low N. Altogether, this study identified key traits involved in high NUE in wheat, facilitating the breeding of new wheat genotypes with enhanced NUE.

Quantitative genetics aims to bridge the genome to phenome gap. The advent of high-throughput genotyping technologies has accelerated the progress of genome to phenome mapping, but a challenge remains in phenotyping. Various high-throughput phenotyping (HTP) platforms have been developed recently to obtain economically important phenotypes in an automated fashion with less human labor and reduced costs. However, the effective way of designing HTP has not been investigated thoroughly. In addition, high-dimensional HTP data bring up a big challenge for statistical analysis by increasing computational demands. A new strategy for modeling high-dimensional HTP data and elucidating the interrelationships among these phenotypes are needed. Previous studies used pedigree-based connectetdness statistics to study the design of phenotyping. The availability of genetic markers provides a new opportunity to evaluate connectedness based on genomic data, which can serve as a means to design HTP. This dissertation first discusses the utility of connectedness spanning in three studies. In the first study, I introduced genomic connectedness and compared it with traditional pedigree-based connectedness. The relationship between genomic connectedness and prediction accuracy based on cross-validation was investigated in the second study. The third study introduced a user-friendly connectedness R package, which provides a suite of functions to evaluate the extent of connectedness. In the last study, I proposed a new statistical approach to model high-dimensional HTP data by leveraging the combination of confirmatory factor analysis and Bayesian network. Collectively, the results from the first three studies suggested the potential usefulness of applying genomic connectedness to design HTP. The statistical approach I introduced in the last study provides a new avenue to model high-dimensional HTP data holistically to further help us understand the interrelationships among phenotypes derived from HTP. Doctor of Philosophy Quantitative genetics aims to bridge the genome to phenome gap. With the advent of genotyping technologies, the genomic information of individuals can be included in a quantitative genetic model. A new challenge is to obtain sufficient and accurate phenotypes in an automated fashion with less human labor and reduced costs. The high-throughput phenotyping (HTP) technologies have emerged recently, opening a new opportunity to address this challenge. However, there is a paucity of research in phenotyping design and modeling high-dimensional HTP data. The main themes of this dissertation are 1) genomic connectedness that could potentially be used as a means to design a phenotyping experiment and 2) a novel statistical approach that aims to handle high-dimensional HTP data. In the first three studies, I first compared genomic connectedness with pedigree-based connectedness. This was followed by investigating the relationship between genomic connectedness and prediction accuracy derived from cross-validation. Additionally, I developed a connectedness R package that implements a variety of connectedness measures. The fourth study investigated a novel statistical approach by leveraging the combination of dimension reduction and graphical models to understand the interrelationships among high-dimensional HTP data.

Most staple food crops were domesticated thousands of years ago through independent processes across different regions of the world. Studies of the history of such crops have been essential to our understanding of plant domestication as a process that started with the collection of wild material and continued with subsequent propagation, cultivation, and selection under human care. Domestication often involves a complex genetic structure with contributions from multiple founder populations, interspecific hybridization, chromosomal introgressions, and polyploidization events that occurred hundreds to thousands of years earlier. Such intricate origins complicate the systematic study of the sources of phenotypic variation. The analysis of recently domesticated, non-traditional, non-model species, such as Sinningia speciosa (Gesneriaceae), can expand the knowledge that we have on phenotypic variation under domestication, and help us to comprehend modern patterns of plant domestication and to broaden our understanding of the general trends. S. speciosa is commonly known as the 'florist's gloxinia', and it has been cultivated for 200 years as an ornamental houseplant. In our genomic study of S. speciosa, we examined an extensive diversity panel consisting of 115 individuals that included different species in the genus, wild representatives, and cultivated accessions, as well as 150 individuals from an F2 segregating population. Our analyses revealed that all of the domesticated varieties are derived from a single founder population that originated in or near the city of Rio de Janeiro in Brazil. We identified two loci associated with domesticated traits (flower symmetry and color) and did not detect any major hybridization or polyploidization events that could have contributed to the rapid increase in phenotypic diversity. Our findings, in conjunction with other features such as a small, low-complexity genome, ease of cultivation, and rapid generation time, makes this species an attractive model for the study of genomic variation under domestication. Basic research on non-model organisms with low economic importance is uncommon but necessary to understand the world from a broader perspective. In such cases, reduced representation approaches like Genotyping-by-Sequencing (GBS) are efficient low-cost alternatives to whole genome resequencing. However, most of these technologies are subject to patent protection, licensing processes, and fees that constrain genomic research for small non-profit research organizations. We have designed a protocol to construct reduced representation libraries from genomic DNA. Our approach, called Targeted Amplification of Scattered Sites (TASS), deviates from the traditional digestion-ligation-amplification process that is the subject of intellectual property that protects most current methods. Instead, TASS relies on 1) targeting and duplicating scattered regions in the genome by annealing and expanding long tail primers with short annealing sites, and 2) amplifying these regions using primers that are complementary to the added overhangs. At the moment GBS is more consistent and delivers more variants than TASS. However, we have established a foundation on which further optimization can produce an accessible, easy to implement, high-throughput genotyping approach. Doctor of Philosophy Most staple food crops were domesticated thousands of years ago through unrelated processes that were initiated across different regions of the world. Studies of the history of such crops have been essential to our understanding of plant domestication, a process that started with the collection of wild material and continued with subsequent propagation and cultivation under human care. Plant domestication has often involved a complex combination of ancestral lineages that encompass multiple populations, crosses with other species, and large DNA reorganizations that occurred hundreds to thousands of years earlier. Such intricate origins make the systematic study of plant domestication very challenging. The analysis of recently domesticated plants such as the 'florist's gloxinia' (Sinningia speciosa), can help us to better understand some of the changes that have occurred during domestication, as well as to comprehend modern patterns of plant domestication and to broaden our understanding of general trends. Florist's gloxinias are ornamental plants that have been cultivated during the last 200 years. In this study we examined 115 specimens, including wild and cultivated types of florist's gloxinias, as well as closely related species in Sinningia. We also constructed and evaluated an artificial population of 150 individuals from the cross of a wild and a cultivated form. Our analyses revealed that all of the domesticated varieties are descendants from a single wild population that originated in or near the city of Rio de Janeiro in Brazil. We also identified two regions of DNA that are responsible for the changes in flower shape and color, and crosses with other species did not introduce such alterations. Our findings, in conjunction with other features such as its small nuclear DNA content, the ease of cultivation indoors, and a rapid generation time, makes the florists' gloxinia an attractive crop to the study the effects of plant domestication. Research on organisms with low economic importance is uncommon but necessary to understand the world from a broader perspective. In such cases, analyzing the entire genetic information that is stored as DNA may be cost-prohibitive. Instead, approaches that sample small portions of DNA from each individual can be utilized. Most of these technologies are currently patented and subject to licensing processes and fees that limit their implementation by small non-profit research organizations. In this study we designed a protocol to sample small portions of DNA, similarly to existing techniques. However, our approach, called Targeted Amplification of Scattered Sites (TASS), employs a sampling process that deviates from the traditional patented procedure that is used in most current methods. At present, TASS is not as consistent and delivers less information than traditional approaches. However, we have established a foundation on which further optimization can produce an accessible and easy to implement technique.

Investigations of potato (Solanum tuberosum) have been hampered by its complicated genetics and high genetic load. This dissertation applies genome reduction techniques to investigate a broad swath of genomic and physiological phenomena. It begins with the presentation and evaluation of a protocol to characterize endoreduplication within potato tubers, demonstrating substantial variation between tissue types and among wild species which may facilitate research into the genesis and growth of these starchy underground stems. Next, we transitioned to explore the distribution and consequences of a specific class of transposable element, Miniature Inverted Transposable Elements (MITEs), showing that they comprise approximately 5% of the potato genome, occur more frequently in genes with stress-related functions, and may be associated with changes, especially decreases, in gene expression. We then combined homology and sparsity based approaches to predict recent MITE activity, identifying five families as especially active. Finally, we expose the gene underlying the potato flower color locus, a homolog of AN2, while showing the effects it exerts on the flavonoid biosynthesis and fruit ripening pathways. This region was shown to be particularly dynamic, replete with MITEs and structural variants which we hypothesize to be the ultimate cause of differences in AN2 expression within the germplasm we examined. While the separate topics of this dissertation are quite disparate, each addresses an important topic in potato genetics, the in-depth study of which is only possible through the utilization of genomic reduction approaches to acquire homozygous genotypes for study and currently available genomic resources. Ph. D. Despite their humble appearance and routine consumption, potatoes have a complex genetic structure and a life cycle capable of both sexual reproduction through flowers, fruit and seed, and asexual reproduction through the tubers which also comprise the edible product. From an agronomic perspective, one of the most important qualities of a potato tuber is size, a feature influenced by genetics and environment. Cell-to-cell variation for the amount of DNA per cell, one component that influences tuber size, is known to occur, yet our ability to measure DNA content in starchy tuber cells has been obscured by debris generated through routine preparation techniques. We present and evaluate a new method for measuring the DNA content of potato tuber cells, which provides reliable results across a range of different potato varieties and species. ‘Jumping genes’ also known as transposons, first reported in maize but now known to occur in most advanced plant and animal species, have been found to comprise ~5% of the recently sequenced potato genome. We show that a particular class of transposons is more likely to occur adjacent or actually in certain types of genes, such as those which confer resistance to disease, where they may have meaningful effects on how those genes operate. We then proceed to predict the current activity of the various families of these jumping genes to understand how they continue to alter the genetic landscape of potato. Finally we identify a particular gene which dictates flower color in potato (purple vs. white). We demonstrate that several transposons occur in some forms of the flower color gene. Originally we hypothesized that transposons were associated with the turning off of the purple flower color form; however, on closer examination, we could express the white flower form in transgenic plants that were originally white-flowered and convert them to have purple flowers, demonstrating that even the white flower form was functional. While the separate topics of this dissertation are quite disparate, each addresses an important topic in potato genetics, the in-depth study of which is only possible through the availability of the special strains of potatoes with reduced chromosome number and the publication of the potato genome.

Abiotic stress is a major constraint for crop productivity worldwide. To better understand the common biological mechanisms of abiotic stress responses in plants, we performed meta-analysis of 652 samples of RNA sequencing (RNA-Seq) data from 43 published abiotic stress experiments in Arabidopsis thaliana. These samples were categorized into eight different abiotic stresses including drought, heat, cold, salt, light and wounding. We developed a multi-step computational pipeline, which performs data downloading, preprocessing, read mapping, read counting and differential expression analyses for RNA-Seq data. We found that 5729 and 5062 genes are induced or repressed by only one type of abiotic stresses. There are only 18 and 12 genes that are induced or repressed by all stresses. The commonly induced genes are related to gene expression regulation by stress hormone abscisic acid. The commonly repressed genes are related to reduced growth and chloroplast activities. We compared stress responsive genes between any two types of stresses and found that heat and cold regulate similar set of genes. We also found that high light affects different set of genes than blue light and red light. Interestingly, ABA regulated genes are different from those regulated by other stresses. Finally, we found that membrane related genes are repressed by ABA, heat, cold and wounding but are up regulated by blue light and red light. The results from this work will be used to further characterize the gene regulatory networks underlying stress responsive genes in plants. Master of Science

Natto and sprout soybeans are produced using small-seeded soybeans and their production is a high value alternative to grow grain soybeans for food in U.S. The development of soybean cultivars with improved natto and sprout quality is crucial for maintaining and increasing the soyfood market. However, there is insufficient information on sprout soybean characteristics. Therefore, the first objective of this study was to evaluate seed and sprout traits as potential selection criteria and study the storage effect on sprout quality. Seeds can be a vehicle for transmission of pathogens capable of causing human illness. That is why, the second objective was to identify seed-borne pathogens on a commercial soybean cultivar and to evaluate different seed decontamination treatments. Finally, seed coat deficiency is an undesirable trait for natto soybean seeds because it causes inferior appearance of the product. Thus, the third objective was to identify quantitative trait loci (QTL) underlying seed coat deficiency (SCD) and associated markers. Results showed that seed size, high-, average- and low-quality sprout percentage, hypocotyl thickness and length and sprout yield are the most important variables for breeding sprout cultivars; and one-year seed storage at room temperature reduced sprout quality. Fusarium, Alternaria and Diaphorte were the most frequent genera isolated from soybean seeds, and 2% calcium hypochlorite and 5% acetic acid were promising seed disinfection treatments. A stable QTL, qSCD20_1, was identified across two years explaining up to 25% of the variation of SCD; and eight molecular markers tightly linked and nearby qSCD20_1 were identified. Information presented will be helpful for sprout and natto soybean cultivar development. Master of Science

Disease resistance is critical in soft red winter wheat (Triticum aestivum L.) cultivars. Leaf rust caused by Puccinia triticina Eriks and stripe rust caused by Puccinia striiformis Westend. f.sp. tritici Eriks. are destructive pathogens of wheat. From 2014 to 2015 phenotypic data was collected at diverse locations for resistance to leaf rust (North Carolina, Texas, and Virginia) and stripe rust (Arkansas, North Carolina, Georgia, Texas, and Virginia) in a Pioneer ‘25R47’ /‘Jamestown’ (P47/JT) population composed of 186 F5:9 recombinant inbred lines (RILs). Analysis of the P47/JT population identified two quantitative trait loci (QTL) for leaf rust resistance on chromosome 5B and two QTL for stripe rust resistance on chromosomes 3B and 6A. Phenotypic variation (%) explained by the putative leaf rust resistance QTL of Jamestown on 5B was as high as 22.1%. Variation explained by the putative stripe rust resistance QTL of Jamestown on 3B and 6A was as high as 11.1 and 14.3%, respectively. Jamestown is postulated to contain gene Lr18. Seedlings of 186 F5:9 recombinant inbred lines from the P47/JT population and 200 F2 seedlings from eight other crosses including Jamestown and/or the Lr18 host differential line RL6009 (Thatcher*6/Africa 43) were screened with P. triticina race TNRJJ. Genetic analysis of the populations was conducted to validate the presence of Lr18 in Jamestown. Results of linkage analysis identified SNP maker IWB41960 linked within 5 cM of gene Lr18 in all three populations. From 2016 to 2017 phenotypic data was collected at diverse locations for resistance to leaf rust (Illinois, North Carolina, and Virginia) in a ‘2013412’ (PI 667644) / VA10W-21 (PI 676295) population (412/21) composed of 157 doubled haploid (DH) lines. The 412/21 DH lines were genotyped via genotyping by sequence (GBS). Analysis of the 412/21 population identified one quantitative trait loci (QTL) region associated with adult plant resistance to leaf rust on chromosome 1B. Phenotypic variation (%) explained by the putative leaf rust resistance QTL of 2013412 on 1B was as high as 40.1%. Kompetitive allele-specific (KASP) markers KASP_S1B_8414614 and KASP_S1B_8566239 were developed as markers for use in marker assisted selection. Ph. D.

Estimated at approximately $43 billion annually, the cultivated soybean Glycine max (L.) Merr., is the second most valuable crop in the United States. Soybeans account for 57% of the world oil-seed production and are utilized as a protein source in products such as animal feed. The value of a soybean crop, measured in seed quality and quantity, is negatively affected by biotic and abiotic stresses. This research is focused on resistance to biotic disease stress in soybean. In particular, we are working on the Phytophthora soja (P. sojae) and Soybean Mosaic Virus (SMV) systems. For each of these diseases, we are working to develop superior soybean germplasm that is resistant to the devastating economic impacts of pathogens. The majority of this research is focused on screening for novel sources of P. sojae resistance with core effectors to identify resistance genes (R-genes) that will be durable under field conditions. Four segregating populations and two recombinant inbred line (RIL) populations have been screened with core effectors. Effector-based screening methods were combined with pathogen-based phenotyping in the form of a mycelium-based trifoliate screening assay. One RIL population has been screened with virulent P. sojae mycelium. Disease phenotyping has generated a preliminary genetic map for resistance in soybean accession PI408132. The identification of novel R-genes will allow for stacking of resistance loci into elite G. max cultivars. The second project covered in this dissertation describes the validation of the SMV resistance gene Rsv3. Utilizing a combination of transient expression and homology modeling; we provide evidence that Glyma14g38533 encodes Rsv3. Ph. D.

In multiple species, genome-wide association (GWA) studies have become an increasingly prevalent method of identifying the quantitative trait loci (QTLs) that underlie complex traits. Despite this, relatively few GWA analyses using high-density genomic markers have been carried out on highly quantitative traits in wheat. We utilized single-nucleotide polymorphism (SNP) data generated via a genotyping-by-sequencing (GBS) protocol to perform GWA on multiple yield-related traits using a panel of 329 soft red winter wheat genotypes grown in four environments. In addition, the SNP data was used to examine linkage disequilibrium and population structure within the testing panel. The results indicated that an alien translocation from the species Triticum timopheevii was responsible for the majority of observed population structure. In addition, a total of 50 significant marker-trait associations were identified. However, a subsequent study cast some doubt upon the reproducibility and reliability of plant QTLs identified via GWA analyses. We used two highly-related panels of different genotypes grown in different sets of environments to attempt to identify highly stable QTLs. No QTLs were shared across panels for any trait, suggesting that QTL-by-environment and QTL-by-genetic background interaction effects are significant, even when testing across many environments. In light of the challenges involved in QTL mapping, prediction of phenotypes using whole-genome marker data is an attractive alternative. However, many evaluations of genomic prediction in crop species have utilized univariate models adapted from animal breeding. These models cannot directly account for genotype-by-environment interaction, and hence are often not suitable for use with lower-heritability traits assessed in multiple environments. We sought to test genomic prediction models capable of more ad-hoc analyses, utilizing highly unbalanced experimental designs consisting of individuals with varying degrees of relatedness. The results suggest that these designs can successfully be used to generate reasonably accurate phenotypic predictions. In addition, multivariate models can dramatically increase predictive accuracy for some traits, though this depends upon the quantity and characteristics of genotype-by-environment interaction. Ph. D.

Soybean and wheat are important crop species due to their significance for human consumption, animal feed, and industrial use. However, increasing global population and worsening climate change have put a major strain on the production system of these crops. Natural disasters such as flooding and drought can severely impact growth and productivity of these crops. In addition, increased application of synthetic nitrogenous fertilizers to meet the global food demand has led to environment related issues. Therefore, with a goal of understanding mechanisms of flooding and drought tolerance in soybean and nitrogen-use-efficiency in wheat, we explored their physiological and transcriptomic regulation. We characterized the fundamental acclimation responses of soybean to flooding and drought and compared the metabolic and transcriptomic regulation during the stresses in a tissue-specific manner. We demonstrated the dynamic reconfiguration of gene expression and metabolism during flooding, drought, and recovery from these stresses. Our study displayed that flooding triggers more dramatic adjustments than drought at the transcriptional level. We also identified that the soybean genome encodes nine members of group VII ERF genes and characterized their responses in leaves and roots under flooding and drought. Based on the expression patterns, it is estimated that two of the nine genes are promising candidate genes regulating tolerance to submergence and drought. In addition, our genome-scale expression analysis discovered commonly induced ERFs and MAPKs across both stresses (flooding and drought) and tissues (leaves and roots), which might play key roles in soybean survival of flooding and drought. In wheat, we evaluated the effect of three different nitrogen rates on yield and its components across four diverse soft red winter wheat genotypes. The cultivar Sisson displayed superior performance in grain yield and nitrogen use efficiency at low nitrogen levels. Our results suggested that improvement of nitrogen use efficiency in low nitrogen environments can be achieved through the selection of three components: grain number/spike, 1000-seed weight, and harvest index. Overall, this study has advanced our understanding of how plants respond to abiotic stresses such as flooding, drought, and nutrient limitation conditions. Ph. D.

Forest tree breeding methods rank among the most efficient ways to increase productivity and quality of forests. With the advent of high-throughput genotyping technology, genome-enabled breeding has started to gain importance and may overcome some weaknesses of traditional tree breeding. Genomic Selection (GS), which involves using genome-wide markers to predict breeding values of individuals in a population, has been proposed for animal and plant breeding programs. GS enables very accurate selection decisions through estimation of genomic estimated breeding values (GEBVs). While the goal of GS is to predict phenotype from genotype, it does not identify the underlying genes that have important roles in a trait. Genome-Wide Association Studies (GWAS) approaches are therefore complementary to GS, enabling identification of these genes, which may be useful for marker-assisted selection in some traits. In this study, we first estimated heritability for several adaptive traits (cold hardiness, dbh, bud flush, height, and bud set) in a population of Populus trichocarpa and for height, diameter, and stem straightness in Pinus taeda. GEBVs accuracies were estimated using a ridge regression–best linear unbiased prediction (rrBLUP) model, and these accuracies were compared with estimated heritabilities. GWAS was also performed for the both imputed and non–imputed data of P. taeda population using TASSEL (Trait Analysis by aSSociation Evolution and Linkage) software, as well as rrBLUP and FFBSKAT (Fast Family-Based Sequence Kernel Association Test) packages in R. Heritabilities ranged from 0.34 to 0.56 for P. trichocarpa and 0.14 to 0.37 for P.taeda. GWAS identified 3244 associations for dbh, 4077 associations for stem straightness, and 5280 SNPs for height (p≤0.05) in TASSEL using the reduced model (marker data only), whereas 2729, 3272 and 3531 associations were found with the full model where we also included population structure as a covariate. FFBSKAT showed a similar number of SNP associations (2989, 3046 and 3058). There was an inflation of SNP associations (~20k) found in rrBLUP, which suggests population structure was not effectively controlled. The GEBVs accuracies ranged from 0.09 and 0.22 for P.trichocarpa and 0.09 to 0.23 for P.taeda using rrBLUP method. Testing the effect of repetation on the accuracy of GEBV for poplar showed that there was no significant difference between the number of cycles. Also, there was no significant difference the accuracy of GEBVs in pine between two different imputation methods, the marker mean value and Beagle software. Master of Science

Temperate tree species exhibit seasonal growth cycling, and the timing of such transition varies with local climate. Under anthropogenic climate change, the local pattern of growth and dormancy in tree populations is expected to become uncoupled with shifting seasonal environmental signals, particularly temperature. Thus, an understanding of the genetic underpinnings of local adaptation is key to predicting the fate of tree populations in the future. In this thesis, we coupled sampling of range-wide natural accessions of P. trichocarpa with adaptive trait phenotyping and genome-wide genotyping to uncover relationships between genotype, phenotype, and environment. We detected strong correlations between adaptive phenotypes, climate, and geography, which suggested climatic selection driving adaptation of these populations to local environments. We subsequently combined genotype-phenotype association tests with sliding window analysis and identified regions strongly associated with these adaptive traits. We also compared adaptive markers identified in two independent GWAS on samples across latitude and altitude transects and found a set of associated variants shared across both transects. We further scanned the genome with three selection tests to identify regions showing evidence of recent positive and divergent selection. By comparing candidate selection regions across altitude and latitude, we detected a set of overlapping regions showing differentiation across gradients of the same climate variables. We validated the functional imortance of these selection regions by combining GWAS and showed that selection regions contain a strong signature of phenotypic associations. We also studied the distribution of deleterious allels across genome and natural populations, and found that deleterious alleles preferentially accumulate in regions of low recombination and hithihking regions. Finally, marginal populations contained more deleterious alleles compared with central populations, which is likely due to ineffective selection in small populations and recent bottlenecks associated with postglacial recolonization. These findings provide new insights into the genomic architecture underlying climatic adaptation and how selection drives adaptive evolution of tree species. Ph. D.

As the global population rises, the demand for food increases which underscores a need for improvement in food security. Disease pressures are a major concern surrounding sustainable agriculture. Static crop populations, containing little to no genetic diversity, are vulnerable to diverse pathogen populations. Wild relatives of crop plants are a reservoir for new disease resistance traits that can be introgressed into cultivated crops. The identification of novel disease resistance is of paramount importance because pathogen co-evolution is not only defeating current resistance genes (R genes) but chemical controls as well. Phytophthora sojae (P. sojae), the causal agent of Phytophthora root and stem rot disease, reduces soybean harvests worldwide. We developed an approach to screen for new R genes that recognize core effectors from P. sojae. We expect R genes identified by these screens to be durable because P. sojae requires core effectors for virulence. We utilized effector-based screening to probe Glycine soja germplasm with core RXLR effectors from P. sojae to search for novel R genes. We developed segregating populations from crosses of P. sojae resistant G. soja germplasm with susceptible G. max cultivar Williams to determine inheritance of potential R genes in germplasm that responded to core effectors. We are using marker assisted breeding to map disease resistance traits in recombinant inbred (RI) lines. To better understand pathosystems, we examined host resistance and susceptibility using bioinformatics. We analyzed the interaction between Arabidopsis thaliana ecotype Col-0 and Hyaloperonospora arabidopsidis isolate Emwa1 using a publicly available RNA time-course experiment. We describe a new algorithm to sort genes into time-point specific clusters using activation and repression parameters. Gene ontology annotations were used to identify defense genes with unique expression profiles, and A. thaliana null mutants for these genes were significantly more susceptible to Emwa1 than wild-type. We plan to use these tools to rapidly identify and guide introgression of durable disease resistance into crop species. Ph. D.

Disease outbreaks caused by oomycetes can be catastrophic. The first part of this dissertation describes development of a system to identify potential new and durable resistance (R) genes against P. sojae in soybean germplasm. We developed a system to screen soybean germplasm for genes that recognize core Phytophthora sojae RXLR effectors that are conserved within the pathogen species and essential for virulence. R genes that recognize these effectors will likely be effective and durable against diverse P. sojae isolates. We developed a system to deliver individual P. sojae effectors by Type III secretion into soybean using the bacterium Pseudomonas, and we screened 12 core effectors on a collection of 30 G. max lines that likely contain new resistance genes against P. sojae. We identified candidate R genes against 10 effectors. Genetic segregation ratios from crosses indicated that three of these genes have a simple inheritance pattern and would be amenable to breeding into elite cultivars. The second part of the dissertation involves use of a model plant-oomycete system to study the genetic basis of susceptibility to oomycete diseases. We compared host transcriptomes from a resistant and a susceptible infection of Arabidopsis thaliana by the downy mildew pathogen Hyaloperonospora arabidopsidis (Hpa). We identified five gene clusters with expression patterns specific to the susceptible interaction. Genes from each cluster were selected and null mutants were tested for altered susceptibility to virulent Hpa. Most A. thaliana null mutants showed enhanced disease susceptibility, suggesting their involvement in pattern-triggered immunity (PTI). A knockout mutant in the AtGcn5 gene was completely resistant to Hpa Emco5 suggesting that the gene/protein is necessary for Hpa to successfully colonize the plant. This study provided new molecular insights into plant-oomycete interaction and revealed a plant gene that could potentially be engineered to provide enhanced resistance to oomycete pathogens. Ph. D.

Phytic acid (PA) accounts for up to 75% of the P in soybean (Glycine max L. Merr.) seeds, but it is indigestible by mono- and agastric animals resulting in economic and environmental detriment. Soybean lines with genetically reduced PA contents have been developed using three distinct mutant alleles at the MIPS1, LPA1, and LPA2 genes resulting in up to a 75% reduction in PA. Low PA (LPA) soymeal-based feeds have been tested on several agricultural species and shown to reduce the P in the animal effluent, but they have not been tested on any aquacultural species. However, LPA soybean lines often exhibit low field emergence making them commercially inviable. The cause of this phenomenon is widely debated with possibilities ranging from increased disease pressure to decreased seedling vigor. The objectives of this research were to 1) enhance field emergence of LPA soybean varieties through pre-planting seed treatments, 2) study the impact of the LPA mutant alleles on agronomic, quality, and seed composition traits, and 3) design a low-error method for studying the effect of LPA soymeal-based feeds on aquatic animals using Pacific White Shrimp (Litoenaeus vannamei). These results describe a variety of agronomic and genetic strategies with which the low field emergence of LPA soybeans can be addressed, reveal a heretofore not reported interaction between the mips1 and lpa2 alleles to further increase the digestibility of soymeal, and a possible method for studying LPA soymeal based feed on aquacultural animals. Master of Science

Switchgrass is a C4 perennial grass that is currently being developed for use as a second generation lignocellulosic biofuel crop. For switchgrass to be fully utilized as a bioenergy crop, large-scale plantings of elite switchgrass germplasm, possibly in monoculture, are likely to occur. This practice may increase the selection pressure on plant pathogens, such as switchgrass rust, which could result in devastating disease epidemics. The identification and deployment of quantitative trait loci (QTLs) and major plant disease resistance genes (R) in switchgrass breeding programs could offer broad spectrum and durable disease resistance in commercial switchgrass cultivars. 'Alamo', a lowland cultivar, is generally resistant to switchgrass rust whereas 'Dacotah', an upland cultivar, is highly susceptible. I hypothesized that major R genes and/or QTLs were contributing to the differences in disease phenotypes of these two cultivars. In this dissertation, bioinformatics and molecular biology approaches were employed to dissect the genetic mechanisms underlying switchgrass rust disease resistance. Novel pseudo-F2 mapping populations were created from a cross derived from 'Alamo' and 'Dacotah'. RNA-sequencing of the pseudo-F2 progenies of 'Alamo' x 'Dacotah' was used to construct a genetic linkage map and to identify potential QTLs correlating with disease resistance. In addition, a homology-based computational method was used to identify 1,011 potential NB-LRR R genes in the switchgrass genome (v 1.1). These potential R genes were characterized for polymorphism and expression differences between 'Alamo' and 'Dacotah'. Moreover, I found that some NB-LRR genes are developmentally regulated in switchgrass. One of the major objectives of switchgrass breeding programs is to develop cultivars with improved feedstock quality; however, changes in the components of the plant cell wall may affect disease resistance. I hypothesized that genetically modified switchgrass plants with altered cell wall components will respond differently than the wild-type to switchgrass rust. Transgenic switchgrass plants overexpressing AtSHN3, a transcription factor with known functions in epicuticular wax accumulation and cell wall deposition, were created. I found that AtSHN3-overexpressing transgenic switchgrass lines were more susceptible than wild-type plants in their response to switchgrass rust. Overall, the results of this dissertation provide a platform for elucidating the molecular mechanisms underlying resistance of switchgrass to switchgrass rust. These findings will help breeders create switchgrass cultivars with improved disease resistance, and will ultimately allow switchgrass to be used for sustainable biomass production. Ph. D.

Soybean is a dominant oilseed cultivated worldwide for its use in multiple sectors such as food and feed industries, animal husbandry, cosmetics and pharmaceutical sectors, and more recently, in production of biodiesel. Increasing demand of soybean, changing environmental conditions, and evolution of pathogens pose challenges to soybean production in limited acreage. Genetic research is the key to ensure the continued growth in soybean production, with enhanced yield and quality, while reducing the losses due to diseases and pests. This research is focused on the understanding of transcriptional regulation of two economically important agronomic traits of soybean: low seed phytic acid and resistance to Soybean mosaic virus (SMV), using the 'transcriptomics' and 'genomics' approaches. The low phytic acid (lpa) soybean is more desirable than conventional soybean, as phytic acid is an anti-nutritional component of seed and is associated with phosphorus pollution. Despite the eco-friendly nature of the lpa soybean, it shows poor emergence, which reduces soybean yield. This research is mainly focused on addressing the impact of lpa-causing mutations on seed development, which is suspected to cause low emergence in lpa soybeans. The differences in transcriptome profiles of developing seeds in lpa and normal phytic acid soybean are revealed and the biological pathways that may potentially be involved in regulation of seed development are suggested. The second research project is focused on Rsv3-type resistance, which is effective against most virulent strains of Soybean mosaic virus. The Rsv3 locus, which maps on to soybean chromosome 14, contains 10 genes including a cluster of coiled coil-nucleotide binding-leucine rich repeat (CC-NB-LRR) protein-encoding genes. This dissertation employed a comparative sequencing approach to narrow down the list of Rsv3 gene candidates to the most promising CC-NB-LRR gene. The evidence provided in this study clearly indicates a single CC-NB-LRR gene as the most promising candidate to deliver Rsv3-type resistance. Ph. D.

The data from genome-wide association studies (GWAS) in humans are still predominantly analyzed using single marker association methods. As an alternative to Single Marker Analysis (SMA), all or subsets of markers can be tested simultaneously. This approach requires a form of Penalized Regression (PR) as the number of SNPs is much larger than the sample size. Here we review PR methods in the context of GWAS, extend them to perform penalty parameter and SNP selection by False Discovery Rate (FDR) control, and assess their performance (including penalties incorporating linkage disequilibrium) in comparison with SMA. PR methods were compared with SMA on realistically simulated GWAS data consisting of genotype data from single and multiple chromosomes and a continuous phenotype and on real data. Based on our comparisons our analytic FDR criterion may currently be the best approach to SNP selection using PR for GWAS. We found that PR with FDR control provides substantially more power than SMA with genome-wide type-I error control but somewhat less power than SMA with Benjamini-Hochberg FDR control. PR controlled the FDR conservatively while SMA-BH may not achieve FDR control in all situations. Differences among PR methods seem quite small when the focus is on variable selection with FDR control. Incorporating LD into PR by adapting penalties developed for covariates measured on graphs can improve power but also generate morel false positives or wider regions for follow-up. We recommend using the Elastic Net with a mixing weight for the Lasso penalty near 0.5 as the best method. Ph. D.

Fusarium Head Blight (FHB) is a devastating fungal disease of wheat (Triticum aestivum L.) and results in significant economic losses due to reductions in grain yield and the accumulation of mycotoxins, such as deoxynivalenol (DON) and nivalenol (NIV). As a result, breeding programs have been working to identify resistance genes in wheat varieties known to be resistant to FHB. Some of the major quantitative trait loci (QTL) for FHB resistance identified to date have been from exotic sources such as 'Sumai3' and the Chinese landrace Wangshuibai, and native resistance has been identified in North American cultivars such as Ernie and Truman which are being used in breeding programs. This study was conducted to characterize and map QTL for resistance to FHB in the soft red winter wheat cultivar Jamestown and to identify tightly linked DNA markers associated with those QTL so that marker-assisted selection (MAS) can be used in pyramiding these and other known QTL into elite backgrounds. Types of resistance assessed in this study include: Type I (resistance to initial infection; incidence), Type II (resistance to spread in wheat spike; severity), and decreases in mycotoxin accumulation (DON) and percentage of Fusarium damaged kernels (FDK). A population composed of 186 F5:7 recombinant-inbred lines (RILs) from the cross Pioneer Brand '25R47'/Jamestown were used to evaluate these traits in six environments (MD, NC, and VA in 2011 and 2012). This study identified a QTL for resistance to DON accumulation and FHB severity on the wheat chromosome 1B. The QTL accounted for 12.7% to 13.3% of the phenotypic variation in DON accumulation and 26.1% of the phenotypic variation in FHB severity. The most diagnostic marker for the QTL on chromosome 1B associated with resistance to FHB severity and DON accumulation is Xwmc500.6 located 7.2 cM from the QTL peak and flanked by markers Xwmc500.7 and Xgwm273.2 (28.2 cM interval). Master of Science

Phytate is the major storage form of phosphorus in seeds of soybeans. Because phytate chelates mineral cations including calcium, iron, and zinc, these mixed salts are often excreted by non-ruminant animals such as humans, swine, poultry, and fish. While this causes iron and zinc deficiencies, phytate is also considered a water pollutant due to the excess phosphorus excreted in animal waste. These negative environmental and nutritional effects, create a need for low phytate soybeans. While several low phytate soybean lines have been developed, a major drawback is the reduced seedling emergence of these lines resulting in low yields. Therefore, understanding the genetic and molecular bases of the low emergence trait in relation to seed phytate content in major crops such as soybean is of great economic importance. This PhD project worked towards the long term goal of developing low phytate soybean cultivars with good seedling emergence and high-yield. This dissertation focused on metabolomic differences between low and normal phytate lines and how these could relate to the low emergence phenotype. The genetic materials used here include four near isogenic lines that differ in mutations in two multi drug resistance-associated proteins (MRPs). Only the line with both mutations was low phytate. The phytate levels, field- and lab-based emergence rates were determined for these lines, their parents and a control line through replicated field experiments for three consecutive years. The emergence rates of the low phytate lines were not always reduced. This showed that the environment the seeds were produced in is highly important, especially when breeding and commercially growing low phytate lines. A protocol was developed for successful metabolomic discrimination of these closely related soybean lines. The polar and non-polar metabolite profiles were determined using ultra performance liquid chromatography mass spectrometry and metabolomic differences between the low and normal phytate lines were identified. The low phytate double mutant did not contain C22 glucose terminated Group A soyasaponins and almost exclusively contained C22 xylose terminated Group A soyasaponins (A4, A5 and A6). Compared to the normal phytate lines, the low phytate soybean line showed a higher concentration of storage lipids (triacylglycerols and diacylglycerols) and certain phospholipids. Ph. D.

Powdery mildew, caused by the pathogen Blumeria graminis (DC) Speer (Syn. Erysiphe graminis DC) f. sp. tritici, is a major disease of wheat (Triticum aestivum L.). Race-specific resistance is easily identified in the field due to its qualitative phenotype and it is easy to incorporate because it is inherited as a single gene. Unfortunately, this type of resistance is easily overcome by the pathogen. Traits associated with quantitative trait loci (QTL) such as adult-plant resistance (APR), have become popular with plant breeders because of their durability over a wide geographic range and time. Due to the quantitative nature of these genes, they are difficult to study requiring multiple assessments of disease development under natural conditions in more than one location over a period of several weeks. Numerous QTL for APR to powdery mildew have been mapped in independent studies in different wheat backgrounds. The wheat cultivar Massey has been the subject of several studies due to its APR to powdery mildew that has remained effective for several decades. However, it has been difficult to identity simple sequence repeat (SSR) markers that are tightly linked to the QTL for APR in Massey. Such markers give breeders an advantage by allowing them to quickly identify and select for traits that would be difficult to distinguish in the field among breeding progeny from several backgrounds. Therefore, identification of tightly linked markers associated with APR to powdery mildew is necessary so that these traits can be selected for reliably in progeny. Master of Science

Potato (Solanum tuberosum L) is an outcrossing, highly heterozygous cultivated in which the elucidation of the genetic basis of quantitative traits, is more complex than in self-pollinated crops. Both a candidate gene approach and a whole genome SNP genotyping analysis were used to assess allelic variation and to identify loci associated with biosynthesis and accumulation of steroidal glycoalkaloids (SGAs). SGAs are secondary metabolites produced in Solanum species as defense against insects and pathogens. Fragments of genomic DNA coding for regions of five SGA biosynthetic candidate genes were amplified, cloned and sequenced [3-hydroxy-3-methylglutaryl coenzyme A reductase 1 and 2 (HMG1, HMG2); 2,3-squalene epoxidase (SQE); solanidine galactosyltransferase (SGT1); and solanidine glucosyltransferase (SGT2)]. A germplasm panel of six wild potato species [Solanum chacoense (chc 80-1), S. commersonii subsp. commersonii, S. demissum, S. sparsipilum, S. spegazzinii, and S. stoloniferum] and a cultivated clone S. tuberosum Group Phureja (phu DH) was used in an allelic variation analysis. A segregating interspecific F2 population phu DH �" chc 80-1 was screened to assess association with SGAs. Sequence diversity analysis showed a tendency of purifying selection and increased frequency of rare alleles in most of the candidate genes. Genes of primary metabolism (HMG1, HMG2 and SQE) had stronger selection constraints than those in secondary metabolism (SGT1 and SGT2). Sequence polymorphism in HMG2, SQE, SGT1 and SGT2 separated either the phu DH clone which produced no SGAs, or chc 80-1, the greatest SGA accumulator, from other accessions in the panel. Segregation analysis of the F2 population revealed that allelic sequences of HMG2 and SGT2 derived from chc 80-1 were significantly associated with the greatest SGA accumulation. In the whole genome analysis, SNP genotyping and cluster analysis based on putative association with SGA accumulation in the germplasm panel, allowed identification of eight informative SNPs that can be used in future studies. In the segregating F2 population, loci located on five pseudochromosomes were associated with SGA synthesis. Loci on pseudochromosomes 1 and 6 explained segregation ratios of synthesis for α-solanine and α-chaconine, the most common SGAs in most potato species. In addition, loci on seven pseudochromosomes were associated with accumulation. New candidate genes, putatively affecting synthesis and accumulation of SGAs, were identified in adjacent genomic regions of significant SNPs. This research demonstrates how the newly available genome sequence of potato and associated biotechnological tools accelerates the identification of genetic factors underling complex traits in a species with a difficult breeding structure. Ph. D.

Phenotypic and genotypic characterization of hulled and hulless winter barley (Hordeum vulgare L.) is necessary for improvement using traditional and molecular breeding techniques. Identification of genomic regions conferring resistance to Fusarium head blight (caused by Fusarium graminearum), leaf rust (caused by Puccinia hordei G. Otth), powdery mildew [caused by Blumeria graminis (DC.) E.O. Speer f. sp. hordei Em. Marchal], net blotch (caused by Pyrenophora teres) and spot blotch [caused by Cochliobolus sativus (Ito & Kuribayashi) Drechsler ex Dastur] will greatly aid in breeding for improved resistance. Determining factors that contribute to yield differences between hulled and hulless genotypes, and identification of markers associated with yield and yield related traits will greatly aid in improvement of hulled and hulless genotypes. The hulled cultivar Nomini, hulless cultivar Eve, and hulless line VA06H-48 were consistently resistant to Fusarium head blight (FHB) and had low deoxynivalenol (DON) accumulation. Screening with molecular markers on chromosomes 2H and 6H for FHB and DON identified quantitative trait loci (QTL) which may confer resistance in Virginia Tech germplasm. Evaluation of hulled and hulless full-sibs from four populations indicated that grain volume weight and protein concentration were significantly (P d 0.05) higher for hulless genotypes, while seedling emergence and grain ash concentration were significantly (P d 0.05) higher for hulled genotypes. In linear regression analysis, none of the assessed traits explained yield variation in all populations and environments. Identification of hulless genotypes having yield potentials similar to those of their hulled sibs should be possible after adjusting for hull weight. A genome wide association study was used to identify chromosome regions governing traits of importance in six-rowed winter barley germplasm and to identify single nucleotide polymorphisms (SNPs) markers for use in a marker-assisted breeding program. Significant SNPs associated with previously described QTL or genes were identified for heading date, test weight, yield, grain protein, polyphenol oxidase activity, and resistance to leaf rust, powdery mildew, net blotch, and spot blotch. Novel QTL also were identified for agronomic, quality, and disease resistance traits. These SNP-trait associations provide the opportunity to directly select for QTL contributing to multiple traits in breeding programs. Ph. D.

Since 2000, stripe rust, caused by the fungal pathogen (Puccinia striiformis Westend. f.sp. tritici Eriks.) has resulted in yield losses of wheat (Triticum aestivium L.) in the United States, that exceeded the combined losses of leaf rust (Puccinia triticina Eriks.) and stem rust (Puccinia gramins Pers.:Pers f. sp. Tritici Eriks. E. Henn.). The objectives of this study are to identify and map adult plant stripe rust resistance quantitative trait loci (QTL) in soft red winter (SRW) wheat that are effective against race PST-100, which is the predominant race of the pathogen in the eastern U.S. Adult plant resistance (APR) was characterized in the resistant wheat lines "USG 3555", VA00W-38, and "Coker 9553". Resistance in each of the lines was evaluated in populations derived from crosses with susceptible lines "Neuse", Pioneer Brand "26R46", and VA01W-21, respectively. On chromosomes 1AS, 4BL, and 7D of USG 3555, QTL were identified that explain on average 12.8, 73.0, and 13.6 percent of the variation for stripe rust infection type (IT), and 13.5, 72.3, and 10.5 percent of the variation for stripe rust severity. A QTL from Neuse was identified on 3A that on average explains 10.9 percent of the variation for IT and 13.0 percent of the variation for severity. On chromosomes 2AS and 4BL of VA00W-38, QTL were identified that on average explain 58.9 and 19.3 percent of the variation for stripe rust IT, and 51.9 and 12.1 percent of the variation for severity. On chromosomes 6BL and 3BL of Pioneer 26R46, QTL were identified that on average explain 8.9 and 2.1 percent of the variation for IT and 11.7 and 3.9 percent of the variation for severity. Coker 9553 possesses the QTL on 4BL that is also present in USG 3555 and VA00W-38. The SSR markers, Xgwm296, Xbarc163, and Xwmc756 were tightly linked to QTL on chromosomes 2AS, 4BL, and 6BL, respectively, and their use and development of additional diagnostic markers will facilitate the incorporation and pyramiding of stripe rust resistance QTL into SRW wheat lines via marker-assisted selection. Ph. D.