Asymmetric subgenome selection and cis-regulatory divergence during cotton domestication

A Retrotransposon Insertion in GhMML3_D12 Is Likely Responsible for the Lintless Locus li3 of Tetraploid Cotton

Cotton (Gossypium) seed fibers can be divided into lint (long) or fuzz (very short). Using fiberless (fuzzless-lintless) mutants, the lint initiation gene Li₃ was identified by map-based cloning. The gene is an R2R3-MYB transcription factor located on chromosome D12 (GhMML3_D12). Sequence analysis revealed that li₃ is a loss-of-function allele containing a retrotransposon insertion in the second exon that completely blocks the gene’s expression. The genetic loci n₂ and n₃ underlying the recessive fuzzless phenotype in Gossypium hirsutum were also mapped. The genomic location of n₃ overlapped with that of the dominant fuzzless locus N₁, and n₃ appeared to be a loss-of-function allele caused by a single nucleotide polymorphism (SNP) mutation in the coding region of GhMML3_A12. The n₂ allele was found to be co-located with li₃ and originated from G. babardense. n₂ and li₃ are possibly the multiple alleles of the GhMML3_D12 gene. Genetic analysis showed that Li₃ and N₃ are a pair of homologs with additive effects for the initiation of fibers (fuzz or lint). In addition, the presence of another locus was speculated, and it appeared to show an inhibitory effect on the expression of GhMML3. These findings provide new information about the genetic factors affecting the initiation of fibers in cotton.

The evolutionary origin and domestication history of goldfish (Carassius auratus)

Goldfish have been subjected to over 1,000 y of intensive domestication and selective breeding. In this report, we describe a high-quality goldfish genome (2n = 100), anchoring 95.75% of contigs into 50 pseudochromosomes. Comparative genomics enabled us to disentangle the two subgenomes that resulted from an ancient hybridization event. Resequencing 185 representative goldfish variants and 16 wild crucian carp revealed the origin of goldfish and identified genomic regions that have been shaped by selective sweeps linked to its domestication. Our comprehensive collection of goldfish varieties enabled us to associate genetic variations with a number of well-known anatomical features, including features that distinguish traditional goldfish clades. Additionally, we identified a tyrosine-protein kinase receptor as a candidate causal gene for the first well-known case of Mendelian inheritance in goldfish-the transparent mutant. The goldfish genome and diversity data offer unique resources to make goldfish a promising model for functional genomics, as well as domestication.

Mapping Regulatory Determinants in Plants

The domestication and improvement of many plant species have frequently involved modulation of transcriptional outputs and continue to offer much promise for targeted trait engineering. The cis-regulatory elements (CREs) controlling these trait-associated transcriptional variants however reside within non-coding regions that are currently poorly annotated in most plant species. This is particularly true in large crop genomes where regulatory regions constitute only a small fraction of the total genomic space. Furthermore, relatively little is known about how CREs function to modulate transcription in plants. Therefore understanding where regulatory regions are located within a genome, what genes they control, and how they are structured are important factors that could be used to guide both traditional and synthetic plant breeding efforts. Here, we describe classic examples of regulatory instances as well as recent advances in plant regulatory genomics. We highlight valuable molecular tools that are enabling large-scale identification of CREs and offering unprecedented insight into how genes are regulated in diverse plant species. We focus on chromatin environment, transcription factor (TF) binding, the role of transposable elements, and the association between regulatory regions and target genes.

Identification, expression, and artificial selection of silkworm epigenetic modification enzymes

Background

Understanding the genetic basis of phenotype variations during domestication and breeding is of great interest. Epigenetics and epigenetic modification enzymes (EMEs) may play a role in phenotypic variations; however, no comprehensive study has been performed to date. Domesticated silkworm (Bombyx mori) may be utilized as a model in determining how EMEs influence domestication traits.

Results

We identified 44 EMEs in the genome of silkworm (Bombyx mori) using homology searching. Phylogenetic analysis showed that genes in a subfamily among different animals were well clustered, and the expression pattern of EMEs is constant among Bombyx mori, Drosophila melanogaster, and Mus musculus. These are most highly expressed in brain, early embryo, and internal genitalia. By gene-related selective sweeping, we identified five BmEMEs under artificial selection during the domestication and breeding of silkworm. Among these selected genes, BmSuv4–20 and BmDNMT2 harbor selective mutations in their upstream regions that alter transcription factor-binding sites. Furthermore, these two genes are expressed higher in the testis and ovary of domesticated silkworm compared to wild silkworms, and correlations between their expression pattern and meiosis of the sperm and ova were observed.

Conclusions

The domestication of silkworm has induced artificial selection on epigenetic modification markers that may have led to phenotypic changes during domestication. We present a novel perspective to understand the genetic basis underlying animal domestication and breeding.

Supplementary information

Supplementary information accompanies this paper at 10.1186/s12864-020-07155-z.

Decoding the plant genome: From epigenome to 3D organization

The linear genome of eukaryotes is partitioned into diverse chromatin states and packaged into a three-dimensional (3D) structure, which has functional implications in DNA replication, DNA repair, and transcriptional regulation. Over the past decades, research on plant functional genomics and epigenomics has made great progress, with thousands of genes cloned and molecular mechanisms of diverse biological processes elucidated. Recently, 3D genome research has gradually attracted great attention of many plant researchers. Herein, we briefly review the progress in genomic and epigenomic research in plants, with a focus on Arabidopsis and rice, and summarize the currently used technologies and advances in plant 3D genome organization studies. We also discuss the relationships between one-dimensional linear genome sequences, epigenomic states, and the 3D chromatin architecture. This review provides basis for future research on plant 3D genomics.

Marchantia TCP transcription factor activity correlates with three-dimensional chromatin structure

Information in the genome is not only encoded within sequence or epigenetic modifications, but is also found in how it folds in three-dimensional space. The formation of self-interacting genomic regions, named topologically associated domains (TADs), is known as a key feature of genome organization beyond the nucleosomal level. However, our understanding of the formation and function of TADs in plants is extremely limited. Here we show that the genome of Marchantia polymorpha, a member of a basal land plant lineage, exhibits TADs with epigenetic features similar to those of higher plants. By analysing various epigenetic marks across Marchantia TADs, we find that these regions generally represent interstitial heterochromatin and their borders are enriched with Marchantia transcription factor TCP1. We also identify a type of TAD that we name 'TCP1-rich TAD', in which genomic regions are highly accessible and are densely bound by TCP1 proteins. Transcription of TCP1 target genes differs on the basis gene location, and those in TCP1-rich TADs clearly show a lower expression level. In tcp1 mutant lines, neither TCP1-bound TAD borders nor TCP1-rich TADs display drastically altered chromatin organization patterns, suggesting that, in Marchantia, TCP1 is dispensable for TAD formation. However, we find that in tcp1 mutants, genes residing in TCP1-rich TADs have a greater extent of expression fold change as opposed to genes that do not belong to these TADs. Our results suggest that, besides standing as spatial chromatin-packing modules, plant TADs function as nuclear microcompartments associated with transcription factor activities.

Genomic signatures and candidate genes of lint yield and fibre quality improvement in Upland cotton in Xinjiang

