Cotton Breeding in Australia: Meeting the Challenges of the 21st Century.

The Commonwealth Scientific and Industrial Research Organisation (CSIRO) cotton breeding program is the sole breeding effort for cotton in Australia, developing high performing cultivars for the local industry which is worth∼AU$3 billion per annum. The program is supported by Cotton Breeding Australia, a Joint Venture between CSIRO and the program's commercial partner, Cotton Seed Distributors Ltd. (CSD). While the Australian industry is the focus, CSIRO cultivars have global impact in North America, South America, and Europe. The program is unique compared with many other public and commercial breeding programs because it focuses on diverse and integrated research with commercial outcomes. It represents the full research pipeline, supporting extensive long-term fundamental molecular research; native and genetically modified (GM) trait development; germplasm enhancement focused on yield and fiber quality improvements; integration of third-party GM traits; all culminating in the release of new commercial cultivars. This review presents evidence of past breeding successes and outlines current breeding efforts, in the areas of yield and fiber quality improvement, as well as the development of germplasm that is resistant to pests, diseases and abiotic stressors. The success of the program is based on the development of superior germplasm largely through field phenotyping, together with strong commercial partnerships with CSD and Bayer CropScience. These relationships assist in having a shared focus and ensuring commercial impact is maintained, while also providing access to markets, traits, and technology. The historical successes, current foci and future requirements of the CSIRO cotton breeding program have been used to develop a framework designed to augment our breeding system for the future. This will focus on utilizing emerging technologies from the genome to phenome, as well as a panomics approach with data management and integration to develop, test and incorporate new technologies into a breeding program. In addition to streamlining the breeding pipeline for increased genetic gain, this technology will increase the speed of trait and marker identification for use in genome editing, genomic selection and molecular assisted breeding, ultimately producing novel germplasm that will meet the coming challenges of the 21st Century.

Correction to: Improved lint yield under field conditions in cotton over-expressing transcription factors regulating fibre initiation.

Due to an unfortunate misunderstanding, an extra middle initial erroneously appeared in the original publication and the full name of the first author should read Shi Ming Liu.

Improved lint yield under field conditions in cotton over-expressing transcription factors regulating fibre initiation.

Only a few transcription factors (TFs) regulating which cells of the ovule epidermis differentiate into lint fibres have been identified in cotton (Gossypium hirsutum L.). In this study, the effect on lint yield and fibre quality of over-expressing three TFs in cotton, GhHD-1, GhMYB25 and GhMYB25Like, and their double and triple combinations, were evaluated in field experiments over two seasons. The expression of single or stacked TFs were all driven either by an ovule-specific promoter, FBP 7, or a constitutive promoter, Stunt 7, in a Coker 315 background. TF type, either singly or in combination, was found to be the most significant factor affecting lint yield. Among 64 transgenic lines tested, seven were higher yielding than null segregant lines in one or both seasons and were all from the sets with single and double over-expressed TF combinations. A reduced yield was associated with the set of triple combinations. The two most stable high yielding lines across the seasons recorded 12-22% higher yields than the nulls, although were not competitive to locally adapted commercial controls. Over-expression of TFs singly or in combination did not significantly alter fibre length and strength, but sometimes increased fibre micronaire. There were positive relationships between lint yield and lint percentage and lint yield and fibre density amongst the transgenic lines. Our preliminary results suggest that manipulating TF expression, either singly or in pairs, can increase the density of fibres initiated on developing seeds and fibre yields under field conditions while maintaining overall fibre quality.

Expansion of MIR482/2118 by a class-II transposable element in cotton.

Some plant microRNA (miRNA) families contain multiple members generating identical or highly similar mature miRNA variants. Mechanisms underlying the expansion of miRNA families remain elusive, although tandem and/or segmental duplications have been proposed. In this study of two tetraploid cottons, Gossypium hirsutum and Gossypium barbadense, and their extant diploid progenitors, Gossypium arboreum and Gossypium raimondii, we investigated the gain and loss of members of the miR482/2118 superfamily, which modulates the expression of nucleotide-binding site leucine-rich repeat (NBS-LRR) disease resistance genes. We found significant expansion of MIR482/2118d in G. barbadense, G. hirsutum and G. raimondii, but not in G. arboreum. Several newly expanded MIR482/2118d loci have mutated to produce different miR482/2118 variants with altered target-gene specificity. Based on detailed analysis of sequences flanking these MIR482/2118 loci, we found that this expansion of MIR482/2118d and its derivatives resulted from an initial capture of an MIR482/2118d by a class-II DNA transposable element (TE) in G. raimondii prior to the tetraploidization event, followed by transposition to new genomic locations in G. barbadense, G. hirsutum and G. raimondii. The 'GosTE' involved in the capture and proliferation of MIR482/2118d and its derivatives belongs to the PIF/Harbinger superfamily, generating a 3-bp target site duplication upon insertion at new locations. All orthologous MIR482/2118 loci in the two diploids were retained in the two tetraploids, but mutation(s) in miR482/2118 were observed across all four species as well as in different cultivars of both G. barbadense and G. hirsutum, suggesting a dynamic co-evolution of miR482/2118 and its NBS-LRR targets. Our results provide fresh insights into the mechanisms contributing to MIRNA proliferation and enrich our knowledge on TEs.

