Development of Active Root System Architecture of Upland Cotton for Improved Sub-surface Water Uptake During Drought Conditions

Progress report for GS23-295

Project Type: Graduate Student
Funds awarded in 2023: $15,900.00
Projected End Date: 08/31/2025
Grant Recipient: Texas Tech University
Region: Southern
State: Texas
Graduate Student:
Major Professor:
Dr. Gunvant Patil
Texas Tech University
Expand All

Project Information

Summary:

Irrigation can be a short-term solution to water scarcity. For many areas of the world, long-term irrigation use can be unsustainable for cotton agriculture. Cotton typically has some drought tolerance; however, the inability to easily phenotype the root system architecture (RSA) has resulted in elite cultivars with reduced growth and plasticity in root systems. To explore cotton’s phenotypic variation of RSA in the germplasm, 660 accessions from around the globe were narrowed down to 200 accessions for RSA phenotyping analysis. Four seed companies each contributed three of their elite cultivars to this research. The accessions were simultaneously grown in the greenhouse and field locations. The greenhouse evaluated baseline RSA under well-watered conditions. The field tested the accessions under drought stress and compared the RSA of seven accessions. The rainfed condition of West Texas during the 2024 summer season allowed sufficient drought stress to affect most of the accessions significantly. A minimum of fifteen accessions had no significant difference between the irrigated and rainfed conditions. Many accessions also had a more robust RSA with higher root length, greater root-to-shoot ratio, and finer average diameter. A pattern of thicker primary root diameter, but finer lateral roots, performed better under drought conditions, along with a higher root length and a higher root-to-shoot ratio, which contributed to performing better under drought conditions. Using the phenotypic data of this study will allow for genetic associations to be completed, and identification of potential genes involved in the RSA will help future research and breeders develop cotton accessions that can integrate the benefits of an increased RSA for drought tolerance. 

Project Objectives:

First Year's Objectives 

1) Identify a subset of Upland cotton accessions (175-200) and elite commercial cultivars (15) by selecting from the 600 accessions grown in the greenhouse.  

2) a) Evaluate the root system architecture (RSA) of the diverse germplasm using image-based root phenotyping platform under controlled greenhouse conditions. The plants will be grown in the greenhouse for about 2 weeks before scanning their RSA with the scanner and software program. b) Select a subset of the 5 most promising and the 5 least promising of the exotic accessions in addition to the control and elite cultivars for contrast. To evaluate the RSA traits, under well-watered and water-deficit treatments in field conditions, five replications will be under irrigated conditions and the other will be rainfed only. 

3) Based on the 1st year’s findings, lines with contrasting RSA will be further explored to identify QTLs or genomic loci using advanced genomic tools (GWAS) and additional field experiments to identify novel genes for optimizing the cotton RSA for enhanced water capture and efficiency. 

Second Year's Objectives

Objective 1: Evaluation of Root System Architecture (RSA) under greenhouse conditions 

Objective 2: Evaluation of selected accessions under field conditions 

Objective 3: Analysis of phenotypes and genotypes through genome-wide association mapping 

These objectives are the foundational elements of a continuous research project to understand our long-term goal: identify novel QTL/genes for optimizing cotton RSA for enhanced water capture and use efficiency. Based on the 2nd year’s findings, encompassing morphological, anatomical, and physiological traits, lines with contrasting RSA will be further explored to identify gene/QTLs using advanced genomics tools and additional field experiments with continued support from TSSC and other sources.   

Research

Materials and methods:

In the first year, the objectives were to determine the cotton varieties for this project;660 cotton accessions were ordered from USDA-GRIN. Those were planted in the greenhouse to grow for seed accession. Those were screened using phenotypes and their origins to reduce to 200 accessions to move to the next objective. The elite commercial cultivars grown in West Texas were collected from various seed industry partners, including BASF, Bayer, Armor, Phytogen, and NexGen. These companies selected varieties they were interested in testing for drought tolerance in this project. Some varieties already have some drought tolerance characteristics. The main goal of this objective is to get a diverse set of genetic variation to compare root systems for the best and the worst RSA, as well as to compare water efficiency for sustainable water use. Traits related to yield, disease resistance, and fiber quality will also be considered.