Xinjiang has been the largest and highest yield cotton production region not only in China, but also in the world. Improvements in Upland cotton cultivars in Xinjiang have occurred via pedigree selection and/or crossing of elite alleles from the former Soviet Union and other cotton producing regions of China. But it is unclear how genomic constitutions from foundation parents have been selected and inherited. Here, we deep-sequenced seven historic foundation parents, comprising four cultivars introduced from the former Soviet Union (108Ф, C1470, 611Б and KK1543) and three from United States and Africa (DPL15, STV2B and UGDM), and re-sequenced sixty-nine Xinjiang modern cultivars. Phylogenetic analysis of more than 2 million high-quality single nucleotide polymorphisms allowed their classification two groups, suggesting that Xinjiang Upland cotton cultivars were not only spawned from 108Ф, C1470, 611Б and KK1543, but also had a close kinship with DPL15, STV2B and UGDM. Notably, identity-by-descent (IBD) tracking demonstrated that the former Soviet Union cultivars have made a huge contribution to modern cultivar improvement in Xinjiang. A total of 156 selective sweeps were identified. Among them, apoptosis-antagonizing transcription factor gene (GhAATF1) and mitochondrial transcription termination factor family protein gene (GhmTERF1) were highly involved in the determination of lint percentage. Additionally, the auxin response factor gene (GhARF3) located in inherited IBD segments from 108Ф and 611Б was highly correlated with fibre quality. These results provide an insight into the genomics of artificial selection for improving cotton production and facilitate next-generation precision breeding of cotton and other crops.

BnaSNPDB: An interactive web portal for the efficient retrieval and analysis of SNPs among 1,007 rapeseed accessions

The rapid development of high-throughput sequencing technology and the decrease in sequencing costs provide valuable resources and great opportunities for researchers to investigate genomic variations across hundreds or even thousands of accessions in the post-genomic era. The management and exploration of these large-scale genomic variations heavily rely on programming and command-line environments, which are challenging and time-consuming for most experimental biologists and plant breeders. Here, we present BnaSNPDB, an interactive web portal with a user-friendly interface that provides multiple analysis modules for retrieving, analyzing, and visualizing single nucleotide polymorphisms among 1,007 accessions of worldwide rapeseed germplasm. It is compatible, extendable, and portable to be easily set up on different operating systems, and can be accessed at http://121.41.229.126:3838/bnasnpdb and http://rapeseed.zju.edu.cn:3838/bnasnpdb. Its whole dataset and code are available at https://github.com/YTLogos/BnaSNPDB. This database is essential for accelerating studies on the functional genomics and screening of the molecular markers of molecular-assisted breeding in rapeseed.

The Genomes of the Allohexaploid Echinochloa crus-galli and Its Progenitors Provide Insights into Polyploidization-Driven Adaptation

The hexaploid species Echinochloa crus-galli is one of the most detrimental weeds in crop fields, especially in rice paddies. Its evolutionary history is similar to that of bread wheat, arising through polyploidization after hybridization between a tetraploid and a diploid species. In this study, we generated and analyzed high-quality genome sequences of diploid (E. haploclada), tetraploid (E. oryzicola), and hexaploid (E. crus-galli) Echinochloa species. Gene family analysis showed a significant loss of disease-resistance genes such as those encoding NB-ARC domain-containing proteins during Echinochloa polyploidization, contrary to their significant expansionduring wheat polyploidization, suggesting that natural selection might favor reduced investment in resistance in this weed to maximize its growth and reproduction. In contrast to the asymmetric patterns of genome evolution observed in wheat and other crops, no significant differences in selection pressure were detected between the subgenomes in E. oryzicola and E. crus-galli. In addition, distinctive differences in subgenome transcriptome dynamics during hexaploidization were observed between E. crus-galli and bread wheat. Collectively, our study documents genomic mechanisms underlying the adaptation of a major agricultural weed during polyploidization. The genomic and transcriptomic resources of three Echinochloa species and new insights into the polyploidization-driven adaptive evolution would be useful for future breeding cereal crops.

Genome-wide association study reveals the genetic basis of fiber quality traits in upland cotton (Gossypium hirsutum L.)

BACKGROUND

Fiber quality is an important economic trait of cotton, and its improvement is a major goal of cotton breeding. To better understand the genetic mechanisms responsible for fiber quality traits, we conducted a genome-wide association study to identify and mine fiber-quality-related quantitative trait loci (QTLs) and genes.

RESULTS

In total, 42 single nucleotide polymorphisms (SNPs) and 31 QTLs were identified as being significantly associated with five fiber quality traits. Twenty-five QTLs were identified in previous studies, and six novel QTLs were firstly identified in this study. In the QTL regions, 822 genes were identified and divided into four clusters based on their expression profiles. We also identified two pleiotropic SNPs. The SNP locus i52359Gb was associated with fiber elongation, strength, length and uniformity, while i11316Gh was associated with fiber strength and length. Moreover, these two SNPs were nonsynonymous and located in genes Gh_D09G2376 and Gh_D06G1908, respectively. RT-qPCR analysis revealed that these two genes were preferentially expressed at one or more stages of cotton fiber development, which was consistent with the RNA-seq data. Thus, Gh_D09G2376 and Gh_D06G1908 may be involved in fiber developmental processes.

CONCLUSIONS

The findings of this study provide insights into the genetic bases of fiber quality traits, and the identified QTLs or genes may be applicable in cotton breeding to improve fiber quality.

Plant and animal chromatin three-dimensional organization: similar structures but different functions

Chromatin is the main carrier of genetic information and is non-randomly distributed within the nucleus. Next-generation sequence-based chromatin conformation capture technologies have enabled us to directly examine its three-dimensional organization at an unprecedented scale and resolution. In the best-studied mammalian models, chromatin folding can be broken down into three hierarchical levels, compartment, domains, and loops, which play important roles in transcriptional regulation. Although similar structures have now been identified in plants, they might not possess exactly the same functions as the mammalian ones. Here, we review recent Hi-C studies in plants, compare plant chromatin structures with their mammalian counterparts, and discuss the differences between plants with different genome sizes.

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

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

Distinct Expression and Methylation Patterns for Genes with Different Fates following a Single Whole-Genome Duplication in Flowering Plants

For most sequenced flowering plants, multiple whole-genome duplications (WGDs) are found. Duplicated genes following WGD often have different fates that can quickly disappear again, be retained for long(er) periods, or subsequently undergo small-scale duplications. However, how different expression, epigenetic regulation, and functional constraints are associated with these different gene fates following a WGD still requires further investigation due to successive WGDs in angiosperms complicating the gene trajectories. In this study, we investigate lotus (Nelumbo nucifera), an angiosperm with a single WGD during the K-pg boundary. Based on improved intraspecific-synteny identification by a chromosome-level assembly, transcriptome, and bisulfite sequencing, we explore not only the fundamental distinctions in genomic features, expression, and methylation patterns of genes with different fates after a WGD but also the factors that shape post-WGD expression divergence and expression bias between duplicates. We found that after a WGD genes that returned to single copies show the highest levels and breadth of expression, gene body methylation, and intron numbers, whereas the long-retained duplicates exhibit the highest degrees of protein-protein interactions and protein lengths and the lowest methylation in gene flanking regions. For those long-retained duplicate pairs, the degree of expression divergence correlates with their sequence divergence, degree in protein-protein interactions, and expression level, whereas their biases in expression level reflecting subgenome dominance are associated with the bias of subgenome fractionation. Overall, our study on the paleopolyploid nature of lotus highlights the impact of different functional constraints on gene fate and duplicate divergence following a single WGD in plant.

Fine-mapping and candidate gene analysis of qFS-Chr. D02, a QTL for fibre strength introgressed from a semi-wild cotton into Gossypium hirsutum