Trichomes at the Base of the Petal Are Regulated by the Same Transcription Factors as Cotton Seed Fibers.

Many polypetalous plants have a constriction at the base of the petal that leaves a small gap that can provide entry into the young flower bud before the reproductive organs are fully developed. In cotton (Gossypium hirsutum L.), this gap is occluded by tufts of short unicellular trichomes superficially resembling the fibers found on cotton seeds. We are just beginning to understand the developmental regulation of the seed fibers and have previously characterized several MIXTA-like MYB transcription factors (TFs) that are critical for correct seed fiber development but know little about the molecular regulation of other types of cotton trichomes. Here, using RNAi or dominant suppression transgenic cotton lines and natural fiber mutants, we investigated the development and regulation of the petal base trichomes. Petal base trichomes and seed trichomes were also examined across several different species within and outside of the Malvoideae. We found that the petal base trichomes are regulated by the same MYB TFs as cotton seed fibers and, since they are more widely distributed across different taxa than the seed fibers, could have preceded them in the evolution of these important textile fibers produced by some cotton species.

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.

GbEXPATR, a species-specific expansin, enhances cotton fibre elongation through cell wall restructuring.

Cotton provides us the most important natural fibre. High fibre quality is the major goal of cotton breeding, and introducing genes conferring longer, finer and stronger fibre from Gossypium barbadense to Gossypium hirsutum is an important breeding strategy. We previously analysed the G. barbadense fibre development mechanism by gene expression profiling and found two homoeologous fibre-specific α-expansins from G. barbadense, GbEXPA2 and GbEXPATR. GbEXPA2 (from the DT genome) is a classical α-expansin, while its homoeolog, GbEXPATR (AT genome), encodes a truncated protein lacking the normal C-terminal polysaccharide-binding domain of other α-expansins and is specifically expressed in G. barbadense. Silencing EXPA in G. hirsutum induced shorter fibres with thicker cell walls. GbEXPA2 overexpression in G. hirsutum had no effect on mature fibre length, but produced fibres with a slightly thicker wall and increased crystalline cellulose content. Interestingly, GbEXPATR overexpression resulted in longer, finer and stronger fibres coupled with significantly thinner cell walls. The longer and thinner fibre was associated with lower expression of a number of secondary wall-associated genes, especially chitinase-like genes, and walls with lower cellulose levels but higher noncellulosic polysaccharides which advocated that a delay in the transition to secondary wall synthesis might be responsible for better fibre. In conclusion, we propose that α-expansins play a critical role in fibre development by loosening the cell wall; furthermore, a truncated form, GbEXPATR, has a more dramatic effect through reorganizing secondary wall synthesis and metabolism and should be a candidate gene for developing G. hirsutum cultivars with superior fibre quality.

Inheritance and QTL Mapping of Leaf Nutrient Concentration in a Cotton Inter-Specific Derived RIL Population.