For the second objective, all the selected lines were screened for root morphology at the seedling stage, 2 weeks after planting. They were screened using a visual scoring and high-throughput phenotyping platform, WinRhizo, Inc. WinRHIZO is an image analysis system specifically designed for measuring root morphological traits. We focused on total root length, diameter, and surface area among other characteristics. According to Kadam et al. (2015), this process is done in Dr. Jagadish’s lab, which has extensive experience working with root morphology and anatomy of various crops. These experiments will help us narrow the lines for in-depth analysis under well-watered and water-deficient conditions.

Washing Set Up - 2023 - Greenhouse WinRHIZO Scanner UsedScanned Image Example

The next phase of objective two is to identify lines with differential response to water deficit conditions regarding an altered RSA. To achieve this goal, we plan on growing a subset of contrasting lines in controlled-environmental greenhouse conditions under well-watered 100% field capacity and water-deficit conditions around 40%-60% field capacity (Kadam et al., 2015). Water deficit stress will be imposed after seedling emergence. Plants under the control treatment will be maintained at 100% field capacity throughout the experiments. After 30 days of stress, plants will be harvested around 50 days after sowing, and traits associated with the RSA, such as shoot biomass, will be measured. Physiology will be collected as detailed in Kadam et al., 2015. The anatomical studies will consider the root diameter, stele diameter, and metaxylem number.

After evaluating RSA and anatomical traits from objective 2, cotton cultivars showing differential or contrasting phenotypes were further evaluated under field conditions. Uprooting cotton plants in the field would have been challenging, so large woven bags filled with field soil were used to simulate the field conditions to overcome this issue. Watering was done manually by adding the recommended amount determined by the University of Georgia.

 

Cotton Weekly Water Requirement
Image Found at: https://site.extension.uga.edu/worthag/2020/05/early-season-irrigation-requirements-for-cotton-production/

The setup of the bags used during the 2023 summer season.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

During the cotton growth, the morphological (shoot) traits were measured. These accessions were grown and intended to be harvested at the boll formation stage (i.e., the peak root biomass accumulation phase). However, a hailstorm destroyed some of the bags, and some were damaged. The time it took to wash these manually was longer than expected. Therefore, this experiment had to be repeated over the winter in the greenhouse with the same field soil and the same bags.

Washing of RSA from Bag Fall 2023 FieldPhoto of the RSA after washing, Fall 2023 Field

 

 

 

 

 

 

 

 

 

 

The repeated experiment was done in the greenhouse because of the winter conditions. With the limited space in the greenhouse and the amount of soil needed to transfer to the greenhouse from the field, the number of genotypes was reduced to 12 instead of the original 24. This was intended to have three replications of 12 rather than fewer replications. Watering was done similarly to manually measuring the field with the graduated cylinder. The volume and weight were calculated for the soil density, and the soil was set to fit within the same mass per volume to get the ideal density and compaction.

A comprehensive analysis of diverse germplasm and elite commercial cultivars grown in West Texas will allow us to identify the functional role of root tissue plasticity in adapting to water-deficit stress. It will be interesting to evaluate whether current commercial cultivars, in addition to deep root systems, possess active shallow roots and whether this helps conserve soil moisture during the vegetative and reproductive stages.

This information will help investigate, understand, and optimize root plasticity in upland cotton and translate this knowledge to develop next-generation cotton cultivars with better RSA under water-deficit conditions to improve sustainable agriculture. Commercial producers can apply this research to their ongoing projects to improve the future elite varieties, enabling farmers to reduce their water consumption and save money while increasing their yields. 

