Heterochromatin evolution in Arachis investigated through genome-wide analysis of repetitive DNA

The Stilbene Synthase Family in Arachis: A Genome-Wide Study and Functional Characterization in Response to Stress

Peanut (Arachis hypogaea) and its wild relatives are among the few species that naturally synthesize resveratrol, a well-known stilbenoid phytoalexin that plays a crucial role in plant defense against biotic and abiotic stresses. Resveratrol has received considerable attention due to its health benefits, such as preventing and treating various human diseases and disorders. Chalcone (CHS) and Stilbene (STS) Synthases are plant-specific type III Polyketide Synthases (PKSs) that share the same substrates and are key branch enzymes in the biosynthesis of flavonoids and stilbenoids, respectively. Although resveratrol accumulation in response to external stimulus has been described in peanut, there are no comprehensive studies of the CHS and STS gene families in the genus Arachis. In the present study, we identified and characterized 6 CHS and 46 STS genes in the tetraploid peanut and an average of 4 CHS and 22 STS genes in three diploid wild species (Arachis duranensis, Arachis ipaënsis and Arachis stenosperma). The CHS and STS gene and protein structures, chromosomal distributions, phylogenetic relationships, conserved amino acid domains, and cis-acting elements in the promoter regions were described for all Arachis species studied. Based on gene expression patterns of wild A. stenosperma STS genes in response to different biotic and abiotic stresses, we selected the candidate AsSTS4 gene, which is strongly induced by ultraviolet (UV) light exposure, for further functional investigation. The AsSTS4 overexpression in peanut hairy roots significantly reduced (47%) root-knot nematode infection, confirming that stilbene synthesis activation in transgenic plants can increase resistance to pathogens. These findings contribute to understanding the role of resveratrol in stress responses in Arachis species and provide the basis for genetic engineering for improved production of valuable secondary metabolites in plants.

The nature and organization of satellite DNAs in Petunia hybrida, related, and ancestral genomes

Introduction

The garden petunia, Petunia hybrida (Solanaceae) is a fertile, diploid, annual hybrid species (2n=14) originating from P. axillaris and P. inflata 200 years ago. To understand the recent evolution of the P. hybrida genome, we examined tandemly repeated or satellite sequences using bioinformatic and molecular cytogenetic analysis.

Methods

Raw reads from available genomic assemblies and survey sequences of P. axillaris N (PaxiN), P. inflata S6, (PinfS6), P. hybrida (PhybR27) and the here sequenced P. parodii S7 (PparS7) were used for graph and k-mer based cluster analysis of TAREAN and RepeatExplorer. Analysis of repeat specific monomer lengths and sequence heterogeneity of the major tandem repeat families with more than 0.01% genome proportion were complemented by fluorescent in situ hybridization (FISH) using consensus sequences as probes to chromosomes of all four species.

Results

Seven repeat families, PSAT1, PSAT3, PSAT4, PSAT5 PSAT6, PSAT7 and PSAT8, shared high consensus sequence similarity and organisation between the four genomes. Additionally, many degenerate copies were present. FISH in P. hybrida and in the three wild petunias confirmed the bioinformatics data and gave corresponding signals on all or some chromosomes. PSAT1 is located at the ends of all chromosomes except the 45S rDNA bearing short arms of chromosomes II and III, and we classify it as a telomere associated sequence (TAS). It is the most abundant satellite repeat with over 300,000 copies, 0.2% of the genomes. PSAT3 and the variant PSAT7 are located adjacent to the centromere or mid-arm of one to three chromosome pairs. PSAT5 has a strong signal at the end of the short arm of chromosome III in P. axillaris and P.inflata, while in P. hybrida additional interstitial sites were present. PSAT6 is located at the centromeres of chromosomes II and III. PSAT4 and PSAT8 were found with only short arrays.

Discussion

These results demonstrate that (i) repeat families occupy distinct niches within chromosomes, (ii) they differ in the copy number, cluster organization and homogenization events, and that (iii) the recent genome hybridization in breeding P. hybrida preserved the chromosomal position of repeats but affected the copy number of repetitive DNA.

Genomic relationships of the polyploid rhizoma peanut (Arachis glabrata Benth.) inferred by genomic in situ hybridization (GISH)