Fibre strength (FS) is an important quality attribute in the modern textile industry, which is genetically controlled by quantitative trait loci (QTLs). Fine-mapping stable QTLs for FS to identify candidate genes would be valuable for uncovering the genetic basis of fibre quality traits in cotton. Here, a single segment introgression line, IL-D2-2, from the cross of (TM-1×TX-1046) reported in our previous studies, was found to have significantly improved FS compared with the recurrent parent TM-1. To fine-map the QTLs of the FS, we further crossed IL-D2-2 with its recurrent parent TM-1 to produce F₂ and F_2:3 populations. QTL analysis and substitution mapping showed qFS-Chr. D02 was anchored into a 550.66 kb-interval between two markers, INTR1027 and JESPR-231. This interval contained 67 genes, among which 27 genes related to cell-wall synthesis were selected to conduct qRT-PCR. The results revealed seven genes were expressed significantly differently during the fibre secondary-wall-thickening stage (10-25 days post-anthesis), three being upregulated and four downregulated in IL-D2-2. Both GH_D02G2269 (UDP-glucosyl transferase 84B1) and GH_D02G2289 (unknown function (DUF869)) with nonsynonymous SNPs in IL-D2-2 had significantly downregulated expression, suggesting they were candidates for qFS-Chr. D02. This research provides information about marker-assisted selection for cotton fibre strength improvement.

Genome-Wide Characterization of DNase I-Hypersensitive Sites and Cold Response Regulatory Landscapes in Grasses

Deep sequencing of DNase-I treated chromatin (DNase-seq) can be used to identify DNase I-hypersensitive sites (DHSs) and facilitates genome-scale mining of de novo cis-regulatory DNA elements. Here, we adapted DNase-seq to generate genome-wide maps of DHSs using control and cold-treated leaf, stem, and root tissues of three widely studied grass species: Brachypodium distachyon, foxtail millet (Setaria italica), and sorghum (Sorghum bicolor). Functional validation demonstrated that 12 of 15 DHSs drove reporter gene expression in transiently transgenic B. distachyon protoplasts. DHSs under both normal and cold treatment substantially differed among tissues and species. Intriguingly, the putative DHS-derived transcription factors (TFs) are largely colocated among tissues and species and include 17 ubiquitous motifs covering all grass taxa and all tissues examined in this study. This feature allowed us to reconstruct a regulatory network that responds to cold stress. Ethylene-responsive TFs SHINE3, ERF2, and ERF9 occurred frequently in cold feedback loops in the tissues examined, pointing to their possible roles in the regulatory network. Overall, we provide experimental annotation of 322,713 DHSs and 93 derived cold-response TF binding motifs in multiple grasses, which could serve as a valuable resource for elucidating the transcriptional networks that function in the cold-stress response and other physiological processes.

Global Role of Crop Genomics in the Face of Climate Change

The development of climate change resilient crops is necessary if we are to meet the challenge of feeding the growing world's population. We must be able to increase food production despite the projected decrease in arable land and unpredictable environmental conditions. This review summarizes the technological and conceptual advances that have the potential to transform plant breeding, help overcome the challenges of climate change, and initiate the next plant breeding revolution. Recent developments in genomics in combination with high-throughput and precision phenotyping facilitate the identification of genes controlling critical agronomic traits. The discovery of these genes can now be paired with genome editing techniques to rapidly develop climate change resilient crops, including plants with better biotic and abiotic stress tolerance and enhanced nutritional value. Utilizing the genetic potential of crop wild relatives (CWRs) enables the domestication of new species and the generation of synthetic polyploids. The high-quality crop plant genome assemblies and annotations provide new, exciting research targets, including long non-coding RNAs (lncRNAs) and cis-regulatory regions. Metagenomic studies give insights into plant-microbiome interactions and guide selection of optimal soils for plant cultivation. Together, all these advances will allow breeders to produce improved, resilient crops in relatively short timeframes meeting the demands of the growing population and changing climate.

Genome-Wide Identification and Expression Profile Analysis of the PHT1 Gene Family in Gossypium hirsutum and Its Two Close Relatives of Subgenome Donor Species

Phosphate transporter (PHT) is responsible for plant phosphorus (P) absorption and transport. PHT1 is a component of the high-affinity phosphate transporter system and plays pivotal roles in P absorption under P starvation conditions. However, in cotton, the number and identity of PHT1 genes that are crucial for P absorption from soil remain unclear. Here, genome-wide identification detected twelve PHT1 genes in Gossypium hirsutum and seven and eight PHT1 genes in two close relatives of the G. hirsutum genome-G. arboreum and G. raimondii, respectively. In addition, under low-phosphate treatment, the expressions of GaPHT1;3, GaPHT1;4, and GaPHT1;5 in roots were upregulated after 3 h of induction, and GhPHT1;3-At, GhPHT1;4-At, GhPHT1;5-At, GhPHT1;3-Dt, GhPHT1;4-Dt, and GhPHT1;5-Dt in the roots began to respond after 1 h of induction. Homologous pairs-GaPHT1;4 and GhPHT1;4-At; GaPHT1;5 and GhPHT1;5-At; GrPHT1;4 and GhPHT1;4-Dt, with GhPHT1;5-Dt and GhPHT1;5-At being syntenic-were all highly expressed in the roots under normal conditions. Among the genes highly expressed in the roots, GhPHT1;4-At, GhPHT1;5-At, GhPHT1;4-Dt and GhPHT1;5-Dt were continuously upregulated by P starvation. Therefore, it is concluded that these four genes might be key genes for P uptake in cotton roots. The results of this study provide insights into the mechanisms of P absorption and transport in cotton.

The miR319-Targeted GhTCP4 Promotes the Transition from Cell Elongation to Wall Thickening in Cotton Fiber

Plant cell growth involves a complex interplay among cell-wall expansion, biosynthesis, and, in specific tissues, secondary cell wall (SCW) deposition, yet the coordination of these processes remains elusive. Cotton fiber cells are developmentally synchronous, highly elongated, and contain nearly pure cellulose when mature. Here, we report that the transcription factor GhTCP4 plays an important role in balancing cotton fiber cell elongation and wall synthesis. During fiber development the expression of miR319 declines while GhTCP4 transcript levels increase, with high levels of the latter promoting SCW deposition. GhTCP4 interacts with a homeobox-containing factor, GhHOX3, and repressing its transcriptional activity. GhTCP4 and GhHOX3 function antagonistically to regulate cell elongation, thereby establishing temporal control of fiber cell transition to the SCW stage. We found that overexpression of GhTCP4A upregulated and accelerated activation of the SCW biosynthetic pathway in fiber cells, as revealed by transcriptome and promoter activity analyses, resulting in shorter fibers with varied lengths and thicker walls. In contrast, GhTCP4 downregulation led to slightly longer fibers and thinner cell walls. The GhHOX3-GhTCP4 complex may represent a general mechanism of cellular development in plants since both are conserved factors in many species, thus providing us a potential molecular tool for the design of fiber traits.

High-density genetic variation maps reveal the correlation between asymmetric interspecific introgressions and improvement of agronomic traits in Upland and Pima cotton varieties developed in Xinjiang, China

The two new world tetraploid cottons, Gossypium hirsutum and Gossypium barbadense, are cultivated worldwide and are characterised by a high yield and superior fibre quality, respectively. Historical genetic introgression has been reported between them; however, the existence of introgression and its genetic effects on agronomic traits remain unclear with regard to independent breeding of G. hirsutum (Upland cotton) and G. barbadense (Pima cotton) elite cultivars. We collected 159 G. hirsutum and 70 G. barbadense cultivars developed in Xinjiang, China, along with 30 semi-wild accessions of G. hirsutum, to perform interspecific introgression tests, intraspecific selection analyses and genome-wide association studies (GWAS) with fibre quality and yield component traits in multiple environments. In total, we identified seven interspecific introgression events and 52 selective sweep loci in G. hirsutum, as well as 17 interspecific introgression events and 19 selective sweep loci in G. barbadense. Correlation tests between agronomic traits and introgressions showed that introgression loci were mutually beneficial for the improvement of fibre quality and yield traits in both species. In addition, the phenotypic effects of four interspecific introgression events could be detected by intraspecific GWAS, with Gb_INT13 significantly improving fibre yield in G. barbadense. The present study describes the landscape of genetic introgression and selection between the two species, and highlights the genetic effects of introgression among populations, which can be used for future improvement of fibre yield and quality in G. barbadense and G. hirsutum, respectively.

