Distribution and differentiation of wild, feral, and cultivated populations of perennial upland cotton (Gossypium hirsutum L.) in Mesoamerica and the Caribbean

Origin and diversity of the wild cottons (Gossypium hirsutum) of Mound Key, Florida

Elucidating genetic diversity within wild forms of modern crops is essential for understanding domestication and the possibilities of wild germplasm utilization. Gossypium hirsutum is a predominant source of natural plant fibers and the most widely cultivated cotton species. Wild forms of G. hirsutum are challenging to distinguish from feral derivatives, and truly wild populations are uncommon. Here we characterize a population from Mound Key Archaeological State Park, Florida using genome-wide SNPs extracted from 25 individuals over three sites. Our results reveal that this population is genetically dissimilar from other known wild, landrace, and domesticated cottons, and likely represents a pocket of previously unrecognized wild genetic diversity. The unexpected level of divergence between the Mound Key population and other wild cotton populations suggests that the species may harbor other remnant and genetically distinct populations that are geographically scattered in suitable habitats throughout the Caribbean. Our work thus has broader conservation genetic implications and suggests that further exploration of natural diversity in this species is warranted.

Gossypium purpurascens genome provides insight into the origin and domestication of upland cotton

INTRODUCTION

Allotetraploid upland cotton (Gossypium hirsutum L.) is native to the Mesoamerican and Caribbean regions, had been improved in the southern United States by the mid-eighteenth century, was then dispersed worldwide. However, a Hainan Island Native Cotton (HIC) has long been grown extensively on Hainan Island, China.

OBJECTIVES

Explore HIC's evolutionary relationship and genomic diversity with other tetraploid cottons, its origin and whether it was used for YAZHOUBU (Yazhou cloth, World Intangible Cultural Heritage) weaving, and the role of structural variations (SVs) in upland cotton domestication.

METHODS

We assembled a high-quality genome of one HIC plant. We performed phylogenetic analysis, divergence time estimation, principal component analysis and population differentiation estimation using cotton assemblies and/or resequencing data. SVs were detected by whole-genome comparison. A F2 population was used for linkage analysis and to study effects of SVs. Buoyancy and salt water tolerance tests for seeds were conducted.

RESULTS

We found that the HIC belongs to G. purpurascens. G. purpurascens is best classified as a primitive race of G. hirsutum. The potential for long range transoceanic dispersal of G. purpurascens seeds was proved. A set of SVs, selective sweep regions between G. hirsutum races and cultivars, and quantitative trait loci (QTLs) of eleven agronomic traits were obtained. SVs, especially large-scale SVs, were found to have important effects on cotton domestication and improvement. Of them, eight large-scale inversions strongly associated with yield and fiber quality have probably undergone artificial selection in domestication.

CONCLUSION

G. purpurascens including HIC is a primitive race of G. hirsutum, probably disperse to Hainan from Central America by floating on ocean currents, may have been partly domesticated, planted and was likely used for YAZHOUBU weaving in Hainan much earlier than the Pre-Columbian period. SV plays an important role in cotton domestication and improvement.

Aphid and caterpillar feeding drive similar patterns of induced defences and resistance to subsequent herbivory in wild cotton

MAIN CONCLUSION

Our results indicate caterpillars and aphids cause similar levels of induced defences and resistance against caterpillars in wild cotton plants. These symmetrical effects are not consistent with patterns predicted by plant defensive signaling crosstalk and call for further work addressing the biochemical mechanisms underpinning these results. Plant-induced responses to attack often mediate interactions between different species of insect herbivores. These effects are predicted to be contingent on the herbivore's feeding guild, whereby prior feeding by insects should negatively impact subsequent feeding by insects of the same guild (induced resistance) but may positively influence insects of a different guild (induced susceptibility) due to interfering crosstalk between plant biochemical pathways specific to each feeding guild. We compared the effects of prior feeding by leaf-chewing caterpillars (Spodoptera frugiperda) vs. sap-sucking aphids (Aphis gossypii) on induced defences in wild cotton (Gossypium hirsutum) and the consequences of these attacks on subsequently feeding caterpillars (S. frugiperda). To this end, we conducted a greenhouse experiment where cotton plants were either left undamaged or first exposed to caterpillar or aphid feeding, and we subsequently placed caterpillars on the plants to assess their performance. We also collected leaves to assess the induction of chemical defences in response to herbivory. We found that prior feeding by both aphids and caterpillars resulted in reductions in consumed leaf area, caterpillar mass gain, and caterpillar survival compared with control plants. Concomitantly, prior aphid and caterpillar herbivory caused similar increases in phenolic compounds (flavonoids and hydroxycinnamic acids) and defensive terpenoids (hemigossypolone) compared with control plants. Overall, these findings indicate that these insects confer a similar mode and level of induced resistance in wild cotton plants, calling for further work addressing the biochemical mechanisms underpinning these effects.

Soil salinization disrupts plant-plant signaling effects on extra-floral nectar induction in wild cotton

