Transposable elements play an important role during cotton genome evolution and fiber cell development

Cotton transposon-related variome reveals roles of transposon-related variations in modern cotton cultivation

INTRODUCTION

Transposon plays a vital role in cotton genome evolution, contributing to the expansion and divergence of genomes within the Gossypium genus. However, knowledge of transposon activity in modern cotton cultivation is limited.

OBJECTIVES

In this study, we aimed to construct transposon-related variome within Gossypium genus and reveal role of transposon-related variations during cotton cultivation. In addition, we try to identify valuable transposon-related variations for cotton breeding.

METHODS

We utilized graphical genome construction to build up the graphical transposon-related variome. Based on the graphical variome, we integrated t-test, eQTL analysis and Mendelian Randomization (MR) to identify valuable transposon activities and elite genes. In addition, a convolutional neural network (CNN) model was constructed to evaluate epigenomic effects of transposon-related variations.

RESULTS

We identified 35,980 transposon activities among 10 cotton genomes, and the diversity of genomic and epigenomic features was observed among 21 transposon categories. The graphical cotton transposon-related variome was constructed, and 9,614 transposon-related variations with plasticity in the modern cotton cohort were used for eQTL, phenotypic t-test and Mendelian Randomization. 128 genes were identified as gene resources improving fiber length and strength simultaneously. 4 genes were selected from 128 genes to construct the elite gene panel whose utility has been validated in a natural cotton cohort and 2 accessions with phenotypic divergence. Based on the eQTL analysis results, we identified transposon-related variations involved in cotton's environmental adaption and human domestication, providing evidence of their role in cotton's adaption-domestication cooperation.

CONCLUSIONS

The cotton transposon-related variome revealed the role of transposon-related variations in modern cotton cultivation, providing genomic resources for cotton molecular breeding.

Pangenome analysis reveals transposon-driven genome evolution in cotton

BACKGROUND

Transposable elements (TEs) have a profound influence on the trajectory of plant evolution, driving genome expansion and catalyzing phenotypic diversification. The pangenome, a comprehensive genetic pool encompassing all variations within a species, serves as an invaluable tool, unaffected by the confounding factors of intraspecific diversity. This allows for a more nuanced exploration of plant TE evolution.

RESULTS

Here, we constructed a pangenome for diploid A-genome cotton using 344 accessions from representative geographical regions, including 223 from China as the main component. We found 511 Mb of non-reference sequences (NRSs) and revealed the presence of 5479 previously undiscovered protein-coding genes. Our comprehensive approach enabled us to decipher the genetic underpinnings of the distinct geographic distributions of cotton. Notably, we identified 3301 presence-absence variations (PAVs) that are closely tied to gene expression patterns within the pangenome, among which 2342 novel expression quantitative trait loci (eQTLs) were found residing in NRSs. Our investigation also unveiled contrasting patterns of transposon proliferation between diploid and tetraploid cotton, with long terminal repeat (LTR) retrotransposons exhibiting a synchronized surge in polyploids. Furthermore, the invasion of LTR retrotransposons from the A subgenome to the D subgenome triggered a substantial expansion of the latter following polyploidization. In addition, we found that TE insertions were responsible for the loss of 36.2% of species-specific genes, as well as the generation of entirely new species-specific genes.

CONCLUSIONS

Our pangenome analyses provide new insights into cotton genomics and subgenome dynamics after polyploidization and demonstrate the power of pangenome approaches for elucidating transposon impacts and genome evolution.

Phased gap-free genome assembly of octoploid cultivated strawberry illustrates the genetic and epigenetic divergence among subgenomes

The genetic and epigenetic mechanisms underlying the coexistence and coordination of the four diverged subgenomes (ABCD) in octoploid strawberries (Fragaria × ananassa) remains poorly understood. In this study, we have assembled a haplotype-phased gap-free octoploid genome for the strawberry, which allowed us to uncover the sequence, structure, and epigenetic divergences among the subgenomes. The diploid progenitors of the octoploid strawberry, apart from subgenome A (Fragaria vesca), have been a subject of public controversy. Phylogenomic analyses revealed a close relationship between diploid species Fragaria iinumae and subgenomes B, C, and D. Subgenome A, closely related to F. vesca, retains the highest number of genes, exhibits the lowest content of transposable elements (TEs), experiences the strongest purifying selection, shows the lowest DNA methylation levels, and displays the highest expression level compared to the other three subgenomes. Transcriptome and DNA methylome analyses revealed that subgenome A-biased genes were enriched in fruit development biological processes. In contrast, although subgenomes B, C, and D contain equivalent amounts of repetitive sequences, they exhibit diverged methylation levels, particularly for TEs located near genes. Taken together, our findings provide valuable insights into the evolutionary patterns of subgenome structure, divergence and epigenetic dynamics in octoploid strawberries, which could be utilized in strawberry genetics and breeding research.