At the end of December 2023, the top 0-12 cm depth of topsoil was moved from the Quaker Research Farm to the Horticulture Greenhouse so field soil could be used for the repeated aboveground experiment to get the root architecture traits for the selected accessions. Accessions were previously selected based on their differences in root phenotypic traits. At the end of the 2023 summer season, many of our replications were destroyed in a hailstorm, and the timing for extraction prevented getting good, replicated data to move forward. The bags were filled with field soil by measuring mass and volume to have as similar a density and compaction as manually allowed (Fig. 1 A&B). 12 Genotypes were used three elite cultivars, TM1, and Coker 312 were used for the reference genotypes, also 3 of the lower and 4 of the higher performing root traits. They were randomly placed for three replications (Fig. 1C).

A) Bag Set up Greenhouse Spring 2025B) Bag Set Up in Greenhouse Spring 2024

C)

Row 1 

Row 2 

Row 3 

189 

141 

184 

107 

70 

105 

99 

132 

197 

193 

66 

181 

70 

105 

197 

141 

66 

181 

189 

184 

132 

107 

99 

193 

105 

181 

107 

10 

70 

197 

141 

11 

66 

184 

193 

12 

99 

132 

189 

 

 

Figure 1: Bag Set up in the Horticultural Greenhouse with three bags across and 12 bags down.  A 5-gallon bucket was used to measure the mass of soil that went into each bag. And the height of the cylinder was measured by marking the inside of the bag. Similar soil density and compaction must be improved to get the same amount of soil mass within a given volume for each bag. B) Representing the three bags across and how the bags look when filled. C) Random design of 3 replications 

 Each replication was planted two days apart to allow time for extraction for 2 days for each replication. Trying to allow them to be at the same growth stage when extracted. Plants were grown for six weeks and were extracted before they started the flowering stage. Many varieties differ in their flowering stages, and it was intended to see their root traits before allocating resources to the reproductive system. For manually watering the plants, a graduated cylinder was used to give them their recommended inches once per week. Above-ground traits were documented at two growth points before extraction. Plant height and health, number of nodes, level of pubescence, anthocyanin present, and number of leaves are some documented traits (Fig. 2). 

 

Selection of Canopy Comparison
Figure 2: Shows the canopy growth for some of the plants. On the left are images taken on March 1, 2024, and on the right are images taken on March 16, 2024. The top row is one sample from accession TTU070. The middle row is a sample from TTU132, and the bottom row is a sample from TTU197.

When extracting the samples, some difficulties arose with the stickiness of the soil, combined with how delicate the roots were. This made collecting the complete root system intact impossible. Nearly all the roots broke off from the primary root. They were manually collected from each bag until no more roots could be found (Fig. 3). The roots were washed by placing them into a bin of water with a mesh over the top, allowing the fine soil to fall away from the root and keeping the roots moist until time for imaging. 

Initial Washing of the Greenhouse Tall Bags
Figure 3: The beginning of extracting the roots from the soil. Gathering and placing the roots onto the backdrop took time to image them for the software. Thirty-six bags were completed within one week. Thank you to the help from Patil Lab, who took the soil back outside while I worked on root extraction.

Each root piece was manually placed on a painted black metal sheet (Fig. 4). They were carefully placed so as not to touch each other to obtain the total root length, surface area, and average diameter. Unfortunately, the correct root angle, number of tips and forks, and tap root depth could not be obtained. 

RSA of four Samples
Figure 4: Extracted roots displayed on the background that can be used in the WinRHIZO software for color analysis to analyze the root length, surface area, and average diameter.

These accessions were intended to be placed in this year’s field experiment to see how their performance was under drought conditions; however, attempts to get enough seeds from some of the landrace varieties (105, 107, 099, 141) in time for this year’s field season were unsuccessful. Only a couple of seeds were collected because of their late flowering due to the assumption of photoperiod sensitivity traits. Therefore, accessions 105, 107, 099, and 141 were planted earlier than the previous year, hopeful that more seeds would be collected for future experimentation. 