Genetic dissection of an allotetraploid interspecific CSSLs guides interspecific genetics and breeding in cotton

Background

The low genetic diversity of Upland cotton limits the potential for genetic improvement. Making full use of the genetic resources of Sea-island cotton will facilitate genetic improvement of widely cultivated Upland cotton varieties. The chromosome segments substitution lines (CSSLs) provide an ideal strategy for mapping quantitative trait loci (QTL) in interspecific hybridization.

Results

In this study, a CSSL population was developed by PCR-based markers assisted selection (MAS), derived from the crossing and backcrossing of Gossypium hirsutum (Gh) and G. barbadense (Gb), firstly. Then, by whole genome re-sequencing, 11,653,661 high-quality single nucleotide polymorphisms (SNPs) were identified which ultimately constructed 1211 recombination chromosome introgression segments from Gb. The sequencing-based physical map provided more accurate introgressions than the PCR-based markers. By exploiting CSSLs with mutant morphological traits, the genes responding for leaf shape and fuzz-less mutation in the Gb were identified. Based on a high-resolution recombination bin map to uncover genetic loci determining the phenotypic variance between Gh and Gb, 64 QTLs were identified for 14 agronomic traits with an interval length of 158 kb to 27 Mb. Surprisingly, multiple alleles of Gb showed extremely high value in enhancing cottonseed oil content (SOC).

Conclusions

This study provides guidance for studying interspecific inheritance, especially breeding researchers, for future studies using the traditional PCR-based molecular markers and high-throughput re-sequencing technology in the study of CSSLs. Available resources include candidate position for controlling cotton quality and quantitative traits, and excellent breeding materials. Collectively, our results provide insights into the genetic effects of Gb alleles on the Gh, and provide guidance for the utilization of Gb alleles in interspecific breeding.

3D genome architecture coordinates trans and cis regulation of differentially expressed ear and tassel genes in maize

BACKGROUND

Maize ears and tassels are two separate types of inflorescence which are initiated by similar developmental processes but gradually develop distinct architectures. However, coordinated trans and cis regulation of differentially expressed genes determining ear and tassel architecture within the 3D genome context is largely unknown.

RESULTS

We identify 56,055 and 52,633 open chromatin regions (OCRs) in developing maize ear and tassel primordia using ATAC-seq and characterize combinatorial epigenome features around these OCRs using ChIP-seq, Bisulfite-seq, and RNA-seq datasets. Our integrative analysis of coordinated epigenetic modification and transcription factor binding to OCRs highlights the cis and trans regulation of differentially expressed genes in ear and tassel controlling inflorescence architecture. We further systematically map chromatin interactions at high-resolution in corresponding tissues using in situ digestion-ligation-only Hi-C (DLO Hi-C). The extensive chromatin loops connecting OCRs and genes provide a 3D view on cis- and trans-regulatory modules responsible for ear- and tassel-specific gene expression. We find that intergenic SNPs tend to locate in distal OCRs, and our chromatin interaction maps provide a potential mechanism for trait-associated intergenic SNPs that may contribute to phenotypic variation by influencing target gene expression through chromatin loops.

CONCLUSIONS

Our comprehensive epigenome annotations and 3D genome maps serve as valuable resource and provide a deep understanding of the complex regulatory mechanisms of genes underlying developmental and morphological diversities between maize ear and tassel.

Genomic Insight into Differentiation and Selection Sweeps in the Improvement of Upland Cotton

Upland cotton is the most economically important fibre crop. The human-mediated selection has resulted in modern upland cultivars with higher yield and better fibre quality. However, changes in genome structure resulted from human-mediated selection are poorly understood. Comparative population genomics offers us tools to dissect the genetic history of domestication and helps to understand the genome-wide effects of human-mediated selection. Hereby, we report a comprehensive assessment of Gossypium hirsutum landraces, obsolete cultivars and modern cultivars based on high throughput genome-wide sequencing of the core set of genotypes. As a result of the genome-wide scan, we identified 93 differential regions and 311 selection sweeps associated with domestication and improvement. Furthermore, we performed genome-wide association studies to identify traits associated with the differential regions and selection sweeps. Our study provides a genetic basis to understand the domestication process in Chinese cotton cultivars. It also provides a comprehensive insight into changes in genome structure due to selection and improvement during the last century. We also identified multiple genome-wide associations (GWAS associations) for fibre yield, quality and other morphological characteristics.

Combined GWAS and eQTL analysis uncovers a genetic regulatory network orchestrating the initiation of secondary cell wall development in cotton

The cotton fibre serves as a valuable experimental system to study cell wall synthesis in plants, but our understanding of the genetic regulation of this process during fibre development remains limited. We performed a genome-wide association study (GWAS) and identified 28 genetic loci associated with fibre quality in allotetraploid cotton. To investigate the regulatory roles of these loci, we sequenced fibre transcriptomes of 251 cotton accessions and identified 15 330 expression quantitative trait loci (eQTL). Analysis of local eQTL and GWAS data prioritised 13 likely causal genes for differential fibre quality in a transcriptome-wide association study (TWAS). Characterisation of distal eQTL revealed unequal genetic regulation patterns between two subgenomes, highlighted by an eQTL hotspot (Hot216) that established a genome-wide genetic network regulating the expression of 962 genes. The primary regulatory role of Hot216, and specifically the gene encoding a KIP-related protein, was found to be the transcriptional regulation of genes responsible for cell wall synthesis, which contributes to fibre length by modulating the developmental transition from rapid cell elongation to secondary cell wall synthesis. This study uncovered the genetic regulation of fibre-cell development and revealed the molecular basis of the temporal modulation of secondary cell wall synthesis during plant cell elongation.

The Genetic Basis of Morphological Diversity in Domesticated Goldfish

Although domesticated goldfish strains exhibit highly diversified phenotypes in morphology, the genetic basis underlying these phenotypes is poorly understood. Here, based on analysis of transposable elements in the allotetraploid goldfish genome, we found that its two subgenomes have evolved asymmetrically since a whole-genome duplication event in the ancestor of goldfish and common carp. We conducted whole-genome sequencing of 27 domesticated goldfish strains and wild goldfish. We identified more than 60 million genetic variations and established a population genetic structure of major goldfish strains. Genome-wide association studies and analysis of strain-specific variants revealed genetic loci associated with several goldfish phenotypes, including dorsal fin loss, long-tail, telescope-eye, albinism, and heart-shaped tail. Our results suggest that accumulated mutations in the asymmetrically evolved subgenomes led to generation of diverse phenotypes in the goldfish domestication history. This study is a key resource for understanding the genetic basis of phenotypic diversity among goldfish strains.

Gossypium Genomics: Trends, Scope, and Utilization for Cotton Improvement

Cotton (Gossypium spp.) is the most important natural fiber crop worldwide. The diversity of Gossypium species also provides an ideal model for investigating evolution and domestication of polyploids. However, the huge and complex cotton genome hinders genomic research. Technical advances in high-throughput sequencing and bioinformatics analysis have now largely overcome these obstacles, bringing about a new era of cotton genomics. Here, we review recent progress in Gossypium genomics based on whole genome sequencing, resequencing, and comparative genomics, which have provided insights about the genomic basis of fiber biogenesis and the landscape of cotton functional genomics. We address current challenges and present multidisciplinary genomics-enabled breeding strategies covering the breadth of high fiber yield, quality, and environmental resilience for future cotton breeding programs.