Developing and deploying cotton cultivars with high nutrient uptake, use efficiency and tolerance to nutrient related soil stresses is desirable to assist sustainable soil management. Genetic variation, heritability, selection response and quantitative trait loci (QTLs) were investigated for five macronutrients (P, K, Ca, Mg, S) and five micronutrients (Fe, Mn, B, Zn, and Cu) in a recombinant inbred line (RIL) population from an inter-specific cross between Gossypium hirsutum cv. Guazuncho 2, and G. barbadense accession VH8-4602. Na and K/Na ratio were also studied as the imbalance between Na and other nutrients is detrimental to cotton growth and development. The concentrations of nutrients were measured for different plant parts of the two parents and for leaf samples of the whole population collected at early to peak flowering in field experiments over two years in a sodic Vertosol soil. Parental contrast was large for most nutrient concentrations in leaves when compared with other plant parts. Segregation for leaf nutrient concentration was observed within the population with transgression for P, K, K/Na ratio and all micronutrients. Genotypic difference was the major factor behind within-population variation for most nutrients, while narrow sense heritability was moderate (0.27 for Mn and Cu, and 0.43 for B). At least one significant QTL was identified for each nutrient except K and more than half of those QTLs were clustered on chromosomes 14, 18 and 22. Selection response was predicted to be low for P and all micronutrients except B, high for K, Na and B, and very high for K/Na ratio. Correlations were more common between macronutrients, Na and K/Na ratio where the nature and strength of the relations varied (r=-0.69 to 0.76). We conclude that there is sufficient genetic diversity between these two tetraploid cotton species that could be exploited to improve cotton nutrient status by introgressing species-unique favourable alleles.

Development of a 63K SNP Array for Cotton and High-Density Mapping of Intraspecific and Interspecific Populations of Gossypium spp.

High-throughput genotyping arrays provide a standardized resource for plant breeding communities that are useful for a breadth of applications including high-density genetic mapping, genome-wide association studies (GWAS), genomic selection (GS), complex trait dissection, and studying patterns of genomic diversity among cultivars and wild accessions. We have developed the CottonSNP63K, an Illumina Infinium array containing assays for 45,104 putative intraspecific single nucleotide polymorphism (SNP) markers for use within the cultivated cotton species Gossypium hirsutum L. and 17,954 putative interspecific SNP markers for use with crosses of other cotton species with G. hirsutum. The SNPs on the array were developed from 13 different discovery sets that represent a diverse range of G. hirsutum germplasm and five other species: G. barbadense L., G. tomentosum Nuttal × Seemann, G. mustelinum Miers × Watt, G. armourianum Kearny, and G. longicalyx J.B. Hutchinson and Lee. The array was validated with 1,156 samples to generate cluster positions to facilitate automated analysis of 38,822 polymorphic markers. Two high-density genetic maps containing a total of 22,829 SNPs were generated for two F2 mapping populations, one intraspecific and one interspecific, and 3,533 SNP markers were co-occurring in both maps. The produced intraspecific genetic map is the first saturated map that associates into 26 linkage groups corresponding to the number of cotton chromosomes for a cross between two G. hirsutum lines. The linkage maps were shown to have high levels of collinearity to the JGI G. raimondii Ulbrich reference genome sequence. The CottonSNP63K array, cluster file and associated marker sequences constitute a major new resource for the global cotton research community.

Sequencing of allotetraploid cotton (Gossypium hirsutum L. acc. TM-1) provides a resource for fiber improvement.

Upland cotton is a model for polyploid crop domestication and transgenic improvement. Here we sequenced the allotetraploid Gossypium hirsutum L. acc. TM-1 genome by integrating whole-genome shotgun reads, bacterial artificial chromosome (BAC)-end sequences and genotype-by-sequencing genetic maps. We assembled and annotated 32,032 A-subgenome genes and 34,402 D-subgenome genes. Structural rearrangements, gene loss, disrupted genes and sequence divergence were more common in the A subgenome than in the D subgenome, suggesting asymmetric evolution. However, no genome-wide expression dominance was found between the subgenomes. Genomic signatures of selection and domestication are associated with positively selected genes (PSGs) for fiber improvement in the A subgenome and for stress tolerance in the D subgenome. This draft genome sequence provides a resource for engineering superior cotton lines.

DNA demethylases target promoter transposable elements to positively regulate stress responsive genes in Arabidopsis.

BACKGROUND: DNA demethylases regulate DNA methylation levels in eukaryotes. Arabidopsis encodes four DNA demethylases, DEMETER (DME), REPRESSOR OF SILENCING 1 (ROS1), DEMETER-LIKE 2 (DML2), and DML3. While DME is involved in maternal specific gene expression during seed development, the biological function of the remaining DNA demethylases remains unclear. RESULTS: We show that ROS1, DML2, and DML3 play a role in fungal disease resistance in Arabidopsis. A triple DNA demethylase mutant, rdd (ros1 dml2 dml3), shows increased susceptibility to the fungal pathogen Fusarium oxysporum. We identify 348 genes differentially expressed in rdd relative to wild type, and a significant proportion of these genes are downregulated in rdd and have functions in stress response, suggesting that DNA demethylases maintain or positively regulate the expression of stress response genes required for F. oxysporum resistance. The rdd-downregulated stress response genes are enriched for short transposable element sequences in their promoters. Many of these transposable elements and their surrounding sequences show localized DNA methylation changes in rdd, and a general reduction in CHH methylation, suggesting that RNA-directed DNA methylation (RdDM), responsible for CHH methylation, may participate in DNA demethylase-mediated regulation of stress response genes. Many of the rdd-downregulated stress response genes are downregulated in the RdDM mutants nrpd1 and nrpe1, and the RdDM mutants nrpe1 and ago4 show enhanced susceptibility to F. oxysporum infection. CONCLUSIONS: Our results suggest that a primary function of DNA demethylases in plants is to regulate the expression of stress response genes by targeting promoter transposable element sequences.

