Genome-wide identification of differentially expressed genes under water deficit stress in upland cotton (Gossypium hirsutum L.)
© Park et al.; licensee BioMed Central Ltd. 2012
Received: 23 September 2011
Accepted: 10 May 2012
Published: 15 June 2012
Cotton is the world’s primary fiber crop and is a major agricultural commodity in over 30 countries. Like many other global commodities, sustainable cotton production is challenged by restricted natural resources. In response to the anticipated increase of agricultural water demand, a major research direction involves developing crops that use less water or that use water more efficiently. In this study, our objective was to identify differentially expressed genes in response to water deficit stress in cotton. A global expression analysis using cDNA-Amplified Fragment Length Polymorphism was conducted to compare root and leaf gene expression profiles from a putative drought resistant cotton cultivar grown under water deficit stressed and well watered field conditions.
We identified a total of 519 differentially expressed transcript derived fragments. Of these, 147 transcript derived fragment sequences were functionally annotated according to their gene ontology. Nearly 70 percent of transcript derived fragments belonged to four major categories: 1) unclassified, 2) stress/defense, 3) metabolism, and 4) gene regulation. We found heat shock protein-related and reactive oxygen species-related transcript derived fragments to be among the major parts of functional pathways induced by water deficit stress. Also, twelve novel transcripts were identified as both water deficit responsive and cotton specific. A subset of differentially expressed transcript derived fragments was verified using reverse transcription-polymerase chain reaction. Differential expression analysis also identified five pairs of duplicated transcript derived fragments in which four pairs responded differentially between each of their two homologues under water deficit stress.
In this study, we detected differentially expressed transcript derived fragments from water deficit stressed root and leaf tissues in tetraploid cotton and provided their gene ontology, functional/biological distribution, and possible roles of gene duplication. This discovery demonstrates complex mechanisms involved with polyploid cotton’s transcriptome response to naturally occurring field water deficit stress. The genes identified in this study will provide candidate targets to manipulate the water use characteristics of cotton at the molecular level.
Predicted world population growth will require more food, feed, and fiber production. Intermittent drought and shortages of water supply negatively impact crop productivity and are predicted to occur more frequently in future agricultural systems [1–3]. In future agricultural systems, an increase in crop productivity must also be accompanied by a reduced environmental impact of production on natural resources . Therefore, to increase productivity while maintaining adequate water resources, it is necessary to understand the mechanisms plants use to cope with water deficit stress . There have been two general efforts towards more sustainable and water use efficient agricultural production. One effort is represented by breeding and genetic engineering [1, 6, 7]. Another effort is focused on the improvement of agricultural practices such as tillage and irrigation systems [8–10]. To successfully pursue the breeding and genetic engineering approach, it is fundamental to understand plant responses to drought stress at the molecular, cellular, and genetics levels. This understanding is critical, because drought tolerance is a quantitative trait influenced by a combination of regulatory pathways [5, 11].
Drought stress impacts a network of plant gene expression mechanisms that lead to the reprogramming of a variety of physiological and metabolic processes in accordance with stress response. Early studies primarily used model plant species to identify a wide spectrum of genes that are involved in different levels of metabolism, signal transduction, osmotic regulation, stress response, and gene regulation [12–14]. More recently, many plant gene expression profiling experiments, using a range of technical approaches, have been conducted on economically important crop plants to determine stress related expression signatures. For example, next generation sequencing technology (NGS) has been applied to provide global gene expression profiles of whole transcriptomes [15, 16]. NGS expression studies are greatly facilitated by the availability of annotated genome sequence (reference or otherwise).
As an alternative to NGS, the cDNA-AFLP technique is a widely adopted gene expression profiling method. This is especially true for crop plants, such as cotton, that are awaiting whole genome- or whole transcriptome- sequence data. cDNA-AFLP profiling has proven to be a robust method to discover differentially expressed genes in a number of plant species [17–19]. With the advantage of no prerequisite reference genome sequence requirement, cDNA-AFLP has been used to understand regulation mechanisms under various biotic and abiotic stress conditions that influence patterns of global gene expression in diverse crop species [20–22].
Cotton, a warm climate crop of tremendous global economic importance, is known to contain intrinsic mechanisms that allow it to partially withstand water deficit- and heat-stresses by its deep and extensive root system and adaptive osmoregulation mechanisms [23, 24]. However, as water deficit stress progressively continues, development and reproduction of cotton is severely affected, and the quality and yield of cotton fiber production is significantly reduced [25, 26]. Therefore, for successful cotton fiber production, it is essential to understand proper management strategies between water demand, irrigation, and plant responses that are seasonally variable .
In cotton, there are several reports on the molecular regulatory mechanisms that orchestrate drought resistance using quantitative trait loci (QTL) mapping or by developing drought stress induced cDNA libraries [28–30]. In addition, microarray gene expression profiling has been used in cotton tissues exposed to variable water deficit stress in a controlled environment growth condition . However, microarray gene expression profiling is limited to ESTs or genes present on the microarray chip and may not provide a complete representation of differentially expressed genes.
Although genome sequencing data of many economically important crop species are widely available as full or high quality draft formats , the progress of cotton genome sequencing has been slow. Genome sequencing of an ancestral diploid genome (D) is near completion [33, 34]. Once complete, the D-genome sequence will be useful as a platform for more comprehensive information on cultivated cotton’s tetraploid genome .
In this study, our objective was to gain comprehensive insight into the water deficit stress-related gene expression profile in cotton grown under field conditions. To accomplish this, we used cDNA-AFLP to identify differentially expressed genes in the roots and leaves of a single cotton cultivar grown in a rain-fed, drought-prone field environment. From this large comparative analysis, fundamental experimental data illustrate complex mechanisms of gene expression in response to water deficit stress. This study should facilitate future efforts to improve the agricultural productivity of cotton and other crops under soil water deficit conditions.
Collection of cDNA-AFLP data sets
Since the cDNA-AFLP procedure uses a two selective nucleotide extension at the 3’ terminus of each primer (See methods), there are 256 available primer combinations. In total, 64 primer combination covers one-quarter of the transcriptome of interest in tissues tested. Although the total gene number of the cotton transcriptome is not clear, it is estimated that the tetraploid cotton genome contains over 70,000 genes . Therefore, we estimate that our expression data represents approximately 17,500 genes.
Sequence analysis and functional annotations
Differentially expressed genes under water deficit stress
Tissue specificity of annotated proteins is described in Additional file 1. Annotated proteins in the stress/defense category covered 18.4% (27 TDFs) of the total TDFs. Ten of these 27 TDFs coded for proteins related to the production of reactive oxygen species (ROS). Examples of ROS-related proteins included glutathione S-transferase (GST), manganase superoxide dismutase (Mn-SOD)-like protein, glutathione peroxidase, mono-dehydroascorbate reductase, phosphoadenosine phosphosulfate reductase (Thioredoxin), aldo/keto reductase, 5’-adenylylsulfate reductase (Glutathione), and superoxide dismutase (SOD) copper chaperone. Six of the 27 TDFs were annotated as heat shock proteins (HSP). Examples of annotated HSP included chloroplast HSP70.1-like proteins, HSP70, DnaJ HSP-like protein, and alpha-crystalline HSP. Also, of the 27 TDFs, two were annotated as alcohol dehydrogenase (ADH) enzymes, two as drought-induced proteins, and seven TDFs with relationships to other stress/defense responses.
Among the 17 gene-regulation categorized TDFs, eight coded for transcription factors and seven for translation-related factors. The signal transduction category (14 TDFs, 9.5%) consisted of four kinase, four phosphatase, and four GTP binding factors. Two TDFs coding for aquaporin proteins (PIP1;3 and PIP1;12) were identified in leaf and root tissues. Five cell-wall biosynthesis-related genes were also identified (Additional file 1).
Figure 3 shows a comparison of up- (Figure 3A) and down-regulated (Figure 3B) TDFs in root and leaf tissues. Compared to leaf tissue, up-regulated TDFs were found more abundantly in the root across the five primary functional categories- unclassified, metabolism, stress/defense, signal transduction, and gene regulation. Down-regulation was more prominent in leaf tissue for stress/defense, signal transduction, and gene regulation. Down-regulation was more prominent in root tissue for metabolism, unclassified, and transport categories.
Validation of expression pattern by RT-PCR
Summary of RT-PCR reactions performed for selected TDFs a
glycosyl transferase family 1 protein
nucleic acid binding protein/cold shock induced
chloroplast heat shock protein 70.1
up L up R
manganese superoxide dismutase-like protein
universal stress protein
phosphoadenosine phosphosulfate reductase family protein
beta-ketoacyl-ACP synthase II
MADS-box protein (AGL84)
translation initiation factor-like protein
RAN2; GTP binding
oligo Peptide transporter, putative
NITRATE TRANSPORTER1/PEPTIDE TRANSPORTER (POT/NAXT) family
Photosystem I reaction center subunit VI
Cytochrome b6-f complex
Gene duplication and water deficit stress
Homologous TDFs with high sequence similarity
Similarity to AD genomes§
Among them, PIP1-related TDFs (99% identity with each other) belonged to a multi-gene aquaporin water channel membrane protein family in cotton . These TDFs showed very high sequence similarity to GhPIP1;3, GhPIP1;4, GhPIP1;12 and GhPIP1;13 (E value < 9E-174 for 06 C03; < 6E-123 for 06E09 as shown by the multiple sequence alignment (Additional file 2)). According to the prediction of functional residues among aquaporin water channel proteins [40, 41], two amino acid substitutions (both Val-Ile conversions) between two TDFs do not affect functional properties. Moreover, these two non-synonymous substitutions belonged to the 4th conserved transmembrane domain, (V/I)GTF(V/I), that is found frequently in cotton EST sequences encoding the aquaporin PIP1 subgroup [39, 42]. This points to the importance of either the Val or Ile at these two sites. Our differential expression revealed that an aquaporin TDF, 06 C03 is up-regulated in root by water deficit stress while a homologous TDF, 06E09 is up-regulated in leaf tissue and down-regulated in root tissue. This finding indicates the possible functional, tissue-specific diversification of duplicated PIP1 genes under water deficit stress. Duplicated TDFs from three other homologous pairs also showed differential expression patterns upon water deficit stress, whereas duplicated TDFs of M6PI showed the same pattern of expression (Table 2). These data indicate the possibility that products of duplicated genes might play roles in a tissue dependent manner and that regulatory elements controlling the expression of duplicated genes evolved along with polyploidization.
Molecular framework under water deficit
Considering that drought-resistant cultivars are not readily available for commercial cotton production, more effort is required to identify or generate traits to withstand soil water deficit stress. In cotton production systems, rain-free periods occurring during reproductive growth are especially damaging. In this study, our objective was to understand a molecular architecture of water deficit stress at the level of gene expression by analyzing differentially expressed transcripts in the field during reproductive growth. Therefore, differentially expressed TDFs identified here can be considered genes functionally active (in case of up-regulation) or inactive (in case of down-regulation) at the stage of adaptation to naturally occurring water deficit stress. Our study revealed more than 500 transcripts with altered gene expression levels. The number of down-regulated genes was 1.5X higher than the number of up-regulated genes. Down-regulated genes occurred more frequently in root tissues. The majority of TDFs with altered expression belonged to functional categories including metabolism, stress/defense, signal transduction, and gene regulation as well as the category unclassified. This result corresponds well with previous findings of other species demonstrating transcriptome alteration under water deficit stress [12, 13].
In terms of cotton, a recent microarray expression profiling experiment conducted with plants in a greenhouse reported a total of 2,106 stress responsive-transcripts . The majority of transcripts showed tissue specific expression and a higher number of stress responsive-transcripts were identified in leaf compared to root. In comparison, our study with field grown plants identified fewer stress responsive-transcripts (519), of which, 215 were leaf specific, 293 were root specific, and 11 were expressed in both leaf and root. The candidate gene targets identified in both studies require future work to determine their potential application to improve cotton water use.
Our expression data suggests the involvement of the ROS related defensive pathway since we identified 10 TDFs related to anti-oxidation mechanisms. The increased level of ROS, driven by water deficit stress, can affect cellular components that are oxidized partially or severely [43, 44]. Therefore, it is critical that plants protect themselves from harmful oxidations with detoxifying mechanisms by using antioxidants and scavenging agents . DST (Drought and Salt Tolerance), a previously unknown zinc finger protein, was found as a negative regulator of drought and salt stress by repressing H2O2 accumulation and stomatal closure in an abscisic acid (ABA)-independent manner . In addition, recent studies on developmental and stress-induced cellular processes suggest that ROS and callose deposition are co-regulated, thereby controlling the cell wall matrix adjacent to the plasmodesmata for intercellular redox signal transduction . In our study, one TDF (17-1A09) encoding a callose synthase-like family protein was isolated as a down-regulated gene and many ROS-related TDFs showed up regulation (seven TDFs) or down-regulation (three TDFs). It is of interest to note that this plasmodesmata-related co-regulation of ROS and callose represents a cellular mechanism in the cotton fiber elongation process leading to water uptake following decreased osmotic potential in elongating fiber cells . These findings highlight a possible interconnection between ROS and callose deposition in the areas adjacent to plasmodesmata that cotton employs in response to water deficit stress.
Previous studies have shown that phosphoethanolamine N-methyltransferase, vacuolar invertase, and aldo/keto reductase (as identified in this study as TDFs 07B12, A05D02, and A04A09, respectively) are involved in water deficit stress-related defense mechanisms [48–50]. When silenced, the decrease of phosphoethanolamine N-methyltransferase produced not only multiple growth defects but also temperature sensitive male sterility in Arabidopsis. The vacuolar invertase was shown to be related to water deficit stress in maize . Aldo/keto reductase catalyzes the detoxification reaction of reactive aldehyde groups generated by abiotic stresses thereby providing plants with stress tolerance .
It is well known in many plant species that ABA acts as a key hormone in the abiotic plant stress response [51, 52]. Upon water deficit stress, stomata in leaves are closed and prevent water loss through transpiration. This process is believed to be regulated by ABA [53, 54]. Interestingly, we did not identify any gene related to biosynthesis or action mechanisms of ABA nor genes involved in the regulation of stomata . However, besides the function of ABA in stomatal closure, it was recently shown in Arabidopsis and tobacco that interactions possibly occur between water deficit-responsive proteins and HSP chaperones [56, 57]. Under water deficit conditions, both HSP70 and HSP90 chaperones were recruited to control stomatal closure, thereby serving as machinery important for stomatal gating. The finding of HSP-coding TDFs also suggests the existence of the HSP related machinery in water deficit stress in cotton. In our study, five HSP-coding TDFs were up-regulated while only one showed down-regulation. Previously, a drought-related alpha-crystalline HSP was identified by differential screening from 10-d drought-stressed G. arboreum cotton seedlings and was effective in providing tetraploid cotton plants with drought stress tolerance when over-expressed [58, 59]. The potential involvement of HSP-related drought tolerance was also suggested recently in a microarray transcriptome analysis in drought-stressed cotton plants .
In accordance with transpiration, molecular mechanisms underlying water uptake and transport are pivotal throughout the plant body. There is accumulating evidence that addresses the relationship between aquaporins and transport of water in various physiological conditions including water deficit stress [60–64]. Aquaporins are frequently identified as water deficit stress responsive genes in diverse plant species [65–67]. In roots of drought stressed chickpea plants, members of the aquaporin gene family appeared up- and down-regulated suggesting that a complex regulation of water status governs plant growth and development through aquaporin activity under water deficit . The identification of two homologous TDFs (06 C03 and 06E09) with similarity to members of the aquaporin water-channel protein family in this study (Table 2 and Additional file 1) indicate the possible involvement of cotton aquaporins under water deficit stressed leaf and root tissues. Recently, the identification of the large cotton aquaporin family illustrated their diverse function . In addition to their fundamental role in intercellular transport of water molecules across the plant body, reports have shown the significance of aquaporins in facilitating leaf CO2 conductivity relevant to plant photosynthetic capacity [69–71]. However, the function of aquaporins in water deficit stress tolerance remains unclear as aquaporin genes have not been identified in water deficit stress QTL studies to date . Therefore, more supportive and quantitative data need to follow.
After sequencing 366 excised differentially expressed fragments, we functionally annotated 147 TDFs following sequence verification. Of the 147, we described a number of water deficit stress-responsive genes functionally relevant to metabolism, signal transduction, gene regulation, and stress/defense mechanisms in root and leaf tissues (Figure 3). In addition, homology searches using BlastX did not classify 24.5% of TDFs due to lack of sequence similarity to known proteins. Among those, twenty-seven unclassified TDFs were differentially expressed in root tissue. The abundance of unclassified TDFs identified in this study provide additional transcriptome coverage not represented in EST populations commonly used in microarray experiments. It was reported that unknown genes such as proteins with obscure functions (POF) cover more than 20% of each new genome sequenced with many being species specific .
In maize seedlings with water shortage, 5 – 11% of genes were differentially expressed across an array of genetically diverse inbred lines. Also, while many of the genes were not repeatedly identified in different maize lines, more than 40% of the cellular pathways were shared across all the lines examined . In our study, a similar percentage of genes (13%) showed expression level changes upon water deficit stress. Since the genotype in this study is believed to be drought tolerant, it would be interesting to determine if biological pathways highlighted in this study (for example, ROS-, or HSP-related defense mechanisms) would appear in common across an array of cotton genotypes. Transcripts involving the HSP-containing functional group were also identified in a microarray based, water deficit stress response gene expression study using another cotton cultivar, FiberMax 989 . The few water deficit stress-related cotton genes identified in the current study could be used in candidate gene-focused association mapping approaches to identify QTL under drought stress. This approach was previously shown to identify SNPs associated with ABA and sugar levels under water deficit .
Implications of gene duplication in cotton abiotic stress response
As evidenced by recent studies, it is becoming clearer that gene duplication has contributed to adaptive evolution and plant diversification that can lead to evolutionary novelty [76–78]. In our study, four pairs of duplicated genes including two PIP1 homologues were regulated spatially under water deficit stress (Table 2). The differential expression of two duplicated PIP1s did not appear to result from differences of deduced amino acid sequences (Additional file 2). Hence, it is of interest to elucidate factors that lead to the differential response of duplicated genes against stress. These responses can be controlled genetically, epigenetically, or by gene specific cis- or trans-elements . Recent studies on genetic and epigenetic responses of plant genomes under environmental stimuli suggest that novel stress-induced genotypes such as methylation-sensitive polymorphisms can contribute to crop diversity that may lead to improvements in productivity . In rice, drought stress-induced changes in genome-wide methylation patterns were evaluated using a methylation sensitive cDNA-AFLP method. This revealed a 12.1% site-specific methylation difference between a drought tolerant line (DK151) and its sensitive progenitor line (IR64) when stressed and showed methylation status was regulated tissue specifically . Thus, water deficit responsive molecular changes presented in our study implicate important mechanisms of each gene or coordinated regulation of genes that should be targets of future study.
In response to field applied water deficit stress, we identified 519 differentially expressed TDFs. This global expression analysis generated sequence information for 147 TDFs and provided gene ontology, distribution of functional categories, and homologous genes. In accordance with data provided here, we highlight possible mechanisms by which cotton responds to water deficit stress. The combination of molecular changes in gene regulation and signal transduction results in increased activity of ROS-related defense mechanisms and HSP-driven protection machinery to provide cellular redox homeostasis and stabilization of functionally/structurally important proteins. A readily available cotton genome sequence will enhance our ability to understand the relationship between gene duplication/polyploidy and the functional, molecular adjustments cotton makes in response to water deficit stress and other environmental stresses.
During 2009, two 8-row plots (6 m row length and 96.5 cm row spacing) of cv. Siokra L-23 were grown at the North Carolina State University Sandhills Research Station near Jackson Springs, NC, USA on a uniformly deep, Candor sand soil with very low water holding capacity. Plots were subjected to irrigated and non-irrigated conditions, respectively, for 4 weeks during reproductive growth. Irrigation was applied supplemental to recorded rainfall weekly using above-ground drip irrigation as described by Campbell and Bauer . Siokra L-23 is generally known to be tolerant to water deficit stress in molecular and physiological levels [83, 84]. All plant measurements and tissue samples were obtained mid-day during the flowering period. On each sample date, to confirm differential plant water status, the water potential of uppermost fully expanded leaves was measured on three representative plants from each plot with a pressure bomb (Model 600, PMS Instrument Company, Albany, OR). From the same plants, leaf and root tissues were harvested and immediately frozen in liquid nitrogen for RNA isolation.
RNA isolation, cDNA synthesis, and cDNA-AFLP
Tissues collected from a single sample date with the greatest water potential difference between well-watered and water deficit stressed plants were used for all nucleic acid analyses. Average water potential of leaf tissue from well-watered and water deficit stressed plants was −1.47 ± 0.13 MPa and −2.58 ± 0.34 MPa, respectively. Total RNA was isolated using the XT buffer system with the addition of chloroform/iso-amyl alcohol extraction and LiCl precipitation steps . Isolation of total RNA and mRNA purification is described previously .
A cDNA template was synthesized as follows. First, poly(A) RNA was purified from 200 μg of total RNA using a microPoly(A) Purist Kit (Ambion). After the treatment of DNase, 1 μg of purified mRNA was used to generate first strand cDNA following the manufacturer’s protocol (SuperScript III 1st strand cDNA synthesis kit, Invitrogen) with no RNase H treatment. Double-stranded cDNA was synthesized using an E.coli DNA polymerase I (NEB) with RNase H (Invitrogen) at 16°C for 1 hour followed by an additional reaction for 1 hour at 22°C with T4 DNA ligase (NEB) treatment. After inactivation and clean-up, double-stranded cDNA template was digested with TaqI and MseI restriction endonucleases sequentially, ligated with two adaptors, and amplified with pre-amplification primer mix (MseI primer: 5’-GATGAGTCCTGAGTAA-3’, TaqI primer: 5’-GTAGACTGCGTACCGA-3’). After dilution (1:150) of pre-amplified DNA samples, selective amplification steps were performed with 64 combinations of eight MseI selective primers (5’-GATGAGTCCTGAGTAANN’-3’) and eight IRDye 700-labeled TaqI selective primers (5’-/5IRD700/GTAGACTGCGTACCGANN’-3’) for all four templates as indicated in the LI-COR cDNA-AFLP protocol. More detailed cDNA-AFLP procedures are described elsewhere [19, 86]. Two separate denaturing polyacryamide gel electrophoresis (PAGE) gels were analyzed for each sample using a LI-COR 4300 DNA Analyzer; one for the preliminary differential expression profiling and the other for gel scanning and band excision. Gel scan images were produced using a LI-COR Odyssey Infrared Imaging System that can detect amplified fragments labeled with IRDye 700. Band excision and recovery of fragments were performed as recommended in the LI-COR manual. To verify technical reproducibility of the cDNA-AFLP reactions, pre-and selective amplification reactions were performed twice independently using three randomly chosen cDNA-AFLP primer sets and the image resolved on a 6.5% PAGE is provided in Additional file 3.
Cloning, sequencing, and sequence analysis
After recovery of individual TDFs from cDNA-AFLP, fragments were re-amplified with primers used for pre-amplification and subcloned into pCR-TOPO cloning vector systems (Invitrogen) following clean-up of PCR reaction products with a Wizard PCR clean-up kit (Promega). Three to four colonies from each band were bi-directionally sequenced. Details about cloning, sequencing, and RT-PCR amplification procedures have been described previously . Sequence information of gene-specific primers used for RT-PCR is available upon request. For the identity of sequenced TDFs, homology-based blast searches from the National Center for Biotechnology Information (NCBI, http://www.ncbi.nlm.nih.gov/) were performed, and predicted protein sequences and ESTs with highest similarity were annotated. Functional annotation of each TDF was also performed following the Gene Ontology tool . When necessary, cotton assembly contig sequences were used for additional blast homology searches with a cotton46 version that is available from the Comparative Evolutionary Genomics of Cotton web site (http://cottonevolution.info/). Alignment and assembly of sequences were processed using VectorNTI (Invitrogen) and ClustalW2 programs (http://www.ebi.ac.uk/Tools/msa/clustalw2/).
Amplified fragment length polymorphism
Reverse transcription-polymerase chain reaction
Transcript derived fragment
Reactive oxygen species
Heat shock protein.
This research project was supported by funding from CRIS No. 6657-21000-006-00D of the U.S. Department of Agriculture and a grant from Cotton Incorporated. Authors acknowledge assistance from Dr. Paxton Payton for the RNA isolation protocol, Drs. Joshua Udall and Martin Wubben for critically reviewing this manuscript, Drs. Michel Claverie and Jean Marc Lacape for the cDNA-AFLP protocol, Fanny Liu for DNA sequencing, and Dr. Greg Constable for supplying cotton seeds. We acknowledge Mr. Jeff Chandler and North Carolina State University Sandhills Research Station staff for providing excellent technical field plot support. We also thank Ernie Strickland, Kendreal Wingate, and Bobby Fisher for their technical assistance. Mention of trade names or commercial products in this publication is solely for the purpose of providing specific information and does not imply recommendation or endorsement by the U.S. Department of Agriculture.
- Chaves MM, Oliveira MM: Mechanisms underlying plant resilience to water deficits: prospects for water-saving agriculture. J Exp Bot. 2004, 55 (407): 2365-2384. 10.1093/jxb/erh269.PubMedView Article
- Fedoroff NV, Battisti DS, Beachy RN, Cooper PJ, Fischhoff DA, Hodges CN, Knauf VC, Lobell D, Mazur BJ, Molden D, et al: Radically rethinking agriculture for the 21st century. Science. 2010, 327 (5967): 833-834. 10.1126/science.1186834.PubMedPubMed CentralView Article
- Long SP, Ort DR: More than taking the heat: crops and global change. Curr Opin Plant Biol. 2010, 13 (3): 241-248.PubMedView Article
- Zhao M, Running SW: Drought-induced reduction in global terrestrial net primary production from 2000 through 2009. Science. 2010, 329 (5994): 940-943. 10.1126/science.1192666.PubMedView Article
- Morison JIL, Baker NR, Mullineaux PM, Davies WJ: Improving water use in crop production. Phil Trans R Soc B. 2008, 363: 639-658. 10.1098/rstb.2007.2175.PubMedPubMed CentralView Article
- Stiller WN, Read JJ, Constable GA, Reid PE: Selection for water use efficiency traits in a cotton breeding program: cultivar differences. Crop Sci. 2005, 45: 1107-1113. 10.2135/cropsci2004.0545.View Article
- Campbell BT, Chee PW, Lubbers E, Bowman DT, Meredith WR, Johnson JM, Fraser DE: Genetic Improvement of the Pee Dee Cotton Germplasm Collection following Seventy Years of Plant Breeding. Crop Sci. 2011, 51: 955-968. 10.2135/cropsci2010.09.0545.View Article
- Bhattarai SP, McHugh AD, Lotz G, Midmore DJ: The response of cotton to subsurface drip and furrow irrigation in a vertisol. Exp Agric. 2006, 42: 29-49. 10.1017/S0014479705003029.View Article
- Truman CC, Strickland TC, Potter TL, Franklin DH, Bosch DD, Bednarz CW: Variable rainfall intensity and tillage effects on runoff, sediment, and carbon losses from a loamy sand under simulated rainfall. J Environ Qual. 2007, 36 (5): 1495-1502. 10.2134/jeq2006.0018.PubMedView Article
- Bauer PJ, Fortnum BA, Frederick JR: Cotton responses to tillage and rotation during the turn of the century drought. Agron J. 2010, 102 (4): 1145-1148. 10.2134/agronj2010.0037.View Article
- Ahuja I, de Vos RC, Bones AM, Hall RD: Plant molecular stress responses face climate change. Trends Plant Sci. 2010, 15 (12): 664-674. 10.1016/j.tplants.2010.08.002.PubMedView Article
- Seki M, Narusaka M, Ishida J, Nanjo T, Fujita M, Oono Y, Kamiya A, Nakajima M, Enju A, Sakurai T, et al: Monitoring the expression profiles of 7000 Arabidopsis genes under drought, cold and high-salinity stresses using a full-length cDNA microarray. Plant J. 2002, 31 (3): 279-292. 10.1046/j.1365-313X.2002.01359.x.PubMedView Article
- Rabbani MA, Maruyama K, Abe H, Khan MA, Katsura K, Ito Y, Yoshiwara K, Seki M, Shinozaki K, Yamaguchi-Shinozaki K: Monitoring expression profiles of rice genes under cold, drought, and high-salinity stresses and abscisic acid application using cDNA microarray and RNA gel-blot analyses. Plant Physiol. 2003, 133 (4): 1755-1767. 10.1104/pp.103.025742.PubMedPubMed CentralView Article
- Umezawa T, Fujita M, Fujita Y, Yamaguchi-Shinozaki K, Shinozaki K: Engineering drought tolerance in plants: discovering and tailoring genes to unlock the future. Curr Opin Biotechnol. 2006, 17 (2): 113-122. 10.1016/j.copbio.2006.02.002.PubMedView Article
- Cnudde F, Hedatale V, de Jong H, Pierson ES, Rainey DY, Zabeau M, Weterings K, Gerats T, Peters JL: Changes in gene expression during male meiosis in Petunia hybrida. Chromosome Res. 2006, 14 (8): 919-932. 10.1007/s10577-006-1099-5.PubMedView Article
- Bellin D, Ferrarini A, Chimento A, Kaiser O, Levenkova N, Bouffard P, Delledonne M: Combining next-generation pyrosequencing with microarray for large scale expression analysis in non-model species. BMC Genomics. 2009, 10: 555-10.1186/1471-2164-10-555.PubMedPubMed CentralView Article
- Bachem CW, van der Hoeven RS, de Bruijn SM, Vreugdenhil D, Zabeau M, Visser RG: Visualization of differential gene expression using a novel method of RNA fingerprinting based on AFLP: analysis of gene expression during potato tuber development. Plant J. 1996, 9 (5): 745-753. 10.1046/j.1365-313X.1996.9050745.x.PubMedView Article
- Akihiro T, Umezawa T, Ueki C, Lobna BM, Mizuno K, Ohta M, Fujimura T: Genome wide cDNA-AFLP analysis of genes rapidly induced by combined sucrose and ABA treatment in rice cultured cells. FEBS Lett. 2006, 580 (25): 5947-5952. 10.1016/j.febslet.2006.09.065.PubMedView Article
- Vuylsteke M, Peleman JD, van Eijk MJ: AFLP-based transcript profiling (cDNA-AFLP) for genome-wide expression analysis. Nat Protoc. 2007, 2 (6): 1399-1413. 10.1038/nprot.2007.174.PubMedView Article
- Durrant WE, Rowland O, Piedras P, Hammond-Kosack KE, Jones JD: cDNA-AFLP reveals a striking overlap in race-specific resistance and wound response gene expression profiles. Plant Cell. 2000, 12 (6): 963-977.PubMedPubMed CentralView Article
- Dubos C, Plomion C: Identification of water-deficit responsive genes in maritime pine (Pinus pinaster Ait.) roots. Plant Mol Biol. 2003, 51 (2): 249-262. 10.1023/A:1021168811590.PubMedView Article
- Wang X, Tang C, Zhang G, Li Y, Wang C, Liu B, Qu Z, Zhao J, Han Q, Huang L, et al: cDNA-AFLP analysis reveals differential gene expression in compatible interaction of wheat challenged with Puccinia striiformis f. sp. tritici. BMC Genomics. 2009, 10: 289-10.1186/1471-2164-10-289.PubMedPubMed CentralView Article
- Ackerson RC, Hebert RR: Osmoregulation in cotton in response to water stress : I. Alterations in photosynthesis, leaf conductance, translocation, and ultrastructure. Plant Physiol. 1981, 67 (3): 484-488. 10.1104/pp.67.3.484.PubMedPubMed CentralView Article
- Grimes DW, Yamada H: Relation of cotton growth and yield to minimum leaf water potential. Crop Sci. 1982, 22: 134-139. 10.2135/cropsci1982.0011183X002200010031x.View Article
- Cothren J: Physiology of the Cotton Plant. In Cotton: Origin, History,Technology, and Production. Edited by Smith CW, Cothren J.: New York, NewYork: Wiley; 1999:207–268.
- Bradow JM, Davidonis GH: Effects of environment on fiber quality. InPhysiology of cotton. Edited by Stewart JM, Oosterhuis DM, Heitholt JJ,Mauney JR.: Heidelberg, Germany: Springer; 2010:229–245.View Article
- Martin EC, Stephens W, Wiedenfeld R, Bittenbender HC, Beasley JPJ, MooreJM, Neibling H, CGallian JJ: Sugars, oil, & fiber crops. In Irrigation ofagricultural crops. 2nd edition. Edited by Lascano RJ, Sojka RE.: Madison,Wisconsin: ASA-CSSA-SSSA; 2007:279–335. vol. 30.
- Saranga Y, Menz M, Jiang CX, Wright RJ, Yakir D, Paterson AH: Genomic dissection of genotype x environment interactions conferring adaptation of cotton to arid conditions. Genome Res. 2001, 11 (12): 1988-1995. 10.1101/gr.157201.PubMedView Article
- Zhang L, Li FG, Liu CL, Zhang CJ, Zhang XY: Construction and analysis of cotton (Gossypium arboreum L.) drought-related cDNA library. BMC Res Notes. 2009, 2: 120-10.1186/1756-0500-2-120.PubMedPubMed CentralView Article
- Levi A, Paterson AH, Cakmak I, Saranga Y: Metabolite and mineral analyses of cotton near-isogenic lines introgressed with QTLs for productivity and drought-related traits. Physiol Plant. 2011, 141 (3): 265-275. 10.1111/j.1399-3054.2010.01438.x.PubMedView Article
- Payton P, Kottapalli KR, Kebede H, Mahan JR, Wright RJ, Allen RD: Examining the drought stress transcriptome in cotton leaf and root tissue. Biotechnol Lett. 2011, 33 (4): 821-828. 10.1007/s10529-010-0499-y.PubMedView Article
- Feuillet C, Leach JE, Rogers J, Schnable PS, Eversole K: Crop genome sequencing: lessons and rationales. Trends Plant Sci. 2011, 16 (2): 77-88. 10.1016/j.tplants.2010.10.005.PubMedView Article
- Taliercio E, Allen RD, Essenberg M, Klueva N, Nguyen H, Patil MA, Payton P, Millena AC, Phillips AL, Pierce ML, et al: Analysis of ESTs from multiple Gossypium hirsutum tissues and identification of SSRs. Genome. 2006, 49 (4): 306-319. 10.1139/G05-115.PubMedView Article
- Udall JA, Swanson JM, Haller K, Rapp RA, Sparks ME, Hatfield J, Yu Y, Wu Y, Dowd C, Arpat AB, et al: A global assembly of cotton ESTs. Genome Res. 2006, 16 (3): 441-450. 10.1101/gr.4602906.PubMedPubMed CentralView Article
- Paterson AH, Rong J, Gingle AR, Chee PW, Dennis ES, Llewellyn D, Dure LS, Haigler C, Myers JO, Peterson DG, et al: Sequencing and Utilization of the Gossypium Genomes. Tropical Plant Biol. 2010, 3: 71-74. 10.1007/s12042-010-9051-4.View Article
- Guo W, Cai C, Wang C, Zhao L, Wang L, Zhang T: A preliminary analysis of genome structure and composition in Gossypium hirsutum. BMC Genomics. 2008, 9: 314-10.1186/1471-2164-9-314.PubMedPubMed CentralView Article
- Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, et al: Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet. 2000, 25 (1): 25-29. 10.1038/75556.PubMedPubMed CentralView Article
- Flagel LE, Wendel JF: Gene duplication and evolutionary novelty in plants. New Phytol. 2009, 183 (3): 557-564. 10.1111/j.1469-8137.2009.02923.x.PubMedView Article
- Park W, Scheffler BE, Bauer PJ, Campbell BT: Identification of the family of aquaporin genes and their expression in upland cotton (Gossypium hirsutum L.). BMC Plant Biol. 2010, 10: 142-10.1186/1471-2229-10-142.PubMedPubMed CentralView Article
- Froger A, Tallur B, Thomas D, Delamarche C: Prediction of functional residues in water channels and related proteins. Protein Sci. 1998, 7 (6): 1458-1468. 10.1002/pro.5560070623.PubMedPubMed CentralView Article
- Maurel C, Verdoucq L, Luu DT, Santoni V: Plant aquaporins: membrane channels with multiple integrated functions. Annu Rev Plant Biol. 2008, 59: 595-624. 10.1146/annurev.arplant.59.032607.092734.PubMedView Article
- Yang Z, Wu Y, Li Y, Ling HQ, Chu C: OsMT1a, a type 1 metallothionein, plays the pivotal role in zinc homeostasis and drought tolerance in rice. Plant Mol Biol. 2009, 70 (1–2): 219-229.PubMedView Article
- Foyer CH, Noctor G: Redox homeostasis and antioxidant signaling: a metabolic interface between stress perception and physiological responses. Plant Cell. 2005, 17 (7): 1866-1875. 10.1105/tpc.105.033589.PubMedPubMed CentralView Article
- Huang XY, Chao DY, Gao JP, Zhu MZ, Shi M, Lin HX: A previously unknown zinc finger protein, DST, regulates drought and salt tolerance in rice via stomatal aperture control. Genes Dev. 2009, 23 (15): 1805-1817. 10.1101/gad.1812409.PubMedPubMed CentralView Article
- Jubany-Mari T, Munne-Bosch S, Alegre L: Redox regulation of water stress responses in field-grown plants. Role of hydrogen peroxide and ascorbate. Plant Physiol Biochem. 2010, 48 (5): 351-358. 10.1016/j.plaphy.2010.01.021.PubMedView Article
- Benitez-Alfonso Y, Jackson D, Maule A: Redox regulation of intercellular transport. Protoplasma. 2011, 248 (1): 131-140. 10.1007/s00709-010-0243-4.PubMedView Article
- Ruan YL, Xu SM, White R, Furbank RT: Genotypic and developmental evidence for the role of plasmodesmatal regulation in cotton fiber elongation mediated by callose turnover. Plant Physiol. 2004, 136 (4): 4104-4113. 10.1104/pp.104.051540.PubMedPubMed CentralView Article
- Kim JY, Mahe A, Brangeon J, Prioul JL: A maize vacuolar invertase, IVR2, is induced by water stress. Organ/tissue specificity and diurnal modulation of expression. Plant Physiol. 2000, 124 (1): 71-84. 10.1104/pp.124.1.71.PubMedPubMed CentralView Article
- Mou Z, Wang X, Fu Z, Dai Y, Han C, Ouyang J, Bao F, Hu Y, Li J: Silencing of phosphoethanolamine N-methyltransferase results in temperature-sensitive male sterility and salt hypersensitivity in Arabidopsis. Plant Cell. 2002, 14 (9): 2031-2043. 10.1105/tpc.001701.PubMedPubMed CentralView Article
- Turoczy Z, Kis P, Torok K, Cserhati M, Lendvai A, Dudits D, Horvath GV: Overproduction of a rice aldo-keto reductase increases oxidative and heat stress tolerance by malondialdehyde and methylglyoxal detoxification. Plant Mol Biol. 2011, 75 (4–5): 399-412.PubMedView Article
- Zhu JK: Salt and drought stress signal transduction in plants. Annu Rev Plant Biol. 2002, 53: 247-273. 10.1146/annurev.arplant.53.091401.143329.PubMedPubMed CentralView Article
- Kim TH, Bohmer M, Hu H, Nishimura N, Schroeder JI: Guard cell signal transduction network: advances in understanding abscisic acid, CO2, and Ca2+ signaling. Annu Rev Plant Biol. 2010, 61: 561-591. 10.1146/annurev-arplant-042809-112226.PubMedPubMed CentralView Article
- Ackerson RC: Synthesis and movement of abscisic Acid in water-stressed cotton leaves. Plant Physiol. 1982, 69 (3): 609-613. 10.1104/pp.69.3.609.PubMedPubMed CentralView Article
- Leung J, Giraudat J: Abscisic Acid Signal Transduction. Annu Rev Plant Physiol Plant Mol Biol. 1998, 49: 199-222. 10.1146/annurev.arplant.49.1.199.PubMedView Article
- Yoo CY, Pence HE, Hasegawa PM, Mickelbart MV: Regulation of Transpiration to Improve Crop Water Use. Critical Reviews in Plant Science. 2009, 28: 410-431. 10.1080/07352680903173175.View Article
- Cho EK, Hong CB: Over-expression of tobacco NtHSP70-1 contributes to drought-stress tolerance in plants. Plant Cell Rep. 2006, 25 (4): 349-358. 10.1007/s00299-005-0093-2.PubMedView Article
- Clement M, Leonhardt N, Droillard MJ, Reiter I, Montillet JL, Genty B, Lauriere C, Nussaume L, Noel LD: The Cytosolic/Nuclear HSC70 and HSP90 Molecular Chaperones Are Important for Stomatal Closure and Modulate Abscisic Acid-Dependent Physiological Responses in Arabidopsis. Plant Physiol. 2011, 156 (3): 1481-1492. 10.1104/pp.111.174425.PubMedPubMed CentralView Article
- Maqbool A, Zahur M, Irfan M, Qaiser U, Rashid B, Husnain T, Riazuddin S: Identification, characterization and expression of drought related a-crystalline heat shock protein gene (GHSP) from Desi Cotton (Gossypium arboreum L.). Crop Sci. 2007, 47: 2437-2444. 10.2135/cropsci2007.03.0120.View Article
- Maqbool A, Abbas W, Rao AQ, Irfan M, Zahur M, Bakhsh A, Riazuddin S, Husnain T: Gossypium arboreum GHSP26 enhances drought tolerance in Gossypium hirsutum. Biotechnol Prog. 2010, 26 (1): 21-25.PubMed
- Jang JY, Kim DG, Kim YO, Kim JS, Kang H: An expression analysis of a gene family encoding plasma membrane aquaporins in response to abiotic stresses in Arabidopsis thaliana. Plant Mol Biol. 2004, 54 (5): 713-725.PubMedView Article
- Lee HK, Cho SK, Son O, Xu Z, Hwang I, Kim WT: Drought stress-induced Rma1H1, a RING membrane-anchor E3 ubiquitin ligase homolog, regulates aquaporin levels via ubiquitination in transgenic Arabidopsis plants. Plant Cell. 2009, 21 (2): 622-641. 10.1105/tpc.108.061994.PubMedPubMed CentralView Article
- Alexandersson E, Danielson JA, Rade J, Moparthi VK, Fontes M, Kjellbom P, Johanson U: Transcriptional regulation of aquaporins in accessions of Arabidopsis in response to drought stress. Plant J. 2010, 61 (4): 650-660. 10.1111/j.1365-313X.2009.04087.x.PubMedView Article
- Rae L: Lao NT. Kavanagh TA: Regulation of multiple aquaporin genes inArabidopsis by a pair of recently duplicated DREB transcription factors.Planta; 2011.
- Shatil-Cohen A, Attia Z, Moshelion M: Bundle-sheath cell regulation of xylem-mesophyll water transport via aquaporins under drought stress: a target of xylem-borne ABA?. Plant J. 2011, 67 (1): 72-80. 10.1111/j.1365-313X.2011.04576.x.PubMedView Article
- Hayano-Kanashiro C, Calderon-Vazquez C, Ibarra-Laclette E, Herrera-Estrella L, Simpson J: Analysis of gene expression and physiological responses in three Mexican maize landraces under drought stress and recovery irrigation. PLoS One. 2009, 4 (10): e7531-10.1371/journal.pone.0007531.PubMedPubMed CentralView Article
- Jain D, Chattopadhyay D: Analysis of gene expression in response to water deficit of chickpea (Cicer arietinum L.) varieties differing in drought tolerance. BMC Plant Biol. 2010, 10: 24-10.1186/1471-2229-10-24.PubMedPubMed CentralView Article
- Lata C, Sahu PP, Prasad M: Comparative transcriptome analysis of differentially expressed genes in foxtail millet (Setaria italica L.) during dehydration stress. Biochem Biophys Res Commun. 2010, 393 (4): 720-727. 10.1016/j.bbrc.2010.02.068.PubMedView Article
- Molina C, Rotter B, Horres R, Udupa SM, Besser B, Bellarmino L, Baum M, Matsumura H, Terauchi R, Kahl G, et al: SuperSAGE: the drought stress-responsive transcriptome of chickpea roots. BMC Genomics. 2008, 9: 553-10.1186/1471-2164-9-553.PubMedPubMed CentralView Article
- Uehlein N, Lovisolo C, Siefritz F, Kaldenhoff R: The tobacco aquaporin NtAQP1 is a membrane CO2 pore with physiological functions. Nature. 2003, 425 (6959): 734-737. 10.1038/nature02027.PubMedView Article
- Flexas J, Ribas-Carbo M, Hanson DT, Bota J, Otto B, Cifre J, McDowell N, Medrano H, Kaldenhoff R: Tobacco aquaporin NtAQP1 is involved in mesophyll conductance to CO2 in vivo. Plant J. 2006, 48 (3): 427-439. 10.1111/j.1365-313X.2006.02879.x.PubMedView Article
- Heckwolf M, Pater D, Hanson DT, Kaldenhoff R: The Arabidopsis thaliana aquaporin AtPIP1;2 is a physiologically relevant CO(2) transport facilitator. Plant J. 2011, 67 (5): 795-804. 10.1111/j.1365-313X.2011.04634.x.PubMedView Article
- Sutka M, Li G, Boudet J, Boursiac Y: Doumas P. Maurel C: Natural variation ofroot hydraulics in Arabidopsis grown in normal and salt stress conditions.Plant Physiol; 2011.
- Gollery M, Harper J, Cushman J, Mittler T, Mittler R: POFs: what we don't know can hurt us. Trends Plant Sci. 2007, 12 (11): 492-496. 10.1016/j.tplants.2007.08.018.PubMedView Article
- Pennisi E: Plant genetics. The blue revolution, drop by drop, gene by gene. Science. 2008, 320 (5873): 171-173. 10.1126/science.320.5873.171.PubMedView Article
- Setter TL, Yan J, Warburton M, Ribaut JM, Xu Y, Sawkins M, Buckler ES, Zhang Z, Gore MA: Genetic association mapping identifies single nucleotide polymorphisms in genes that affect abscisic acid levels in maize floral tissues during drought. J Exp Bot. 2011, 62 (2): 701-716. 10.1093/jxb/erq308.PubMedPubMed CentralView Article
- Des Marais DL, Rausher MD: Escape from adaptive conflict after duplication in an anthocyanin pathway gene. Nature. 2008, 454 (7205): 762-765.PubMed
- Buggs RJ, Zhang L, Miles N, Tate JA, Gao L, Wei W, Schnable PS, Barbazuk WB, Soltis PS, Soltis DE: Transcriptomic shock generates evolutionary novelty in a newly formed, natural allopolyploid plant. Curr Biol. 2011, 21 (7): 551-556. 10.1016/j.cub.2011.02.016.PubMedView Article
- Dong S, Adams KL: Differential contributions to the transcriptome of duplicated genes in response to abiotic stresses in natural and synthetic polyploids. New Phytol. 2011, 190 (4): 1045-1057. 10.1111/j.1469-8137.2011.03650.x.PubMedView Article
- Udall JA, Swanson JM, Nettleton D, Percifield RJ, Wendel JF: A novel approach for characterizing expression levels of genes duplicated by polyploidy. Genetics. 2006, 173 (3): 1823-1827. 10.1534/genetics.106.058271.PubMedPubMed CentralView Article
- Lukens LN, Zhan S: The plant genome's methylation status and response to stress: implications for plant improvement. Curr Opin Plant Biol. 2007, 10 (3): 317-322. 10.1016/j.pbi.2007.04.012.PubMedView Article
- Wang WS, Pan YJ, Zhao XQ, Dwivedi D, Zhu LH, Ali J, Fu BY, Li ZK: Drought-induced site-specific DNA methylation and its association with drought tolerance in rice (Oryza sativa L.). J Exp Bot. 2011, 62 (6): 1951-1960. 10.1093/jxb/erq391.PubMedPubMed CentralView Article
- Campbell BT, Bauer PJ: Genetic Variation for Yield and Fiber Quality Response to Supplemental Irrigation within the Pee Dee Upland Cotton Germplasm Collection. Crop Sci. 2007, 47: 589-597.View Article
- Nepomuceno AL, Oosterhuisb DM, Stewartb JM: Physiological responses of cotton leaves and roots to water deficit induced by polyethylene glycol. Environ Exp Bot. 1998, 40 (1): 29-41. 10.1016/S0098-8472(98)00018-5.View Article
- Voloudakis AE, Kosmas SA, Tsakas S, Eliopoulos E, Loukas M, Kosmidou K: Expression of selected drought-related genes and physiological response of Greek cotton varieties. FuncPlant Biol. 2002, 29 (10): 1237-1245.
- Wan CY, Wilkins TA: A modified hot borate method significantly enhances the yield of high-quality RNA from cotton (Gossypium hirsutum L.). Anal Biochem. 1994, 223 (1): 7-12. 10.1006/abio.1994.1538.PubMedView Article
- Bachem CW, Oomen RJ, Visser RG: Transcript imaging with cDNA-AFLP: a step-by-step protocol. Plant Mol Biol Rep. 1998, 16: 157-173. 10.1023/A:1007468801806.View Article
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.