Plant-plant interactions via volatile organic compounds (VOCs) have received much attention, but how abiotic stresses affect these interactions is poorly understood. We tested the effect of VOCs exposure from damaged conspecifics on the production of extra-floral nectar (EFN) in wild cotton plants (Gossypium hirsutum), a coastal species in northern Yucatan (Mexico), and whether soil salinization affected these responses. We placed plants in mesh cages, and within each cage assigned plants as emitters or receivers. We exposed emitters to either ambient or augmented soil salinity to simulate a salinity shock, and within each group subjected half of the emitters to no damage or artificial leaf damage with caterpillar regurgitant. Damage increased the emission of sesquiterpenes and aromatic compounds under ambient but not under augmented salinity. Correspondingly, exposure to VOCs from damaged emitters had effect on receiver EFN induction, but this effect was contingent on salinization. Receivers produced more EFN in response to damage after being exposed to VOCs from damaged emitters when the latter were grown under ambient salinity, but not when they were subjected to salinization. These results suggest complex effects of abiotic factors on VOC-mediated plant interactions.

First Report on the Genetic Diversity of Populations of Gossypium barbadense L. and Gossypium hirsutun L. in the Amazonian Native Communities, Cusco-Peru

The genus Gossypium has important ethnobotanical and economic value for Amazonian Native Communities (A.N.C.). However, little research has been undertaken on the distribution and genetic diversity of cotton populations maintained in the Peruvian rainforest. This work aims to present the first report on the genetic diversity of Gossypium spp. populations in the A.N.C. of the province of La Convención, Cusco-Peru. The methodology was based on exploring, collecting, identifying, and characterizing the Gossypium populations present in the A.N.C. Twenty-six descriptors were evaluated (9 quantitative and 17 qualitative), and with this information, distribution, correlation, and principal component (PC) analyses were carried out. As a result, plants of two species [G. barbadense L. (44 samples) and G. hirsutum L. (19 samples)], one variety [G. barbadense var. brasiliensis (75 samples)], and three previously unidentified variations (9 samples) were identified. Altogether, 147 samples were collected. G. barbadense var. brasiliensis, which was always found in association with other economic crops within an altitude range of 338 to 1086 m, was the most predominant (51%), distributed in eleven A.N.C. and always in small plots (up to 2 ha). G. barbadense L. was cultivated between 397 and 1137 m of altitude in eight A.N.C. in plots of up to 3 ha in marginal lands. G. hirsutum L., with a smaller distribution (13%), was found between 334 and 497 m of altitude in only three communities; this species is cultivated in marginal areas throughout the year. The variability found for the first two PCs when considering the quantitative and qualitative descriptors was high (74.7%) and moderate (48.2%), respectively. When combining all the descriptors, the analysis showed that the first two PCs accounted for 51.8% of the total variability of the data. The PCs of the two types of data and their combination confirmed that the three populations found were grouped. The nine undefined samples were close to or intermediate between the described ones, showing that these samples may be the result of spontaneous crosses; as such, these samples need to be better evaluated with other tools for further definition. The information obtained shows that in the A.N.C. of Cusco-Peru, there is variability conserved by the inhabitants, who have been able to maintain and use these genotypes, even from their Amazonian indigenous ancestry, and the environment has been able to generate variability among the species, as will be highlighted in future works.

Ecology of Echinops giganteus A. Rich. in Sub-Saharan Africa: Distribution, Ecoclimatic Niche, and Phytosociology

Echinops giganteus A. Rich. is an aromatic and medicinal plant of the Asteraceae family exploited in Cameroon under the access and benefit sharing (ABS) standard. Despite its importance, little information exists on the ecology of E. giganteus. The aim of the present study was to contribute to a better understanding of its ecology for sustainable management in the Western Highlands of Cameroon. Occurrence data as well as stationary ecological information were collected in the field and from different databases. Bioclimatic data were extracted from the WorldClim database and processed using DIVA-GIS and Maxent software. The Braun-Blanquet quadrat method was used for the phytosociological study. Results showed that the distribution of E. giganteus in its wild state is restricted to sub-Saharan Africa. This distribution is likely conditioned by altitude (1000 m–2600 m), light, temperature, and rainfall. The bioclimatic variables that best explained this distribution were the mean annual temperature (Bio1: 38.8%) and the precipitation of the coldest quarter (Bio19: 24.9%), and their favorable ranges were between 2°C–32°C and 300 mm–1800 mm, respectively. E. giganteus is a heliophilic plant that prefers well-drained substrates and would not have a requirement for organic matter. The floristic analysis of the E. giganteus community identified 68 plant species in 59 genera and 28 botanical families, with the most represented family being the Asteraceae (49%). The average species richness per quadrat was 8 species, dominated by herbs. Species consistent with E. giganteus were Aspilia africana (Pers.) C. D. Adams and Imperata cylindrica (L.) P. Beauv. Chamaephytes and Phanerophytes were dominant among the biological types, while the phytogeographic types were dominated by Pantropical species (38.23%). The most represented diaspore types and modes of dissemination were pogonochores (35.85%) and anemochores (55.38%).

Perennial Cotton Ratoon Cultivation: A Sustainable Method for Cotton Production and Breeding

