Application of an Improved 2-Dimensional High-Throughput Soybean Root Phenotyping Platform to Identify Novel Genetic Variants Regulating Root Architecture Traits

Nutrient-efficient root system architecture (RSA) is becoming an important breeding objective for generating crop varieties with improved nutrient and water acquisition efficiency. Genetic variants shaping soybean RSA is key in improving nutrient and water acquisition. Here, we report on the use of an improved 2-dimensional high-throughput root phenotyping platform that minimizes background noise by imaging pouch-grown root systems submerged in water. We also developed a background image cleaning Python pipeline that computationally removes images of small pieces of debris and filter paper fibers, which can be erroneously quantified as root tips. This platform was used to phenotype root traits in 286 soybean lines genotyped with 5.4 million single-nucleotide polymorphisms. There was a substantially higher correlation in manually counted number of root tips with computationally quantified root tips (95% correlation), when the background was cleaned of nonroot materials compared to root images without the background corrected (79%). Improvements in our RSA phenotyping pipeline significantly reduced overestimation of the root traits influenced by the number of root tips. Genome-wide association studies conducted on the root phenotypic data and quantitative gene expression analysis of candidate genes resulted in the identification of 3 putative positive regulators of root system depth, total root length and surface area, and root system volume and surface area of thicker roots (DOF1-like zinc finger transcription factor, protein of unknown function, and C2H2 zinc finger protein). We also identified a putative negative regulator (gibberellin 20 oxidase 3) of the total number of lateral roots.

The transcription factors, STOP1 and TCP20, are required for root system architecture alterations in response to nitrate deficiency

Nitrate distribution in soils is often heterogeneous. Plants have adapted to this by modifying their root system architecture (RSA). Previous studies showed that NITRATE-TRANSPORTER1.1 (NRT1.1), which also transports auxin, helps inhibit lateral root primordia (LRP) emergence in nitrate-poor patches, by preferentially transporting auxin away from the LRP. In this study, we identified the regulatory system for this response involving the transcription factor (TF), SENSITIVE-TO-PROTON-RHIZOTOXICITY1 (STOP1), which is accumulated in the nuclei of LRP cells under nitrate deficiency and directly regulates Arabidopsis NRT1.1 expression. Mutations in STOP1 mimic the root phenotype of the loss-of-function NRT1.1 mutant under nitrate deficiency, compared to wild-type plants, including increased LR growth and higher DR5promoter activity (i.e., higher LRP auxin signaling/activity). Nitrate deficiency-induced LR growth inhibition was almost completely reversed when STOP1 and the TF, TEOSINTE-BRANCHED1,-CYCLOIDEA,-PCF-DOMAIN-FAMILY-PROTEIN20 (TCP20), a known activator of NRT1.1 expression, were both mutated. Thus, the STOP1-TCP20 system is required for activation of NRT1.1 expression under nitrate deficiency, leading to reduced LR growth in nitrate-poor regions. We found this STOP1-mediated system is more active as growth media becomes more acidic, which correlates with reductions in soil nitrate as the soil pH becomes more acidic. STOP1 has been shown to be involved in RSA modifications in response to phosphate deficiency and increased potassium uptake, hence, our findings indicate that root growth regulation in response to low availability of the major fertilizer nutrients, nitrogen, phosphorus and potassium, all involve STOP1, which may allow plants to maintain appropriate root growth under the complex and varying soil distribution of nutrients.

Comparative transmission genetics of introgressed chromatin in reciprocal advanced backcross populations in Gossypium (cotton) polyploids