A comprehensive overview of cotton genomics, biotechnology and molecular biological studies

Cotton is an irreplaceable economic crop currently domesticated in the human world for its extremely elongated fiber cells specialized in seed epidermis, which makes it of high research and application value. To date, numerous research on cotton has navigated various aspects, from multi-genome assembly, genome editing, mechanism of fiber development, metabolite biosynthesis, and analysis to genetic breeding. Genomic and 3D genomic studies reveal the origin of cotton species and the spatiotemporal asymmetric chromatin structure in fibers. Mature multiple genome editing systems, such as CRISPR/Cas9, Cas12 (Cpf1) and cytidine base editing (CBE), have been widely used in the study of candidate genes affecting fiber development. Based on this, the cotton fiber cell development network has been preliminarily drawn. Among them, the MYB-bHLH-WDR (MBW) transcription factor complex and IAA and BR signaling pathway regulate the initiation; various plant hormones, including ethylene, mediated regulatory network and membrane protein overlap fine-regulate elongation. Multistage transcription factors targeting CesA 4, 7, and 8 specifically dominate the whole process of secondary cell wall thickening. And fluorescently labeled cytoskeletal proteins can observe real-time dynamic changes in fiber development. Furthermore, research on the synthesis of cotton secondary metabolite gossypol, resistance to diseases and insect pests, plant architecture regulation, and seed oil utilization are all conducive to finding more high-quality breeding-related genes and subsequently facilitating the cultivation of better cotton varieties. This review summarizes the paramount research achievements in cotton molecular biology over the last few decades from the above aspects, thereby enabling us to conduct a status review on the current studies of cotton and provide strong theoretical support for the future direction.

Functional characterization of TBL genes revealed the role of GhTBL7 and GhTBL58 in cotton fiber elongation

TBL (Trichome Birefringence Like) gene family members are involved in trichome initiation and xylan acetylation in several plant species. In our research, we identified 102 TBLs from G. hirsutum. The phylogenetic tree classified TBL genes into five groups. Collinearity analysis of TBL genes indicated 136 paralogous gene pairs in G. hirsutum. Gene duplication indicated that WGD or segmental duplication contributed to the GhTBL gene family expansion. Promoter cis-elements of GhTBLs were related to growth and development, seed-specific regulation, light, and stress responses. GhTBL genes (GhTBL7, GhTBL15, GhTBL21, GhTBL25, GhTBL45, GhTBL54, GhTBL67, GhTBL72, and GhTBL77) exhibited upregulated response under exposure to cold, heat, NaCl, and PEG. GhTBL genes exhibited high expression during fiber development stages. Two GhTBL genes (GhTBL7 and GhTBL58) showed differential expression at 10 DPA fiber, as 10 DPA is a fast fiber elongation stage and fiber elongation is a very important stage of cotton fiber development. Subcellular localization of GhTBL7 and GhTBL58 revealed that these genes reside inside the cell membrane. Promoter GUS activity of GhTBL7 and GhTBL58 exhibited deep staining in roots. To further validate the role of these genes in cotton fiber elongation, we silenced these genes and observed a significant reduction in the fiber length at 10 DPA. In conclusion, the functional study of cell membrane-associated genes (GhTBL7 and GhTBL58) showed deep staining in root tissues and potential function during cotton fiber elongation at 10 DPA fiber.

Gene B5 in Cotton Confers High and Broad Resistance to Bacterial Blight and Conditions High Amounts of Sesquiterpenoid Phytoalexins

Bacterial blight resistance gene B₅ has received little attention since it was first described in 1950. A near-isogenic line (NIL) of Gossypium hirsutum cotton, AcB₅, was generated in an otherwise bacterial-blight-susceptible 'Acala 44' background. The introgressed locus B₅ in AcB₅ conferred strong and broad-spectrum resistance to bacterial blight. Segregation patterns of test crosses under Oklahoma field conditions indicated that AcB₅ is likely homozygous for resistance at two loci with partial dominance gene action. In controlled-environment conditions, two of the four copies of B₅ were required for effective resistance. Contrary to expectations of gene-for-gene theory, AcB₅ conferred high resistance toward isogenic strains of Xanthomonas citri subsp. malvacearum carrying cloned avirulence genes avrB4, avrb7, avrBIn, avrB101, and avrB102, respectively, and weaker resistance toward the strain carrying cloned avrb6. The hypothesis that each B gene, in the absence of a polygenic complex, triggers sesquiterpenoid phytoalexin production was tested by measurement of cadalene and lacinilene phytoalexins during resistant responses in five NILs carrying different B genes, four other lines carrying multiple resistance genes, as well as susceptible Ac44E. Phytoalexin production was an obvious, but variable, response in all nine resistant lines. AcB₅ accumulated an order of magnitude more of all four phytoalexins than any of the other resistant NILs. Its total levels were comparable to those detected in OK1.2, a highly resistant line that possesses several B genes in a polygenic background.