Cotton production is challenged by high costs with multiple management and material inputs including seed, pesticide, and fertilizer application. The production costs can be decreased and profits can be increased by developing efficient crop management strategies, including perennial cotton ratoon cultivation. This review focuses on the role of ratoon cultivation in cotton productivity and breeding. In areas that are frost-free throughout the year, when the soil temperature is suitable for cotton growth in spring, the buds of survived plants begin to sprout, and so their flowering and fruiting periods are approximately 4-6 weeks earlier than those of sown cotton. Due to the absence of frost damage, the ratoon cotton continues to grow, and the renewed plants can offer a higher yield than cotton sown in the following season. Moreover, ratoon cultivation from the last crop without sowing can help conserve seeds, reduce labor inputs, and reduce soil and water loss. In this review, the preservation of perennial cotton germplasm resources, the classification and genome assignment of perennial species in the cotton gene pools, and effective strategies for the collection, preservation, identification, and utilization of perennial cotton germplasms are discussed. Ratoon cultivation is the main driver of cotton production and breeding, especially to maintain male sterility for the utilization and fixation of heterosis. Ratoon cultivation of cotton is worth adopting because it has succeeded in Brazil, China, and India. Therefore, taking advantages of the warm environment to exploit the indeterminant growth habit of perennial cotton for breeding would be an efficiency-increasing, cost-saving, and eco-friendly approach in frost-free regions. In the future, more attention should be given to ratooning perennial cotton for breeding male-sterile lines.

Potential Distribution of Wild Host Plants of the Boll Weevil (Anthonomus grandis) in the United States and Mexico

The boll weevil (Anthonomus grandis Boheman) reproduces on a reported 13 species of wild host plants in North America, two in the United States and 12 in Mexico. The distributions of these plants are of economic importance to pest management and provide insight into the evolutionary history and origin of the BW. However, detailed information regarding the distributions of many of these species is lacking. In this article, we present distribution models for all of the reported significant BW host plants from Mexico and the United States using spatial distribution modelling software. Host plant distributions were divided into two groups: “eastern” and “western.” In Mexico, Hampea nutricia along the Gulf Coast was the most important of the eastern group, and the wild cottons, Gossypium aridum and Gossypium thurberi were most important in the western group. Other species of Hampea, Gossypium, and Cienfuegosia rosei have relatively restricted distributions and are of apparent minimal economic importance. Cienfuegosia drummondii is the only truly wild host in the southern United States, east of New Mexico. Factors determining potential distributions were variable and indicated that species were present in five vegetation types. Ecological and economic considerations of host plant distributions are discussed, as well as threats to host plant conservation.

The witchweed Striga gesnerioides and the cultivated cowpea: A geographical and historical analysis of their West African distribution points to the prevalence of agro-ecological factors and the parasite's multilocal evolution potential

The increasing severity of Striga gesnerioides attacks on cowpea across West Africa has been related to its prolificity, seed mobility and longevity, and adaptation to aridity, in a context of agricultural intensification. To understand this fast extension, we analyzed (1) the distributions of the crop and the witchweed with ecological niche modeling and multivariate climate analysis, and (2) the chronological information available from collections and the literature. The ecoclimatic envelope of S. gesnerioides attacks on cowpea is the same as on wild hosts. Consistently, the modeled distribution of cowpea infestations is closely similar to the simple superposition of the parasite model (involving all hosts) and the crop model. Striga gesnerioides infestations are restricted to the driest component of the cultivated cowpea ecoclimatic niche, corresponding to the Sahelian and Sudano-Sahelian belts and the Dahomey gap. Thus, the parasite distribution, determined by its own requirements, does not constrain cowpea cultivation under Guinean climates. The spatial and temporal distributions of S. gesnerioides field infestations are consistent with an earlier impact on cowpea production in eastern West Africa, related itself to a similar trend in cowpea cultivation intensification from Niger, Nigeria and Benin to Burkina Faso and Ghana. Mali and Senegal were affected later, and literature reports of Senegalese strains of S. gesnerioides from the wild developing virulence on cowpea offer a model for the diffusion of witchweed parasitism by multilocal evolution, through host-driven selection, instead of epidemic diffusion. A contrario, in Côte d’Ivoire, cowpea is much less widespread, so the parasite has remained confined to the wild compartment. Thus, both historical and ecogeographic analyses refute the vision of S. gesnerioides as an invader. Instead, they point to the increasing importance and intensification of the crop, and the consequent loss of biodiversity, as the main drivers of the extension and diversification of its crop-specific strains.

Eco-Geography of Feral Cotton: A Missing Piece in the Puzzle of Gene Flow Dynamics Among Members of Gossypium hirsutum Primary Gene Pool

Mexico is the center of origin and genetic diversity of upland cotton (Gossypium hirsutum L.), the most important source of natural fiber in the world. Currently, wild and domesticated populations (including genetically modified [GM] varieties) occur in this country and gene flow among them has shaped the species’ genetic diversity and structure, setting a complex and challenging scenario for its conservation. Moreover, recent gene flow from GM cultivars to wild Mexican cotton populations has been reported since 2011. In situ conservation of G. hirsutum requires knowledge about the extent of its geographic distribution, both wild and domesticated, as well as the possible routes and mechanisms that contribute to gene flow between the members of the species wild-to-domesticated continuum (i.e., the primary gene pool). However, little is known about the distribution of feral populations that could facilitate gene flow by acting as bridges. In this study, we analyzed the potential distribution of feral cotton based on an ecological niche modeling approach and discussed its implications in the light of the distribution of wild and domesticated cotton. Then, we examined the processes that could be leading to the escape of seeds from the cultivated fields. Our results indicate that the climatic suitability of feral plants in the environmental and geographic space is broad and overlaps with areas of wild cotton habitat and crop fields, suggesting a region that could bridge cultivated cotton and its wild relatives by allowing gene flow between them. This study provides information for management efforts focused on the conservation of wild populations, native landraces, and non-GM domesticated cotton at its center of origin and genetic diversity.

