Asymmetric subgenome selection and cis-regulatory divergence during cotton domestication

Extensive intraspecific gene order and gene structural variations in upland cotton cultivars

Multiple cotton genomes (diploid and tetraploid) have been assembled. However, genomic variations between cultivars of allotetraploid upland cotton (Gossypium hirsutum L.), the most widely planted cotton species in the world, remain unexplored. Here, we use single-molecule long read and Hi-C sequencing technologies to assemble genomes of the two upland cotton cultivars TM-1 and zhongmiansuo24 (ZM24). Comparisons among TM-1 and ZM24 assemblies and the genomes of the diploid ancestors reveal a large amount of genetic variations. Among them, the top three longest structural variations are located on chromosome A08 of the tetraploid upland cotton, which account for ~30% total length of this chromosome. Haplotype analyses of the mapping population derived from these two cultivars and the germplasm panel show suppressed recombination rates in this region. This study provides additional genomic resources for the community, and the identified genetic variations, especially the reduced meiotic recombination on chromosome A08, will help future breeding.

Dissection of the genetic variation and candidate genes of lint percentage by a genome-wide association study in upland cotton

A genome-wide associated study identified six novel QTLs for lint percentage. Two candidate genes underlying this trait were also detected. Increasing lint percentage (LP) is a core goal of cotton breeding. To better understand the genetic basis of LP, a genome-wide association study (GWAS) was conducted using 276 upland cotton accessions planted in multiple environments and genotyped with a CottonSNP63K array. After filtering, 10,660 high-quality single-nucleotide polymorphisms (SNPs) were retained. Population structure, principal component and neighbor-joining phylogenetic tree analyses divided the accessions into two subpopulations. These results along with linkage disequilibrium decay indicated accessions were not highly structured and exhibited weak relatedness. GWAS uncovered 23 polymorphic SNPs and 15 QTLs significantly associated with LP, with six new QTLs identified. Two candidate genes, Gh_D05G0313 and Gh_D05G1124, both contained one significant SNP, highly expressed during ovule and fiber development stages, implying that the two genes may act as the most promising regulators of LP. Furthermore, the phenotypic value of LP was found to be positively correlated with the number of favorable SNP alleles. These favorable alleles for LP identified in the study may be useful for improving lint yield.

Diversity of a wall-associated kinase gene in wild and cultivated barley

Domestication of barley and other cereals was accompanied by an increase in seed size which has been ascribed to human selection, large seeds being preferred by early farmers or favoured by cultivation practices such as deep sowing. An alternative suggestion is that the increase in seed size was an indirect consequence of selection for plants with more vigorous growth. To begin to address the latter hypothesis we studied the diversity of HvWAK1, a wall-associated kinase gene involved in root proliferation, in 220 wild barley accessions and 200 domesticated landraces. A 3655-bp sequence comprising the gene and upstream region contained 69 single nucleotide polymorphisms (SNPs), one indel and four short tandem repeats. A network of 50 haplotypes revealed a complex evolutionary relationship, but with landraces largely restricted to two parts of the topology. SNPs in the HvWAK1 coding region resulted in nonsynonymous substitutions at nine positions in the translation product, but none of these changes were predicted to have a significant effect on the protein structure. In contrast, the region upstream of the coding sequence contained five SNPs that were invariant in the domesticated population, fixation of these SNPs decreasing the likelihood that the upstream of a pair of TATA boxes and transcription start sites would be used to promote transcription of HvWAK1. The sequence diversity therefore suggests that the cis-regulatory region of HvWAK1 might have been subject to selection during barley domestication. The extent of root proliferation has been linked with traits such as above-ground biomass, so selection for particular cis-regulatory variants of HvWAK1 would be consistent with the hypothesis that seed size increases during domestication were the indirect consequence of selection for plants with increased growth vigour.

Genome-wide comparative analysis of RNA-binding Glycine-rich protein family genes between Gossypium arboreum and Gossypium raimondii

RB-GRP (RNA-binding Glycine-rich protein gene) family belongs to the fourth subfamily of the GRP (Glycine-rich protein gene) superfamily, which plays a great role in plant growth and development, as well as in abiotic stresses response, while it has not been identified in cotton. Here, we identified 37 and 32 RB-GRPs from two cotton species (Gossypium arboreum and Gossypium raimondii, respectively), which were divided into four distinct subfamilies based on the presence of additional motifs and the arrangement of the glycine repeats. The distribution of RB-GRPs was nonrandom and uneven among the chromosomes both in two cotton species. The expansion of RB-GRP gene family between two cultivars was mainly attributed to segmental and tandem duplication events indicated by synteny analysis, and the tandem duplicated genes were mapped into homologous collinear blocks, indicated that they shared a common ancestral gene in both species. Furthermore, most RB-GRPs in two cotton species undergone stronger negative selective pressure by evolutionary analysis of RB-GRP orthologous genes. Meanwhile, RB-GRPs participated in different abiotic stresses (Abscisic acid, salt and Polyethylene glycol) responses and tissues at different developmental stages between two cotton species were showed by gene expression analysis. This research would provide insight into the evolution and function of the RB-GRPs in Gossypium species.

Genome-wide analysis of genetic variations between dominant and recessive NILs of glanded and glandless cottons

Cotton is an important economic crop in worldwide. It produces fiber for the textile industry and provides cottonseeds with high-quality protein and oil. However, the presence of gossypol limits the utilization of cottonseed. Two pairs of cotton near isogenic lines (NILs) with different pigment glands, i.e., Coker 312 vs Coker 312 W and CCRI12 vs CCRI12W, exhibit different gossypol contents. The glandless traits of Coker 312 W and CCRI12W are controlled by recessive and dominant genes, respectively. However, knowledge regarding the genomic variations in the NILs is limited. Therefore, the NILs genomes were resequenced and the sequencing depths were greater than 34×. Compared with the TM-1 genome, numerous SNPs, Indels, SVs, and CNVs were discovered. KEGG pathway analysis revealed that genes with SNPs and Indels from the recessive NILs and genes with Indels from the dominant NILs shared only one enriched pathway, i.e., the sesquiterpenoid and triterpenoid biosynthesis pathway, which is relevant to gossypol biosynthesis. Expression analysis revealed that key genes with variations that participate in the gossypol biosynthesis and pigment gland formation pathways had different expression patterns among the dominant, recessive glandless and glanded plants. The expression levels in the glanded organs were higher than those in their NILs. Altogether, our results provide deeper insight into cotton NILs with different pigment glands.

Building Transcription Factor Binding Site Models to Understand Gene Regulation in Plants

Transcription factors (TFs) are key cellular components that control gene expression. They recognize specific DNA sequences, the TF binding sites (TFBSs), and thus are targeted to specific regions of the genome where they can recruit transcriptional co-factors and/or chromatin regulators to fine-tune spatiotemporal gene regulation. Therefore, the identification of TFBSs in genomic sequences and their subsequent quantitative modeling is of crucial importance for understanding and predicting gene expression. Here, we review how TFBSs can be determined experimentally, how the TFBS models can be constructed in silico, and how they can be optimized by taking into account features such as position interdependence within TFBSs, DNA shape, and/or by introducing state-of-the-art computational algorithms such as deep learning methods. In addition, we discuss the integration of context variables into the TFBS modeling, including nucleosome positioning, chromatin states, methylation patterns, 3D genome architectures, and TF cooperative binding, in order to better predict TF binding under cellular contexts. Finally, we explore the possibilities of combining the optimized TFBS model with technological advances, such as targeted TFBS perturbation by CRISPR, to better understand gene regulation, evolution, and plant diversity.

Whole genome sequencing reveals rare off-target mutations and considerable inherent genetic or/and somaclonal variations in CRISPR/Cas9-edited cotton plants