Genetic and DNA methylation changes in cotton (Gossypium) genotypes and tissues.

In plants, epigenetic regulation is important in normal development and in modulating some agronomic traits. The potential contribution of DNA methylation mediated gene regulation to phenotypic diversity and development in cotton was investigated between cotton genotypes and various tissues. DNA methylation diversity, genetic diversity, and changes in methylation context were investigated using methylation-sensitive amplified polymorphism (MSAP) assays including a methylation insensitive enzyme (BsiSI), and the total DNA methylation level was measured by high-performance liquid chromatography (HPLC). DNA methylation diversity was greater than the genetic diversity in the selected cotton genotypes and significantly different levels of DNA methylation were identified between tissues, including fibre. The higher DNA methylation diversity (CHG methylation being more diverse than CG methylation) in cotton genotypes suggest epigenetic regulation may be important for cotton, and the change in DNA methylation between fibre and other tissues hints that some genes may be epigenetically regulated for fibre development. The novel approach using BsiSI allowed direct comparison between genetic and epigenetic diversity, and also measured CC methylation level that cannot be detected by conventional MSAP.

Pectin methylesterase and pectin remodelling differ in the fibre walls of two gossypium species with very different fibre properties.

Pectin, a major component of the primary cell walls of dicot plants, is synthesized in Golgi, secreted into the wall as methylesters and subsequently de-esterified by pectin methylesterase (PME). Pectin remodelling by PMEs is known to be important in regulating cell expansion in plants, but has been poorly studied in cotton. In this study, genome-wide analysis showed that PMEs are a large multi-gene family (81 genes) in diploid cotton (Gossypium raimondii), an expansion over the 66 in Arabidopsis and suggests the evolution of new functions in cotton. Relatively few PME genes are expressed highly in fibres based on EST abundance and the five most abundant in fibres were cloned and sequenced from two cotton species. Their significant sequence differences and their stage-specific expression in fibres within a species suggest sub-specialisation during fibre development. We determined the transcript abundance of the five fibre PMEs, total PME enzyme activity, pectin content and extent of de-methylesterification of the pectin in fibre walls of the two cotton species over the first 25-30 days of fibre growth. There was a higher transcript abundance of fibre-PMEs and a higher total PME enzyme activity in G. barbadense (Gb) than in G. hirsutum (Gh) fibres, particularly during late fibre elongation. Total pectin was high, but de-esterified pectin was low during fibre elongation (5-12 dpa) in both Gh and Gb. De-esterified pectin levels rose thereafter when total PME activity increased and this occurred earlier in Gb fibres resulting in a lower degree of esterification in Gb fibres between 17 and 22 dpa. Gb fibres are finer and longer than those of Gh, so differences in pectin remodelling during the transition to wall thickening may be an important factor in influencing final fibre diameter and length, two key quality attributes of cotton fibres.

Epidermal cell differentiation in cotton mediated by the homeodomain leucine zipper gene, GhHD-1.

Gossypium hirsutum L. (cotton) fibres are specialized trichomes a few centimetres in length that grow from the seed coat. Few genes directly involved in the differentiation of these epidermal cells have been identified. These include GhMYB25-like and GhMYB25, two related MYB transcription factors that regulate fibre cell initiation and expansion. We have also identified a putative homeodomain leucine zipper (HD-ZIP) transcription factor, GhHD-1, expressed in trichomes and early fibres that might play a role in cotton fibre initiation. Here, we characterize GhHD-1 homoeologues from tetraploid G. hirsutum and show, using reporter constructs and quantitative real-time PCR (qRT-PCR), that they are expressed predominantly in epidermal tissues during early fibre development, and in other tissues bearing epidermal trichomes. Silencing of GhHD-1 reduced trichome formation and delayed the timing of fibre initiation. Constitutive overexpression of GhHD-1 increased the number of fibres initiating on the seed, but did not affect leaf trichomes. Expression of GhHD-1 in cotton silenced for different fibre MYBs suggest that in ovules it acts downstream of GhMYB25-like, but is unaffected in GhMYB25- or GhMYB109-silenced plants. Microarray analysis of silencing and overexpression lines of GhHD-1 indicated that it potentially regulates the levels of ethylene and reactive oxidation species (ROS) through a WRKY transcription factor and calcium-signalling pathway genes to activate downstream genes necessary for cell expansion and elongation.