Parallel and Intertwining Threads of Domestication in Allopolyploid Cotton

The two cultivated allopolyploid cottons, Gossypium hirsutum and Gossypium barbadense, represent a remarkable example of parallel independent domestication, both involving dramatic morphological transformations under selection from wild perennial plants to annualized row crops. Deep resequencing of 643 newly sampled accessions spanning the wild-to-domesticated continuum of both species, and their allopolyploid relatives, are combined with existing data to resolve species relationships and elucidate multiple aspects of their parallel domestication. It is confirmed that wild G. hirsutum and G. barbadense were initially domesticated in the Yucatan Peninsula and NW South America, respectively, and subsequently spread under domestication over 4000-8000 years to encompass most of the American tropics. A robust phylogenomic analysis of infraspecific relationships in each species is presented, quantify genetic diversity in both, and describe genetic bottlenecks associated with domestication and subsequent diffusion. As these species became sympatric over the last several millennia, pervasive genome-wide bidirectional introgression occurred, often with striking asymmetries involving the two co-resident genomes of these allopolyploids. Diversity scans revealed genomic regions and genes unknowingly targeted during domestication and additional subgenomic asymmetries. These analyses provide a comprehensive depiction of the origin, divergence, and adaptation of cotton, and serve as a rich resource for cotton improvement.

Use of hydraulic traits for modeling genotype-specific acclimation in cotton under drought

Understanding the genetic and physiological basis of abiotic stress tolerance under field conditions is key to varietal crop improvement in the face of climate variability. Here, we investigate dynamic physiological responses to water stress in silico and their relationships to genotypic variation in hydraulic traits of cotton (Gossypium hirsutum), an economically important species for renewable textile fiber production. In conjunction with an ecophysiological process-based model, heterogeneous data (plant hydraulic traits, spatially-distributed soil texture, soil water content and canopy temperature) were used to examine hydraulic characteristics of cotton, evaluate their consequences on whole plant performance under drought, and explore potential genotype × environment effects. Cotton was found to have R-shaped hydraulic vulnerability curves (VCs), which were consistent under drought stress initiated at flowering. Stem VCs, expressed as percent loss of conductivity, differed across genotypes, whereas root VCs did not. Simulation results demonstrated how plant physiological stress can depend on the interaction between soil properties and irrigation management, which in turn affect genotypic rankings of transpiration in a time-dependent manner. Our study shows how a process-based modeling framework can be used to link genotypic variation in hydraulic traits to differential acclimating behaviors under drought.

Conservation and Divergence in Duplicated Fiber Coexpression Networks Accompanying Domestication of the Polyploid Gossypium hirsutum L

Gossypium hirsutum L. (Upland cotton) has an evolutionary history involving inter-genomic hybridization, polyploidization, and subsequent domestication. We analyzed the developmental dynamics of the cotton fiber transcriptome accompanying domestication using gene coexpression networks for both joint and homoeologous networks. Remarkably, most genes exhibited expression for at least one homoeolog, confirming previous reports of widespread gene usage in cotton fibers. Most coexpression modules comprising the joint network are preserved in each subgenomic network and are enriched for similar biological processes, showing a general preservation of network modular structure for the two co-resident genomes in the polyploid. Interestingly, only one fifth of homoeologs co-occur in the same module when separated, despite similar modular structures between the joint and homoeologous networks. These results suggest that the genome-wide divergence between homoeologous genes is sufficient to separate their co-expression profiles at the intermodular level, despite conservation of intramodular relationships within each subgenome. Most modules exhibit D-homoeolog expression bias, although specific modules do exhibit A-homoeolog bias. Comparisons between wild and domesticated coexpression networks revealed a much tighter and denser network structure in domesticated fiber, as evidenced by its fewer modules, 13-fold increase in the number of development-related module member genes, and the poor preservation of the wild network topology. These results demonstrate the amazing complexity that underlies the domestication of cotton fiber.

Genetic structure, gene flow pattern, and association analysis of superior germplasm resources in domesticated upland cotton (Gossypium hirsutum L.)