The CRISPR/Cas9 system has been extensively applied for crop improvement. However, our understanding of Cas9 specificity is very limited in Cas9-edited plants. To identify on- and off-target mutation in an edited crop, we described whole genome sequencing (WGS) of 14 Cas9-edited cotton plants targeted to three genes, and three negative (Ne) control and three wild-type (WT) plants. In total, 4188-6404 unique single-nucleotide polymorphisms (SNPs) and 312-745 insertions/deletions (indels) were detected in 14 Cas9-edited plants compared to WT, negative and cotton reference genome sequences. Since the majority of these variations lack a protospacer-adjacent motif (PAM), we demonstrated that the most variations following Cas9-edited are due either to somaclonal variation or/and pre-existing/inherent variation from maternal plants, but not off-target effects. Of a total of 4413 potential off-target sites (allowing ≤5 mismatches within the 20-bp sgRNA and 3-bp PAM sequences), the WGS data revealed that only four are bona fide off-target indel mutations, validated by Sanger sequencing. Moreover, inherent genetic variation of WT can generate novel off-target sites and destroy PAMs, which suggested great care should be taken to design sgRNA for the minimizing of off-target effect. These findings suggested that CRISPR/Cas9 system is highly specific for cotton plants.

Live-cell imaging of the cytoskeleton in elongating cotton fibres

Cotton (Gossypium hirsutum) fibres consist of single cells that grow in a highly polarized manner, assumed to be controlled by the cytoskeleton^1-3. However, how the cytoskeletal organization and dynamics underpin fibre development remains unexplored. Moreover, it is unclear whether cotton fibres expand via tip growth or diffuse growth^2-4. We generated stable transgenic cotton plants expressing fluorescent markers of the actin and microtubule cytoskeleton. Live-cell imaging revealed that elongating cotton fibres assemble a cortical filamentous actin network that extends along the cell axis to finally form actin strands with closed loops in the tapered fibre tip. Analyses of F-actin network properties indicate that cotton fibres have a unique actin organization that blends features of both diffuse and tip growth modes. Interestingly, typical actin organization and endosomal vesicle aggregation found in tip-growing cell apices were not observed in fibre tips. Instead, endomembrane compartments were evenly distributed along the elongating fibre cells and moved bi-directionally along the fibre shank to the fibre tip. Moreover, plus-end tracked microtubules transversely encircled elongating fibre shanks, reminiscent of diffusely growing cells. Collectively, our findings indicate that cotton fibres elongate via a unique tip-biased diffuse growth mode.

Deep Transcriptome Analysis Reveals Reactive Oxygen Species (ROS) Network Evolution, Response to Abiotic Stress, and Regulation of Fiber Development in Cotton

Reactive oxygen species (ROS) are important molecules in the plant, which are involved in many biological processes, including fiber development and adaptation to abiotic stress in cotton. We carried out transcription analysis to determine the evolution of the ROS genes and analyzed their expression levels in various tissues of cotton plant under abiotic stress conditions. There were 515, 260, and 261 genes of ROS network that were identified in Gossypium hirsutum (AD₁ genome), G. arboreum (A genome), and G. raimondii (D genome), respectively. The ROS network genes were found to be distributed in all the cotton chromosomes, but with a tendency of aggregating on either the lower or upper arms of the chromosomes. Moreover, all the cotton ROS network genes were grouped into 17 families as per the phylogenetic tress analysis. A total of 243 gene pairs were orthologous in G. arboreum and G. raimondii. There were 240 gene pairs that were orthologous in G. arboreum, G. raimondii, and G. hirsutum. The synonymous substitution value (K_s) peaks of orthologous gene pairs between the At subgenome and the A progenitor genome (G. arboreum), D subgenome and D progenitor genome (G. raimondii) were 0.004 and 0.015, respectively. The K_s peaks of ROS network orthologous gene pairs between the two progenitor genomes (A and D genomes) and two subgenomes (At and Dt subgenome) were 0.045. The majority of K_a/K_s value of orthologous gene pairs between the A, D genomes and two subgenomes of TM-1 were lower than 1.0. RNA seq. analysis and RT-qPCR validation, showed that, CSD1,2,3,5,6; FSD1,2; MSD1,2; APX3,11; FRO5.6; and RBOH6 played a major role in fiber development while CSD1, APX1, APX2, MDAR1, GPX4-6-7, FER2, RBOH6, RBOH11, and FRO5 were integral for enhancing salt stress in cotton. ROS network-mediated signal pathway enhances the mechanism of fiber development and regulation of abiotic stress in Gossypium. This study will enhance the understanding of ROS network and form the basic foundation in exploring the mechanism of ROS network-involving the fiber development and regulation of abiotic stress in cotton.

Dynamic control of enhancer activity drives stage-specific gene expression during flower morphogenesis

Enhancers are critical for developmental stage-specific gene expression, but their dynamic regulation in plants remains poorly understood. Here we compare genome-wide localization of H3K27ac, chromatin accessibility and transcriptomic changes during flower development in Arabidopsis. H3K27ac prevalently marks promoter-proximal regions, suggesting that H3K27ac is not a hallmark for enhancers in Arabidopsis. We provide computational and experimental evidence to confirm that distal DNase І hypersensitive sites are predictive of enhancers. The predicted enhancers are highly stage-specific across flower development, significantly associated with SNPs for flowering-related phenotypes, and conserved across crucifer species. Through the integration of genome-wide transcription factor (TF) binding datasets, we find that floral master regulators and stage-specific TFs are largely enriched at developmentally dynamic enhancers. Finally, we show that enhancer clusters and intronic enhancers significantly associate with stage-specific gene regulation by floral master TFs. Our study provides insights into the functional flexibility of enhancers during plant development, as well as hints to annotate plant enhancers.

Enhancer elements can control spatial and temporal patterns of gene expression. Here the authors profile chromatin accessibility, histone modifications and gene expression during Arabidopsis flower development providing evidence for sets of distal enhancers acting in a highly stage-specific manner.

Identifying favorable alleles for improving key agronomic traits in upland cotton

BACKGROUND

Gossypium hirsutum L. is grown worldwide and is the largest source of natural fiber crop. We focus on exploring the favorable alleles (FAs) for upland cotton varieties improvement, and further understanding the history of accessions selection and acumination of favorable allele during breeding.

RESULTS

The genetic basis of phenotypic variation has been studied. But the accumulation of favorable alleles in cotton breeding history in unknown, and potential favorable alleles to enhance key agronomic traits in the future cotton varieties have not yet been identified. Therefore, 419 upland cotton accessions were screened, representing a diversity of phenotypic variations of 7362 G. hirsutum, and 15 major traits were investigated in 6 environments. These accessions were categorized into 3 periods (early, medium, and modern) according to breeding history. All accessions were divided into two major groups using 299 polymorphic microsatellite markers: G1 (high fiber yield and quality, late maturity) and G2 (low fiber yield and quality, early maturity). The proportion of G1 genotype gradually increased from early to modern breeding periods. Furthermore, 21 markers (71 alleles) were significantly associated (-log P > 4) with 15 agronomic traits in multiple environments. Seventeen alleles were identified as FAs; these alleles accumulated more in the modern period than in other periods, consistent with their phenotypic variation trends in breeding history. Our results demonstrate that the favorable alleles accumulated through breeding effects, especially for common favorable alleles. However, the potential elite accessions could be rapidly screened by rare favorable alleles.

CONCLUSION

In our study, genetic variation and genome-wide associations for 419 upland cotton accessions were analyzed. Two favorable allele types were identified during three breeding periods, providing important information for yield/quality improvement of upland cotton germplasm.