Characterization of centromeric DNA of Gossypium anomalum reveals sequence-independent enrichment dynamics of centromeric repeats

Centromeres in eukaryotes are composed of highly repetitive DNAs, which evolve rapidly and are thought to achieve a favorable structure in mature centromeres. However, how the centromeric repeat evolves into an adaptive structure is largely unknown. We characterized the centromeric sequences of Gossypium anomalum through chromatin immunoprecipitation against CENH3 antibodies. We revealed that the G. anomalum centromeres contained only retrotransposon-like repeats but were depleted in long arrays of satellites. These retrotransposon-like centromeric repeats were present in the African-Asian and Australian lineage species, suggesting that they might have arisen in the common ancestor of these diploid species. Intriguingly, we observed a substantial increase and decrease in copy numbers among African-Asian and Australian lineages, respectively, for the retrotransposon-derived centromeric repeats without apparent structure or sequence variation in cotton. This result indicates that the sequence content is not a decisive aspect of the adaptive evolution of centromeric repeats or at least retrotransposon-like centromeric repeats. In addition, two active genes with potential roles in gametogenesis or flowering were identified in CENH3 nucleosome-binding regions. Our results provide new insights into the constitution of centromeric repetitive DNA and the adaptive evolution of centromeric repeats in plants.

Annotation of the Turnera subulata (Passifloraceae) Draft Genome Reveals the S-Locus Evolved after the Divergence of Turneroideae from Passifloroideae in a Stepwise Manner

A majority of Turnera species (Passifloraceae) exhibit distyly, a reproductive system involving both self-incompatibility and reciprocal herkogamy. This system differs from self-incompatibility in Passiflora species. The genetic basis of distyly in Turnera is a supergene, restricted to the S-morph, and containing three S-genes. How supergenes and distyly evolved in Turnera, and the other Angiosperm families exhibiting distyly remain largely unknown. Unraveling the evolutionary origins in Turnera requires the generation of genomic resources and extensive phylogenetic analyses. Here, we present the annotated draft genome of the S-morph of distylous Turnera subulata. Our annotation allowed for phylogenetic analyses of the three S-genes' families across 56 plant species ranging from non-seed plants to eudicots. In addition to the phylogenetic analysis, we identified the three S-genes' closest paralogs in two species of Passiflora. Our analyses suggest that the S-locus evolved after the divergence of Passiflora and Turnera. Finally, to provide insights into the neofunctionalization of the S-genes, we compared expression patterns of the S-genes with close paralogs in Arabidopsis and Populus trichocarpa. The annotation of the T. subulata genome will provide a useful resource for future comparative work. Additionally, this work has provided insights into the convergent nature of distyly and the origin of supergenes.

Polyploidy before and after domestication of crop species

Recent advances in the genomics of polyploid species answer some of the long-standing questions about the role of polyploidy in crop species. Here, we summarize the current literature to reexamine scenarios in which polyploidy played a role both before and after domestication. The prevalence of polyploidy can help to explain environmental robustness in agroecosystems. This review also clarifies the molecular basis of some agriculturally advantageous traits of polyploid crops, including yield increments in polyploid cotton via subfunctionalization, modification of a separated sexuality to selfing in polyploid persimmon via neofunctionalization, and transition to a selfing system via nonfunctionalization combined with epistatic interaction between duplicated S-loci. The rapid progress in genomics and genetics is discussed along with how this will facilitate functional studies of understudied polyploid crop species.

Function deficiency of GhOMT1 causes anthocyanidins over-accumulation and diversifies fibre colours in cotton (Gossypium hirsutum)

Naturally coloured cotton (NCC) fibres need little or no dyeing process in textile industry to low-carbon emission and are environment-friendly. Proanthocyanidins (PAs) and their derivatives were considered as the main components causing fibre coloration and made NCCs very popular and healthy, but the monotonous fibre colours greatly limit the wide application of NCCs. Here a G. hirsutum empurpled mutant (HS2) caused by T-DNA insertion is found to enhance the anthocyanidins biosynthesis and accumulate anthocyanidins in the whole plant. HPLC and LC/MS-ESI analysis confirmed the anthocyanidins methylation and peonidin, petunidin and malvidin formation are blocked. The deficiency of GhOMT1 in HS2 was associated with the activation of the anthocyanidin biosynthesis and the altered components of anthocyanidins. The transcripts of key genes in anthocyanidin biosynthesis pathway are significantly up-regulated in HS2, while transcripts of the genes for transport and decoration were at similar levels as in WT. To investigate the potential mechanism of GhOMT1 deficiency in cotton fibre coloration, HS2 mutant was crossed with NCCs. Surprisingly, offsprings of HS2 and NCCs enhanced PAs biosynthesis and increased PAs levels in their fibres from the accumulated anthocyanidins through up-regulated GhANR and GhLAR. As expected, multiple novel lines with improved fibre colours including orange red and navy blue were produced in their generations. Based on this work, a new strategy for breeding diversified NCCs was brought out by promoting PA biosynthesis. This work will help shed light on mechanisms of PA biosynthesis and bring out potential molecular breeding strategy to increase PA levels in NCCs.