The rhizoma peanut (Arachis glabrata Benth., section Rhizomatosae) is a tetraploid perennial legume. Although several A. glabrata cultivars have been developed as forage and ornamental turf, the origin and genomic constitution of this species are still unknown. In this study, we evaluated the affinity between the genomes of A. glabrata and the probable diploid donors of the sections Rhizomatosae, Arachis, Erectoides and Procumbentes by genomic in situ hybridization (GISH). Single GISH analyses detected that species of the sections Erectoides (E2 subgenome) and Procumbentes (E3 subgenome) were the diploid species with the highest degree of genomic affinity with A. glabrata. Based on single GISH experiments and DNA sequence similarity, three species -A. duranensis, A. paraguariensis subsp. capibarensis, and A. rigonii-, which showed the most uniform and brightest hybridization patterns and lowest genetic distance, were selected as probes for double GISH experiments. Double GISH experiments showed that A. glabrata is constituted by four identical or very similar chromosome complements. In these assays, A. paraguariensis subsp. capibarensis showed the highest brightness onto A. glabrata chromosomes. Thus, our results support the autopolyploid origin of A. glabrata and show that the species with E2 subgenome are the most probable ancestors of this polyploid legume forage.

Comparative repeatome analysis reveals new evidence on genome evolution in wild diploid Arachis (Fabaceae) species

Opposing changes in the abundance of satellite DNA and long terminal repeat (LTR) retroelements are the main contributors to the variation in genome size and heterochromatin amount in Arachis diploids. The South American genus Arachis (Fabaceae) comprises 83 species organized in nine taxonomic sections. Among them, section Arachis is characterized by species with a wide genome and karyotype diversity. Such diversity is determined mainly by the amount and composition of repetitive DNA. Here we performed computational analysis on low coverage genome sequencing to infer the dynamics of changes in major repeat families that led to the differentiation of genomes in diploid species (x = 10) of genus Arachis, focusing on section Arachis. Estimated repeat content ranged from 62.50 to 71.68% of the genomes. Species with different genome composition tended to have different landscapes of repeated sequences. Athila family retrotransposons were the most abundant and variable lineage among Arachis repeatomes, with peaks of transpositional activity inferred at different times in the evolution of the species. Satellite DNAs (satDNAs) were less abundant, but differentially represented among species. High rates of evolution of an AT-rich superfamily of satDNAs led to the differential accumulation of heterochromatin in Arachis genomes. The relationship between genome size variation and the repetitive content is complex. However, largest genomes presented a higher accumulation of LTR elements and lower contents of satDNAs. In contrast, species with lowest genome sizes tended to accumulate satDNAs in detriment of LTR elements. Phylogenetic analysis based on repetitive DNA supported the genome arrangement of section Arachis. Altogether, our results provide the most comprehensive picture on the repeatome dynamics that led to the genome differentiation of Arachis species.

Chromosome-level genome of Camellia lanceoleosa provides a valuable resource for understanding genome evolution and self-incompatibility

The section Oleifera (Theaceae) has attracted attention for the high levels of unsaturated fatty acids found in its seeds. Here, we report the chromosome-scale genome of the sect. Oleifera using diploid wild Camellia lanceoleosa with a final size of 3.00 Gb and an N50 scaffold size of 186.43 Mb. Repetitive sequences accounted for 80.63% and were distributed unevenly across the genome. Camellia lanceoleosa underwent a whole-genome duplication event approximately 65 million years ago (65 Mya), prior to the divergence of C. lanceoleosa and Camellia sinensis (approx. 6-7 Mya). Syntenic comparisons of these two species elucidated the genomic rearrangement, appearing to be driven in part by the activity of transposable elements. The expanded and positively selected genes in C. lanceoleosa were significantly enriched in oil biosynthesis, and the expansion of homomeric acetyl-coenzyme A carboxylase (ACCase) genes and the seed-biased expression of genes encoding heteromeric ACCase, diacylglycerol acyltransferase, glyceraldehyde-3-phosphate dehydrogenase and stearoyl-ACP desaturase could be of primary importance for the high oil and oleic acid content found in C. lanceoleosa. Theanine and catechins were present in the leaves of C. lanceoleosa. However, caffeine can not be dectected in the leaves but was abundant in the seeds and roots. The functional and transcriptional divergence of genes encoding SAM-dependent N-methyltransferases may be associated with caffeine accumulation and distribution. Gene expression profiles, structural composition and chromosomal location suggest that the late-acting self-incompatibility of C. lanceoleosa is likely to have favoured a novel mechanism co-occurring with gametophytic self-incompatibility. This study provides valuable resources for quantitative and qualitative improvements and genome assembly of polyploid plants in sect. Oleifera.