Resequencing core accessions of a pedigree identifies derivation of genomic segments and key agronomic trait loci during cotton improvement

Upland cotton (Gossypium hirsutum) is the world's largest source of natural fibre and dominates the global textile industry. Hybrid cotton varieties exhibit strong heterosis that confers high fibre yields, yet the genome-wide effects of artificial selection that have influenced Upland cotton during its breeding history are poorly understood. Here, we resequenced Upland cotton genomes and constructed a variation map of an intact breeding pedigree comprising seven elite and 19 backbone parents. Compared to wild accessions, the 26 pedigree accessions underwent strong artificial selection during domestication that has resulted in reduced genetic diversity but stronger linkage disequilibrium and higher extents of selective sweeps. In contrast to the backbone parents, the elite parents have acquired significantly improved agronomic traits, with an especially pronounced increase in the lint percentage. Notably, identify by descent (IBD) tracking revealed that the elite parents inherited abundant beneficial trait segments and loci from the backbone parents and our combined analyses led to the identification of a core genomic segment which was inherited in the elite lines from the parents Zhong 7263 and Ejing 1 and that was strongly associated with lint percentage. Additionally, SNP correlation analysis of this core segment showed that a non-synonymous SNP (A-to-G) site in a gene encoding the cell wall-associated receptor-like kinase 3 (GhWAKL3) protein was highly correlated with increased lint percentage. Our results substantially increase the valuable genomics resources available for future genetic and functional genomics studies of cotton and reveal insights that will facilitate yield increases in the molecular breeding of cotton.

Suppression of tryptophan synthase activates cotton immunity by triggering cell death via promoting SA synthesis

Primary metabolism plays an important role in plant growth and development, however the relationship between primary metabolism and the adaptive immune response is largely unknown. Here, we employed RNA interference (RNAi), virus-induced gene silencing (VIGS) technology, phytohormone profiling, genetic studies, and transcriptome and metabolome analysis to investigate the function of the tryptophan synthesis pathway in the resistance of cotton to V. dahliae. We found that knock-down of GbTSA1 (Tryptophan Synthase α) and GbTSB1 (tryptophan synthase β) induced a spontaneous cell death phenotype in a salicylic acid (SA)-dependent manner and enhanced resistance to V. dahliae in cotton plants. Metabolome analysis showed that indole and indolic metabolites were highly accumulated in GbTSA1- or GbTSB1-silenced plants. Transcriptomic analysis showed that exogenous indole promotes the expression levels of genes involved in SA synthesis and the defense response. Similarly, indole application strongly enhanced cotton resistance to V. dahliae. These results suggested that metabolic intermediates in the Trp synthesis pathway may be a signal to activate SA synthesis. These results also provided a strategy to elicit plant defense responses by the application of indole.

3D genome organization: a role for phase separation and loop extrusion?

In eukaryotes, genomic information is encoded in chromosomes, which occupy distinct territories within the nucleus. Inside these territories, chromosomes are folded in a hierarchical set of topological structures, called compartments, topologically associated domains and loops. Phase separation and loop extrusion are the mechanisms indicated to mediate the 3D organization of the genome, and gene activity and epigenetic marks determine the activity level of the formed chromatin domains. The main difference between plants and animals may be the absence of canonical insulator elements in plants. Comparison across plant species indicates that the identification of chromatin domains is affected by genome size, gene density, and the linear distribution of genes and transposable elements.

Population-level approaches reveal novel aspects of lignin biosynthesis, content, composition and structure

Population-level studies enabled by high-throughput phenotyping have revealed significant variation in lignin characteristics including content, S:G:H ratio, inter-unit linkage distributions, and molecular weights across multiple plant species. Coupled with genome-wide association mapping studies (GWAS) targeted at linking genetic mutations to phenotype, significant progress has been made in associating putative causal mutations to variation in lignin characteristics. Despite this progress, there are few examples, in which these associations have been molecularly validated to provide new insights into the genetic regulation of lignin biosynthesis. Given a recent report of a GWAS-discovered 5-enolpyruvylshikimate-3-phosphate (EPSP) synthase moonlighting as a transcriptional regulator of lignin biosynthesis, the potential to bridge scientific disciplines in order to uncover hidden elements of lignin biosynthesis has been demonstrated, offering a path to alter lignin characteristics via genetic manipulation in order to expedite lignin valorization. To maximize this potential, however, there is a crucial need for (1) broader surveys of naturally varying diverse plant populations and (2) analytical platforms that can resolve subtle properties at fine chemical and biological scales.

QTL analysis for yield and fibre quality traits using three sets of introgression lines developed from three Gossypium hirsutum race stocks

Upland cotton (Gossypium hirsutum L.) race stocks may possess desirable traits for the genetic improvement of cotton. Quantitative trait locus (QTL) analysis can assist in uncovering new alleles from unadapted race stocks. In this study, three sets of chromosome segment introgression lines (ILs) were developed from three backcrosses (BC₃) between three race stocks, G. hirsutum races latifolium accs. TX-34 and TX-48 and punctatum acc. TX-114, as donor parents and Texas Marker-1 (TM-1) as the recurrent parent. Based on a total of 452 polymorphic simple sequence repeat (SSR) markers in BC₃F₂ genotyping, 149, 150 and 184 ILs were obtained from TM-1 × TX-34, TM-1 × TX-48 and TM-1 × TX-114, respectively. The average introgressed chromosomal segment length was 12.7 cM, and the total genetic distance was 3268 cM covering approximately 73.4% of the Upland cotton genome. The BC₃F₂, BC₃F_2:3 and BC₃F_2:4 progeny, which produced the ILs, were evaluated for yield and fibre quality traits. A total of 128 QTLs were detected, each of which explained 1.6-13.0% of the phenotypic variation. Thirty-five common QTLs related to eight traits were detected. Six QTL clusters were found on five chromosomes. Thirty-eight QTLs were previously unreported, and they may be footprints of cotton domestication. Domestication or artificial selection by humans successfully eliminated most unfavourable QTLs (21/38); however, some favourable QTLs (17/38) are not present in modern cultivars, demonstrating the importance of race stocks for improving cotton cultivars. The 26 elite ILs developed could be used to improve the yield and fibre quality components simultaneously. These results provide information on desirable QTLs for cotton improvement.

Whole-genome resequencing reveals Brassica napus origin and genetic loci involved in its improvement

Brassica napus (2n = 4x = 38, AACC) is an important allopolyploid crop derived from interspecific crosses between Brassica rapa (2n = 2x = 20, AA) and Brassica oleracea (2n = 2x = 18, CC). However, no truly wild B. napus populations are known; its origin and improvement processes remain unclear. Here, we resequence 588 B. napus accessions. We uncover that the A subgenome may evolve from the ancestor of European turnip and the C subgenome may evolve from the common ancestor of kohlrabi, cauliflower, broccoli, and Chinese kale. Additionally, winter oilseed may be the original form of B. napus. Subgenome-specific selection of defense-response genes has contributed to environmental adaptation after formation of the species, whereas asymmetrical subgenomic selection has led to ecotype change. By integrating genome-wide association studies, selection signals, and transcriptome analyses, we identify genes associated with improved stress tolerance, oil content, seed quality, and ecotype improvement. They are candidates for further functional characterization and genetic improvement of B. napus.

Brassica napus is a globally important oil crop, but the origin of the allotetraploid genome and its improvement process are largely unknown. Here, the authors take a population genetic approach to resolve its origin and evolutionary history, and identify candidate genes related to important agricultural traits.

TaZIM-A1 negatively regulates flowering time in common wheat (Triticum aestivum L.)