An Overview of Cotton Gland Development and Its Transcriptional Regulation

Cotton refers to species in the genus Gossypium that bear spinnable seed coat fibers. A total of 50 species in the genus Gossypium have been described to date. Of these, only four species, viz. Gossypium, hirsutum, G. barbadense, G. arboretum, and G. herbaceum are cultivated; the rest are wild. The black dot-like structures on the surfaces of cotton organs or tissues, such as the leaves, stem, calyx, bracts, and boll surface, are called gossypol glands or pigment glands, which store terpenoid aldehydes, including gossypol. The cotton (Gossypium hirsutum) pigment gland is a distinctive structure that stores gossypol and its derivatives. It provides an ideal system for studying cell differentiation and organogenesis. However, only a few genes involved in the process of gland formation have been identified to date, and the molecular mechanisms underlying gland initiation remain unclear. The terpenoid aldehydes in the lysigenous glands of Gossypium species are important secondary phytoalexins (with gossypol being the most important) and one of the main defenses of plants against pests and diseases. Here, we review recent research on the development of gossypol glands in Gossypium species, the regulation of the terpenoid aldehyde biosynthesis pathway, discoveries from genetic engineering studies, and future research directions.

Recent Advances and Future Perspectives in Cotton Research

Cotton is not only the world's most important natural fiber crop, but it is also an ideal system in which to study genome evolution, polyploidization, and cell elongation. With the assembly of five different cotton genomes, a cotton-specific whole-genome duplication with an allopolyploidization process that combined the A- and D-genomes became evident. All existing A-genomes seemed to originate from the A₀-genome as a common ancestor, and several transposable element bursts contributed to A-genome size expansion and speciation. The ethylene production pathway is shown to regulate fiber elongation. A tip-biased diffuse growth mode and several regulatory mechanisms, including plant hormones, transcription factors, and epigenetic modifications, are involved in fiber development. Finally, we describe the involvement of the gossypol biosynthetic pathway in the manipulation of herbivorous insects, the role of GoPGF in gland formation, and host-induced gene silencing for pest and disease control. These new genes, modules, and pathways will accelerate the genetic improvement of cotton.

Full-length annotation with multistrategy RNA-seq uncovers transcriptional regulation of lncRNAs in cotton

Long noncoding RNAs (lncRNAs) are crucial factors during plant development and environmental responses. To build an accurate atlas of lncRNAs in the diploid cotton Gossypium arboreum, we combined Isoform-sequencing, strand-specific RNA-seq (ssRNA-seq), and cap analysis gene expression (CAGE-seq) with PolyA-seq and compiled a pipeline named plant full-length lncRNA to integrate multi-strategy RNA-seq data. In total, 9,240 lncRNAs from 21 tissue samples were identified. 4,405 and 4,805 lncRNA transcripts were supported by CAGE-seq and PolyA-seq, respectively, among which 6.7% and 7.2% had multiple transcription start sites (TSSs) and transcription termination sites (TTSs). We revealed that alternative usage of TSS and TTS of lncRNAs occurs pervasively during plant growth. Besides, we uncovered that many lncRNAs act in cis to regulate adjacent protein-coding genes (PCGs). It was especially interesting to observe 64 cases wherein the lncRNAs were involved in the TSS alternative usage of PCGs. We identified lncRNAs that are coexpressed with ovule- and fiber development-associated PCGs, or linked to GWAS single-nucleotide polymorphisms. We mapped the genome-wide binding sites of two lncRNAs with chromatin isolation by RNA purification sequencing. We also validated the transcriptional regulatory role of lnc-Ga13g0352 via virus-induced gene suppression assay, indicating that this lncRNA might act as a dual-functional regulator that either activates or inhibits the transcription of target genes.

Analysis and review of trichomes in plants

BACKGROUND

Trichomes play a key role in the development of plants and exist in a wide variety of species.

RESULTS

In this paper, it was reviewed that the structure and morphology characteristics of trichomes, alongside the biological functions and classical regulatory mechanisms of trichome development in plants. The environment factors, hormones, transcription factor, non-coding RNA, etc., play important roles in regulating the initialization, branching, growth, and development of trichomes. In addition, it was further investigated the atypical regulation mechanism in a non-model plant, found that regulating the growth and development of tea (Camellia sinensis) trichome is mainly affected by hormones and the novel regulation factors.