Resequencing 93 accessions of coffee unveils independent and parallel selection during Coffea species divergence

Coffea arabica, C. canephora and C. excelsa, with differentiated morphological traits and distinct agro-climatic conditions, compose the majority of the global coffee plantation. To comprehensively understand their genetic diversity and divergence for future genetic improvement requires high-density markers. Here, we sequenced 93 accessions encompassing these three Coffea species, uncovering 15,367,960 single-nucleotide polymorphisms (SNPs). These SNPs are unequally distributed across different genomic regions and gene families, with two disease-resistant gene families showing the highest SNP density, suggesting strong balancing selection. Meanwhile, the allotetraploid C. arabica exhibits greater nucleotide diversity, followed by C. canephora and C. excelsa. Population divergence (F_ST), population stratification and phylogeny all support strong divergence among species, with C. arabica and its parental species C. canephora being closer genetically. Scanning of genomic islands with elevated F_ST and structure-disruptive SNPs contributing to species divergence revealed that most of the selected genes in each lineage are independent, with a few being selected in parallel for two or three species, such as genes in root hair cell development, flavonols accumulation and disease-resistant genes. Moreover, some of the SNPs associated with coffee lipids exhibit significantly biased allele frequency among species, being valuable for interspecific breeding. Overall, our study not only uncovers the key population genomic patterns among species but also contributes a substantial genomic resource for coffee breeding.

Genome sequence of Gossypium herbaceum and genome updates of Gossypium arboreum and Gossypium hirsutum provide insights into cotton A-genome evolution

Upon assembling the first Gossypium herbaceum (A1) genome and substantially improving the existing Gossypium arboreum (A2) and Gossypium hirsutum ((AD)1) genomes, we showed that all existing A-genomes may have originated from a common ancestor, referred to here as A0, which was more phylogenetically related to A1 than A2. Further, allotetraploid formation was shown to have preceded the speciation of A1 and A2. Both A-genomes evolved independently, with no ancestor-progeny relationship. Gaussian probability density function analysis indicates that several long-terminal-repeat bursts that occurred from 5.7 million years ago to less than 0.61 million years ago contributed compellingly to A-genome size expansion, speciation and evolution. Abundant species-specific structural variations in genic regions changed the expression of many important genes, which may have led to fiber cell improvement in (AD)1. Our findings resolve existing controversial concepts surrounding A-genome origins and provide valuable genomic resources for cotton genetic improvement.

Stepwise selection on homeologous PRR genes controlling flowering and maturity during soybean domestication

Adaptive changes in plant phenology are often considered to be a feature of the so-called 'domestication syndrome' that distinguishes modern crops from their wild progenitors, but little detailed evidence supports this idea. In soybean, a major legume crop, flowering time variation is well characterized within domesticated germplasm and is critical for modern production, but its importance during domestication is unclear. Here, we identify sequential contributions of two homeologous pseudo-response-regulator genes, Tof12 and Tof11, to ancient flowering time adaptation, and demonstrate that they act via LHY homologs to promote expression of the legume-specific E1 gene and delay flowering under long photoperiods. We show that Tof12-dependent acceleration of maturity accompanied a reduction in dormancy and seed dispersal during soybean domestication, possibly predisposing the incipient crop to latitudinal expansion. Better understanding of this early phase of crop evolution will help to identify functional variation lost during domestication and exploit its potential for future crop improvement.

Chromatin domains in space and their functional implications

Genome organization displays functional compartmentalization. Many factors, including epigenetic modifications, transcription factors, chromatin remodelers, and RNAs, shape chromatin domains and the three-dimensional genome organization. Various types of chromatin domains with distinct epigenetic and spatial features exhibit different transcriptional activities. As part of the efforts to better understand plant functional genomics, over the past a few years, spatial distribution patterns of plant chromatin domains have been brought to light. In this review, we discuss chromatin domains associated with the nuclear periphery and the nucleolus, as well as chromatin domains staying in proximity and showing physical interactions. The functional implication of these domains is discussed, with a particular focus on the transcriptional regulation and replication timing. Finally, from a biophysical point of view, we discuss potential roles of liquid-liquid phase separation in plant nuclei in the genesis and maintenance of spatial chromatin domains.

Identification and evolution of gene regulatory networks: insights from comparative studies in plants

The availability of genome sequences, genome-wide assays of transcription factor binding, and accessible chromatin maps have unveiled gene regulatory landscapes in plants. This understanding has ushered in comparative gene regulatory network studies that assess network rewiring between species, across time, and between biological tissues. Comparisons of cis-regulatory elements across the plant kingdom have uncovered examples of conserved sequences, but also of divergence, indicating that selective pressures can vary in different plant families. Transcription factor duplication, followed by spatiotemporal expression divergence of the duplicates, also appears to be a key mechanism of network evolution. Here, we review recent literature describing the regulation of gene expression in plants, and how comparative studies provide insights into how these regulatory interactions change and lead to gene regulatory network rewiring.

Identification of candidate genes for key fibre-related QTLs and derivation of favourable alleles in Gossypium hirsutum recombinant inbred lines with G. barbadense introgressions

Fine mapping QTLs and identifying candidate genes for cotton fibre-quality and yield traits would be beneficial to cotton breeding. Here, we constructed a high-density genetic map by specific-locus amplified fragment sequencing (SLAF-seq) to identify QTLs associated with fibre-quality and yield traits using 239 recombinant inbred lines (RILs), which was developed from LMY22 (a high-yield Gossypium hirsutumL. cultivar) × LY343 (a superior fibre-quality germplasm with G. barbadenseL. introgressions). The genetic map spanned 3426.57 cM, including 3556 SLAF-based SNPs and 199 SSR marker loci. A total of 104 QTLs, including 67 QTLs for fibre quality and 37 QTLs for yield traits, were identified with phenotypic data collected from 7 environments. Among these, 66 QTLs were co-located in 19 QTL clusters on 12 chromosomes, and 24 QTLs were detected in three or more environments and determined to be stable. We also investigated the genomic components of LY343 and their contributions to fibre-related traits by deep sequencing the whole genome of LY343, and we found that genomic components from G. hirsutum races (which entered LY343 via its G. barbadense parent) contributed more favourable alleles than those from G. barbadense. We further identified six putative candidate genes for stable QTLs, including Gh_A03G1147 (GhPEL6), Gh_D07G1598 (GhCSLC6) and Gh_D13G1921 (GhTBL5) for fibre-length QTLs and Gh_D03G0919 (GhCOBL4), Gh_D09G1659 (GhMYB4) and Gh_D09G1690 (GhMYB85) for lint-percentage QTLs. Our results provide comprehensive insight into the genetic basis of the formation of fibre-related traits and would be helpful for cloning fibre-development-related genes as well as for marker-assisted genetic improvement in cotton.

Chromatin Organization in Early Land Plants Reveals an Ancestral Association between H3K27me3, Transposons, and Constitutive Heterochromatin