Flowering time is a critical determinant of regional adaptation for crops and has strong effects on crop yields. Here, we report that TaZIM-A1, an atypical GATA-like transcription factor, is a negative regulator of flowering in wheat. TaZIM-A1 possessed weak transcriptional repression activity, with its CCT domain functioning as the major inhibitory region. TaZIM-A1 expression exhibited a typical circadian oscillation pattern under various light regimes. Overexpression of TaZIM-A1 caused a delay in flowering time and a decrease in thousand-kernel weight (TKW) in wheat under long-day conditions. Moreover, TaZIM-A1 directly bound to the promoters of TaCO-1 and TaFT-1 and downregulated their expression. Sequence analysis of a collection of common wheat cultivars identified three and two haplotypes for TaZIM-A1 and TaZIM-B1, respectively. Association analysis revealed that TaZIM-A1-HapI/-HapIII and TaZIM-B1-HapI have undergone strong positive selection during modern wheat breeding, likely due to their association with earlier heading and higher TKW. Diagnostic markers were developed for these haplotypes that can be used for wheat cultivar improvement, via marker-assisted breeding.

Genomic analyses of an extensive collection of wild and cultivated accessions provide new insights into peach breeding history

BACKGROUND

Human selection has a long history of transforming crop genomes. Peach (Prunus persica) has undergone more than 5000 years of domestication that led to remarkable changes in a series of agronomically important traits, but genetic bases underlying these changes and the effects of artificial selection on genomic diversity are not well understood.

RESULTS

Here, we report a comprehensive analysis of peach evolution based on genome sequences of 480 wild and cultivated accessions. By focusing on a set of quantitative trait loci (QTLs), we provide evidence supporting that distinct phases of domestication and improvement have led to an increase in fruit size and taste and extended its geographic distribution. Fruit size was predominantly selected during domestication, and selection for large fruits has led to the loss of genetic diversity in several fruit weight QTLs. In contrast, fruit taste-related QTLs were successively selected for by domestication and improvement, with more QTLs selected for during improvement. Genome-wide association studies of 11 agronomic traits suggest a set of candidate genes controlling these traits and potential markers for molecular breeding. Candidate loci for genes that contributed to the adaption to low-chill regions were identified. Furthermore, the genomic bases of divergent selection for fruit texture and local breeding for different flavors between Asian and European/North American cultivars were also determined.

CONCLUSIONS

Our results elucidate the genetic basis of peach evolution and provide new resources for future genomics-guided peach breeding.

Exploring and Exploiting Pan-genomics for Crop Improvement

Genetic variation ranging from single-nucleotide polymorphisms to large structural variants (SVs) can cause variation of gene content among individuals within the same species. There is an increasing appreciation that a single reference genome is insufficient to capture the full landscape of genetic diversity of a species. Pan-genome analysis offers a platform to evaluate the genetic diversity of a species via investigation of its entire genome repertoire. Although a recent wave of pan-genomic studies has shed new light on crop diversity and improvement using advanced sequencing technology, the potential applications of crop pan-genomics in crop improvement are yet to be fully exploited. In this review, we highlight the progress achieved in understanding crop pan-genomics, discuss biological activities that cause SVs, review important agronomical traits affected by SVs, and present our perspective on the application of pan-genomics in crop improvement.

Decrease of gene expression diversity during domestication of animals and plants

BACKGROUND

The genetic mechanisms underlying the domestication of animals and plants have been of great interest to biologists since Darwin. To date, little is known about the global pattern of gene expression changes during domestication.

RESULTS

We generated and collected transcriptome data for seven pairs of domestic animals and plants including dog, silkworm, chicken, rice, cotton, soybean and maize and their wild progenitors and compared the expression profiles between the domestic and wild species. Intriguingly, although the number of expressed genes varied little, the domestic species generally exhibited lower gene expression diversity than did the wild species, and this lower diversity was observed for both domestic plants and different kinds of domestic animals including insect, bird and mammal in the whole-genome gene set (WGGS), candidate selected gene set (CSGS) and non-CSGS, with CSGS exhibiting a higher degree of decreased expression diversity. Moreover, different from previous reports which found 2 to 4% of genes were selected by human, we identified 6892 candidate selected genes accounting for 7.57% of the whole-genome genes in rice and revealed that fewer than 8% of the whole-genome genes had been affected by domestication.

CONCLUSIONS

Our results showed that domestication affected the pattern of variation in gene expression throughout the genome and generally decreased the expression diversity across species, and this decrease may have been associated with decreased genetic diversity. This pattern might have profound effects on the phenotypic and physiological changes of domestic animals and plants and provide insights into the genetic mechanisms at the transcriptome level other than decreased genetic diversity and increased linkage disequilibrium underpinning artificial selection.

Fiber Quality Improvement in Upland Cotton (Gossypium hirsutum L.): Quantitative Trait Loci Mapping and Marker Assisted Selection Application

Genetic improvement in fiber quality is one of the main challenges for cotton breeders. Fiber quality traits are controlled by multiple genes and are classified as complex quantitative traits, with a negative relationship with yield potential, so the genetic gain is low in traditional genetic improvement by phenotypic selection. The availability of Gossypium genomic sequences facilitates the development of high-throughput molecular markers, quantitative trait loci (QTL) fine mapping and gene identification, which helps us to validate candidate genes and to use marker assisted selection (MAS) on fiber quality in breeding programs. Based on developments of high density linkage maps, QTLs fine mapping, marker selection and omics, we have performed trait dissection on fiber quality traits in diverse populations of upland cotton. QTL mapping combined with multi-omics approaches such as, RNA sequencing datasets to identify differentially expressed genes have benefited the improvement of fiber quality. In this review, we discuss the application of molecular markers, QTL mapping and MAS for fiber quality improvement in upland cotton.

Identification of Loci and Candidate Genes Responsible for Fiber Length in Upland Cotton (Gossypium hirsutum L.) via Association Mapping and Linkage Analyses

Fiber length (FL) is an important fiber quality trait in cotton. Although many fiber quality quantitative trait loci (QTL) responsible for FL have been identified, most cannot be applied to breeding programs, mainly due to unstable environments or large confidence intervals. In this study, we combined a genome-wide association study (GWAS) and linkage mapping to identify and validate high-quality QTLs responsible for FL. For the GWAS, we developed 93,250 high-quality single-nucleotide polymorphism (SNP) markers based on 355 accessions, and the FL was measured in eight different environments. For the linkage mapping, we constructed an F ₂ population from two extreme accessions. The high-density linkage maps spanned 3,848.29 cM, with an average marker interval of 1.41 cM. In total, 14 and 13 QTLs were identified in the association and linkage mapping analyses, respectively. Most importantly, a major QTL on chromosome D03 identified in both populations explained more than 10% of the phenotypic variation (PV). Furthermore, we found that a sucrose synthesis-related gene (Gh_D03G1338) was associated with FL in this QTL region. The RNA-seq data showed that Gh_D03G1338 was highly expressed during the fiber development stage, and the qRT-PCR analysis showed significant expression differences between the long fiber and short fiber varieties. These results suggest that Gh_D03G1338 may determine cotton fiber elongation by regulating the synthesis of sucrose. Favorable QTLs and candidate genes should be useful for increasing fiber quality in cotton breeding.

Crop genome-wide association study: a harvest of biological relevance

With the advent of rapid genotyping and next-generation sequencing technologies, genome-wide association study (GWAS) has become a routine strategy for decoding genotype-phenotype associations in many species. More than 1000 such studies over the last decade have revealed substantial genotype-phenotype associations in crops and provided unparalleled opportunities to probe functional genomics. Beyond the many 'hits' obtained, this review summarizes recent efforts to increase our understanding of the genetic architecture of complex traits by focusing on non-main effects including epistasis, pleiotropy, and phenotypic plasticity. We also discuss how these achievements and the remaining gaps in our knowledge will guide future studies. Synthetic association is highlighted as leading to false causality, which is prevalent but largely underestimated. Furthermore, validation evidence is appealing for future GWAS, especially in the context of emerging genome-editing technologies.