Introgression is a potential source of valuable genetic variation and interspecific introgression lines are important resources for plant breeders to access novel alleles. Experimental advanced-generation backcross populations contain individuals with genomic compositions similar to those resulting from natural interspecific hybridization and provide opportunities to study the nature and transmission pattern of donor chromatin in recipient genomes. Here, we analyze transmission of donor chromatin in reciprocal backcrosses between G. hirsutum and G. barbadense. Across the genome, recurrent backcrossing in both backgrounds yielded donor chromatin at slightly higher frequencies than the Mendelian expectation in BC5F1 plants, while the average frequency of donor alleles in BC5F2 segregating families was less than expected. In the two subgenomes of polyploid cotton, the rate of donor chromatin introgression was similar. Although donor chromatin was tolerated over much of the recipient genomes, 21 regions recalcitrant to donor alleles were identified. Only limited correspondence is observed between the recalcitrant regions in the two backgrounds, suggesting the effect of species background on introgression of donor segments. Genetic breakdown was progressive, with floral abscission and seed inviability ongoing during backcrossing cycles. Regions of either high or low introgression tended to be in terminal chromosomal regions that are generally rich in both genes and crossover events, with long stretches around the centromere having limited crossover activity resulting in relatively constant low introgression frequencies. Constraints on fixation and selection of donor alleles highlights the challenges of utilizing introgression breeding in crop improvement.

Iterative image segmentation of plant roots for high-throughput phenotyping

Accurate segmentation of root system architecture (RSA) from 2D images is an important step in studying phenotypic traits of root systems. Various approaches to image segmentation exist but many of them are not well suited to the thin and reticulated structures characteristic of root systems. The findings presented here describe an approach to RSA segmentation that takes advantage of the inherent structural properties of the root system, a segmentation network architecture we call ITErRoot. We have also generated a novel 2D root image dataset which utilizes an annotation tool developed for producing high quality ground truth segmentation of root systems. Our approach makes use of an iterative neural network architecture to leverage the thin and highly branched properties of root systems for accurate segmentation. Rigorous analysis of model properties was carried out to obtain a high-quality model for 2D root segmentation. Results show a significant improvement over other recent approaches to root segmentation. Validation results show that the model generalizes to plant species with fine and highly branched RSA’s, and performs particularly well in the presence of non-root objects.

Improved Upland Cotton Germplasm for Multiple Fiber Traits Mediated by Transferring and Pyramiding Novel Alleles From Ethyl Methanesulfonate-Generated Mutant Lines Into Elite Genotypes

Ethyl methanesulfonate (EMS) mutagenesis offers important advantages for improving crops, such as cotton, with limited diversity in elite gene pools. EMS-induced point mutations are less frequently associated with deleterious traits than alleles from wild or exotic germplasm. From 157 mutant lines that have significantly improved fiber properties, we focused on nine mutant lines here. A total of eight populations were developed by crossing mutant lines in different combinations into GA230 (GA2004230) background. Multiple lines in each population were significantly improved for the fiber trait that distinguished the donor parent(s), demonstrating that an elite breeding line (GA230) could be improved for fiber qualities using the mutant lines. Genotypes improved for multiple fiber traits of interest suggesting that allele pyramiding is possible. Compared to midparent values, individual progeny in the population conferred fiber quality improvements of as much as 31.7% (in population O) for micronaire (MIC), 16.1% (in population P) for length, 22.4% (in population K) for strength, 4.1% (in population Q) for uniformity, 45.8% (in population N) for elongation, and 13.9% (in population O) for lint percentage (lint%). While further testing for stability of the phenotype and estimation of yield potential is necessary, mutation breeding shows promise as an approach to reduce the problem of the genetic bottleneck of upland cotton. The populations developed here may also contribute to identifying candidate genes and causal mutations for fiber quality improvement.

Identification of small effect quantitative trait loci of plant architectural, flowering, and early maturity traits in reciprocal interspecific introgression population in cotton

Plant architecture, flowering time and maturity traits are important determinants of yield and fiber quality of cotton. Genetic dissection of loci determining these yield and quality components is complicated by numerous loci with alleles conferring small differences. Therefore, mapping populations segregating for smaller numbers and sizes of introgressed segments is expected to facilitate dissection of these complex quantitative traits. At an advanced stage in the development of reciprocal advanced backcross populations from crosses between elite Gossypium hirsutum cultivar 'Acala Maxxa' (GH) and G. barbadense 'Pima S6' (GB), we undertook mapping of plant architectural traits, flowering time and maturity. A total of 284 BC₄F₁ and BC₄F₂ progeny rows, 120 in GH and 164 in GB background, were evaluated for phenotype, with only 4 and 3 (of 7) traits showing significant differences among progenies. Genotyping by sequencing yielded 3,186 and 3,026 SNPs, respectively, that revealed a total of 27 QTLs in GH background and 22 in GB, for plant height, days to flowering, residual flowering at maturity and maturity. More than of 90% QTLs identified in both backgrounds had small effects (%PV < 10), supporting the merit of this population structure to reduce background noise and small effect QTLs. Germplasm developed in this study may serve as potential pre-breeding material to develop improved cotton cultivars.