Shortly after removing the tall bags from the greenhouse, another experiment was conducted to collect the RSA phenotypes for the GWAS analysis. Two hundred twelve accessions, 12 elite cultivars, 200 various accessions from obsolete cultivars, landraces, breeding lines, and TM1 and Coker 312, were selected to represent cotton germplasm from around the world and West Texas (Fig. 5). These accessions were selected based on their ability to flower in this area and their similarities and differences in visual phenotypes. 

A map of the germplasm's origins.
Figure 5: The countries of origin of the seeds ordered from the GRIN database and the companies that supplied the elite cultivar seeds.

The bags used were 2 gallons filled with the Quaker Research Farm soil. Three replications of each accession were randomly distributed throughout the four tables in the same greenhouse room. Each day, 64 bags were completed. Five seeds were placed in each of the 636 bags. After germination, the weakest seedlings were thinned out, leaving only one plant per bag. This experiment attempted to keep the conditions optimal for well-watered conditions for all accessions because there was not enough room or time for another 636 bags for drought conditions. Each bag was saturated with ~ 600 ml and left to dry for 24 hours to reach field capacity before planting the seeds. Each set of 64 bags was one day offset from the previous 64 bags to allow them to grow the same amount of time and be extracted on the same days after planting (DAP). Within 11 days, all bags were planted (Fig. 6A). Each bag received 200 ml of water on 3, 8, and 15 DAP as the greenhouse temperature would reach 110 degrees during the day. Water was added using a graduated cylinder to ensure the same amount was added to each bag. Plants were extracted on 17 DAP (Fig. 6B). Washing the roots was done by first soaking the bags, allowing them to become overly saturated to reduce the stickiness of the soil to the roots. The bags were cut and ripped away as one held onto the plant's stem. Then, working the way down the tap root, try to loosen as many roots as possible from the soil and gather any roots that broke off. Samples were stored in Ziplock baggies with 50% ethanol and placed into a 4° C room until they could be scanned (Fig. 6C). 

Set up, washing, and storage of the cotton.
Figure 6: A) The overview of all 636 bags. B) Washing was carried out; it took approximately 15 minutes per bag to collect all the roots, preferably all intact, but many fine roots broke off and were collected by hand. C) All roots and shoots were stored with their labels in 50% ethanol to prevent decomposition.

Samples were scanned the same way as the previous year's sample with the WinRHIZO scanner and software by placing an acrylic tray with some water to ease the roots' separation. Plastic tweezers were used to pry apart the roots to prevent all bunching up in one section. However, some roots were so prolific that they made separation challenging. Each sample took around 7 minutes to scan, and while one was being scanned, another sample was being prepped and separated.  

The same 212 accessions grown in the greenhouse were also planted in the Quaker Research Farm in a 16-row plot with sub-surface drip irrigation that I could manually turn on or off for eight rows at a time (Fig. 7). For the first eight rows, I left the irrigation off for the entire season and only received water from the rain. In the second 8 rows, I supplemented the weeks with no rain with irrigation. Aimed for about 1” of water per week. The irrigated plot received 5 weeks of irrigation more than the rainfed plants.

Field Overview Summer 2024
Figure 7: Overview of the field, several weeks after planting and just after manually weeding. Eight rows have 212 accessions and will be irrigated ~1 inch per week. Weeks that rained, irrigation was skipped. Another eight rows with the same 212 accessions were watered only by rainfall for the entire season.

Figure 8 shows the difference in canopy growth between the two plots in September. Irrigation was stopped at the beginning of September. A sensor from AquaSpy was placed in the rainfed and irrigated fields to take measurements until the end of the season. It measured every 15-minute interval with moisture levels, salinity, and temperature. Figure 8 focuses on the moisture levels collected by both sensors. The sensors were placed at the beginning of July; irrigation began four weeks after planting. 