Genetic Basis of Fiber Improvement and Decreased Stress Tolerance in Cultivated Versus Semi-Domesticated Upland Cotton

Crop domestication from wild ancestors has resulted in the wide adaptation coupled with improved yield and quality traits. However, the genetic basis of many domesticated characteristics remains to be explored. Upland cotton (Gossypium hirsutum) is the most important tetraploid cotton species, accounting for about 90% of world cotton commerce. Here, we reveal the effects of domestication on fiber and stress traits through comprehensive analyses of semi-domesticated races and cultivated cotton accessions. A total of 416 cotton accessions were genotyped, and a decrease in genetic diversity from races to landraces and modern cultivars was detected. Furthermore, 71 domestication selective sweeps (DSS) and 14 improvement selective sweeps (ISS) were identified, with the Dt sub-genome experiencing stronger selection than the At sub-genome during the both selection types. The more expressed genes and a delay in the expression peak of genes related to secondary cell wall (SCW) development in modern cultivars compared to semi-domesticated cotton races, may have contributed to long fibers in these plants. However, down-regulation of genes related to stress response was responsible for decreasing stress tolerance in modern cultivars. We further experimentally confirmed that silencing of PR1 and WRKY20, genes that showed higher expression in the semi-domesticated races, drastically compromised cotton resistance to V. dahliae. Our results reveal fiber improvement and decreased stress tolerance as a result of the domestication of modern upland cotton cultivars.

Nulliplex-branch, a TERMINAL FLOWER 1 ortholog, controls plant growth habit in cotton

Nulliplex-branch (nb) mutants in cotton display a specific architecture. The gene responsible for the nb phenotype was identified, and its modulation mode was further studied. Plant architecture is an important agronomic factor influencing various traits such as yield and variety adaptability in crop plants. Cotton (Gossypium) simultaneously displays monopodial and sympodial growth. Nulliplex-branch (nb) mutants showing determinate sympodial shoots have been reported in both G. hirsutum (Ghnb) and G. barbadense (Gbnb). In this study, the gene responsible for the nb phenotype was identified. GhNB and GbNB were found to be allelic loci and are TERMINAL FLOWER 1 orthologs on the Dt subgenome, though the At copies remain native. Sequencing and association analyses identified four (Gh-nb1-Gh-nb4) and one (Gb-nb1) type of point mutation in the coding sequences of Ghnb and Gbnb, respectively. The NB gene was mainly expressed in the root and shoot apex, and expression rhythms were also observed in these tissues, suggesting that the expression of the NB gene could be regulated by photoperiod. Constitutive overexpression of GhNB suppresses the differentiation of the reproductive shoots. Knockout of both copies of GhNB caused the main and lateral shoots to terminate in flowers, which is a more determinate architecture than that of the nb mutants and implies that its function might be dosage dependent. A protein lipid overlay assay indicated that the amino acid substitutions in Gh-nb1 and Gb-nb1 weaken the ligand-binding activity of the NB protein in vitro. These findings suggest that the NB gene plays crucial roles in regulating the determinacy of shoots, and the modulation of this gene should constitute an effective crop improvement approach through adjusting the growth habit of cotton.

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.

The cotton HD-Zip transcription factor GhHB12 regulates flowering time and plant architecture via the GhmiR157-GhSPL pathway

Domestication converts perennial and photoperiodic ancestral cotton to day-neutral cotton varieties, and the selection of short-season cotton varieties is one of the major objectives of cotton breeding. However, little is known about the mechanism of flowering time in cotton. Here, we report a cotton HD-ZIP I-class transcription factor (GhHB12) specifically expressed in axillary buds, which antagonisticlly interacts with GhSPL10/13 to repress the expression of GhFT, GhFUL, and GhSOC1, resulting in bushy architecture and delayed flowering under long-day conditions. We found that GhHB12-mediated ancestral upland cotton phenotypes (bushy architecture and delayed flowering) could be rescued under short-day conditions. We showed that overexpressing of GhrSPL10 partially rescues the bushy architecture and delayed flowering phenotypes, while overexpression of GhmiR157 reinforced these phenotypes in GhHB12-overexpressing plants. This study defines a regulatory module which regulates cotton architecture, phase transition and could be applied in the breeding of early maturing cotton varieties.

Xin He et al. present a characterization of GhHB12, a HD-ZIP family transcription factor expressed in upland cotton axillary buds. They show that GhHB12 regulates flowering time, plant architecture and phase transition via a regulatory module that could be harnessed to improve cotton for mechanical harvesting.

Genetic variation of dynamic fiber elongation and developmental quantitative trait locus mapping of fiber length in upland cotton (Gossypium hirsutum L.)

Background

In upland cotton (Gossypium hirsutum L.), genotypes with the same mature fiber length (FL) might possess different genes and exhibit differential expression of genes related to fiber elongation at different fiber developmental stages. However, there is a lack of information on the genetic variation influencing fiber length and its quantitative trait loci (QTLs) during the fiber elongation stage. In this study, a subset of upland cotton accessions was selected based on a previous GWAS conducted in China and grown in multiple environments to determine the dynamic fiber length at 10, 15, 20, and 25 days post-anthesis (DPA) and maturity. The germplasm lines were genotyped with the Cotton 63 K Illumina single-nucleotide polymorphism (SNP) array for GWAS.

Results

A total of 25, 38, 57, 89 and 88 SNPs showed significant correlations with fiber length at 10, 15, 20 and 25 DPA and maturity, respectively. In addition, 60 more promising SNPs were detected in at least two tests and two FL developmental time points, and 20 SNPs were located within the confidence intervals of QTLs identified in previous studies. The fastest fiber-length growth rates were obtained at 10 to 15 DPA in 69 upland cotton lines and at 15 to 20 DPA in 14 upland cotton accessions, and 10 SNPs showed significant correlations with the fiber-length growth rate. A combined transcriptome and qRT-PCR analysis revealed that two genes (D10G1008 and D13G2037) showed differential expression between two long-fiber genotypes and two short-fiber genotypes.

Conclusions

This study provides important new insights into the genetic basis of the time-dependent fiber-length trait and reveals candidate SNPs and genes for improving fiber length in upland cotton.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5309-2) contains supplementary material, which is available to authorized users.

A high density SLAF-SNP genetic map and QTL detection for fibre quality traits in Gossypium hirsutum

Background

Upland Cotton (Gossypium hirsutum) is a very important cash crop known for its high quality natural fiber. Recent advances in sequencing technologies provide powerful tools with which to explore the cotton genome for single nucleotide polymorphism marker identification and high density genetic map construction toward more reliable quantitative trait locus mapping.

Results

In the present study, a RIL population was developed by crossing a Chinese high fiber quality cultivar (Yumian 1) and an American high fiber quality line (CA3084), with distinct genetic backgrounds. Specific locus amplified fragment sequencing (SLAF-seq) technology was used to discover SNPs, and a genetic map containing 6254 SNPs was constructed, covering 3141.72 cM with an average distance of 0.5 cM between markers. A total of 95 QTL were detected for fiber quality traits in three environments, explaining 5.5-24.6% of the phenotypic variance. Fifty-five QTL found in multiple environments were considered stable QTL. Nine of the stable QTL were found in all three environments. We identified 14 QTL clusters on 13 chromosomes, each containing one or more stable QTL.

Conclusion