Gene flow patterns and the genetic structure of domesticated crops like cotton are not well understood. Furthermore, marker-assisted breeding of cotton has lagged far behind that of other major crops because the loci associated with cotton traits such as fiber yield and quality have scarcely been identified. In this study, we used 19 microsatellites to first determine the population genetic structure and patterns of gene flow of superior germplasm resources in upland cotton. We then used association analysis to identify which markers were associated with 15 agronomic traits (including ten yield and five fiber quality traits). The results showed that the upland cotton accessions have low levels of genetic diversity (polymorphism information content = 0.427), although extensive gene flow occurred among different ecological and geographic regions. Bayesian clustering analysis indicated that the cotton resources used in this study did not belong to obvious geographic populations, which may be the consequence of a single source of domestication followed by frequent genetic introgression mediated by human transference. A total of 82 maker-trait associations were examined in association analysis and the related ratios for phenotypic variations ranged from 3.04% to 47.14%. Interestingly, nine SSR markers were detected in more than one environmental condition. In addition, 14 SSR markers were co-associated with two or more different traits. It was noteworthy that NAU4860 and NAU5077 markers detected at least in two environments were simultaneously associated with three fiber quality traits (uniformity index, specific breaking strength and micronaire value). In conclusion, these findings provide new insights into the population structure and genetic exchange pattern of cultivated cotton accessions. The quantitative trait loci of domesticated cotton identified will also be very useful for improvement of yield and fiber quality of cotton in molecular breeding programs.

Genetic Analysis of the Transition from Wild to Domesticated Cotton (Gossypium hirsutum L.)

The evolution and domestication of cotton is of great interest from both economic and evolutionary standpoints. Although many genetic and genomic resources have been generated for cotton, the genetic underpinnings of the transition from wild to domesticated cotton remain poorly known. Here we generated an intraspecific QTL mapping population specifically targeting domesticated cotton phenotypes. We used 466 F₂ individuals derived from an intraspecific cross between the wild Gossypium hirsutum var. yucatanense (TX2094) and the elite cultivar G. hirsutum cv. Acala Maxxa, in two environments, to identify 120 QTL associated with phenotypic changes under domestication. While the number of QTL recovered in each subpopulation was similar, only 22 QTL were considered coincident (i.e., shared) between the two locations, eight of which shared peak markers. Although approximately half of QTL were located in the A-subgenome, many key fiber QTL were detected in the D-subgenome, which was derived from a species with unspinnable fiber. We found that many QTL are environment-specific, with few shared between the two environments, indicating that QTL associated with G. hirsutum domestication are genomically clustered but environmentally labile. Possible candidate genes were recovered and are discussed in the context of the phenotype. We conclude that the evolutionary forces that shape intraspecific divergence and domestication in cotton are complex, and that phenotypic transformations likely involved multiple interacting and environmentally responsive factors.

Introgression Leads to Genomic Divergence and Responsible for Important Traits in Upland Cotton

Understanding the genetic diversity and population structure of germplasms is essential when selecting parents for crop breeding. The genomic changes that occurred during the domestication and improvement of Upland cotton (Gossypium hirsutum) remains poorly understood. Besides, the available genetic resources from cotton cultivars are limited. By applying restriction site-associated DNA marker sequencing (RAD-seq) technology to 582 tetraploid cotton accessions, we confirmed distinct genomic regions on chromosomes A06 and A08 in Upland cotton cultivar subgroups. Based on the pedigree, reported QTLs, introgression analyses, and genome-wide association study (GWAS), we suggest that these divergent regions might have resulted from the introgression of exotic lineages of G. hirsutum landraces and their wild relatives. These regions were the typical genomic signatures that might be responsible for maturity and fiber quality on chromosome A06 and chromosome A08, respectively. Moreover, these genomic regions are located in the putative pericentromeric regions, implying that their application will be challenging. In the study, based on high-density SNP markers, we reported two genomic signatures on chromosomes A06 and A08, which might originate from the introgression events in the Upland cotton population. Our study provides new insights for understanding the impact of historic introgressions on population divergence and important agronomic traits of modern Upland cotton cultivars.

Genotyping-by-Sequencing of Gossypium hirsutum Races and Cultivars Uncovers Novel Patterns of Genetic Relationships and Domestication Footprints

Determining the genetic rearrangement and domestication footprints in Gossypium hirsutum cultivars and primitive race genotypes are essential for effective gene conservation efforts and the development of advanced breeding molecular markers for marker-assisted breeding. In this study, 94 accessions representing the 7 primitive races of G hirsutum, along with 9 G hirsutum and 12 Gossypium barbadense cultivated accessions were evaluated. The genotyping-by-sequencing (GBS) approach was employed and 146 558 single nucleotide polymorphisms (SNP) were generated. Distinct SNP signatures were identified through the combination of selection scans and association analyses. Phylogenetic analyses were also conducted, and we concluded that the Latifolium, Richmondi, and Marie-Galante race accessions were more genetically related to the G hirsutum cultivars and tend to cluster together. Fifty-four outlier SNP loci were identified by selection-scan analysis, and 3 SNPs were located in genes related to the processes of plant responding to stress conditions and confirmed through further genome-wide signals of marker-phenotype association analysis, which indicate a clear selection signature for such trait. These results identified useful candidate gene locus for cotton breeding programs.

Compensation to simulated insect leaf herbivory in wild cotton (Gossypium hirsutum): responses to multiple levels of damage and associated traits