Satellite DNA probes of Alstroemeria longistaminea (Alstroemeriaceae) paint the heterochromatin and the B chromosome, reveal a G-like banding pattern, and point to a strong structural karyotype conservation

Alstroemeria species present a well-conserved and asymmetric karyotype. The genus is divided into a Chilean clade, rich in heterochromatin, and a Brazilian clade, poor in heterochromatin. We investigated the distribution of the main repetitive sequences in the chromosomes of the Brazilian species A. longistaminea (2n = 16 + 0-6B) aiming to evaluate the role played by these sequences on the structural organization of the karyotype. In situ hybridization of the three most abundant retrotransposons, corresponding to ~ 45% of the genome, was uniformly distributed. Three satellite DNA sequences, representing near half of the whole satellite fraction (1.93% of the genome), were mainly concentrated on the heterochromatin and one of them painted the whole B chromosome. Noteworthy, some satellites were located on euchromatin, either dispersed or concentrated in clusters along the chromosomes, revealing a G-band-like pattern. The two satellites that presented more C-band- and G-band-like labeling were also hybridized in situ in two other Alstroemeria species. They revealed astonishing similar patterns of distribution, indicating an unusually structural karyotype conservation among Brazilian species.

An 82 bp tandem repeat family typical of 3' non-coding end of Gypsy/TAT LTR retrotransposons is conserved in Coffea spp. pericentromeres

Coffea spp. chromosomes are very small and accumulate a variety of repetitive DNA families around the centromeres. However, the proximal regions of Coffea chromosomes remain poorly understood, especially regarding the nature and organisation of the sequences. Taking advantage of the genome sequences of C. arabica (2n = 44), C. canephora, and C. eugenioides (C. arabica progenitors with 2n = 22) and good coverage genome sequencing of dozens of other wild Coffea spp., repetitive DNA sequences were identified, and the genomes were compared to decipher particularities of pericentromeric structures. The searches revealed a short tandem repeat (82 bp length) typical of Gypsy/TAT LTR retrotransposons, named Coffea_sat11. This repeat organises clusters with fragments of other transposable elements, comprising regions of non-coding RNA production. Cytogenomic analyses showed that Coffea_sat11 extends from the pericentromeres towards the middle of the chromosomal arms. This arrangement was observed in the allotetraploid C. arabica chromosomes, as well as in its progenitors. This study improves our understanding of the role of the Gypsy/TAT LTR retrotransposon lineage in the organisation of Coffea pericentromeres, as well as the conservation of Coffea_sat11 within the genus. The relationships between fragments of other transposable elements and the functional aspects of these sequences on the pericentromere chromatin were also evaluated. Highlights: A scattered short tandem repeat, typical of Gypsy/TAT LTR retrotransposons, associated with several fragments of other transposable elements, accumulates in the pericentromeres of Coffea chromosomes. This arrangement is preserved in all clades of the genus and appears to have a strong regulatory role in the organisation of chromatin around centromeres.

How diverse is heterochromatin in the Caesalpinia group? Cytogenomic characterization of Erythrostemon hughesii Gagnon & G.P. Lewis (Leguminosae: Caesalpinioideae)