Treatment Comparison
Figure 8: A) Overview of the rainfed field. B) Overview of the irrigated field. C) Soil moisture measurements over time at varying soil depths from AquaSpy sensor technology.

By the peaks at 4 inches depth, it is possible to determine when it rained during the months of July through October. There was no rain during the hottest months of the year, from August through September, contributing to the drought stress for the cotton during this time. The difference in water between the treatments through August and September is a clear separation at most depths. Looking at the 8 to 12-inch depth where the drip irrigation emitters are located, there are two steady points in the month of July. But a slow decrease in August, which could suggest that the 1 inch of irrigation was not enough to bring it back up to saturation levels, but most likely helped to prevent the plants from being too stressed during August. A repeated experiment next summer, with time to do physiological measurements, would be helpful to know the stress of the leaves between the two treatments.

Leaf samples were collected from all the accessions and sent to a company, SGS Technologies, to do the genotype sequencing by running an 18K SNParray, we received SNPs for around 17K markers with around 10K that were homologous with unique markers and the other markers were either heterozygous and/or the marker matched multiple locations in the genome—using R version 4.4.2 and RStudio package called GAPIT (Wang and Zhang, 2021).

 

Research results and discussion:

The field experiment included a field of 550 Upland Cotton varieties, for which leaf samples were collected to genotype the germplasm for eventual association mapping. The genotyping is planned to take place later in 2024. There were 24 genotypes selected for the above ground bag experiment. Only 2 replications could be completed with the lateness of the season, getting them planted, and the time and labor it took to fill the bags. The images of the roots were uploaded into the WinRHIZO software for the analysis of the traits with the ruler sticks set for image calibration. The results from the end of October, with the extraction of the bags, show some differences in the genotypes and amount of root length (Fig. 9). However, none can be statistically significant with the limited replications and missing data, making it impossible to validate these results statistically. A repeated experiment will need to be done. Only two replications could be done with the time it took to fill and the summer heat. The two replications also vary greatly because of the time difference in extracting the roots. Also, 12 samples were damaged in the hailstorm, which could account for some variation. The plants were also past their boll development in many of the plants because of the time it took to extract them. This could also have had an impact on the roots that were present. A lot was learned from this experiment; hopefully, the repeated experiment will not repeat the same mistakes.

Field 2023 Total Root Length Results
Figure 9: The two replications' total root length was averaged to show a quick comparison.

The next experiment was done in the greenhouse over the winter break and beginning of January 2024. The same bags were used with fewer accessions, but with three replications instead of two. The results in Figure 10 contain the total root length and the surface area of the shoots for each of the 12 accessions. The three elite cultivars have yellow stars to distinguish them from the other varieties. TTU193 is the Texas Marker 1 (TM1), and TTU197 is the Coker 312. These findings follow a similar pattern to the previous year’s greenhouse experiment in that the elite cultivars seemed to do poorer than TM1 and Coker312, and the landraces (132,099, 105, and 107) are continuing to outperform in root length compared to the elite cultivars and TM1 and Coker312. A couple of surprises for this are that TTU141 and TTU070 were expected to be lower compared to other accessions when looking at the previous year’s greenhouse data. However, it was observed that these accessions withstood some of the greenhouse’s insect pressures better than the elite cultivars did and could have had an advantage in their growth. Some of these landraces' more prolific root growth could potentially contribute to drought tolerance when placed in field conditions. To test this hypothesis, these accessions went on to objective two where field testing with drought stress was added. 

Greenhouse February 2024 RSA Results
Figure 10: The bars in tan represent the Total Root Length results for each sample. The green bars represent the shoot’s total surface area for each accession. The yellow stars signify the elite cultivars.

All accessions' total root length (TRL) shows a normally distributed pattern with a mean of around 570 cm of root length for 17 DAP. The range varies from 150 to around 900 cm, suggesting a large variation in the germplasm for total root length (Fig. 9). The colors represent the country of origin from which the accession seed was collected. Each color seems randomly spread out between high, middle, and low; suggesting that not just one origin has the highest root length.  