Root Adaptation via Common Genetic Factors Conditioning Tolerance to Multiple Stresses for Crops Cultivated on Acidic Tropical Soils

Crop tolerance to multiple abiotic stresses has long been pursued as a Holy Grail in plant breeding efforts that target crop adaptation to tropical soils. On tropical, acidic soils, aluminum (Al) toxicity, low phosphorus (P) availability and drought stress are the major limitations to yield stability. Molecular breeding based on a small suite of pleiotropic genes, particularly those with moderate to major phenotypic effects, could help circumvent the need for complex breeding designs and large population sizes aimed at selecting transgressive progeny accumulating favorable alleles controlling polygenic traits. The underlying question is twofold: do common tolerance mechanisms to Al toxicity, P deficiency and drought exist? And if they do, will they be useful in a plant breeding program that targets stress-prone environments. The selective environments in tropical regions are such that multiple, co-existing regulatory networks may drive the fixation of either distinctly different or a smaller number of pleiotropic abiotic stress tolerance genes. Recent studies suggest that genes contributing to crop adaptation to acidic soils, such as the major Arabidopsis Al tolerance protein, AtALMT1, which encodes an aluminum-activated root malate transporter, may influence both Al tolerance and P acquisition via changes in root system morphology and architecture. However, trans-acting elements such as transcription factors (TFs) may be the best option for pleiotropic control of multiple abiotic stress genes, due to their small and often multiple binding sequences in the genome. One such example is the C2H2-type zinc finger, AtSTOP1, which is a transcriptional regulator of a number of Arabidopsis Al tolerance genes, including AtMATE and AtALMT1, and has been shown to activate AtALMT1, not only in response to Al but also low soil P. The large WRKY family of transcription factors are also known to affect a broad spectrum of phenotypes, some of which are related to acidic soil abiotic stress responses. Hence, we focus here on signaling proteins such as TFs and protein kinases to identify, from the literature, evidence for unifying regulatory networks controlling Al tolerance, P efficiency and, also possibly drought tolerance. Particular emphasis will be given to modification of root system morphology and architecture, which could be an important physiological "hub" leading to crop adaptation to multiple soil-based abiotic stress factors.

The Ligon lintless -2 Short Fiber Mutation Is Located within a Terminal Deletion of Chromosome 18 in Cotton

Extreme elongation distinguishes about one-fourth of cotton (Gossypium sp.) seed epidermal cells as "lint" fibers, useful for the textile industry, from "fuzz" fibers (<5 mm). Ligon lintless-2 (Li ₂ ), a dominant mutation that results in no lint fiber but normal fuzz fiber, offers insight into pathways and mechanisms that differentiate spinnable cotton from its progenitors. A genetic map developed using 1,545 F2 plants showed that marker CISP15 was 0.4 cM from Li ₂ , and "dominant" simple sequence repeat (SSR) markers (i.e. with null alleles in the Li ₂ genotype) SSR7 and SSR18 showed complete linkage with Li ₂ Nonrandom distribution of markers with null alleles suggests that the Li ₂ phenotype results from a 176- to 221-kb deletion of the terminal region of chromosome 18 that may have been masked in prior pooled-sample mapping strategies. The deletion includes 10 genes with putative roles in fiber development. Two Glycosyltransferase Family 1 genes showed striking expression differences during elongation of wild-type versus Li ₂ fiber, and virus-induced silencing of these genes in the wild type induced Li ₂ -like phenotypes. Further, at least 7 of the 10 putative fiber development genes in the deletion region showed higher expression in the wild type than in Li ₂ mutants during fiber development stages, suggesting coordinated regulation of processes in cell wall development and cell elongation, consistent with the hypothesis that some fiber-related quantitative trait loci comprise closely spaced groups of functionally diverse but coordinately regulated genes.