A high-density genetic map of Gossypium hirsutum developed by using specific locus amplified fragment sequencing technology provides detailed mapping of fiber quality QTL, and identification of ‘stable QTL’ found in multiple environments. A marker-rich genetic map provides a foundation for fine mapping, candidate gene identification and marker-assisted selection of favorable alleles at stable QTL in breeding programs.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5294-5) contains supplementary material, which is available to authorized users.

Reference genome sequences of two cultivated allotetraploid cottons, Gossypium hirsutum and Gossypium barbadense

Allotetraploid cotton species (Gossypium hirsutum and Gossypium barbadense) have long been cultivated worldwide for natural renewable textile fibers. The draft genome sequences of both species are available but they are highly fragmented and incomplete^1-4. Here we report reference-grade genome assemblies and annotations for G. hirsutum accession Texas Marker-1 (TM-1) and G. barbadense accession 3-79 by integrating single-molecule real-time sequencing, BioNano optical mapping and high-throughput chromosome conformation capture techniques. Compared with previous assembled draft genomes^1,3, these genome sequences show considerable improvements in contiguity and completeness for regions with high content of repeats such as centromeres. Comparative genomics analyses identify extensive structural variations that probably occurred after polyploidization, highlighted by large paracentric/pericentric inversions in 14 chromosomes. We constructed an introgression line population to introduce favorable chromosome segments from G. barbadense to G. hirsutum, allowing us to identify 13 quantitative trait loci associated with superior fiber quality. These resources will accelerate evolutionary and functional genomic studies in cotton and inform future breeding programs for fiber improvement.

Coping with drought: stress and adaptive mechanisms, and management through cultural and molecular alternatives in cotton as vital constituents for plant stress resilience and fitness

Increased levels of greenhouse gases in the atmosphere and associated climatic variability is primarily responsible for inducing heat waves, flooding and drought stress. Among these, water scarcity is a major limitation to crop productivity. Water stress can severely reduce crop yield and both the severity and duration of the stress are critical. Water availability is a key driver for sustainable cotton production and its limitations can adversely affect physiological and biochemical processes of plants, leading towards lint yield reduction. Adaptation of crop husbandry techniques suitable for cotton crop requires a sound understanding of environmental factors, influencing cotton lint yield and fiber quality. Various defense mechanisms e.g. maintenance of membrane stability, carbon fixation rate, hormone regulation, generation of antioxidants and induction of stress proteins have been found play a vital role in plant survival under moisture stress. Plant molecular breeding plays a functional role to ascertain superior genes for important traits and can offer breeder ready markers for developing ideotypes. This review highlights drought-induced damage to cotton plants at structural, physiological and molecular levels. It also discusses the opportunities for increasing drought tolerance in cotton either through modern gene editing technology like clustered regularly interspaced short palindromic repeat (CRISPR/Cas9), zinc finger nuclease, molecular breeding as well as through crop management, such as use of appropriate fertilization, growth regulator application and soil amendments.

LncRNAs in polyploid cotton interspecific hybrids are derived from transposon neofunctionalization

BACKGROUND

Interspecific hybridization and whole genome duplication are driving forces of genomic and organism diversification. But the effect of interspecific hybridization and whole genome duplication on the non-coding portion of the genome in particular remains largely unknown. In this study, we examine the profile of long non-coding RNAs (lncRNAs), comparing them with that of coding genes in allotetraploid cotton (Gossypium hirsutum), its putative diploid ancestors (G. arboreum; G. raimondii), and an F1 hybrid (G. arboreum × G. raimondii, AD).

RESULTS

We find that most lncRNAs (80%) that were allelic expressed in the allotetraploid genome. Moreover, the genome shock of hybridization reprograms the non-coding transcriptome in the F1 hybrid. Interestingly, the activated lncRNAs are predominantly transcribed from demethylated TE regions, especially from long interspersed nuclear elements (LINEs). The DNA methylation dynamics in the interspecies hybridization are predominantly associated with the drastic expression variation of lncRNAs. Similar trends of lncRNA bursting are also observed in the progress of polyploidization. Additionally, we find that a representative novel lncRNA XLOC_409583 activated after polyploidization from a LINE in the A subgenome of allotetraploid cotton was involved in control of cotton seedling height.

CONCLUSION

Our results reveal that the processes of hybridization and polyploidization enable the neofunctionalization of lncRNA transcripts, acting as important sources of increased plasticity for plants.

The Rise of Cotton Genomics

Cotton is the most important fiber crop, but its polyploid genome is complex and difficult to improve. Recent advances in obtaining high-quality genome sequences of cultivated tetraploid and diploid cotton species may now overcome this hurdle, as well as provide new insights on important fiber quality and disease resistance genes.

Genome-Wide Association Study Identifies Candidate Genes Related to Seed Oil Composition and Protein Content in Gossypium hirsutum L

Cotton (Gossypium spp.) is a leading natural fiber crop and an important source of vegetable protein and oil for humans and livestock. To investigate the genetic architecture of seed nutrients in upland cotton, a genome-wide association study (GWAS) was conducted in a panel of 196 germplasm resources under three environments using a CottonSNP80K chip of 77,774 loci. Relatively high genetic diversity (average gene diversity being 0.331) and phenotypic variation (coefficient of variation, CV, exceeding 3.9%) were detected in this panel. Correlation analysis revealed that the well-documented negative association between seed protein (PR) and oil may be to some extent attributable to the negative correlation between oleic acid (OA) and PR. Linkage disequilibrium (LD) was unevenly distributed among chromosomes and subgenomes. It ranged from 0.10-0.20 Mb (Chr19) to 5.65-5.75 Mb (Chr25) among the chromosomes and the range of Dt-subgenomes LD decay distances was smaller than At-subgenomes. This panel was divided into two subpopulations based on the information of 41,815 polymorphic single-nucleotide polymorphism (SNP) markers. The mixed linear model considering both Q-matrix and K-matrix [MLM(Q+K)] was employed to estimate the association between the SNP markers and the seed nutrients, considering the false positives caused by population structure and the kinship. A total of 47 SNP markers and 28 candidate quantitative trait loci (QTLs) regions were found to be significantly associated with seven cottonseed nutrients, including protein, total fatty acid, and five main fatty acid compositions. In addition, the candidate genes in these regions were analyzed, which included three genes, Gh_D12G1161, Gh_D12G1162, and Gh_D12G1165 that were most likely involved in the control of cottonseed protein concentration. These results improved our understanding of the genetic control of cottonseed nutrients and provided potential molecular tools to develop cultivars with high protein and improved fatty acid compositions in cotton breeding programs through marker-assisted selection.

Genome-Wide Analysis Elucidates the Role of CONSTANS-like Genes in Stress Responses of Cotton

The CONSTANS (CO)-like gene family has been well studied for its role in the regulation of plant flowering time. However, their role remains poorly understood in cotton. To better understand the possible roles of CO-like in cotton, we performed a comprehensive genome-wide analysis of CO-like genes in cotton. Phylogenetic tree analysis showed that CO-like genes naturally clustered into three groups. Segmental duplication and whole genome duplication (WGD), which occurred before polyploidy, were important contributors to its expansion within the At (“t” indicates tetraploid) and Dt subgenomes, particularly in Group III. Long-terminal repeat retroelements were identified as the main transposable elements accompanying 18 genes. The genotype of GhCOL12_Dt displayed low diversity; it was a candidate involved in domestication. Selection pressure analyses indicated that relaxed purifying selection might have provided the main impetus during the evolution of CO-like genes in upland cotton. In addition, the high expression in the torus and calycle indicated that CO-like genes might affect flowering. The genes from Group II, and those from Group III involved in segmental duplication or WGD, might play important roles in response to drought and salt stress. Overall, this comprehensive genome-wide study of the CO-like gene family would facilitate further detailed studies in cotton.