Identifying the mechanisms of compensation to insect herbivory remains a major challenge in plant biology and evolutionary ecology. Most previous studies have addressed plant compensatory responses to one or two levels of insect herbivory, and the underlying traits mediating such responses remain elusive in many cases. We evaluated responses associated with compensation to multiple intensities of leaf damage (0% control, 10%, 25%, 50%, 75% of leaf area removed) by means of mechanical removal of foliar tissue and application of a caterpillar (Spodoptera exigua) oral secretions in 3-month-old wild cotton plants (Gossypium hirsutum). Four weeks post-treatment, we measured plant growth and multiple traits associated with compensation, namely: changes in above- and belowground, biomass and the concentration of nutrients (nitrogen and phosphorus) and non-structural carbon reserves (starch and soluble sugars) in roots, stems and leaves. We found that wild cotton fully compensated in terms of growth and biomass allocation when leaf damage was low (10%), whereas moderate (25%) to high leaf damage in some cases led to under-compensation. Nonetheless, high levels of leaf removal (50% and 75%) in most cases did not cause further reductions in height and allocation to leaf and stem biomass relative to low and moderate damage. There were significant positive effects of leaf damage on P concentration in leaves and stems, but not roots, as well as a negative effect on soluble sugars in roots. These results indicate that wild cotton fully compensated for a low level of leaf damage but under-compensated under moderate to high leaf damage, but can nonetheless sustain growth despite increasing losses to herbivory. Such responses were possibly mediated by a re-allocation of carbohydrate reserves from roots to shoots.

Bottom-up control of geographic variation in insect herbivory on wild cotton (Gossypium hirsutum) by plant defenses and climate

PREMISE

The occurrence and amount of herbivory are shaped by bottom-up forces, primarily plant traits (e.g., defenses), and by abiotic factors. Addressing these concurrent effects in a spatial context has been useful in efforts to understand the mechanisms governing variation in plant-herbivore interactions. Still, few studies have evaluated the simultaneous influence of multiple sources of bottom-up variation on spatial variation in herbivory.

METHODS

We tested to what extent chemical (phenolics, production of terpenoid glands) and physical (pubescence) defensive plant traits and climatic factors are associated with variation in herbivory by leaf-chewing insects across populations of wild cotton (Gossypium hirsutum).

RESULTS

We found substantial population variation in cotton leaf defenses and insect leaf herbivory. Leaf pubescence, but not gossypol gland density or phenolic content, was significantly negatively associated with herbivory by leaf-chewing insects. In addition, there were direct effects of climate on defenses and herbivory, with leaf pubescence increasing toward drier conditions and leaf damage increasing toward wetter and cooler conditions. There was no evidence, however, of indirect effects (via plant defenses) of climate on herbivory.

CONCLUSIONS

These results suggest that spatial variation in insect herbivory on wild G. hirsutum is predominantly driven by concurrent and independent influences of population variation in leaf pubescence and climatic factors.

Abundant Small Genetic Alterations after Upland Cotton Domestication

Domestication has long been recognized as the most direct and effective way to intentionally influence morphological and physiological phenotypes in plants and animals. Consequently, understanding how small genetic alterations contribute to domestication is of considerable importance. In this study, we resequenced the genome of the wild upland cotton variety Gossypium hirsutum var. yucatanense, the putative wild ancestor of cultivated upland cotton, and then compared single nucleotide polymorphism (SNP) and short insertion and deletion (InDel) variations of the genome with the cultivated accession (TM-1) of G. hirsutum. We found approximately 6.6 million SNPs and 0.7 million InDels between the two genomes. Most of the small genetic variations were anchored in the noncoding regions. With regard to potential coding genes, we found 24,035 genes with nonsynonymous SNPs. Interestingly, 2603 genes in domesticated cotton are found that have changed the positions of stop codons or shifted reading frames from that in G. hirsutum var. yucatanense. This suggests that domestication may have been selected for mutations that restored gene function or that wild cotton has undergone a number of gene inactivation events since its divergence from cultivated cotton. The former scenario seems most likely due to the intense selective pressure applied during the domestication process. These results demonstrate that, within a relatively short period of time, the cotton genome has been readjusted through small genetic changes. The current study provides useful clues for seeking interesting genes for cotton improvement.

Evaluation of the Impact of Genetically Modified Cotton After 20 Years of Cultivation in Mexico

For more than 20 years cotton has been the most widely sown genetically modified (GM) crop in Mexico. Its cultivation has fulfilled all requirements and has gone through the different regulatory stages. During the last 20 years, both research-institutions and biotech-companies have generated scientific and technical information regarding GM cotton cultivation in Mexico. In this work, we collected data in order to analyze the environmental and agronomic effects of the use of GM cotton in Mexico. In 1996, the introduction of Bt cotton made it possible to reactivate this crop, which in previous years was greatly reduced due to pest problems, production costs and environmental concerns. Bt cotton is a widely accepted tool for cotton producers and has proven to be efficient for the control of lepidopteran pests. The economic benefits of its use are variable, and depend on factors such as the international cotton-prices and other costs associated with its inputs. So far, the management strategies used to prevent development of insect resistance to GM cotton has been successful, and there are no reports of insect resistance development to Bt cotton in Mexico. In addition, no effects have been observed on non-target organisms. For herbicide tolerant cotton, the prevention of herbicide resistance has also been successful since unlike other countries, the onset of resistance weeds is still slow, apparently due to cultural practices and rotation of different herbicides. Environmental benefits have been achieved with a reduction in chemical insecticide applications and the subsequent decrease in primary pest populations, so that the inclusion of other technologies—e.g., use of non-Bt cotton- can be explored. Nevertheless, control measures need to be implemented during transport of the bolls and fiber to prevent dispersal of volunteer plants and subsequent gene flow to wild relatives distributed outside the GM cotton growing areas. It is still necessary to implement national research programs, so that biotechnology and plant breeding advances can be used in the development of cotton varieties adapted to the Mexican particular environmental conditions and to control insect pests of regional importance.

