The cotton genus (Gossypium) provides a superior system for the study of diversification, genome evolution, polyploidization, and human-mediated selection. To gain insight into phenotypic diversification in cotton seeds, we conducted coexpression network analysis of developing seeds from diploid and allopolyploid cotton species and explored network properties. Key network modules and functional associations were identified related to seed oil content and seed weight. We compared species-specific networks to reveal topological changes, including rewired edges and differentially coexpressed genes, associated with speciation, polyploidy, and cotton domestication. Network comparisons among species indicate that topologies are altered in addition to gene expression profiles, indicating that changes in transcriptomic coexpression relationships play a role in the developmental architecture of cotton seed development. The global network topology of allopolyploids, especially for domesticated G. hirsutum, resembles the network of the A-genome diploid more than that of the D-genome parent, despite its D-like phenotype in oil content. Expression modifications associated with allopolyploidy include coexpression level dominance and transgressive expression, suggesting that the transcriptomic architecture in polyploids is to some extent a modular combination of that of its progenitor genomes. Among allopolyploids, intermodular relationships are more preserved between two different wild allopolyploid species than they are between wild and domesticated forms of a cultivated cotton, and regulatory connections of oil synthesis-related pathways are denser and more closely clustered in domesticated vs. wild G. hirsutum. These results demonstrate substantial modification of genic coexpression under domestication. Our work demonstrates how network inference informs our understanding of the transcriptomic architecture of phenotypic variation associated with temporal scales ranging from thousands (domestication) to millions (speciation) of years, and by polyploidy.

Cotton ranks among the world's important oilseed crops, yet relative to other oilseeds there are few studies of oil-related biosynthetic and regulatory pathways. We present global transcriptome analyses of cotton seed development using RNA-seq and four developmental time-points. Because Upland cotton (Gossypium hirsutum L.) is an allopolyploid containing two genomes (A/D), we partitioned expression into the individual contributions of each homeologous gene copy. Data were explored with respect to genic and subgenomic patterns of expression, globally and with respect to seed pathways and networks. The most dynamic period of transcriptome change is from 20-30 d postanthesis (DPA), with about 20% of genes showing homeolog expression bias. Co-expression analysis shows largely congruent homeolog networks, but also homeolog-specific divergence. Functional enrichment tests show that flavonoid biosynthesis and lipid related genes were significantly represented early and later in seed development, respectively. An involvement of new features in oil biosynthesis was found, like the contribution of DGAT3 (diacylglycerol acyltransferase) to the total triglyceride expression pool. Also, catechin-based and epicatechin-based proanthocyanidin expression are reciprocally biased with respect to homeolog usage. This study provides the first temporal analysis of duplicated gene expression in cotton seed and a resource for understanding new aspects of oil and flavonoid biosynthetic processes.

The growth habit of lateral shoots (also termed “branching habit”) is an important descriptive and agronomic character of peanut. Yet, both the inheritance of branching habit and the genetic mechanism that controls it in this crop remain unclear. In addition, the low degree of polymorphism among cultivated peanut varieties hinders fine-mapping of this and other traits in non-homozygous genetic structures. Here, we combined high-throughput sequencing with a well-defined genetic system to study these issues in peanut. Initially, segregating F2 populations derived from a reciprocal cross between very closely related Virginia-type peanut cultivars with spreading and bunch growth habits were examined. The spreading/bunch trait was shown to be controlled by a single gene with no cytoplasmic effect. That gene was named Bunch1 and was significantly correlated with pod yield per plant, time to maturation and the ratio of “dead-end” pods. Subsequently, bulked segregant analysis was performed on 52 completely bunch, and 47 completely spreading F3 families. In order to facilitate the process of SNP detection and candidate-gene analysis, the transcriptome was used instead of genomic DNA. Young leaves were sampled and bulked. Reads from Illumina sequencing were aligned against the peanut reference transcriptome and the diploid genomes. Inter-varietal SNPs were detected, scored and quality-filtered. Thirty-four candidate SNPs were found to have a bulk frequency ratio value >10 and 6 of those SNPs were found to be located in the genomic region of linkage group B5. Three best hits from that over-represented region were further analyzed in the segregating population. The trait locus was found to be located in a ~1.1 Mbp segment between markers M875 (B5:145,553,897; 1.9 cM) and M255 (B5:146,649,943; 2.25 cM). The method was validated using a population of recombinant inbreed lines of the same cross and a new DNA SNP-array. This study demonstrates the relatively straight-forward utilization of bulk segregant analysis for trait fine-mapping in the low polymeric and heterozygous germplasm of cultivated peanut and provides a baseline for candidate gene discovery and map-based cloning of Bunch1.

Temperature compensation, expressed as the ability to maintain clock characteristics (mainly period) in face of temperature changes, that is, robustness, is considered a key feature of circadian clock systems. In this study, we explore the genetic basis for lack of robustness, that is, plasticity, of circadian clock as reflected by photosynthesis rhythmicity. The clock rhythmicity of a new wild barley reciprocal doubled haploid population was analysed with a high temporal resolution of pulsed amplitude modulation of chlorophyll fluorescence under optimal (22°C) and high (32°C) temperature. This comparison between two environments pointed to the prevalence of clock acceleration under heat. Genotyping by sequencing of doubled haploid lines indicated a rich recombination landscape with minor fixation (less than 8%) for one of the parental alleles. Quantitative genetic analysis included genotype by environment interactions and binary‐threshold models. Variation in the circadian rhythm plasticity phenotypes, expressed as change (delta) of period and amplitude under two temperatures, was associated with maternal organelle genome (the plasmotype), as well as with several nuclear loci. This first reported rhythmicity driven by nuclear loci and plasmotype with few identified variants, paves the way for studying impact of cytonuclear variations on clock robustness and on plant adaptation to changing environments.

No abstract