Overexpression of a potato sucrose synthase gene in cotton accelerates leaf expansion, reduces seed abortion, and enhances fiber production.

Sucrose synthase (Sus) is a key enzyme in the breakdown of sucrose and is considered a biochemical marker for sink strength, especially in crop species, based on mutational and gene suppression studies. It remains elusive, however, whether, or to what extent, increase in Sus activity may enhance sink development. We aimed to address this question by expressing a potato Sus gene in cotton where Sus expression has been previously shown to be critical for normal seed and fiber development. Segregation analyses at T1 generation followed by studies in homozygous progeny lines revealed that increased Sus activity in cotton (1) enhanced leaf expansion with the effect evident from young leaves emerging from shoot apex; (2) improved early seed development, which reduced seed abortion, hence enhanced seed set, and (3) promoted fiber elongation. In young leaves of Sus overexpressing lines, fructose concentrations were significantly increased whereas, in elongating fibers, both fructose and glucose levels were increased. Since hexoses contribute little to osmolality in leaves, in contrast to developing fibers, it is concluded that high Sus activity promotes leaf development independently of osmotic regulation, probably through sugar signaling. The analyses also showed that doubling the Sus activity in 0-d cotton seeds increased their fresh weight by about 30%. However, further increase in Sus activity did not lead to any further increase in seed weight, indicating an upper limit for the Sus overexpression effect. Finally, based on the observed additive effect on fiber yield from increased fiber length and seed number, a new strategy is proposed to increase cotton fiber yield by improving seed development as a whole, rather than solely focusing on manipulating fiber growth.

GhMYB25-like: a key factor in early cotton fibre development.

MYB transcription factors have been implicated in regulation of the development of ovule epidermal cells into the elongated seed fibres of cotton. An R2R3 MYB, GhMYB25-like, identified from its reduced expression in a fibreless mutant of cotton (Xu142 fl), is here shown to play a key role in the very early stages of fibre cell differentiation. A GhMYB25-like promoter-GUS construct was expressed predominantly in the epidermal layers of cotton ovules before anthesis (-3days post-anthesis, dpa), increasing in expression in 0-dpa ovules, primarily in those epidermal cells expanding into fibres, and then in elongating fibres at +3dpa, declining thereafter. This was consistent with GhMYB25-like transcript abundance during fibre development. RNA interference suppression of GhMYB25-like resulted in cotton plants with fibreless seeds, but normal trichomes elsewhere, phenocopying the Xu142 fl mutant. Like Xu142 fl these plants had reduced expression of the fibre-expressed MYBs, GhMYB25 and GhMYB109, indicating that GhMYB25-like is upstream from those MYBs. This hierarchy was supported by the absence of any change in transcript level of GhMYB25-like in GhMYB25- and GhMYB109-silenced transgenic lines. Transgenic cotton with an additional copy of the native gene had elevated expression of GhMYB25-like in ovules, but no obvious increase in fibre initials, suggesting that there may be other factors that interact with GhMYB25-like to differentiate epidermal cells into fibre cells.

Comparisons of early transcriptome responses to low-oxygen environments in three dicotyledonous plant species.

Waterlogging is a serious impediment to crop productivity worldwide which acts to reduce oxygen levels in the rhizosphere due to the low diffusion rate of molecular oxygen in water. Plants respond to low oxygen through rapid and specific changes at both the transcriptional and translational levels. Transcriptional changes to low-oxygen (hypoxia) stress have been studied in a number of plant species using whole genome microarrays. Using transcriptome data from root tissue from early time points (4-5 h) from cotton (Gossypium hirsutum), Arabidopsis and gray poplar (Populus x canescens), we have identified a core set of orthologous genes that responded to hypoxia in similar ways between species, and others that showed species specific responses . Responses to hypoxia were most similar between Arabidopsis and cotton, while the waterlogging tolerant poplar species exhibited some significant differences.