Role of GhHDA5 in H3K9 deacetylation and fiber initiation in Gossypium hirsutum

Cotton fibers are single-celled trichomes that initiate from the epidermal cells of the ovules at or before anthesis. Here, we identified that the histone deacetylase (HDAC) activity is essential for proper cotton fiber initiation. We further identified 15 HDACs homoeologs in each of the A- and D-subgenomes of Gossypium hirsutum. Few of these HDAC homoeologs expressed preferentially during the early stages of fiber development [-1, 0 and 6 days post-anthesis (DPA)]. Among them, GhHDA5 expressed significantly at the time of fiber initiation (-1 and 0 DPA). The in vitro assay for HDAC activity indicated that GhHDA5 primarily deacetylates H3K9 acetylation marks. Moreover, the reduced expression of GhHDA5 also suppresses fiber initiation and lint yield in the RNA interference (RNAi) lines. The 0 DPA ovules of GhHDA5RNAi lines also showed alterations in reactive oxygen species homeostasis and elevated autophagic cell death in the developing fibers. The differentially expressed genes (DEGs) identified through RNA-seq of RNAi line (DEP12) and their pathway analysis showed that GhHDA5 modulates expression of many stress and development-related genes involved in fiber development. The reduced expression of GhHDA5 in the RNAi lines also resulted in H3K9 hyper-acetylation on the promoter region of few DEGs assessed by chromatin immunoprecipitation assay. The positively co-expressed genes with GhHDA5 showed cumulative higher expression during fiber initiation, and gene ontology annotation suggests their involvement in fiber development. Furthermore, the predicted protein interaction network in the positively co-expressed genes indicates HDA5 modulates fiber initiation-specific gene expression through a complex involving reported repressors.

Demography and its effects on genomic variation in crop domestication

Over two thousand plant species have been modified morphologically through cultivation and human use. Here, we review three aspects of crop domestication that are currently undergoing marked revisions, due to analytical advancements and their application to whole genome resequencing (WGS) data. We begin by discussing the duration and demographic history of domestication. There has been debate as to whether domestication occurred quickly or slowly. The latter is tentatively supported both by fossil data and application of WGS data to sequentially Markovian coalescent methods that infer the history of effective population size. This history suggests the possibility of extended human impacts on domesticated lineages prior to their purposeful cultivation. We also make the point that demographic history matters, because it shapes patterns and levels of extant genetic diversity. We illustrate this point by discussing the evolutionary processes that contribute to the empirical observation that most crops examined to date have more putatively deleterious alleles than their wild relatives. These deleterious alleles may contribute to genetic load within crops and may be fitting targets for crop improvement. Finally, the same demographic factors are likely to shape the spectrum of structural variants (SVs) within crops. SVs are known to underlie many of the phenotypic changes associated with domestication and crop improvement, but we currently lack sufficient knowledge about the mechanisms that create SVs, their rates of origin, their population frequencies and their phenotypic effects.

Three-dimensional chromatin packing and positioning of plant genomes

Information and function of a genome are not only decorated with epigenetic marks in the linear DNA sequence but also in their non-random spatial organization in the nucleus. Recent research has revealed that three-dimensional (3D) chromatin organization is highly correlated with the functionality of the genome, contributing to many cellular processes. Driven by the improvements in chromatin conformation capture methods and visualization techniques, the past decade has been an exciting period for the study of plants' 3D genome structures, and our knowledge in this area has been substantially advanced. This Review describes our current understanding of plant chromatin organization and positioning beyond the nucleosomal level, and discusses future directions.

Comprehensive analysis of NAC transcription factors uncovers their roles during fiber development and stress response in cotton

BACKGROUND

Transcription factors operate as important switches of transcription networks, and NAC (NAM, ATAF, and CUC) transcription factors are a plant-specific family involved in multiple biological processes. However, this gene family has not been systematically characterized in cotton.

RESULTS

Here we identify a large number of genes with conservative NAC domains in four cotton species, with 147 found in Gossypium arboreum, 149 in G. raimondii, 267 in G. barbadense and 283 in G. hirsutum. Predicted membrane-bound NAC genes were also identified. Phylogenetic analysis showed that cotton NAC proteins clustered into seven subfamilies and homologous protein pairs showed similar characteristics. Evolutionary property analysis revealed that purifying selection of NAC genes occurred between diploid and polyploid cotton species, and variation analysis showed GhNAC genes may have been subjected to selection and domestication. NAC proteins showed extensive transactivation and this was dependent on the C-terminus. Some development and stress related cis-elements were enriched in the promoters of GhNAC genes. Comprehensive expression analysis indicated that 38 GhNAC genes were candidates for involvement in fiber development, and 120 in stress responses. Gene co-expression network analysis revealed relationships between fiber-associated NAC genes and secondary cell wall (SCW) biosynthesis genes.

CONCLUSIONS

NAC genes were identified in diploid and tetraploid cotton, revealing new insights into their evolution, variation and homology relationships. Transcriptome analysis and co-expression network indicated roles for GhNAC genes in cotton fiber development and stress response, and NAC genes may prove useful in molecular breeding programmes.

Asymmetric expression patterns reveal a strong maternal effect and dosage compensation in polyploid hybrid fish

Background

Hybridization and polyploidization are regarded as the major driving forces in plant speciation, diversification, and ecological adaptation. Our knowledge regarding the mechanisms of duplicated-gene regulation following genomic merging or doubling is primarily derived from plants and is sparse for vertebrates.

Results

We successfully obtained an F1 generation (including allodiploid hybrids and triploid hybrids) from female Megalobrama amblycephala Yih (BSB, 2n = 48) × male Xenocypri davidi Bleeker (YB, 2n = 48). The duplicated-gene expression patterns of the two types of hybrids were explored using RNA-Seq data. In total, 5.44 × 10⁸ (69.32 GB) clean reads and 499,631 assembled unigenes were obtained from the testis transcriptomes. The sequence similarity analysis of 4265 orthologs revealed that the merged genomes were dominantly expressed in different ploidy hybrids. The differentially expressed genes in the two types of hybrids were asymmetric compared with those in both parents. Furthermore, the genome-wide expression level dominance (ELD) was biased toward the maternal BSB genome in both the allodiploid and triploid hybrids. In addition, the dosage-compensation mechanisms that reduced the triploid expression levels to the diploid state were determined in the triploid hybrids.

Conclusions

Our results indicate that divergent genomes undergo strong interactions and domination in allopolyploid offspring. Genomic merger has a greater effect on the gene-expression patterns than genomic doubling. The various expression mechanisms (including maternal effect and dosage compensation) in different ploidy hybrids suggest that the initial genomic merger and doubling play important roles in polyploidy adaptation and evolution.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4883-7) contains supplementary material, which is available to authorized users.

Entering the Next Dimension: Plant Genomes in 3D

After linear sequences of genomes and epigenomic landscape data, the 3D organization of chromatin in the nucleus is the next level to be explored. Different organisms present a general hierarchical organization, with chromosome territories at the top. Chromatin interaction maps, obtained by chromosome conformation capture (3C)-based methodologies, for eight plant species reveal commonalities, but also differences, among them and with animals. The smallest structures, found in high-resolution maps of the Arabidopsis genome, are single genes. Epigenetic marks (histone modification and DNA methylation), transcriptional activity, and chromatin interaction appear to be correlated, and whether structure is the cause or consequence of the function of interacting regions is being actively investigated.

Genome-wide Hi-C analysis reveals extensive hierarchical chromatin interactions in rice