CONCLUSIONS

This review further displayed the complex and differential regulatory networks in trichome initiation and development, provided a reference for basic and applied research on trichomes in plants.

Evolutionary dynamics of transposable elements and satellite DNAs in polyploid Spartina species

Repeated sequences and polyploidy play a central role in plant genome dynamics. Here, we analyze the evolutionary dynamics of repeats in tetraploid and hexaploid Spartina species that diverged during the last 10 million years within the Chloridoideae, one of the poorest investigated grass lineages. From high-throughput genome sequencing, we annotated Spartina repeats and determined what sequence types account for the genome size variation among species. We examined whether differential genome size evolution correlated with ploidy levels and phylogenetic relationships. We also examined the tempo of repeat sequence dynamics associated with allopatric speciation over the last 3-6 million years between hexaploid species that diverged on the American and European Atlantic coasts and tetraploid species from North and South America. The tetraploid S. spartinae, whose phylogenetic placement has been debated, exhibits a similar repeat content as hexaploid species, suggesting common ancestry. Genome expansion or contraction resulting from repeat dynamics seems to be explained mostly by the contrasting divergence times between species, rather than by genome changes triggered by ploidy level change per se. One 370 bp satellite may be exhibiting 'meiotic drive' and driving chromosome evolution in S. alterniflora. Our results provide crucial insights for investigating the genetic and epigenetic consequences of such differential repeat dynamics on the ecology and distribution of the meso- and neopolyploid Spartina species.

RES transformation for biosynthesis and detoxification

The reactive electrophilic species (RES), typically the molecules bearing α,β-unsaturated carbonyl group, are widespread in living organisms and notoriously known for their damaging effects. Many of the mycotoxins released from phytopathogenic fungi are RES and their contamination to cereals threatens food safety worldwide. However, due to their high reactivity, RES are also used by host organisms to synthesize specific metabolites. The evolutionary conserved glyoxalase (GLX) system scavenges the cytotoxic α-oxoaldehydes that bear RES groups, which cause host disorders and diseases. In cotton, a specialized enzyme derived from glyoxalase I (GLXI) through gene duplications and named as specialized GLXI (SPG), acts as a distinct type of aromatase in the gossypol pathway to transform the RES intermediates into the phenolic products. In this review, we briefly introduce the research progress in understanding the RES, especially the RES-type mycotoxins, the GLX system and SPG, and discuss their application potential in detoxification and synthetic biology.

Identification of a genome-specific repetitive element in the Gossypium D genome

The activity of genome-specific repetitive sequences is the main cause of genome variation between Gossypium A and D genomes. Through comparative analysis of the two genomes, we retrieved a repetitive element termed ICRd motif, which appears frequently in the diploid Gossypium raimondii (D₅) genome but rarely in the diploid Gossypium arboreum (A₂) genome. We further explored the existence of the ICRd motif in chromosomes of G. raimondii, G. arboreum, and two tetraploid (AADD) cotton species, Gossypium hirsutum and Gossypium barbadense, by fluorescence in situ hybridization (FISH), and observed that the ICRd motif exists in the D₅ and D-subgenomes but not in the A₂ and A-subgenomes. The ICRd motif comprises two components, a variable tandem repeat (TR) region and a conservative sequence (CS). The two constituents each have hundreds of repeats that evenly distribute across 13 chromosomes of the D₅genome. The ICRd motif (and its repeats) was revealed as the common conservative region harbored by ancient Long Terminal Repeat Retrotransposons. Identification and investigation of the ICRd motif promotes the study of A and D genome differences, facilitates research on Gossypium genome evolution, and provides assistance to subgenome identification and genome assembling.

Identification of a genome-specific repetitive element in the Gossypium D genome

The activity of genome-specific repetitive sequences is the main cause of genome variation between Gossypium A and D genomes. Through comparative analysis of the two genomes, we retrieved a repetitive element termed ICRd motif, which appears frequently in the diploid Gossypium raimondii (D₅) genome but rarely in the diploid Gossypium arboreum (A₂) genome. We further explored the existence of the ICRd motif in chromosomes of G. raimondii, G. arboreum, and two tetraploid (AADD) cotton species, Gossypium hirsutum and Gossypium barbadense, by fluorescence in situ hybridization (FISH), and observed that the ICRd motif exists in the D₅ and D-subgenomes but not in the A₂ and A-subgenomes. The ICRd motif comprises two components, a variable tandem repeat (TR) region and a conservative sequence (CS). The two constituents each have hundreds of repeats that evenly distribute across 13 chromosomes of the D₅genome. The ICRd motif (and its repeats) was revealed as the common conservative region harbored by ancient Long Terminal Repeat Retrotransposons. Identification and investigation of the ICRd motif promotes the study of A and D genome differences, facilitates research on Gossypium genome evolution, and provides assistance to subgenome identification and genome assembling.

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 Role of Sequence Duplication in Transcriptional Regulation and Genome Evolution