The total root length of 636 samples labeled by origin
Figure 11: Total root length for each accession shows some with considerable variation within a genotype and others with slight variation within a genotype. Overall, there is a wide variation between the genotypes.

Looking at the correlation between the RSA traits, most of them are positively correlated. The exceptions include the root to shoot ratio with several traits that can be negatively correlated and average diameter with the root length and average diameter with the root to shoot ratio with no correlation (Fig. 12). This is most likely because of the age of the seedlings being only two weeks old, most of the diameter is the same regardless of length or how tall the plant is. A good overall picture of each accession’s RSA can be determined by comparing these three traits.

Correlation Matrix of the RSA traits from May 2024
Figure 12: Correlation matrix of the RSA traits of the 212 accessions. Direction of the diagonal and color describe a positive or negative correlation between the traits. If it is purple, it is negatively correlated. If it is green, then it has a positive correlation. The shading of the color and thickness of the oval indicates how strong the relationship is between the two traits. The traits are as follows:
shootsurfarea = Shoot Surface Area (cm2)
shootprojarea = Shoot Projected Area (cm2)
shootbiomass = Shoot Dry Biomass (g)
totalbiomass = Total Dry Biomass (g) (root + shoot)
rootbiomass = Root Dry Biomass (g)
totalrootlen = Total Root Length (cm)
rootprojarea = Root Projected Area (cm2)
rootsurfarea = Root Surface Area (cm2)
avgdiam = Average Diameter (cm)
roottoshootratio = Root to Shoot Ratio (root biomass / shoot biomass)

When comparing the accessions with the most and least total root length, they may or may not also be found at the most and least root to shoot ratio or average diameter (Fig. 13). Accessions may have different strategies for drought stress response, by only looking at accessions really good at prolific root length other traits may be missed. Selecting some of the best and worst from each trait may be good to explore in additional experimentation with drought stress. Some elite varieties are found in the bottom ten of these traits. In the case of the average diameter the smallest ten is accession 80 which is also the number one spot for root length. Will 80 do well in drought conditions? To find out the same 212 accessions will be evaluated in objective 2.

Histograms of the three RSA traits that appear to have importance
Figure 13: Histograms of the total root length (a,b), root to shoot ratio (c,d), and average diameter (e,f) to show the distribution of the germplasm. Below the histograms are the top and bottom ten accessions showing contrasting phenotypes. The point is the mean of the three replications for those accessions.

Some of the 12 accessions with more than 10 plants were selected to be pulled for just the crown of root systems to see if there was much difference in their root systems from rainfed to irrigated, and also to compare with each other. One sample that showed potential drought tolerance had an improved root system in the rainfed condition compared to the irrigated condition, while two accessions that showed drought tolerance had a reduced root length but a similar or increased root diameter. This could suggest that their drought tolerance came from above canopy physiological changes and/or the change in the diameter of the root could have helped with the drought stress response (Fig. 14). These three seemingly different strategies may support the finding that there are multiple strategies of drought stress tolerance and varying ways of changing the physiology and anatomy in response to that stress. Could there be a way to combine some of those varying strategies to improve or optimize the drought tolerance? 

Two-way Anova for the 7 accessions uprooted Field 2024
Figure 14: A) Accession TTU041 shows a more robust root system in the rainfed condition, a similar height, and a greater yield than its irrigated counterpart. B) Accession TTU184 appears to have a reduced root system in the rainfed condition when it comes to surface area, but the average diameter of the root is increased. The performance of the plant height was the same, but the number of bolls was increased in the rainfed condition. C) Accession TTU640 had a significant decrease in root surface area but a similar reduction in diameter and plant height, but yield performed about the same.

An ANOVA was performed to see if the treatments were significant for root traits, plant height, and number of bolls (Fig 15). 