Genome packaging by nucleosomes is a hallmark of eukaryotes. Histones and the pathways that deposit, remove, and read histone modifications are deeply conserved. Yet, we lack information regarding chromatin landscapes in extant representatives of ancestors of the main groups of eukaryotes, and our knowledge of the evolution of chromatin-related processes is limited. We used the bryophyte Marchantia polymorpha, which diverged from vascular plants circa 400 mya, to obtain a whole chromosome genome assembly and explore the chromatin landscape and three-dimensional genome organization in an early diverging land plant lineage. Based on genomic profiles of ten chromatin marks, we conclude that the relationship between active marks and gene expression is conserved across land plants. In contrast, we observed distinctive features of transposons and other repetitive sequences in Marchantia compared with flowering plants. Silenced transposons and repeats did not accumulate around centromeres. Although a large fraction of constitutive heterochromatin was marked by H3K9 methylation as in flowering plants, a significant proportion of transposons were marked by H3K27me3, which is otherwise dedicated to the transcriptional repression of protein-coding genes in flowering plants. Chromatin compartmentalization analyses of Hi-C data revealed that repressed B compartments were densely decorated with H3K27me3 but not H3K9 or DNA methylation as reported in flowering plants. We conclude that, in early plants, H3K27me3 played an essential role in heterochromatin function, suggesting an ancestral role of this mark in transposon silencing.

Open-Bud Duplicate Loci Are Identified as MML10s, Orthologs of MIXTA-Like Genes on Homologous Chromosomes of Allotetraploid Cotton

The open-bud (ob) mutants in cotton display abnormal flower buds with the stigma and upper anthers exposed before blooming. This characteristic is potentially useful for the efficient production of hybrid seeds. The recessive inheritance pattern of the ob phenotype in allotetraploid cotton is determined by duplicated recessive loci (ob1ob1ob2ob2). In this study, ob1, which is a MIXTA-like MYB gene on chromosome D13 (MML10_Dt), was identified by map-based cloning. In Gossypium barbadense (Gb) acc. 3-79, a single nucleotide polymorphism (SNP) (G/A) at the splice site of the first intron and an 8-bp deletion in the third exon of MML10_Dt were found, which are the causative mutations at the ob1 loci. A 1783-bp deletion that leads to the loss of the third exon and accounts for the causal variation at the ob2 loci was found in MML10_At of Gossypium hirsutum (Gh) acc. TM-1. The ob phenotype results from the combination of these two loss-of-function loci. Genotyping assays showed that the ob1 and ob2 loci appeared after the formation of allotetraploid cotton and were specific for Gb and Gh, respectively. All Gb lines and most Gh cultivars carry the single corresponding mutant alleles. Genome-wide transcriptome analysis showed that some of the MYB genes and genes related to cell wall biogenesis, trichome differentiation, cytokinin signal transduction, and cell division were repressed in the ob mutants, which may lead to suppression of petal growth. These findings should be of value for breeding superior ob lines in cotton.

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

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

Comparison of Arachis monticola with Diploid and Cultivated Tetraploid Genomes Reveals Asymmetric Subgenome Evolution and Improvement of Peanut

Like many important crops, peanut is a polyploid that underwent polyploidization, evolution, and domestication. The wild allotetraploid peanut species Arachis monticola (A. monticola) is an important and unique link from the wild diploid species to cultivated tetraploid species in the Arachis lineage. However, little is known about A. monticola and its role in the evolution and domestication of this important crop. A fully annotated sequence of ≈2.6 Gb A. monticola genome and comparative genomics of the Arachis species is reported. Genomic reconstruction of 17 wild diploids from AA, BB, EE, KK, and CC groups and 30 tetraploids demonstrates a monophyletic origin of A and B subgenomes in allotetraploid peanuts. The wild and cultivated tetraploids undergo asymmetric subgenome evolution, including homoeologous exchanges, homoeolog expression bias, and structural variation (SV), leading to subgenome functional divergence during peanut domestication. Significantly, SV-associated homoeologs tend to show expression bias and correlation with pod size increase from diploids to wild and cultivated tetraploids. Moreover, genomic analysis of disease resistance genes shows the unique alleles present in the wild peanut can be introduced into breeding programs to improve some resistance traits in the cultivated peanuts. These genomic resources are valuable for studying polyploid genome evolution, domestication, and improvement of peanut production and resistance.

Tissue-specific Hi-C analyses of rice, foxtail millet and maize suggest non-canonical function of plant chromatin domains

Chromatins are not randomly packaged in the nucleus and their organization plays important roles in transcription regulation, which is best studied in the mammalian models. Using in situ Hi-C, we have compared the 3D chromatin architectures of rice mesophyll and endosperm, foxtail millet bundle sheath and mesophyll, and maize bundle sheath, mesophyll and endosperm tissues. We found that their global A/B compartment partitions are stable across tissues, while local A/B compartment has tissue-specific dynamic associated with differential gene expression. Plant domains are largely stable across tissues, while new domain border formations are often associated with transcriptional activation in the region. Genes inside plant domains are not conserved across species, and lack significant co-expression behavior unlike those in mammalian TADs. Although we only observed chromatin loops between gene islands in the large genomes, the maize loop gene pairs' syntenic orthologs have shorter physical distances in small genome monocots, suggesting that loops instead of domains might have conserved biological function. Our study showed that plants' chromatin features might not have conserved biological functions as the mammalian ones.

Identifying statistically significant chromatin contacts from Hi-C data with FitHiC2

Fit-Hi-C is a programming application to compute statistical confidence estimates for Hi-C contact maps to identify significant chromatin contacts. By fitting a monotonically non-increasing spline, Fit-Hi-C captures the relationship between genomic distance and contact probability without any parametric assumption. The spline fit together with the correction of contact probabilities with respect to bin- or locus-specific biases accounts for previously characterized covariates impacting Hi-C contact counts. Fit-Hi-C is best applied for the study of mid-range (e.g., 20 kb-2 Mb for human genome) intra-chromosomal contacts; however, with the latest reimplementation, named FitHiC2, it is possible to perform genome-wide analysis for high-resolution Hi-C data, including all intra-chromosomal distances and inter-chromosomal contacts. FitHiC2 also offers a merging filter module, which eliminates indirect/bystander interactions, leading to significant reduction in the number of reported contacts without sacrificing recovery of key loops such as those between convergent CTCF binding sites. Here, we describe how to apply the FitHiC2 protocol to three use cases: (i) 5-kb resolution Hi-C data of chromosome 5 from GM12878 (a human lymphoblastoid cell line), (ii) 40-kb resolution whole-genome Hi-C data from IMR90 (human lung fibroblast), and (iii) budding yeast whole-genome Hi-C data at a single restriction cut site (EcoRI) resolution. The procedure takes ~12 h with preprocessing when all use cases are run sequentially (~4 h when run parallel). With the recent improvements in its implementation, FitHiC2 (8 processors and 16 GB memory) is also scalable to genome-wide analysis of the highest resolution (1 kb) Hi-C data available to date (~48 h with 32 GB peak memory). FitHiC2 is available through Bioconda, GitHub and the Python Package Index.

Genome-wide association reveals genetic variation of lint yield components under salty field conditions in cotton (Gossypium hirsutum L.)

BACKGROUND

Salinity is one of the most significant environmental factors limiting the productivity of cotton. However, the key genetic components responsible for the reduction in cotton yield in saline-alkali soils are still unclear.

RESULTS