Sequence duplication is nowadays recognized as an important mechanism that underlies the evolution of eukaryote genomes, being indeed one of the most powerful strategies for the generation of adaptive diversity by modulating transcriptional activity. The evolutionary novelties simultaneously associated with sequence duplication and differential gene expression can be collectively referred to as duplication-mediated transcriptional regulation. In the last years, evidence has emerged supporting the idea that sequence duplication and functionalization represent important evolutionary strategies acting at the genome level, and both coding and non-coding sequences have been found to be targets of such events. Moreover, it has been proposed that deleterious effects of sequence duplication might be potentially silenced by endogenous cell machinery (i.e., RNA interference, epigenetic repressive marks, etc). Along these lines, our aim is to highlight the role of sequence duplication on transcriptional activity and the importance of both in genome evolution.

Analysis of the chromatin binding affinity of retrotransposases reveals novel roles in diploid and tetraploid cotton

LTR-retrotransposable elements are major components of diploid (Gossypium arboreum) and tetraploid (Gossypium hirsutum) cotton genomes that have undergone dramatic increases in copy number during the course of evolution. However, little is known about the biological functions of LTR-retrotransposable elements in cotton. Here, we show that a copia-like LTR-retrotransposable element has maintained considerable activity in both G. arboreum and G. hirsutum. We identified two functional domains of the retrotransposon and analyzed their expression levels in various cotton tissues, including leaves, ovules, and germinating seeds. ChIP-qPCR (chromatin immunoprecipitation followed by quantitative PCR), using a copia-specific antibody, established that copia-like proteins primarily bind to the first exons of several protein-coding genes in cotton cells. This finding suggests that retrotransposons play a novel, important role in regulating the transcriptional activities of protein-coding genes with various biological activities.

Comparative effect of allopolyploidy on transposable element composition and gene expression between Gossypium hirsutum and its two diploid progenitors

An allopolyploidization event formed allotetraploid Gossypium species from an A-genome diploid species and a D-genome diploid species. To explore the responses of transposable elements (TEs) to allopolyploidy, we assembled parallel TE datasets from G. hirsutum, G. arboreum and G. raimondii and analyzed the TE types and the effects of TEs on orthologous gene expression in the three Gossypium genomes. Gypsy was the most abundant TE type and most TEs were located ∼500 bp from genes in all three genomes. In G. hirsutum, 35.6% of genes harbored TE insertions, whereas insertions were more frequent in G. arboreum and G. raimondii. G. hirsutum had the highest proportion of uniquely matching 24-nt small interfering RNAs (siRNAs) that targeted TEs. TEs, particularly those targeted by 24-nt siRNAs, were associated with reduced gene expression, but the effect of TEs on orthologous gene expression varied substantially among species. Orthologous gene expression levels in G. hirsutum were intermediate between those of G. arboreum and G. raimondii, which did not experience TE expansion or reduction resulting from allopolyploidization. This study underscores the diversity of TEs co-opted by host genes and provides insights into the roles of TEs in regulating gene expression in Gossypium.

Genome-wide analysis of BES1 genes in Gossypium revealed their evolutionary conserved roles in brassinosteroid signaling

Brassinosteroids (BRs), which are essential phytohormones for plant growth and development, are important for cotton fiber development. Additionally, BES1 transcription factors are critical for BR signal transduction. However, cotton BES1 family genes have not been comprehensively characterized. In this study, we identified 11 BES1 genes in G. arboreum, 11 in G. raimondii, 16 in G. barbadense, and 22 in G. hirsutum. The BES1 sequences were significantly conserved in the Arabidopsis thaliana, rice, and upland cotton genomes. A total of 94 BES1 genes from 10 different plant species were divided into three clades according to the neighbor-joining and minimum-evolution methods. Moreover, the exon/intron patterns and motif distributions were highly conserved among the A. thaliana and cotton BES1 genes. The collinearity among the orthologs from the At and Dt subgenomes was estimated. Segmental duplications in the At and Dt subgenomes were primarily responsible for the expansion of the cotton BES1 gene family. Of the GhBES1 genes, GhBES1.4_At/Dt exhibited BL-induced expression and was predominantly expressed in fibers. Furthermore, Col-0/mGhBES1.4_At plants produced curled leaves with long and bent petioles. These transgenic plants also exhibited decreased hypocotyl sensitivity to brassinazole and constitutive BR induced/repressed gene expression patterns. The constitutive BR responses of the plants overexpressing mGhBES1.4_At were similar to those of the bes1-D mutant.