Cytogenomic characterization of Erythrostemon hughesii reveals a heterogeneity of repeats in its subtelomeric heterochromatin. Comparative analyses with other Caesalpinia group species reveal a significant reduction in the abundance of Ty3-gypsy/Chromovirus Tekay retrotransposons during its evolution. In numerically stable karyotypes, repetitive DNA variability is one of the main causes of genome and chromosome variation and evolution. Species from the Caesalpinia group (Leguminosae) are karyotypically characterized by 2n = 24, with small chromosomes and highly variable CMA⁺ heterochromatin banding patterns that correlate with environmental variables. Erythrostemon hughesii differs from other species of the group examined to date for having subtelomeric CMA⁺ bands; this contrasts with most species in the group which have proximal bands. Here we analyse the repeatome of E. hughesii using genome skimming and chromosomal mapping approaches to characterize the identity of the most abundant repetitive elements and their physical location. The repetitive fraction of E. hughesii comprises 28.73% of the genome. The most abundant elements were retrotransposons (RT) with long terminal repeats (LTR-RT; 9.76%) and satellite DNAs (7.83%). Within the LTR-RTs, the most abundant lineages were: Ty1/copia-Ale (1%), Ty3/gypsy CRM (0.88%) and Ty3/gypsy Athila (0.75%). Using fluorescent in situ hybridization four satellite DNAs and several LTR-RT elements were shown to be present in most subtelomeric CMA⁺ bands. These results highlight how the repeatome in E. hughesii, a species from Oaxaca state in Mexico, is clearly distinct from Northeast Brazilian species of the Caesalpinia group, mainly due to its high diversity of repeats in its subtelomeric heterochromatic bands and low amount of LTR-RT Ty3/gypsy-Tekay elements. Comparative sequence analysis of Tekay elements from different species is congruent with a clade-specific origin of this LTR-RT after the divergence of the Caesalpinia group. We hypothesize that repeat-rich heterochromatin may play a role in leading to faster genomic divergence between individuals, increasing speciation and diversification.

Chromosome diversity in species of the genus Arachis, revealed by FISH and CMA/DAPI banding, and inferences about their karyotype differentiation

The species of the genus Arachis (Leguminosae) are ordered into nine sections. The assignment of genome types in this genus has been based on cross-compatibility analysis and molecular cytogenetic studies. The latter has also allowed karyotypically establishing well-defined genomes and reassigning the genome of several species. However, most of these studies have been focused mainly on the sections Arachis and Rhizomatosae. To increase the knowledge about the chromosome diversity of the whole genus, here we performed a detailed karyotype characterization of representative species of most of the sections and genomes of Arachis. This characterization included chromosome morphology, CMA/DAPI chromosome banding, and chromosome marker localization (rDNAloci and one satDNA sequence) by fluorescent in situ hybridization (FISH). Based on the data obtained and other previously published data, we established the karyotype similarities by cluster analysis and defined eleven karyotype groups. The grouping was partly coincident with the traditional genome assignment, except for some groups and some individual species. Karyotype similarities among some genomes were also found. The main characteristics of each karyotype group of Arachis were summarized. Together, our results provide information that may be beneficial for future cytogenetic and evolutionary studies, and also contribute to the identification of interspecific hybrids.

Genome-wide analysis of microsatellite and sex-linked marker identification in Gleditsia sinensis

BACKGROUND

Gleditsia sinensis Lam. (Leguminosae), a dioecious perennial arbor, demonstrates important medicinal properties and economic value. These properties can be harnessed depending on the sex of the plant. However, the sex of the plants is difficult to identify accurately through morphological methods before the flowering.

RESULTS

We used bulked segregant analysis to screen sex-specific simple sequence repeat (SSR) markers in G. sinensis. Five male and five female plants were pooled to form the male and female bulks, respectively, and subjected to whole-genome sequencing. After high-throughput sequencing, 5,350,359 sequences were obtained, in which 2,065,210 SSRs were searched. Among them, the number of duplicated SSRs was the highest. The male plants could reach 857,874, which accounted for 60.86% of the total number of male plants. The female plants could reach 1,447,603, which accounted for 56.25% of the total model of the female plants. Among all the nucleotide repeat types, the A/T-rich motif was the most abundant. A total of 309,516 female strain-specific SSRs were selected by clustering. After designing the primers, the male and female gene pools were amplified, and five pairs of primers (i.e., 27, 34, 36, 39, and 41) were found to amplify the differential bands in the male and female gene pools. Using the five pairs of primers, we performed PCR verification on 10 individuals of known sex, which constructed the gene pool. The female plants amplified a single fragment of lengths (i.e., 186, 305, 266, 203, and 260 bp) and no male plant strip, thereby completing the identification of the male and female sexes of the G. sinensis.