Here, we evaluated three main components of lint yield, single boll weight (SBW), lint percentage (LP) and boll number per plant (BNPP), across 316 G. hirsutum accessions under four salt conditions over two years. Phenotypic analysis indicated that LP was unchanged under different salt conditions, however BNPP decreased significantly and SBW increased slightly under high salt conditions. Based on 57,413 high-quality single nucleotide polymorphisms (SNPs) and genome-wide association study (GWAS) analysis, a total of 42, 91 and 25 stable quantitative trait loci (QTLs) were identified for SBW, LP and BNPP, respectively. Phenotypic and QTL analysis suggested that there was little correlation among the three traits. For LP, 8 stable QTLs were detected simultaneously in four different salt conditions, while fewer repeated QTLs for SBW or BNPP were identified. Gene Ontology (GO) analysis indicated that their regulatory mechanisms were also quite different. Via transcriptome profile data, we detected that 10 genes from the 8 stable LP QTLs were predominantly expressed during fiber development. Further, haplotype analyses found that a MYB gene (GhMYB103), with the two SNP variations in cis-regulatory and coding regions, was significantly correlated with lint percentage, implying a crucial role in lint yield. We also identified that 40 candidate genes from BNPP QTLs were salt-inducible. Genes related to carbohydrate metabolism and cell structure maintenance were rich in plants grown in high salt conditions, while genes related to ion transport were active in plants grown in low salt conditions, implying different regulatory mechanisms for BNPP at high and low salt conditions.

CONCLUSIONS

This study provides a foundation for elucidating cotton salt tolerance mechanisms and contributes gene resources for developing upland cotton varieties with high yields and salt stress tolerance.

Molecular Evolution and Expansion of the KUP Family in the Allopolyploid Cotton Species Gossypium hirsutum and Gossypium barbadense

The comprehensive analysis of gene family evolution will elucidate the origin and evolution of gene families. The K+ uptake (KUP) gene family plays important roles in K+ uptake and transport, plant growth and development, and abiotic stress responses. However, the current understanding of the KUP family in cotton is limited. In this study, 51 and 53 KUPs were identified in Gossypium barbadense and Gossypium hirsutum, respectively. These KUPs were divided into five KUP subfamilies, with subfamily 2 containing three groups. Different subfamilies had different member numbers, conserved motifs, gene structures, regulatory elements, and gene expansion and loss rates. A paleohexaploidization event caused the expansion of GhKUP and GbKUP in cotton, and duplication events in G. hirsutum and G. barbadense have happened in a common ancestor of Gossypium. Meanwhile, the KUP members of the two allopolyploid subgenomes of G. hirsutum and G. barbadense exhibited unequal gene proportions, gene structural diversity, uneven chromosomal distributions, asymmetric expansion rates, and biased gene loss rates. In addition, the KUP families of G. hirsutum and G. barbadense displayed evolutionary conservation and divergence. Taken together, these results illustrated the molecular evolution and expansion of the KUP family in allopolyploid cotton species.

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

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

Revealing Distinctions in Genetic Diversity and Adaptive Evolution Between Two Varieties of Camellia sinensis by Whole-Genome Resequencing

Camellia sinensis var. sinensis (CSS) and C. sinensis var. assamica (CSA) are the two most economically important tea varieties. They have different characteristics and geographical distribution. Their genetic diversity and differentiation are unclear. Here, we identified 18,903,625 single nucleotide polymorphisms (SNPs) and 7,314,133 insertion-deletion mutations (indels) by whole-genome resequencing of 30 cultivated and three wild related species. Population structure and phylogenetic tree analyses divided the cultivated accessions into CSS and CSA containing 6,440,419 and 6,176,510 unique variations, respectively. The CSS subgroup possessed higher genetic diversity and was enriched for rare alleles. The CSA subgroup had more non-synonymous mutations and might have experienced a greater degree of balancing selection. The evolution rate (dN/dS) and KEGG enrichment indicated that genes involved in the synthesis and metabolism of flavor substances were positively selected in both CSS and CSA subpopulations. However, there are extensive genome differentiation regions (2959 bins and approximately 148 M in size) between the two subgroups. Compared with CSA (141 selected regions containing 124 genes), the CSS subgroup (830 selected regions containing 687 genes) displayed more selection regions potentially related to environmental adaptability. Fifty-three pairs of polymorphic indel markers were developed. Some markers were located in hormone-related genes with distinct alleles in the two cultivated subgroups. These identified variations and selected regions provide clues for the differentiation and adaptive evolution of tea varieties. The newly developed indel markers will be valuable in further genetic research on tea plants.

Genome-wide quantitative trait loci reveal the genetic basis of cotton fibre quality and yield-related traits in a Gossypium hirsutum recombinant inbred line population

Cotton is widely cultivated globally because it provides natural fibre for the textile industry and human use. To identify quantitative trait loci (QTLs)/genes associated with fibre quality and yield, a recombinant inbred line (RIL) population was developed in upland cotton. A consensus map covering the whole genome was constructed with three types of markers (8295 markers, 5197.17 centimorgans (cM)). Six fibre yield and quality traits were evaluated in 17 environments, and 983 QTLs were identified, 198 of which were stable and mainly distributed on chromosomes 4, 6, 7, 13, 21 and 25. Thirty-seven QTL clusters were identified, in which 92.8% of paired traits with significant medium or high positive correlations had the same QTL additive effect directions, and all of the paired traits with significant medium or high negative correlations had opposite additive effect directions. In total, 1297 genes were discovered in the QTL clusters, 414 of which were expressed in two RNA-Seq data sets. Many genes were discovered, 23 of which were promising candidates. Six important QTL clusters that included both fibre quality and yield traits were identified with opposite additive effect directions, and those on chromosome 13 (qClu-chr13-2) could increase fibre quality but reduce yield; this result was validated in a natural population using three markers. These data could provide information about the genetic basis of cotton fibre quality and yield and help cotton breeders to improve fibre quality and yield simultaneously.

Insights into Drought Stress Signaling in Plants and the Molecular Genetic Basis of Cotton Drought Tolerance

Drought stress restricts plant growth and development by altering metabolic activity and biological functions. However, plants have evolved several cellular and molecular mechanisms to overcome drought stress. Drought tolerance is a multiplex trait involving the activation of signaling mechanisms and differentially expressed molecular responses. Broadly, drought tolerance comprises two steps: stress sensing/signaling and activation of various parallel stress responses (including physiological, molecular, and biochemical mechanisms) in plants. At the cellular level, drought induces oxidative stress by overproduction of reactive oxygen species (ROS), ultimately causing the cell membrane to rupture and stimulating various stress signaling pathways (ROS, mitogen-activated-protein-kinase, Ca2+, and hormone-mediated signaling). Drought-induced transcription factors activation and abscisic acid concentration co-ordinate the stress signaling and responses in cotton. The key responses against drought stress, are root development, stomatal closure, photosynthesis, hormone production, and ROS scavenging. The genetic basis, quantitative trait loci and genes of cotton drought tolerance are presented as examples of genetic resources in plants. Sustainable genetic improvements could be achieved through functional genomic approaches and genome modification techniques such as the CRISPR/Cas9 system aid the characterization of genes, sorted out from stress-related candidate single nucleotide polymorphisms, quantitative trait loci, and genes. Exploration of the genetic basis for superior candidate genes linked to stress physiology can be facilitated by integrated functional genomic approaches. We propose a third-generation sequencing approach coupled with genome-wide studies and functional genomic tools, including a comparative sequenced data (transcriptomics, proteomics, and epigenomic) analysis, which offer a platform to identify and characterize novel genes. This will provide information for better understanding the complex stress cellular biology of plants.

Unraveling cis and trans regulatory evolution during cotton domestication

Cis and trans regulatory divergence underlies phenotypic and evolutionary diversification. Relatively little is understood about the complexity of regulatory evolution accompanying crop domestication, particularly for polyploid plants. Here, we compare the fiber transcriptomes between wild and domesticated cotton (Gossypium hirsutum) and their reciprocal F1 hybrids, revealing genome-wide (~15%) and often compensatory cis and trans regulatory changes under divergence and domestication. The high level of trans evolution (54%-64%) observed is likely enabled by genomic redundancy following polyploidy. Our results reveal that regulatory variation is significantly associated with sequence evolution, inheritance of parental expression patterns, co-expression gene network properties, and genomic loci responsible for domestication traits. With respect to regulatory evolution, the two subgenomes of allotetraploid cotton are often uncoupled. Overall, our work underscores the complexity of regulatory evolution during fiber domestication and may facilitate new approaches for improving cotton and other polyploid plants.