De novo assembly of a Chinese soybean genome

Soybean was domesticated in China and has become one of the most important oilseed crops. Due to bottlenecks in their introduction and dissemination, soybeans from different geographic areas exhibit extensive genetic diversity. Asia is the largest soybean market; therefore, a high-quality soybean reference genome from this area is critical for soybean research and breeding. Here, we report the de novo assembly and sequence analysis of a Chinese soybean genome for "Zhonghuang 13" by a combination of SMRT, Hi-C and optical mapping data. The assembled genome size is 1.025 Gb with a contig N50 of 3.46 Mb and a scaffold N50 of 51.87 Mb. Comparisons between this genome and the previously reported reference genome (cv. Williams 82) uncovered more than 250,000 structure variations. A total of 52,051 protein coding genes and 36,429 transposable elements were annotated for this genome, and a gene co-expression network including 39,967 genes was also established. This high quality Chinese soybean genome and its sequence analysis will provide valuable information for soybean improvement in the future.

Horizontal transfers of LTR retrotransposons in seven species of Rosales

Horizontal transposable element transfer (HTT) events have occurred among a large number of species and play important roles in the composition and evolution of eukaryotic genomes. HTTs are also regarded as effective forces in promoting genomic variation and biological innovation. In the present study, HTT events were identified and analyzed in seven sequenced species of Rosales using bioinformatics methods by comparing sequence conservation and K_a/K_s value of reverse transcriptase (RT) with 20 conserved genes, estimating the dating of HTTs, and analyzing the phylogenetic relationships. Seven HTT events involving long terminal repeat (LTR) retrotransposons, two HTTs between Morus notabilis and Ziziphus jujuba, and five between Malus domestica and Pyrus bretschneideri were identified. Further analysis revealed that these LTR retrotransposons had functional structures, and the copy insertion times were lower than the dating of HTTs, particularly in Mn.Zj.1 and Md.Pb.3. Altogether, the results demonstrate that LTR retrotransposons still have potential transposition activity in host genomes. These results indicate that HTT events are another strategy for exchanging genetic material among species and are important for the evolution of genomes.

Genome-Wide Survey and Comparative Analysis of Long Terminal Repeat (LTR) Retrotransposon Families in Four Gossypium Species

Long terminal repeat (LTR) retrotransposon is the most abundant DNA component and is largely responsible for plant genome size variation. Although it has been studied in plant species, very limited data is available for cotton, the most important fiber and texture crop. In this study, we performed a comprehensive analysis of LTR retrotransposon families across four cotton species. In tetraploid Gossypium species, LTR retrotransposon families from the progenitor D genome had more copies in D-subgenome, and families from the progenitor A genome had more copies in A-subgenome. Some LTR retrotransposon families that insert after polyploid formation may still distribute the majority of its copies in one of the subgenomes. The data also shows that families of 10~200 copies are abundant and they have a great influence on the Gossypium genome size; on the contrary, a small number of high copy LTR retrotransposon families have less contribution to the genome size. Kimura distance distribution indicates that high copy number family is not a recent outbreak, and there is no obvious relationship between family copy number and the period of evolution. Further analysis reveals that each LTR retrotransposon family may have their own distribution characteristics in cotton.

Cotton functional genomics reveals global insight into genome evolution and fiber development

Due to the economic value of natural textile fiber, cotton has attracted much research attention, which has led to the publication of two diploid genomes and two tetraploid genomes. These big data facilitate functional genomic study in cotton, and allow researchers to investigate cotton genome structure, gene expression, and protein function on the global scale using high-throughput methods. In this review, we summarized recent studies of cotton genomes. Population genomic analyses revealed the domestication history of cultivated upland cotton and the roles of transposable elements in cotton genome evolution. Alternative splicing of cotton transcriptomes was evaluated genome-widely. Several important gene families like MYC, NAC, Sus and GhPLDα1 were systematically identified and classified based on genetic structure and biological function. High-throughput proteomics also unraveled the key functional proteins correlated with fiber development. Functional genomic studies have provided unprecedented insights into global-scale methods for cotton research.

A D-genome-originated Ty1/Copia-type retrotransposon family expanded significantly in tetraploid cottons