One-way ANOVA for differences in treatment with RSA, growth, and yield traits.
Figure 15: One-Way ANOVA (α = 0.05) for the treatment effect on the following RSA and field phenotypes: Analyzed Region Area (cm2) (a), Total Root Length (cm) (b), Root Surface Area (cm2) (c), Average Diameter (cm) (d), Plant Height (cm) (e), Number of Bolls (f). Seven accessions and three replications were included in this comparison.

To get an overall look at the differences in growth and yield, a principal component analysis was conducted to see if any clusters or samples showed interesting patterns. Starting from the middle of the plot, as one goes to the left, the plant height decreases. As one goes to the right of the plot, the number of bolls, plant height, and weight of the plant increase. Closer to the bottom of the PCA grid is an increased boll weight (Fig. 16). There is a clear distinction between the irrigated (green) samples and the rainfed (purple) samples. Nearly cut vertically through the center. However, some purple samples perform just as well as the higher functioning irrigated samples. And some green samples perform low, close to the middle of the drought-stressed plants. If plants undergoing intensive drought stress can still have the same growth and yield as other irrigated plants, that could suggest having some drought stress tolerance.

PCA results from Field summer 2024
Figure 16: Principal Component Analysis (PCA) of plant height(cm) and number of bolls. The length of the arrow explains how strongly that trait represents the distribution. PC 1 (Number of Bolls) explains almost 74.67% of the variation, and PC 2 (Plant Height) explains 16.22% of the variation. Those points that cluster closer together indicate their similarities, while those farther away have more differences from the rest. Each green dot represents one accession in the irrigated condition. Each purple dot represents one accession in the rainfed condition.

To get an idea of some of the accessions that performed similarly in both the irrigated and drought conditions, the means of the boll number and plant height were compared, and those with the greatest slope showed the greatest difference. Negative slopes suggest that drought hindered yield, while positive slope suggests that the accessions performed better in the drought condition compared to the irrigated condition. A small or near horizontal line suggests that they performed nearly the same in both conditions (Fig. 17). 

 

Accessions that may be drought tolerant vs those that may be drought sensitive.
Figure 17: Each accession was averaged for plant height and number of bolls and grouped by their treatment. Then, they were plotted on the graph to see if their performance increased, decreased, or stayed the same. With all 212 accessions on a plot, it was impossible to see any distinct patterns. To narrow it down, the 15 accessions with the most significant difference in means are on the top row (a,b), and the 15 accessions with the least difference in means between the two treatments are on the bottom row (c,d).

These results together can work to select high and low performance with RSA traits and field performance under drought stress conditions. This selection of contrasting lines can then continue to the next level of experimentation with increased replications, added drought stress in the greenhouse, and another year of field experimentation. In addition to all the previous parameters, physiological data will be critical to measure in the upcoming experimentation. Physiological measurements like stomatal conductance, chlorophyll levels, and relative water content of the leaves give a greater understanding of the condition of the plant and how the physiology plays a role in drought tolerance, along with what the accession maintains.  

Objective 3: Analysis of phenotypes and genotypes through genome-wide association mapping 

By collecting leaf samples of all the accessions and sending them to a company, SGS Technologies, to do the genotype sequencing by running an 18K SNParray, we received SNPs for around 17K markers with around 10K that were homologous with unique markers and the other markers were either heterozygous and/or the marker matched multiple locations in the genome. Using R version 4.4.2 and RStudio package called GAPIT (Wang and Zhang, 2021). Several traits had markers that showed significant QTL peaks (Fig. 19).  