A Genome-Wide Association Study Revealed Key SNPs/Genes Associated With Salinity Stress Tolerance In Upland Cotton

Millions of hectares of land are too saline to produce economically valuable crop yields. Salt tolerance in cotton is an imperative approach for improvement in response to ever-increasing soil salinization. Little is known about the genetic basis of salt tolerance in cotton at the seedling stage. To address this issue, a genome-wide association study (GWAS) was conducted on a core collection of a genetically diverse population of upland cotton (Gossypium hirsutum L.) comprising of 419 accessions, representing various geographic origins, including China, USA, Pakistan, the former Soviet Union, Chad, Australia, Brazil, Mexico, Sudan, and Uganda. Phenotypic evaluation of 7 traits under control (0 mM) and treatment (150 mM) NaCl conditions depicted the presence of broad natural variation in the studied population. The association study was carried out with the efficient mixed-model association eXpedited software package. A total of 17,264 single-nucleotide polymorphisms (SNPs) associated with different salinity stress tolerance related traits were found. Twenty-three candidate SNPs related to salinity stress-related traits were selected. Final key SNPs were selected based on the r² value with nearby SNPs in a linkage disequilibrium (LD) block. Twenty putative candidate genes surrounding SNPs, A10_95330133 and D10_61258588, associated with leaf relative water content, RWC_150, and leaf fresh weight, FW_150, were identified, respectively. We further validated the expression patterns of twelve candidate genes with qRT-PCR, which revealed different expression levels in salt-tolerant and salt-sensitive genotypes. The results of our GWAS provide useful knowledge about the genetic control of salt tolerance at the seedling stage, which could assist in elucidating the genetic and molecular mechanisms of salinity stress tolerance in cotton plants.

Genome-wide identification and characterization of TALE superfamily genes in cotton reveals their functions in regulating secondary cell wall biosynthesis

BACKGROUND

Cotton fiber length and strength are both key traits of fiber quality, and fiber strength (FS) is tightly correlated with secondary cell wall (SCW) biosynthesis. The three-amino-acid-loop-extension (TALE) superclass homeoproteins are involved in regulating diverse biological processes in plants, and some TALE members has been identified to play a key role in regulating SCW formation. However, little is known about the functions of TALE members in cotton (Gossypium spp.).

RESULTS

In the present study, based on gene homology, 46, 47, 88 and 94 TALE superfamily genes were identified in G. arboreum, G. raimondii, G. barbadense and G. hirsutum, respectively. Phylogenetic and evolutionary analysis showed the evolutionary conservation of two cotton TALE families (including BEL1-like and KNOX families). Gene structure analysis also indicated the conservation of GhTALE members under selection. The analysis of promoter cis-elements and expression patterns suggested potential transcriptional regulation functions in fiber SCW biosynthesis and responses to some phytohormones for GhTALE proteins. Genome-wide analysis of colocalization of TALE transcription factors with SCW-related QTLs revealed that some BEL1-like genes and KNAT7 homologs may participate in the regulation of cotton fiber strength formation. Overexpression of GhKNAT7-A03 and GhBLH6-A13 significantly inhibited the synthesis of lignocellulose in interfascicular fibers of Arabidopsis. Yeast two-hybrid (Y2H) experiments showed extensive heteromeric interactions between GhKNAT7 homologs and some GhBEL1-like proteins. Yeast one-hybrid (Y1H) experiments identified the upstream GhMYB46 binding sites in the promoter region of GhTALE members and defined the downstream genes that can be directly bound and regulated by GhTALE heterodimers.

CONCLUSION

We comprehensively identified TALE superfamily genes in cotton. Some GhTALE members are predominantly expressed during the cotton fiber SCW thicking stage, and may genetically correlated with the formation of FS. Class II KNOX member GhKNAT7 can interact with some GhBEL1-like members to form the heterodimers to regulate the downstream targets, and this regulatory relationship is partially conserved with Arabidopsis. In summary, this study provides important clues for further elucidating the functions of TALE genes in regulating cotton growth and development, especially in the fiber SCW biosynthesis network, and it also contributes genetic resources to the improvement of cotton fiber quality.

Multi-strategic RNA-seq analysis reveals a high-resolution transcriptional landscape in cotton

Cotton is an important natural fiber crop, however, its comprehensive and high-resolution gene map is lacking. Here we integrate four complementary high-throughput techniques, including Pacbio long read Iso-seq, strand-specific RNA-seq, CAGE-seq, and PolyA-seq, to systematically explore the transcription landscape across 16 tissues or different organ types in Gossypium arboreum. We devise a computational pipeline, named IGIA, to reconstruct accurate gene structures from the integrated data. Our results reveal a dynamic and diverse transcriptional map in cotton: tissue-specific gene expression, alternative usage of TSSs and polyadenylation sites, hotspot of alternative splicing, and transcriptional read-through. These regulated events affect many genes in various aspects such as gain or loss of functional RNA motifs and protein domains, fine-tuning of DNA binding activity, and co-regulation for genes in the same complex or pathway. The methods and findings provide valuable resources for further functional genomic studies such as understanding natural SNP variations for plant community.

The allotetraploid origin and asymmetrical genome evolution of the common carp Cyprinus carpio

Common carp (Cyprinus carpio) is an allotetraploid species derived from recent whole genome duplication and provides a model to study polyploid genome evolution in vertebrates. Here, we generate three chromosome-level reference genomes of C. carpio and compare to related diploid Cyprinid genomes. We identify a Barbinae lineage as potential diploid progenitor of C. carpio and then divide the allotetraploid genome into two subgenomes marked by a distinct genome similarity to the diploid progenitor. We estimate that the two diploid progenitors diverged around 23 Mya and merged around 12.4 Mya based on the divergence rates of homoeologous genes and transposable elements in two subgenomes. No extensive gene losses are observed in either subgenome. Instead, we find gene expression bias across surveyed tissues such that subgenome B is more dominant in homoeologous expression. CG methylation in promoter regions may play an important role in altering gene expression in allotetraploid C. carpio.

The common carp is derived from recent whole genome duplication and represents a model for polyploid genome evolution, rare in vertebrates. Here, the authors generate and analyse chromosome-level reference genomes for common carp, and describe subgenome gene expression changes.

Evolutionary transcriptomics reveals the origins of olives and the genomic changes associated with their domestication

The olive (Olea europaea L. subsp. europaea) is one of the oldest and most socio-economically important cultivated perennial crop in the Mediterranean region. Yet, its origins are still under debate and the genetic bases of the phenotypic changes associated with its domestication are unknown. We generated RNA-sequencing data for 68 wild and cultivated olive trees to study the genetic diversity and structure both at the transcription and sequence levels. To localize putative genes or expression pathways targeted by artificial selection during domestication, we employed a two-step approach in which we identified differentially expressed genes and screened the transcriptome for signatures of selection. Our analyses support a major domestication event in the eastern part of the Mediterranean basin followed by dispersion towards the West and subsequent admixture with western wild olives. While we found large changes in gene expression when comparing cultivated and wild olives, we found no major signature of selection on coding variants and weak signals primarily affected transcription factors. Our results indicated that the domestication of olives resulted in only moderate genomic consequences and that the domestication syndrome is mainly related to changes in gene expression, consistent with its evolutionary history and life history traits.