Genetic Analysis of Gossypium Fiber Quality Traits in Reciprocal Advanced Backcross Populations

In mapping populations segregating for many loci, the large amount of variation among genotypes often masks small-effect quantitative trait loci (QTL). This problem can be reduced by development of populations with fewer chromosome segments segregating. Here, we report early QTL detection in reciprocal advanced backcross populations from crosses between elite Gossypium hirsutum L. 'Acala Maxxa' (GH) and G. barbadense L. 'Pima S6' (GB). A total of 297 BCF and BCF progeny rows-127 segregating for GB chromosome segments in GH background and 170 segregating for GH chromosome segments in GB background-were evaluated in three environments. Totals of 3186 and 3026 polymorphic single-nucleotide polymorphisms (SNPs) in GH and GB backgrounds, respectively, were identified and used for trait mapping. Small-effect QTL (<10% variance explained) made up 87 and 100% of QTL in GH and GB backgrounds, respectively. In both species, favorable alleles were found with effects being masked or neutralized by unfavorable alleles, with greater scope for improvement of GH than GB by introgressive breeding. A total of three stable QTL-two in GH background for fiber elongation (ELO) and micronaire (MIC) and one in GB background for upper-half mean length (UHM)-were identified in two out of three environments. Curiously, only four QTL-three for UHM and one for ELO-showed the expected opposite effects in reciprocal backgrounds, perhaps reflecting the combined consequences of epistasis, small phenotypic effects, and low coverage of some genomic regions. Along with new information for marker-assisted breeding, this study adds to knowledge that can be used to unravel complex genetic networks governing fiber quality traits.

Genetic Analysis of Gossypium Fiber Quality Traits in Reciprocal Advanced Backcross Populations

In mapping populations segregating for many loci, the large amount of variation among genotypes often masks small-effect quantitative trait loci (QTL). This problem can be reduced by development of populations with fewer chromosome segments segregating. Here, we report early QTL detection in reciprocal advanced backcross populations from crosses between elite Gossypium hirsutum L. 'Acala Maxxa' (GH) and G. barbadense L. 'Pima S6' (GB). A total of 297 BCF and BCF progeny rows-127 segregating for GB chromosome segments in GH background and 170 segregating for GH chromosome segments in GB background-were evaluated in three environments. Totals of 3186 and 3026 polymorphic single-nucleotide polymorphisms (SNPs) in GH and GB backgrounds, respectively, were identified and used for trait mapping. Small-effect QTL (<10% variance explained) made up 87 and 100% of QTL in GH and GB backgrounds, respectively. In both species, favorable alleles were found with effects being masked or neutralized by unfavorable alleles, with greater scope for improvement of GH than GB by introgressive breeding. A total of three stable QTL-two in GH background for fiber elongation (ELO) and micronaire (MIC) and one in GB background for upper-half mean length (UHM)-were identified in two out of three environments. Curiously, only four QTL-three for UHM and one for ELO-showed the expected opposite effects in reciprocal backgrounds, perhaps reflecting the combined consequences of epistasis, small phenotypic effects, and low coverage of some genomic regions. Along with new information for marker-assisted breeding, this study adds to knowledge that can be used to unravel complex genetic networks governing fiber quality traits.

Segregation distortion and genome-wide digenic interactions affect transmission of introgressed chromatin from wild cotton species