Highly structured genetic diversity of Bixa orellana var. urucurana, the wild ancestor of annatto, in Brazilian Amazonia

Annatto (Bixa orellana L.) is a tropical American crop, commercially valuable due to its application in the food and cosmetics industries as a natural dye. The wild ancestor of cultivated annatto is B. orellana var. urucurana. Although never cultivated, this variety occurs in open forests and anthropogenic landscapes, and is always associated with riparian environments. In this study, we evaluated the genetic diversity and structure of B. orellana var. urucurana populations in Brazilian Amazonia using 16 microsatellite loci. We used Ecological Niche Modeling (ENM) to characterize the potential geographical range of this variety in northern South America. We analyzed 170 samples from 10 municipalities in the states of Rondônia, Pará and Roraima. A total of 194 alleles was observed, with an average of 12.1 alleles per locus. Higher levels of expected (H_E) than observed (H_O) heterozygosities were found for all populations. Bayesian analysis, Neighbor-Joining dendrograms and PCAs suggest the existence of three strongly structured groups of populations. A strong and positive correlation between genetic and geographic distances was found, suggesting that genetic differentiation might be caused by geographic isolation. From species distribution modelling, we detected that South Rondônia, Madre di Dios River basin, Llanos de Mojos, Llanos de Orinoco and eastern Ecuador are highly suitable areas for wild annatto to occur, providing additional targets for future exploration and conservation. Climatic adaptation analyses revealed strong differentiation among populations, suggesting that precipitation plays a key role in wild annatto’s current and potential distribution patterns.

Agriculture sows pests: how crop domestication, host shifts, and agricultural intensification can create insect pests from herbivores

We argue that agriculture as practiced creates pests. We use three examples (Corn leafhopper, Dalbulus maidis; Western corn rootworm, Diabrotica virgifera virgifera; Cotton fleahopper, Pseudatomoscelis seriatus) to illustrate: firstly, how since its origins, agriculture has proven conducive to transforming selected herbivores into pests, particularly through crop domestication and spread, and agricultural intensification, and; secondly, that the herbivores that became pests were among those hosted by crop wild relatives or associates, and were pre-adapted either as whole species or component subpopulations. Two of our examples, Corn leafhopper and Western corn rootworm, illustrate how following a host shift to a domesticated host, emergent pests 'hopped' onto crops and rode expansion waves to spread far beyond the geographic ranges of their wild hosts. Western corn rootworm exemplifies how an herbivore-tolerant crop was left vulnerable when it was bred for yield and protected with insecticides. Cotton fleahopper illustrates how removing preferred wild host plants from landscapes and replacing them with crops, allows herbivores with flexible host preferences to reach pest-level populations. We conclude by arguing that in the new geological epoch we face, the Anthropocene, we can improve agriculture by looking to our past to identify and avoid missteps of early and recent farmers.

Independent introductions and admixtures have contributed to adaptation of European maize and its American counterparts

Through the local selection of landraces, humans have guided the adaptation of crops to a vast range of climatic and ecological conditions. This is particularly true of maize, which was domesticated in a restricted area of Mexico but now displays one of the broadest cultivated ranges worldwide. Here, we sequenced 67 genomes with an average sequencing depth of 18x to document routes of introduction, admixture and selective history of European maize and its American counterparts. To avoid the confounding effects of recent breeding, we targeted germplasm (lines) directly derived from landraces. Among our lines, we discovered 22,294,769 SNPs and between 0.9% to 4.1% residual heterozygosity. Using a segmentation method, we identified 6,978 segments of unexpectedly high rate of heterozygosity. These segments point to genes potentially involved in inbreeding depression, and to a lesser extent to the presence of structural variants. Genetic structuring and inferences of historical splits revealed 5 genetic groups and two independent European introductions, with modest bottleneck signatures. Our results further revealed admixtures between distinct sources that have contributed to the establishment of 3 groups at intermediate latitudes in North America and Europe. We combined differentiation- and diversity-based statistics to identify both genes and gene networks displaying strong signals of selection. These include genes/gene networks involved in flowering time, drought and cold tolerance, plant defense and starch properties. Overall, our results provide novel insights into the evolutionary history of European maize and highlight a major role of admixture in environmental adaptation, paralleling recent findings in humans.

Genomic insights into divergence and dual domestication of cultivated allotetraploid cottons

Background

Cotton has been cultivated and used to make fabrics for at least 7000 years. Two allotetraploid species of great commercial importance, Gossypium hirsutum and Gossypium barbadense, were domesticated after polyploidization and are cultivated worldwide. Although the overall genetic diversity between these two cultivated species has been studied with limited accessions, their population structure and genetic variations remain largely unknown.