Retrotransposons comprise of a major fraction of higher plant genomes, and their proliferation and elimination have profound effects on genome evolution and gene functions as well. Previously we found a D-genome-originated Ty1/Copia-type LTR (DOCL) retrotransposon in the chromosome A08 of upland cotton. To further characterize the DOCL retrotransposon family, a total of 342 DOCL retrotransposons were identified in the sequenced cotton genomes, including 73, 157, and 112 from Gossypium raimondii, G. hirsutum, and G. barbadense, respectively. According to phylogenetic analysis, the DOCL family was divided into nine groups (G1-G9), among which five groups (G1-G4 and G9, including 292 members) were proliferated after the formation of tetraploid cottons. It was found that the majority of DOCL retrotransposons (especially those in G2, G3 and G9) inserted in non-allelic loci in G. hirsutum and G. barbadense, suggesting that their proliferations were relatively independent in different tetraploid cottons. Furthermore, DOCL retrotransposons inserted in coding regions largely eliminated expression of the targeted genes in G. hirsutum or G. barbadense. Our data suggested that recent proliferation of retrotransposon families like DOCL was one of important evolutionary forces driving diversification and evolution of tetraploid cottons.

Global analysis of DNA methylation in young (J1) and senescent (J2) Gossypium hirsutum L. cotyledons by MeDIP-Seq

DNA methylation is an important epigenetic modification regulating gene expression, genomic imprinting, transposon silencing and chromatin structure in plants and plays an important role in leaf senescence. However, the DNA methylation pattern during Gossypium hirsutum L. cotyledon senescence is poorly understood. In this study, global DNA methylation patterns were compared between two cotyledon development stages, young (J1) and senescence (J2), using methylated DNA immunoprecipitation (MeDIP-Seq). Methylated cytosine occurred mostly in repeat elements, especially LTR/Gypsy in both J1 and J2. When comparing J1 against J2, there were 1222 down-methylated genes and 623 up-methylated genes. Methylated genes were significantly enriched in carbohydrate metabolism, biosynthesis of other secondary metabolites and amino acid metabolism pathways. The global DNA methylation level decreased from J1 to J2, especially in gene promoters, transcriptional termination regions and regions around CpG islands. We further investigated the expression patterns of 9 DNA methyltransferase-associated genes and 2 DNA demethyltransferase-associated genes from young to senescent cotyledons, which were down-regulated during cotyledon development. In this paper, we first reported that senescent cotton cotyledons exhibited lower DNA methylation levels, primarily due to decreased DNA methyltransferase activity and which also play important role in regulating secondary metabolite process.

Genome-wide identification and characterization of phospholipase C gene family in cotton (Gossypium spp.)

Phospholipase C (PLC) are important regulatory enzymes involved in several lipid and Ca²⁺-dependent signaling pathways. Previous studies have elucidated the versatile roles of PLC genes in growth, development and stress responses of many plants, however, the systematic analyses of PLC genes in the important fiber-producing plant, cotton, are still deficient. In this study, through genome-wide survey, we identified twelve phosphatidylinositol-specific PLC (PI-PLC) and nine non-specific PLC (NPC) genes in the allotetraploid upland cotton Gossypium hirsutum and nine PI-PLC and six NPC genes in two diploid cotton G. arboretum and G.raimondii, respectively. The PI-PLC and NPC genes of G. hirsutum showed close phylogenetic relationship with their homologous genes in the diploid cottons and Arabidopsis. Segmental and tandem duplication contributed greatly to the formation of the gene family. Expression profiling indicated that few of the PLC genes are constitutely expressed, whereas most of the PLC genes are preferentially expressed in specific tissues and abiotic stress conditions. Promoter analyses further implied that the expression of these PLC genes might be regulated by MYB transcription factors and different phytohormones. These results not only suggest an important role of phospholipase C members in cotton plant development and abiotic stress response but also provide good candidate targets for future molecular breeding of superior cotton cultivars.

Epigenetic regulation and epigenomic landscape in rice

Epigenetic regulation has been implicated in the control of complex agronomic traits in rice (Oryza sativa), a staple food crop and model monocot plant. Recent advances in high-throughput sequencing and the moderately complex genome of rice have made it possible to study epigenetic regulation in rice on a genome-wide scale. This review discusses recent advances in our understanding of epigenetic regulation in rice, with an emphasis on the roles of key epigenetic regulators, the epigenomic landscape, epigenetic variation, transposon repression, and plant development.

Recent advances in the research for the homolog of breast cancer associated gene AtROW1 in higher plants

BARD1 (BRCA1 associated RING domain protein 1), as an important animal tumor suppressor gene associated with many kinds of cancers, has been intensively studied for decades. Surprisingly, homolog of BARD1 was found in plants and it was renamed AtROW1 (repressor of Wuschel-1) according to its extremely important function with regard to plant stem cell homeostasis. Although great advances have been made in human BARD1, the function of this animal tumor-suppressor like gene in plant is not well studied and need to be further elucidated. Here, we review and summarize past and present work regarding this protein. Apart from its previously proposed role in DNA repair, recently it is found essential for shoot and root stem cell development and differentiation in plants. The study of AtROW1 in plant may provide an ideal model for further elucidating the functional mechanism of BARD1 in mammals.