This study reports transmission genetics of chromosomal segments into Gossypium hirsutum from its most distant euploid relative, Gossypium mustelinum . Mutilocus interactions and structural rearrangements affect introgression and segregation of donor chromatin. Wild allotetraploid relatives of cotton are a rich source of genetic diversity that can be used in genetic improvement, but linkage drag and non-Mendelian transmission genetics are prevalent in interspecific crosses. These problems necessitate knowledge of transmission patterns of chromatin from wild donor species in cultivated recipient species. From an interspecific cross, Gossypium hirsutum × Gossypium mustelinum, we studied G. mustelinum (the most distant tetraploid relative of Upland cotton) allele retention in 35 BC₃F₁ plants and segregation patterns in BC₃F₂ populations totaling 3202 individuals, using 216 DNA marker loci. The average retention of donor alleles across BC₃F₁ plants was higher than expected and the average frequency of G. mustelinum alleles in BC₃F₂ segregating families was less than expected. Despite surprisingly high retention of G. mustelinum alleles in BC₃F₁, 46 genomic regions showed no introgression. Regions on chromosomes 3 and 15 lacking introgression were closely associated with possible small inversions previously reported. Nonlinear two-locus interactions are abundant among loci with single-locus segregation distortion, and among loci originating from one of the two subgenomes. Comparison of the present results with those of prior studies indicates different permeability of Upland cotton for donor chromatin from different allotetraploid relatives. Different contributions of subgenomes to two-locus interactions suggest different fates of subgenomes in the evolution of allotetraploid cottons. Transmission genetics of G. hirsutum × G. mustelinum crosses reveals allelic interactions, constraints on fixation and selection of donor alleles, and challenges with retention of introgressed chromatin for crop improvement.

Application of genotyping by sequencing technology to a variety of crop breeding programs

Since the Arabidopsis genome was completed, draft sequences or pseudomolecules have been published for more than 100 plant genomes including green algae, in large part due to advances in sequencing technologies. Advanced DNA sequencing technologies have also conferred new opportunities for high-throughput low-cost crop genotyping, based on single-nucleotide polymorphisms (SNPs). However, a recurring complication in crop genotyping that differs from other taxa is a higher level of DNA sequence duplication, noting that all angiosperms are thought to have polyploidy in their evolutionary history. In the current article, we briefly review current genotyping methods using next-generation sequencing (NGS) technologies. We also explore case studies of genotyping-by-sequencing (GBS) applications to several crops differing in genome size, organization and breeding system (paleopolyploids, neo-allopolyploids, neo-autopolyploids). GBS typically shows good results when it is applied to an inbred diploid species with a well-established reference genome. However, we have also made some progress toward GBS of outcrossing species lacking reference genomes and of polyploid populations, which still need much improvement. Regardless of some limitations, low-cost and multiplexed genotyping offered by GBS will be beneficial to breed superior cultivars in many crop species.

Alleles conferring improved fiber quality from EMS mutagenesis of elite cotton genotypes

Genetic improvements for many fiber traits are obtained by mutagenesis of elite cottons, mitigating genetic uniformity in this inbred polyploid by contributing novel alleles important to ongoing crop improvement. The elite gene pool of cotton (Gossypium spp.) has less diversity than those of most other major crops, making identification of novel alleles important to ongoing crop improvement. A total of 3,164 M5 lines resulting from ethyl methanesulfonate (EMS) mutagenesis of two G. hirsutum breeding lines, TAM 94L-25 and Acala 1517-99, were characterized for basic components of fiber quality and selected yield components. Across all measured traits, the ranges of phenotypic values among the mutant lines were consistently larger than could be explained by chance (5.27-10.1 for TAM 94 L-25 and 5.29-7.94 standard deviations for Acala 1517-99-derived lines). Multi-year replicated studies confirmed a genetic basis for these differences, showing significant correlations between lines across years and environments. A subset of 157 lines selected for superior fiber qualities, including fiber elongation (22 lines), length (22), lint percent (17), fineness (23), Rd value (21), strength (19), uniformity (21) and multiple attributes in a selection index (26) were compared to 55 control lines in replicated trials in both Texas and Georgia. For all traits, mutant lines showing substantial and statistically significant improvements over control lines were found, in most cases from each of the two genetic backgrounds. This indicates that genetic improvements for a wide range of fiber traits may be obtained from mutagenesis of elite cottons. Indeed, lines selected for one fiber trait sometimes conferred additional attributes, suggesting pleiotropic effects of some mutations and offering multiple benefits for the incorporation of some alleles into mainstream breeding programs.