CONCLUSIONS

This study provides accurate sex identification strategies between female and male plants, thus improving the utilization rate of G. sinensis resources.

Chromosome-level genome assembly, annotation and evolutionary analysis of the ornamental plant Asparagus setaceus

Asparagus setaceus is a popular ornamental plant cultivated in tropical and subtropical regions globally. Here, we constructed a chromosome-scale reference genome of A. setaceus to facilitate the investigation of its genome characteristics and evolution. Using a combination of Nanopore long reads, Illumina short reads, 10× Genomics linked reads, and Hi-C data, we generated a high-quality genome assembly of A. setaceus covering 710.15 Mb, accounting for 98.63% of the estimated genome size. A total of 96.85% of the sequences were anchored to ten superscaffolds corresponding to the ten chromosomes. The genome of A. setaceus was predicted to contain 28,410 genes, 25,649 (90.28%) of which were functionally annotated. A total of 65.59% of the genome was occupied by repetitive sequences, among which long terminal repeats were predominant (42.51% of the whole genome). Evolutionary analysis revealed an estimated divergence time of A. setaceus from its close relative A. officinalis of ~9.66 million years ago, and A. setaceus underwent two rounds of whole-genome duplication. In addition, 762 specific gene families, 96 positively selected genes, and 76 resistance (R) genes were detected and functionally predicted in A. setaceus. These findings provide new knowledge about the characteristics and evolution of the A. setaceus genome, and will facilitate comparative genetic and genomic research on the genus Asparagus.

Evolutionary convergence or homology? Comparative cytogenomics of Caesalpinia group species (Leguminosae) reveals diversification in the pericentromeric heterochromatic composition

We demonstrated by cytogenomic analysis that the proximal heterochromatin of the Northeast Brazilian species of Caesalpinia group is enriched with phylogenetically conserved Ty3/Gypsy-Tekay RT, but diverge in the presence of Ty3/Gypsy-Athila RT and satDNA. The Caesalpinia Group includes 225 species and 27 monophyletic genera of which four occur in Northeastern Brazil: Erythrostemon (1 sp.), Cenostigma (7 spp.), Libidibia (1 sp.), and Paubrasilia (1 sp.). The last three genera are placed in different clades in the Caesalpinia Group phylogeny, and yet they are characterized by having a numerically stable karyotype 2n = 24 (16 M+8A) and GC-rich heterochromatic bands (chromomycin A₃ positive/CMA⁺ bands) in the proximal chromosome regions. To characterize the composition of their heterochromatin and test for the homology of these chromosomal regions, genomic DNA was extracted from Cenostigma microphyllum, Libidibia ferrea, and Paubrasilia echinata, and sequenced at low coverage using the Illumina platform. The genomic repetitive fractions were characterized using a Galaxy/RepeatExplorer-Elixir platform. The most abundant elements of each genome were chromosomally located by fluorescent in situ hybridization (FISH) and compared to the CMA⁺ heterochromatin distribution. The repetitive fraction of the genomes of C. microphyllum, L. ferrea, and P. echinata were estimated to be 41.70%, 38.44%, and 72.51%, respectively. Ty3/Gypsy retrotransposons (RT), specifically the Tekay lineage, were the most abundant repeats in each of the three genomes. FISH mapping revealed species-specific patterns for the Tekay elements in the proximal regions of the chromosomes, co-localized with CMA⁺ bands. Other species-specific patterns were observed, e.g., for the Ty3/Gypsy RT Athila elements which were found in all the proximal heterochromatin of L. ferrea or restricted to the acrocentric chromosomes of C. microphyllum. This Athila labeling co-localized with satellite DNAs (satDNAs). Although the Caesalpinia Group diverged around 55 Mya, our results suggest an ancestral colonization of Tekay RT in the proximal heterochromatin. Thus, the present-day composition of the pericentromeric heterochromatin in these Northeast Brazilian species is a combination of the maintenance of an ancestral Tekay distribution with a species-specific accumulation of other repeats.