The non-random spatial packing of chromosomes in the nucleus plays a critical role in orchestrating gene expression and genome function. Here, we present a Hi-C analysis of the chromatin interaction patterns in rice (Oryza sativa L.) at hierarchical architectural levels. We confirm that rice chromosomes occupy their own territories with certain preferential inter-chromosomal associations. Moderate compartment delimitation and extensive TADs (Topologically Associated Domains) were determined to be associated with heterogeneous genomic compositions and epigenetic marks in the rice genome. We found subtle features including chromatin loops, gene loops, and off-/near-diagonal intensive interaction regions. Gene chromatin loops associated with H3K27me3 could be positively involved in gene expression. In addition to insulated enhancing effects for neighbor gene expression, the identified rice gene loops could bi-directionally (+/-) affect the expression of looped genes themselves. Finally, web-interleaved off-diagonal IHIs/KEEs (Interactive Heterochromatic Islands or KNOT ENGAGED ELEMENTs) could trap transposable elements (TEs) via the enrichment of silencing epigenetic marks. In parallel, the near-diagonal FIREs (Frequently Interacting Regions) could positively affect the expression of involved genes. Our results suggest that the chromatin packing pattern in rice is generally similar to that in Arabidopsis thaliana but with clear differences at specific structural levels. We conclude that genomic composition, epigenetic modification, and transcriptional activity could act in combination to shape global and local chromatin packing in rice. Our results confirm recent observations in rice and A. thaliana but also provide additional insights into the patterns and features of chromatin organization in higher plants.

Long noncoding RNAs involve in resistance to Verticillium dahliae, a fungal disease in cotton

Long noncoding RNAs (lncRNAs) have several known functions in plant development, but their possible roles in responding to plant disease remain largely unresolved. In this study, we described a comprehensive disease-responding lncRNA profiles in defence against a cotton fungal disease Verticillium dahliae. We further revealed the conserved and specific characters of disease-responding process between two cotton species. Conservatively for two cotton species, we found the expression dominance of induced lncRNAs in the Dt subgenome, indicating a biased induction pattern in the co-existing subgenomes of allotetraploid cotton. Comparative analysis of lncRNA expression and their proposed functions in resistant Gossypium barbadense cv. '7124' versus susceptible Gossypium hirsutum cv. 'YZ1' revealed their distinct disease response mechanisms. Species-specific (LS) lncRNAs containing more SNPs displayed a fiercer inducing level postinfection than the species-conserved (core) lncRNAs. Gene Ontology enrichment of LS lncRNAs and core lncRNAs indicates distinct roles in the process of biotic stimulus. Further functional analysis showed that two core lncRNAs, GhlncNAT-ANX2- and GhlncNAT-RLP7-silenced seedlings, displayed an enhanced resistance towards V. dahliae and Botrytis cinerea, possibly associated with the increased expression of LOX1 and LOX2. This study represents the first characterization of lncRNAs involved in resistance to fungal disease and provides new clues to elucidate cotton disease response mechanism.

Resequencing of 243 diploid cotton accessions based on an updated A genome identifies the genetic basis of key agronomic traits

The ancestors of Gossypium arboreum and Gossypium herbaceum provided the A subgenome for the modern cultivated allotetraploid cotton. Here, we upgraded the G. arboreum genome assembly by integrating different technologies. We resequenced 243 G. arboreum and G. herbaceum accessions to generate a map of genome variations and found that they are equally diverged from Gossypium raimondii. Independent analysis suggested that Chinese G. arboreum originated in South China and was subsequently introduced to the Yangtze and Yellow River regions. Most accessions with domestication-related traits experienced geographic isolation. Genome-wide association study (GWAS) identified 98 significant peak associations for 11 agronomically important traits in G. arboreum. A nonsynonymous substitution (cysteine-to-arginine substitution) of GaKASIII seems to confer substantial fatty acid composition (C16:0 and C16:1) changes in cotton seeds. Resistance to fusarium wilt disease is associated with activation of GaGSTF9 expression. Our work represents a major step toward understanding the evolution of the A genome of cotton.

Resequencing a core collection of upland cotton identifies genomic variation and loci influencing fiber quality and yield

Upland cotton is the most important natural-fiber crop. The genomic variation of diverse germplasms and alleles underpinning fiber quality and yield should be extensively explored. Here, we resequenced a core collection comprising 419 accessions with 6.55-fold coverage depth and identified approximately 3.66 million SNPs for evaluating the genomic variation. We performed phenotyping across 12 environments and conducted genome-wide association study of 13 fiber-related traits. 7,383 unique SNPs were significantly associated with these traits and were located within or near 4,820 genes; more associated loci were detected for fiber quality than fiber yield, and more fiber genes were detected in the D than the A subgenome. Several previously undescribed causal genes for days to flowering, fiber length, and fiber strength were identified. Phenotypic selection for these traits increased the frequency of elite alleles during domestication and breeding. These results provide targets for molecular selection and genetic manipulation in cotton improvement.

Core cis-element variation confers subgenome-biased expression of a transcription factor that functions in cotton fiber elongation

Cotton cultivars have evolved to produce extensive, long, seed-born fibers important for the textile industry, but we know little about the molecular mechanism underlying spinnable fiber formation. Here, we report how PACLOBUTRAZOL RESISTANCE 1 (PRE1) in cotton, which encodes a basic helix-loop-helix (bHLH) transcription factor, is a target gene of spinnable fiber evolution. Differential expression of homoeologous genes in polyploids is thought to be important to plant adaptation and novel phenotypes. PRE1 expression is specific to cotton fiber cells, upregulated during their rapid elongation stage and A-homoeologous biased in allotetraploid cultivars. Transgenic studies demonstrated that PRE1 is a positive regulator of fiber elongation. We determined that the natural variation of the canonical TATA-box, a regulatory element commonly found in many eukaryotic core promoters, is necessary for subgenome-biased PRE1 expression, representing a mechanism underlying the selection of homoeologous genes. Thus, variations in the promoter of the cell elongation regulator gene PRE1 have contributed to spinnable fiber formation in cotton. Overexpression of GhPRE1 in transgenic cotton yields longer fibers with improved quality parameters, indicating that this bHLH gene is useful for improving cotton fiber quality.

Gene retention, fractionation and subgenome differences in polyploid plants

All natural plant species are evolved from ancient polyploids. Polyloidization plays an important role in plant genome evolution, species divergence and crop domestication. We review how the pattern of polyploidy within the plant phylogenetic tree has engendered hypotheses involving mass extinctions, lag-times following polyploidy, and epochs of asexuality. Polyploidization has happened repeatedly in plant evolution and, we conclude, is important for crop domestication. Once duplicated, the effect of purifying selection on any one duplicated gene is relaxed, permitting duplicate gene and regulatory element loss (fractionation). We review the general topic of fractionation, and how some gene categories are retained more than others. Several explanations, including neofunctionalization, subfunctionalization and gene product dosage balance, have been shown to influence gene content over time. For allopolyploids, genetic differences between parental lines immediately manifest as subgenome dominance in the wide-hybrid, and persist and propagate for tens of millions of years. While epigenetic modifications are certainly involved in genome dominance, it has been difficult to determine which came first, the chromatin marks being measured or gene expression. Data support the conclusion that genome dominance and heterosis are antagonistic and mechanically entangled; both happen immediately in the synthetic wide-cross hybrid. Also operating in this hybrid are mechanisms of 'paralogue interference'. We present a foundation model to explain gene expression and vigour in a wide hybrid/new allotetraploid. This Review concludes that some mechanisms operate immediately at the wide-hybrid, and other mechanisms begin their operations later. Direct interaction of new paralogous genes, as measured using high-resolution chromatin conformation capture, should inform future research and single cell transcriptome sequencing should help achieve specificity while studying gene sub- and neo-functionalization.