Results

We resequence the genomes of 147 cotton accessions, including diverse wild relatives, landraces, and modern cultivars, and construct a comprehensive variation map to provide genomic insights into the divergence and dual domestication of these two important cultivated tetraploid cotton species. Phylogenetic analysis shows two divergent groups for G. hirsutum and G. barbadense, suggesting a dual domestication processes in tetraploid cottons. In spite of the strong genetic divergence, a small number of interspecific reciprocal introgression events are found between these species and the introgression pattern is significantly biased towards the gene flow from G. hirsutum into G. barbadense. We identify selective sweeps, some of which are associated with relatively highly expressed genes for fiber development and seed germination.

Conclusions

We report a comprehensive analysis of the evolution and domestication history of allotetraploid cottons based on the whole genomic variation between G. hirsutum and G. barbadense and between wild accessions and modern cultivars. These results provide genomic bases for improving cotton production and for further evolution analysis of polyploid crops.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-017-1167-5) contains supplementary material, which is available to authorized users.

Diversity analysis of cotton (Gossypium hirsutum L.) germplasm using the CottonSNP63K Array

BACKGROUND

Cotton germplasm resources contain beneficial alleles that can be exploited to develop germplasm adapted to emerging environmental and climate conditions. Accessions and lines have traditionally been characterized based on phenotypes, but phenotypic profiles are limited by the cost, time, and space required to make visual observations and measurements. With advances in molecular genetic methods, genotypic profiles are increasingly able to identify differences among accessions due to the larger number of genetic markers that can be measured. A combination of both methods would greatly enhance our ability to characterize germplasm resources. Recent efforts have culminated in the identification of sufficient SNP markers to establish high-throughput genotyping systems, such as the CottonSNP63K array, which enables a researcher to efficiently analyze large numbers of SNP markers and obtain highly repeatable results. In the current investigation, we have utilized the SNP array for analyzing genetic diversity primarily among cotton cultivars, making comparisons to SSR-based phylogenetic analyses, and identifying loci associated with seed nutritional traits.

RESULTS

The SNP markers distinctly separated G. hirsutum from other Gossypium species and distinguished the wild from cultivated types of G. hirsutum. The markers also efficiently discerned differences among cultivars, which was the primary goal when designing the CottonSNP63K array. Population structure within the genus compared favorably with previous results obtained using SSR markers, and an association study identified loci linked to factors that affect cottonseed protein content.

CONCLUSIONS

Our results provide a large genome-wide variation data set for primarily cultivated cotton. Thousands of SNPs in representative cotton genotypes provide an opportunity to finely discriminate among cultivated cotton from around the world. The SNPs will be relevant as dense markers of genome variation for association mapping approaches aimed at correlating molecular polymorphisms with variation in phenotypic traits, as well as for molecular breeding approaches in cotton.

Genome-wide divergence, haplotype distribution and population demographic histories for Gossypium hirsutum and Gossypium barbadense as revealed by genome-anchored SNPs

Use of 10,129 singleton SNPs of known genomic location in tetraploid cotton provided unique opportunities to characterize genome-wide diversity among 440 Gossypium hirsutum and 219 G. barbadense cultivars and landrace accessions of widespread origin. Using the SNPs distributed genome-wide, we examined genetic diversity, haplotype distribution and linkage disequilibrium patterns in the G. hirsutum and G. barbadense genomes to clarify population demographic history. Diversity and identity-by-state analyses have revealed little sharing of alleles between the two cultivated allotetraploid genomes, with a few exceptions that indicated sporadic gene flow. We found a high number of new alleles, representing increased nucleotide diversity, on chromosomes 1 and 2 in cultivated G. hirsutum as compared with low nucleotide diversity on these chromosomes in landrace G. hirsutum. In contrast, G. barbadense chromosomes showed negative Tajima's D on several chromosomes for both cultivated and landrace types, which indicate that speciation of G. barbadense itself, might have occurred with relatively narrow genetic diversity. The presence of conserved linkage disequilibrium (LD) blocks and haplotypes between G. hirsutum and G. barbadense provides strong evidence for comparable patterns of evolution in their domestication processes. Our study illustrates the potential use of population genetic techniques to identify genomic regions for domestication.

DNA Sequence Evolution and Rare Homoeologous Conversion in Tetraploid Cotton

Allotetraploid cotton species are a vital source of spinnable fiber for textiles. The polyploid nature of the cotton genome raises many evolutionary questions as to the relationships between duplicated genomes. We describe the evolution of the cotton genome (SNPs and structural variants) with the greatly improved resolution of 34 deeply re-sequenced genomes. We also explore the evolution of homoeologous regions in the AT- and DT-genomes and especially the phenomenon of conversion between genomes. We did not find any compelling evidence for homoeologous conversion between genomes. These findings are very different from other recent reports of frequent conversion events between genomes. We also identified several distinct regions of the genome that have been introgressed between G. hirsutum and G. barbadense, which presumably resulted from breeding efforts targeting associated beneficial alleles. Finally, the genotypic data resulting from this study provides access to a wealth of diversity sorely needed in the narrow germplasm of cotton cultivars.