Phenogenetic tree, PCA, and Population Structure of the SNP data.
Figure 19: A) Distance-based neighbor-joining (NJ) tree from the vcf2poptree software. Showing 280 Upland cotton accessions from the 17,500 SNPs that were identified using the cottonSNP18K array. The coloring scheme is based on the origin in the GRIN database of where the seed was collected. Those labeled had contrasting performances in their RSA traits. B) PCA of the SNP data using PLINK analysis. C) This population structure used the ADMIXTURE software to determine the K clusters, and the elbow package in R calculated the optimal number of clusters for plotting. The colored clusters correspond to a 70% membership threshold.
Manhattan Plots Genome View of multiple traits.
Figure 20: Manhattan plots for the traits that showed significant markers: Total Root Length, Root to Shoot Ratio, and the Number of Bolls for rainfed and irrigated conditions.

Because this data is very recent, more time is needed to analyze and process it all before conclusions can be drawn about these potential QTLs.   

Participation Summary

Educational & Outreach Activities

3 Webinars / talks / presentations
1 Workshop field days

Participation Summary:

100 Farmers participated
300 Ag professionals participated
Education/outreach description:

Presented my research at the 10th Annual Water College in January 2024. Farmers, sponsors, and educators were present at the event of over 200 hundred people. 

Lamb, M. L. (2024, January 24). Optimizing cotton's root system to improve response to water-deficit stress [PowerPoint slides]. Institute of Genomics for Crop Abiotic Stress Tolerance, Texas Tech University. https://texastechuniversity-my.sharepoint.com/:p:/g/personal/miclamb_ttu_edu/Ede1MWpL2SxNu8MjQb-bWEYBcdg5L-_jzAUcg3wI-n_O8Q?e=dWC8qk&nav=eyJzSWQiOjI1NiwiY0lkIjoxMDk4NTcyMjJ9

Presented at the Agricultural Water Sustainability Summit, Lubbock, Texas, on August 9, 2024, a few hundred farmers, ranchers, and agricultural professionals attended the conference. My presentation was on the last day, towards the end of the conference, so only 12 people were present for my presentation. 

I presented at the Beltwide National Cotton Conference on January 16, 2025, at the Center New Orleans, Louisianna. Although thousands of farmers/ranchers and agricultural professionals participated, around 50 were in attendance for my presentation. 

Presented my field research at the Quaker Farm Field Day August 2024 with around 25 farmers, or Agricultural professionals 

Project Outcomes

Project outcomes:

Suppose we can figure out more about the roots, how they are related to shoot traits, and the genetic relationships with the roots. In that case, It will save future generations time, labor, and money to develop their crops with a better understanding of the Root System Architecture. Because it is tough to view the roots to know what traits are being bred for and what root traits will be beneficial or not in different abiotic stresses, this project's goal is to get additional understanding of methods of how to assess the RSA and also do it in an economically inexpensive way to better our understanding of the genetic mechanisms of the RSA.

The data suggest that the cotton germplasm has a significant variation in RSA traits, and various accessions are more tolerant of drought than others. There could be strategies that include altering the root system during drought stress to help the plant endure the periods of drought without suffering as much yield loss. By comparing the results of this study, we can select accessions that may have desirable RSA traits that can be further studied and understood to be applied to future studies and breeding programs. Further analysis of the GWAS results to further the understanding of potential genes that play a role in those traits will also help determine how to optimize the traits for drought tolerance. Accession TTU 132, TTU 184, TTU 80, and a few others that show promise of drought tolerance and those that show drought sensitivity can be narrowed down into more focused experimentation that can lead to even greater understanding of RSA's role in drought tolerance for cotton.  

Knowledge Gained:

We are still learning as we continue these few months since September 2023. I have been learning how much work is required to learn something unknown and also how much knowledge that comes from so many other people before is needed to build upon each other. I am excited for the next year of learning and hope to come up with new information to better serve society and agricultural partners. 

Recommendations:

It is recommended that additional greenhouse testing with selected accessions do additional drought testing and repeated year experiments to further validate results. It is also recommended that additional investigation into the GWAS results and further study of any QTLs and candidate genes to implement into experiments for the following year. 

Any opinions, findings, conclusions, or recommendations expressed in this publication are those of the author(s) and should not be construed to represent any official USDA or U.S. Government determination or policy.