The landscape of transposable elements and satellite DNAs in the genome of a dioecious plant spinach (Spinacia oleracea L.)

The landscape of transposable element distribution in the genome of Neotropical fish Apareiodon sp. (Characiformes: Parodontidae)

Transposable elements (TEs) are widely present in eukaryotic genomes, where they can contribute to genome size and functional modifications. As new genomes are sequenced and annotated, more studies can be conducted regarding TE content, distribution, and genome evolution. TEs are extensively diversified in fish genomes resulting in an important role in genome and chromosome evolution. However, curated TE libraries are still scarce in non-model organisms, making it difficult to evaluate TE's impact on genomic modifications thoroughly. Here, we aimed to obtain a curated TE library from the neotropical fish Apareiodon sp. genome. The prospection and curation of the TE library resulted in 244 families from 18 superfamilies of DNA transposons and retrotransposons, which comprise about 10% of the genome, with most insertions fitting in one or a few families. A greater diversity of retrotransposon families is present, especially for Ty3 superfamily. Despite the greater number of retrotransposon families, DNA transposons are the most abundant in the genome, with 37% of all TE insertions belonging to the Tc1-Mariner superfamily. Complete TE copies were observed for almost all superfamilies, with most of the sequences on the Tc1-Mariner group. DNA transposons and SINEs presented older insertions in the genome, followed by LINEs and LTR retrotransposons. TE genome density is highest in the cs25 scaffold, and enriched for Helitron elements. With these data, allied to previous studies on chromosome evolution, we suggest that cs25 bears the W chromosome specific region of the Apareiodon sp. genome, with the presence of significant amount of Helitron insertions.

A dominant role of transcriptional regulation during the evolution of C4 photosynthesis in Flaveria species

C4 photosynthesis exemplifies convergent evolution of complex traits. Herein, we construct chromosome-scale genome assemblies and perform multi-omics analysis for five Flaveria species, which represent evolutionary stages from C3 to C4 photosynthesis. Chromosome-scale genome sequence analyses reveal a gradual increase in genome size during the evolution of C4 photosynthesis attributed to the expansion of transposable elements. Systematic annotation of genes encoding C4 enzymes and transporters identify additional copies of three C4 enzyme genes through retrotranspositions in C4 species. C4 genes exhibit elevated mRNA and protein abundances, reduced protein-to-RNA ratios, and comparable translation efficiencies in C4 species, highlighting a critical role of transcriptional regulation in C4 evolution. Furthermore, we observe an increased abundance of ethylene response factor (ERF) transcription factors and cognate cis-regulatory elements associated with C4 genes regulation. Altogether, our study provides valuable genomic resources for the Flaveria genus and sheds lights on evolutionary and regulatory mechanisms underlying C4 photosynthesis.

Chromosome distribution of four LTR retrotransposons and 18 S rDNA in coffea eugenioides

Repetitive sequences are recognized for their roles in plant genome organization and function. Mobile elements are notable repeatome sequences due to their intrinsic mutagenic potential, which is related to the origin of adaptive novelties. Understanding the genomic organization and dynamics of the repeatome is fundamental to enlighten their role in plant genome evolution. We aimed to map and assemble the first karyogram for a Coffea species with a closer look at mobile elements. Four LTR-retrotransposons (LTR-RTs) and the 18S rDNA of Coffea eugenioides, a diploid progenitor of the allotetraploid Coffea arabica, were unprecedently mapped in prometaphase/metaphase chromosomes and interphase nuclei. The LTR-RTs included three Ty1/Copia (Bianca, TAR and Tork) and one Ty3/Gypsy (Athila) identified based on homology searches. The four LTR-RTs were mainly distributed in a clustered pattern throughout different portions of the 2n = 22 chromosomes. Athila showed the most intense fluorescence signals and co-located with the secondary constriction of chromosome 3. In addition, the 18S rDNA was mapped in the distal portions of the short arms of chromosome pairs 3 and 5. The obstacles related to obtaining high-quality chromosomes from Coffea species have long been hampering the cytogenomics, which associates in silico analysis with the in situ mapping. Thus, we hope that the results presented here enlighten not only the composition, but also the distribution of mobile elements in the C. eugenioides genome, providing background for further cytogenomic investigations regarding Coffea repeatome.

Combined metabolomic and transcriptomic analysis reveals the key genes for triterpenoid biosynthesis in Cyclocarya paliurus

BACKGROUND

Cyclocarya paliurus is a high-value tree, and it contains a variety of bioactive secondary metabolites which have broad application prospects in medicine, food and health care. Triterpenoids can improve the bioactive function of C. paliurus health tea and also improve the efficacy of health care tea.

RESULTS

The results of this study showed that there were 69 kinds were terpenoids, and triterpenoids accounted for more than 80%. We excavated 5 kinds of triterpenoid metabolites with high content and significant difference dynamics, namely, corosolic acid, asiatic acid, maslinic acid, ursolic acid and oleanolic acid. The co-expression analysis identified CYP71D8 and CYP716A15 co-expressed with β-AS may generate oleanane type triterpenoids by modifying β-amyrin, while CYP71AN24 and CYP98A2 co-expressed with LUS may play a key role in lupine type triterpenoids biosynthesis. MYB,Whirly,WRKY and bHLH families, which showed strong correlation with function genes, may play an important role in the regulation of P450 and OSC expression. A total of 20 modules were identified by WGCNA analysis, and CYP71AU50 and CYP716A15 in tan and orange modules may play a major role in the synthesis of oleanolic acid, ursolic acid and asiatic acid, while CYP82D47 in lightcyan 1 module may be the hub gene for the biosynthesis of corosolic acid and maslinic acid.

CONCLUSIONS

Our findings mined candidate genes closely related to triterpenoid synthesis in C. paliurus. The results of this paper can provide scientific reference for breeding high-content triterpenoid varieties of C. paliurus.

Comparative analysis of transposable elements dynamics in fish with different sex chromosome systems

Transposable elements (TEs) are widespread genomic components with substantial roles in genome evolution and sex chromosome differentiation. In this study, we compared the TE composition of three closely related fish with different sex chromosome systems: Megaleporinus elongatus (Z1Z1Z2Z2/Z1W1Z2W2), Megaleporinus macrocephalus (ZZ/ZW) (both with highly differentiated W sex chromosomes), and Leporinus friderici (without heteromorphic sex chromosomes). We created custom TE libraries for each species using clustering methods and manual annotation and prediction, and we predicted TE temporal dynamics through divergence-based analysis. The TE abundance ranged from 16% to 21% in the three mobilomes, with L. friderici having the lowest overall. Despite the recent amplification of TEs in all three species, we observed differing expansion activities, particularly between the two genera. Both Megaleporinus recently experienced high retrotransposon activity, with a reduction in DNA TEs, which could have implications in sex chromosome composition. In contrast, L. friderici showed the opposite pattern. Therefore, despite having similar TE compositions, Megaleporinus and Leporinus exhibit distinct TE histories that likely evolved after their separation, highlighting a rapid TE expansion over short evolutionary periods.

Evolution of the spinach sex-linked region within a rarely recombining pericentromeric region

Sex chromosomes have evolved independently in many different plant lineages. Here, we describe reference genomes for spinach (Spinacia oleracea) X and Y haplotypes by sequencing homozygous XX females and YY males. The long arm of 185-Mb chromosome 4 carries a 13-Mb X-linked region (XLR) and 24.1-Mb Y-linked region (YLR), of which 10 Mb is Y specific. We describe evidence that this reflects insertions of autosomal sequences creating a "Y duplication region" or "YDR" whose presence probably directly reduces genetic recombination in the immediately flanking regions, although both the X and Y sex-linked regions are within a large pericentromeric region of chromosome 4 that recombines rarely in meiosis of both sexes. Sequence divergence estimates using synonymous sites indicate that YDR genes started diverging from their likely autosomal progenitors about 3 MYA, around the time when the flanking YLR stopped recombining with the XLR. These flanking regions have a higher density of repetitive sequences in the YY than the XX assembly and include slightly more pseudogenes compared with the XLR, and the YLR has lost about 11% of the ancestral genes, suggesting some degeneration. Insertion of a male-determining factor would have caused Y linkage across the entire pericentromeric region, creating physically small, highly recombining, terminal pseudoautosomal regions. These findings provide a broader understanding of the origin of sex chromosomes in spinach.

Effect of Growth Stages and Lactic Acid Fermentation on Anti-Nutrients and Nutritional Attributes of Spinach (Spinacia oleracea)

Spinach (Spinacia oleracea) is a winter-season green, leafy vegetable grown all over the world, belonging to the family Amaranthus, sub-family Chenopodiaceae. Spinach is a low-caloric food and an enormous source of micronutrients, e.g., calcium, folates, zinc, retinol, iron, ascorbic acid and magnesium. Contrarily, it also contains a variety of anti-nutritional factors, e.g., alkaloids, phytates, saponins, oxalates, tannins and many other natural toxicants which may hinder nutrient-absorption. This study was aimed at investigating the effect of fermentation on improving the nutrient-delivering potential of spinach and mitigating its burden of antinutrients and toxicants at three growth stages: the 1st growth stage as baby leaves, the 2nd growth stage at the coarse stage, and the 3rd growth stage at maturation. The results revealed the significant (p < 0.05) effect of fermentation on increasing the protein and fiber content of spinach powder from 2.53 to 3.53% and 19.33 to 22.03%, respectively, and on reducing total carbohydrate content from 52.92 to 40.52%; the effect was consistent in all three growth stages. A significant decline in alkaloids (6.45 to 2.20 mg/100 g), oxalates (0.07 mg/100 g to 0.02 mg/100 g), phytates (1.97 to 0.43 mg/100 g) and glucosinolates (201 to 10.50 µmol/g) was observed as a result of fermentation using Lactiplantibacillus plantarum. Fermentation had no impact on total phenolic content and the antioxidant potential of spinach, as evaluated using 2,2-diphenyl-1-picrylhydrazyl (DPPH) and ferric-reducing antioxidant power (FRAP) assays. This study proposes fermentation as a safer bioprocess for improving the nutrient-delivering potential of spinach, and suggests processed powders made from spinach as a cost-effective complement to existing plant proteins.

Oligo-FISH of Populus simonii Pachytene Chromosomes Improves Karyotyping and Genome Assembly

Poplar was one of the first woody species whose individual chromosomes could be identified using chromosome specific painting probes. Nevertheless, high-resolution karyotype construction remains a challenge. Here, we developed a karyotype based on the meiotic pachytene chromosome of Populus simonii which is a Chinese native species with many excellent traits. This karyotype was anchored by oligonucleotide (oligo)-based chromosome specific painting probes, a centromere-specific repeat (Ps34), ribosomal DNA, and telomeric DNA. We updated the known karyotype formula for P. simonii to 2n = 2x = 38 = 26m + 8st + 4t and the karyotype was 2C. The fluorescence in situ hybridization (FISH) results revealed some errors in the current P. simonii genome assembly. The 45S rDNA loci were located at the end of the short arm of chromosomes 8 and 14 by FISH. However, they were assembled on pseudochromosomes 8 and 15. In addition, the Ps34 loci were distributed in every centromere of the P. simonii chromosome in the FISH results, but they were only found to be present in pseudochromosomes 1, 3, 6, 10, 16, 17, 18, and 19. Our results reveal that pachytene chromosomes oligo-FISH is a powerful tool for constructing high-resolution karyotypes and improving the quality of genome assembly.

Characterization of repetitive DNA on the genome of the marsh rat Holochilus nanus (Cricetidae: Sigmodontinae)

Repetitive DNA are sequences repeated hundreds or thousands of times and an abundant part of eukaryotic genomes. SatDNA represents the majority of the repetitive sequences, followed by transposable elements. The species Holochilus nanus (HNA) belongs to the rodent tribe Oryzomyini, the most taxonomically diverse of Sigmodontinae subfamily. Cytogenetic studies on Oryzomyini reflect such diversity by revealing an exceptional range of karyotype variability. However, little is known about the repetitive DNA content and its involvement in chromosomal diversification of these species. In the search for a more detailed understanding about the composition of repetitive DNA on the genome of HNA and other species of Oryzomyini, we employed a combination of bioinformatic, cytogenetic and molecular techniques to characterize the repetitive DNA content of these species. RepeatExplorer analysis showed that almost half of repetitive content of HNA genome are composed by Long Terminal Repeats and a less significant portion are composed by Short Interspersed Nuclear Elements and Long Interspersed Nuclear Elements. RepeatMasker showed that more than 30% of HNA genome are composed by repetitive sequences, with two main waves of repetitive element insertion. It was also possible to identify a satellite DNA sequence present in the centromeric region of Oryzomyini species, and a repetitive sequence enriched on the long arm of HNA X chromosome. Also, comparative analysis between HNA genome with and without B chromosome did not evidence any repeat element enriched on the supernumerary, suggesting that B chromosome of HNA is composed by a fraction of repeats from all the genome.

Development and applications of a collection of single copy gene-based cytogenetic DNA markers in garden asparagus

Garden asparagus (Asparagus officinalis, 2n = 2x = 20 chromosomes) is an important dioecious vegetable crop and a model species for studying sex chromosome formation and evolution. However, few molecular cytogenetic studies on garden asparagus have been reported because of its small metaphase chromosomes, the scarcity of distinguished cytogenetic markers, and the high content of repetitive sequences. In this study, a set of single copy genes free of repetitive sequences with sizes ranging from 4.3 kb to 8.2 kb were screened and used as probes for fluorescence in situ hybridization (FISH) to identify individual chromosomes of garden asparagus. The chromosome-specific signal distribution patterns of these probes enabled the distinguishment of each pair of chromosomes. The sequence assembly and cytogenetic map were successfully integrated, and the results confirmed that the chromosome 1 representing the sex chromosome in the genome assembly is chromosome 5 in the karyotype analysis. The cytogenetic identification of the male-specific region of the Y chromosome (MSY) was implemented using a mixed probe derived from a number of MSY-specific single copy sequences. In addition, the chromosome orthologous relationship between garden asparagus (A1-A10, karyotypic analysis) and its hermaphrodite close relative, A. setaceus (B1-B10, karyotypic analysis), was analyzed using this collection of chromosome-specific cytological markers. The results showed that B3 is the ortholog of sex chromosome A5 and thus may represent the ancestral autosome of the current sex chromosome in garden asparagus. Chromosomes B5, B4, B1, B8, B7, and B9 are the orthologs of A2, A3, A4, A7, A8, and A10, respectively. The chromosome identification, cytogenetic recognition of MSY, and the orthologous relationship analysis between garden asparagus and A. setaceus are valuable for the further investigation of the sex chromosome emergence and evolutionary mechanism of garden asparagus and genome structure evolution in the Asparagus genus.

Impact of LTR-Retrotransposons on Genome Structure, Evolution, and Function in Curcurbitaceae Species

Long terminal repeat (LTR)-retrotransposons (LTR-RTs) comprise a major portion of many plant genomes and may exert a profound impact on genome structure, function, and evolution. Although many studies have focused on these elements in an individual species, their dynamics on a family level remains elusive. Here, we investigated the abundance, evolutionary dynamics, and impact on associated genes of LTR-RTs in 16 species in an economically important plant family, Cucurbitaceae. Results showed that full-length LTR-RT numbers and LTR-RT content varied greatly among different species, and they were highly correlated with genome size. Most of the full-length LTR-RTs were amplified after the speciation event, reflecting the ongoing rapid evolution of these genomes. LTR-RTs highly contributed to genome size variation via species-specific distinct proliferations. The Angela and Tekay lineages with a greater evolutionary age were amplified in Trichosanthes anguina, whereas a recent activity burst of Reina and another ancient round of Tekay activity burst were examined in Sechium edule. In addition, Tekay and Retand lineages belonging to the Gypsy superfamily underwent a recent burst in Gynostemma pentaphyllum. Detailed investigation of genes with intronic and promoter LTR-RT insertion showed diverse functions, but the term of metabolism was enriched in most species. Further gene expression analysis in G.pentaphyllum revealed that the LTR-RTs within introns suppress the corresponding gene expression, whereas the LTR-RTs within promoters exert a complex influence on the downstream gene expression, with the main function of promoting gene expression. This study provides novel insights into the organization, evolution, and function of LTR-RTs in Cucurbitaceae genomes.

The Genomics of Plant Satellite DNA

The twenty-first century began with a certain indifference to the research of satellite DNA (satDNA). Neither genome sequencing projects were able to accurately encompass the study of satDNA nor classic methodologies were able to go further in undertaking a better comprehensive study of the whole set of satDNA sequences of a genome. Nonetheless, knowledge of satDNA has progressively advanced during this century with the advent of new analytical techniques. The enormous advantages that genome-wide approaches have brought to its analysis have now stimulated a renewed interest in the study of satDNA. At this point, we can look back and try to assess more accurately many of the key questions that were left unsolved in the past about this enigmatic and important component of the genome. I review here the understanding gathered on plant satDNAs over the last few decades with an eye on the near future.

Using carrot centromeric repeats to study karyotype relationships in the genus Daucus (Apiaceae)

Background

In the course of evolution, chromosomes undergo evolutionary changes; thus, karyotypes may differ considerably among groups of organisms, even within closely related taxa. The genus Daucus seems to be a promising model for exploring the dynamics of karyotype evolution. It comprises some 40 wild species and the cultivated carrot, a crop of great economic significance. However, Daucus species are very diverse morphologically and genetically, and despite extensive research, the taxonomic and phylogenetic relationships between them have still not been fully resolved. Although several molecular cytogenetic studies have been conducted to investigate the chromosomal structure and karyotype evolution of carrot and other Daucus species, detailed karyomorphological research has been limited to carrot and only a few wild species. Therefore, to better understand the karyotype relationships within Daucus, we (1) explored the chromosomal distribution of carrot centromeric repeats (CentDc) in 34 accessions of Daucus and related species by means of fluorescence in situ hybridization (FISH) and (2) performed detailed karyomorphological analysis in 16 of them.

Results

We determined the genomic organization of CentDc in 26 accessions of Daucus (belonging to both Daucus I and II subclades) and one accession of closely related species. The CentDc repeats were present in the centromeric regions of all chromosomes of 20 accessions (representing 11 taxa). In the other Daucus taxa, the number of chromosome pairs with CentDc signals varied depending on the species, yet their centromeric localization was conserved. In addition, precise chromosome measurements performed in 16 accessions showed the inter- and intraspecific karyological relationships among them.

Conclusions

The presence of the CentDc repeats in the genomes of taxa belonging to both Daucus subclades and one outgroup species indicated the ancestral status of the repeat. The results of our study provide useful information for further evolutionary, cytotaxonomic, and phylogenetic research on the genus Daucus and may contribute to a better understanding of the dynamic evolution of centromeric satellites in plants.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07853-2.

Full-length LTR retroelements in Capsicum annuum revealed a few species-specific family bursts with insertional preferences

Capsicum annuum is a species that has undergone an expansion of the size of its genome caused mainly by the amplification of repetitive DNA sequences, including mobile genetic elements. Based on information obtained from sequencing the genome of pepper, the estimated fraction of retroelements is approximately 81%, and previous results revealed an important contribution of lineages derived from Gypsy superfamily. However, the dynamics of the retroelements in the C. annuum genome is poorly understood. In this way, the present work seeks to investigate the phylogenetic diversity and genomic abundance of the families of autonomous (complete and intact) LTR retroelements from C. annuum and inspect their distribution along its chromosomes. In total, we identified 1151 structurally full-length retroelements (340 Copia; 811 Gypsy) grouped in 124 phylogenetic families in the base of their retrotranscriptase. All the evolutive lineages of LTR retroelements identified in plants were present in pepper; however, three of them comprise 83% of the entire LTR retroelements population, the lineages Athila, Del/Tekay, and Ale/Retrofit. From them, only three families represent 70.8% of the total number of the identified retroelements. A massive family-specific wave of amplification of two of them occurred in the last 0.5 Mya (GypsyCa_16; CopiaCa_01), whereas the third is more ancient and occurred 3.0 Mya (GypsyCa_13). Fluorescent in situ hybridization performed with family and lineage-specific probes revealed contrasting patterns of chromosomal affinity. Our results provide a database of the populations LTR retroelements specific to C. annuum genome. The most abundant families were analyzed according to chromosome insertional preferences, suppling useful tools to the design of retroelement-based markers specific to the species.

Multiple independent recombinations led to hermaphroditism in grapevine

We studied the grape sex-determining region (SDR) in 12 Vitis genomes and demonstrated its conservation across 556 genotypes including 193 accessions from 47 world-wide wild grapevine species and 363 accessions of cultivated grapevine. Although the grape SDR is recombination free in all wild species, we found two distinct hermaphrodite (H) haplotypes (H1 and H2) among the cultivated grapevines, both chimeras of male (M) and female (f) haplotypes. The two independent recombinations carry different genetic signatures which long predate the domestication of grapevine, suggesting independent evolutions of this trait in wild European grapevine gene pools prior to human domestication.

Hermaphroditic (perfect) flowers were a key trait in grapevine domestication, enabling a drastic increase in yields due to the efficiency of self-pollination in the domesticated grapevine (Vitis vinifera L. ssp. vinifera). In contrast, all extant wild Vitis species are dioecious, each plant having only male or female flowers. In this study, we identified the male (M) and female (f) haplotypes of the sex-determining region (SDR) in the wild grapevine species V. cinerea and confirmed the boundaries of the SDR. We also demonstrated that the SDR and its boundaries are precisely conserved across the Vitis genus using shotgun resequencing data of 556 wild and domesticated accessions from North America, East Asia, and Europe. A high linkage disequilibrium was found at the SDR in all wild grape species, while different recombination signatures were observed along the hermaphrodite (H) haplotype of 363 cultivated accessions, revealing two distinct H haplotypes, named H1 and H2. To further examine the H2 haplotype, we sequenced the genome of two grapevine cultivars, 'Riesling' and 'Chardonnay'. By reconstructing the first two H2 haplotypes, we estimated the divergence time between H1 and H2 haplotypes at ∼6 million years ago, which predates the domestication of grapevine (∼8,000 y ago). Our findings emphasize the important role of recombination suppression in maintaining dioecy in wild grape species and lend additional support to the hypothesis that at least two independent recombination events led to the reversion to hermaphroditism in grapevine.

Nutritional characterization and food value addition properties of dehydrated spinach powder

This study aimed at investigating the physicochemical and bread-making features of dehydrated spinach. Physicochemical composition of spinach powder was compared with wheat flour and the effect of spinach powder supplementation on the nutritional composition, dough rheology, and quality attributes of chapatti were assessed. The results suggested spinach powder to be holding 8.2% crude fiber, 19.2% protein, 1,304 mg/100g calcium, and 40.4 mg/100g iron. Spinach powder indicated significantly increased values for hygroscopicity, swelling power, and water solubility index values, that is, 6.4%, 7.1 g/g, and 4.2%, respectively, when compared with wheat flour. Supplementation of spinach powder in wheat flour at 20% substitution level significantly reduced dough development properties including water absorption, dough stability, and peak dough development time. Color measurements of baked chapatti indicated a significant reduction in L*, a*, and chroma values with increasing the level of spinach powder supplementation; however, sensory profiling confirmed that supplementation of spinach powder at 7.5% had an optimum effect on the overall acceptability of the baked product. The results further suggested that replacing wheat flour with spinach powder (5%-7.5%, w/w) in baked products could be a viable dietary approach to enhance the optimum supply of micronutrients and to combat micronutrient deficiencies among various population segments.

Genome-Wide Analysis of Transposable Elements and Satellite DNAs in Spinacia Species to Shed Light on Their Roles in Sex Chromosome Evolution

Sex chromosome evolution has mostly been studied in species with heteromorphic sex chromosomes. The Spinacia genus serves as an ideal model for investigating evolutionary mechanisms underlying the transition from homomorphic to heteromorphic sex chromosomes. Among evolutionary factors, repetitive sequences play multiple roles in sex chromosome evolution while their forces have not been fully explored in Spinacia species. Here, we identified major repetitive sequence classes in male and female genomes of Spinacia species and their ancestral relative sugar beet to elucidate the evolutionary processes of sex chromosome evolution using next-generation sequencing (NGS) data. Comparative analysis revealed that the repeat elements of Spinacia species are considerably higher than of sugar beet, especially the Ty3/Gypsy and Ty1/Copia retrotransposons. The long terminal repeat retroelements (LTR) Angela, Athila, and Ogre may be accounted for the higher proportion of repeats in the spinach genome. Comparison of the repeats proportion between female and male genomes of three Spinacia species indicated the different representation in Spinacia tetrandra samples but not in the S. oleracea or S. turkestanica samples. From these results, we speculated that emergence of repetitive DNA sequences may correlate the formation of sex chromosome and the transition from homomorphic sex chromosomes to heteromorphic sex chromosomes as heteromorphic sex chromosomes exclusively existed in Spinacia tetrandra. Three novel sugar beet-specific satellites were identified and confirmed by fluorescence in situ hybridization (FISH); six out of eight new spinach-specific satellites were mapped to the short arm of sex chromosomes. A total of 141 copies of SolSat01-171-s were found in the sex determination region (SDR). Thus, the accumulation of satellite DNA on the short arm of chromosome 1 may be involved in the sex chromosome evolution in Spinacia species. Our study provides a fundamental resource for understanding repeat sequences in Spinacia species and their roles in sex chromosome evolution.

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.

DNA methylation and its effects on gene expression during primary to secondary growth in poplar stems

Background

As an important epigenetic mark, 5-methylcytosine (5mC) methylation is involved in many DNA-dependent biological processes and plays a role during development and differentiation of multicellular organisms. However, there is still a lack of knowledge about the dynamic aspects and the roles of global 5mC methylation in wood formation in tree trunks. In this study, we not only scrutinized single-base resolution methylomes of primary stems (PS), transitional stems (TS), and secondary stems (SS) of Populus trichocarpa using a high-throughput bisulfite sequencing technique, but also analyzed the effects of 5mC methylation on the expression of genes involved in wood formation.

Results

The overall average percentages of CG, CHG, and CHH methylation in poplar stems were ~ 53.6%, ~ 37.7%, and ~ 8.5%, respectively, and the differences of 5mC in genome-wide CG/CHG/CHH contexts among PS, TS, and SS were statistically significant (p < 0.05). The evident differences in CG, CHG, and CHH methylation contexts among 2 kb proximal promoters, gene bodies, and 2 kb downstream regions were observed among PS, TS, and SS. Further analysis revealed a perceptible global correlation between 5mC methylation levels of gene bodies and transcript levels but failed to reveal a correlation between 5mC methylation levels of proximal promoter regions and transcript levels. We identified 653 and 858 DMGs and 4978 and 4780 DEGs in PS vs TS and TS vs SS comparisons, respectively. Only 113 genes of 653 DMGs and 4978 DEGs, and 114 genes of 858 DMGs and 4780 DEG were common. Counterparts of some of these common genes in other species, including Arabidopsis thaliana, are known to be involved in secondary cell wall biosynthesis and hormone signaling. This indicates that methylation may directly modulate wood formation genes and indirectly attune hormone signaling genes, which in turn impact wood formation.

Conclusions

DNA methylation only marginally affects pathway genes or regulators involved in wood formation, suggesting that further studies of wood formation should lean towards the indirect effects of methylation. The information and data we provide here will be instrumental for understanding the roles of methylation in wood formation in tree species.

Characterization of Two Satellite DNA Families in the Genome of the Oomycete Plant Pathogen Phytophthora parasitica

Satellite DNA is a class of repetitive sequences that are organized in long arrays of tandemly repeated units in most eukaryotes. Long considered as selfish DNA, satellite sequences are now proposed to contribute to genome integrity. Despite their potential impact on the architecture and evolution of the genome, satellite DNAs have not been investigated in oomycetes due to the paucity of genomic data and the difficulty of assembling highly conserved satellite arrays. Yet gaining knowledge on the structure and evolution of genomes of oomycete pathogens is crucial to understanding the mechanisms underlying adaptation to their environment and to proposing efficient disease control strategies. A de novo assembly of the genome of Phytophthora parasitica, an important oomycete plant pathogen, led to the identification of several families of tandemly repeated sequences varying in size, copy number, and sequence conservation. Among them, two abundant families, designated as PpSat1 and PpSat2, displayed typical features of satellite DNA and were collectively designated as PpSat. These two satellite families differ by their length, sequence, organization, genomic environment, and evolutionary dynamics. PpSat1, but not PpSat2, presented homologs among oomycetes. This observation, as well as the characterization of transcripts of PpSat families, suggested that these satellite DNA families likely play a conserved role within this important group of pathogens.

Sex Chromosome Evolution: So Many Exceptions to the Rules

Genomic analysis of many nonmodel species has uncovered an incredible diversity of sex chromosome systems, making it possible to empirically test the rich body of evolutionary theory that describes each stage of sex chromosome evolution. Classic theory predicts that sex chromosomes originate from a pair of homologous autosomes and recombination between them is suppressed via inversions to resolve sexual conflict. The resulting degradation of the Y chromosome gene content creates the need for dosage compensation in the heterogametic sex. Sex chromosome theory also implies a linear process, starting from sex chromosome origin and progressing to heteromorphism. Despite many convergent genomic patterns exhibited by independently evolved sex chromosome systems, and many case studies supporting these theoretical predictions, emerging data provide numerous interesting exceptions to these long-standing theories, and suggest that the remarkable diversity of sex chromosomes is matched by a similar diversity in their evolution. For example, it is clear that sex chromosome pairs are not always derived from homologous autosomes. In addition, both the cause and the mechanism of recombination suppression between sex chromosome pairs remain unclear, and it may be that the spread of recombination suppression is a more gradual process than previously thought. It is also clear that dosage compensation can be achieved in many ways, and displays a range of efficacy in different systems. Finally, the remarkable turnover of sex chromosomes in many systems, as well as variation in the rate of sex chromosome divergence, suggest that assumptions about the inevitable linearity of sex chromosome evolution are not always empirically supported, and the drivers of the birth-death cycle of sex chromosome evolution remain to be elucidated. Here, we concentrate on how the diversity in sex chromosomes across taxa highlights an equal diversity in each stage of sex chromosome evolution.

A Chromosome-Scale Assembly of the Garden Orach (Atriplex hortensis L.) Genome Using Oxford Nanopore Sequencing

Atriplex hortensis (2n = 2x = 18, 1C genome size ∼1.1 gigabases), also known as garden orach and mountain-spinach, is a highly nutritious, broadleaf annual of the Amaranthaceae-Chenopodiaceae alliance (Chenopodiaceae sensu stricto, subfam. Chenopodioideae) that has spread in cultivation from its native primary domestication area in Eurasia to other temperate and subtropical regions worldwide. Atriplex L. is a highly complex but, as understood now, a monophyletic group of mainly halophytic and/or xerophytic plants, of which A. hortensis has been a vegetable of minor importance in some areas of Eurasia (from Central Asia to the Mediterranean) at least since antiquity. Nonetheless, it is a crop with tremendous nutritional potential due primarily to its exceptional leaf and seed protein quantities (approaching 30%) and quality (high levels of lysine). Although there is some literature describing the taxonomy and production of A. hortensis, there is a general lack of genetic and genomic data that would otherwise help elucidate the genetic variation, phylogenetic positioning, and future potential of the species. Here, we report the assembly of the first high-quality, chromosome-scale reference genome for A. hortensis cv. "Golden." Long-read data from Oxford Nanopore's MinION DNA sequencer was assembled with the program Canu and polished with Illumina short reads. Contigs were scaffolded to chromosome scale using chromatin-proximity maps (Hi-C) yielding a final assembly containing 1,325 scaffolds with a N50 of 98.9 Mb - with 94.7% of the assembly represented in the nine largest, chromosome-scale scaffolds. Sixty-six percent of the genome was classified as highly repetitive DNA, with the most common repetitive elements being Gypsy-(32%) and Copia-like (11%) long-terminal repeats. The annotation was completed using MAKER which identified 37,083 gene models and 2,555 tRNA genes. Completeness of the genome, assessed using the Benchmarking Universal Single Copy Orthologs (BUSCO) metric, identified 97.5% of the conserved orthologs as complete, with only 2.2% being duplicated, reflecting the diploid nature of A. hortensis. A resequencing panel of 21 wild, unimproved and cultivated A. hortensis accessions revealed three distinct populations with little variation within subpopulations. These resources provide vital information to better understand A. hortensis and facilitate future study.

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.

Whole genome sequencing and comparative genomic analysis of oleaginous red yeast Sporobolomyces pararoseus NGR identifies candidate genes for biotechnological potential and ballistospores-shooting

BACKGROUND

Sporobolomyces pararoseus is regarded as an oleaginous red yeast, which synthesizes numerous valuable compounds with wide industrial usages. This species hold biotechnological interests in biodiesel, food and cosmetics industries. Moreover, the ballistospores-shooting promotes the colonizing of S. pararoseus in most terrestrial and marine ecosystems. However, very little is known about the basic genomic features of S. pararoseus. To assess the biotechnological potential and ballistospores-shooting mechanism of S. pararoseus on genome-scale, the whole genome sequencing was performed by next-generation sequencing technology.

RESULTS

Here, we used Illumina Hiseq platform to firstly assemble S. pararoseus genome into 20.9 Mb containing 54 scaffolds and 5963 predicted genes with a N50 length of 2,038,020 bp and GC content of 47.59%. Genome completeness (BUSCO alignment: 95.4%) and RNA-seq analysis (expressed genes: 98.68%) indicated the high-quality features of the current genome. Through the annotation information of the genome, we screened many key genes involved in carotenoids, lipids, carbohydrate metabolism and signal transduction pathways. A phylogenetic assessment suggested that the evolutionary trajectory of the order Sporidiobolales species was evolved from genus Sporobolomyces to Rhodotorula through the mediator Rhodosporidiobolus. Compared to the lacking ballistospores Rhodotorula toruloides and Saccharomyces cerevisiae, we found genes enriched for spore germination and sugar metabolism. These genes might be responsible for the ballistospores-shooting in S. pararoseus NGR.

CONCLUSION

These results greatly advance our understanding of S. pararoseus NGR in biotechnological potential and ballistospores-shooting, which help further research of genetic manipulation, metabolic engineering as well as its evolutionary direction.

Development and Deployment of High-Throughput Retrotransposon-Based Markers Reveal Genetic Diversity and Population Structure of Asian Bamboo

Bamboo, a non-timber grass species, known for exceptionally fast growth is a commercially viable crop. Long terminal repeat (LTR) retrotransposons, the main class I mobile genetic elements in plant genomes, are highly abundant (46%) in bamboo, contributing to genome diversity. They play significant roles in the regulation of gene expression, chromosome size and structure as well as in genome integrity. Due to their random insertion behavior, interspaces of retrotransposons can vary significantly among bamboo genotypes. Capitalizing this feature, inter-retrotransposon amplified polymorphism (IRAP) is a high-throughput marker system to study the genetic diversity of plant species. To date, there are no transposon based markers reported from the bamboo genome and particularly using IRAP markers on genetic diversity. Phyllostachys genus of Asian bamboo is the largest of the Bambusoideae subfamily, with great economic importance. We report structure-based analysis of bamboo genome for the LTR-retrotransposon superfamilies, Ty3-gypsy and Ty1-copia, which revealed a total of 98,850 retrotransposons with intact LTR sequences at both the ends. Grouped into 64,281 clusters/scaffold using CD-HIT-EST software, only 13 clusters of retroelements were found with more than 30 LTR sequences and with at least one copy having all intact protein domains such as gag and polyprotein. A total of 16 IRAP primers were synthesized, based on the high copy numbers of conserved LTR sequences. A study using these IRAP markers on genetic diversity and population structure of 58 Asian bamboo accessions belonging to the genus Phyllostachys revealed 3340 amplicons with an average of 98% polymorphism. The bamboo accessions were collected from nine different provinces of China, as well as from Italy and America. A three phased approach using hierarchical clustering, principal components and a model based population structure divided the bamboo accessions into four sub-populations, PhSP1, PhSP2, PhSP3 and PhSP4. All the three analyses produced significant sub-population wise consensus. Further, all the sub-populations revealed admixture of alleles. The analysis of molecular variance (AMOVA) among the sub-populations revealed high intra-population genetic variation (75%) than inter-population. The results suggest that Phyllostachys bamboos are not well evolutionarily diversified, although geographic speciation could have occurred at a limited level. This study highlights the usability of IRAP markers in determining the inter-species variability of Asian bamboos.

Cytogenetic and genomic organization analyses of chloroplast DNA invasions in the nuclear genome of Asparagus officinalis L. provides signatures of evolutionary complexity and informativity in sex chromosome evolution

BACKGROUND

The transfer of chloroplast DNA into nuclear genome is a common process in plants. These transfers form nuclear integrants of plastid DNAs (NUPTs), which are thought to be driving forces in genome evolution, including sex chromosome evolution. In this study, NUPTs in the genome of a dioecious plant Asparagus officinalis L. were systematically analyzed, in order to investigate the characteristics of NUPTs in the nuclear genome and the relationship between NUPTs and sex chromosome evolution in this species.

RESULTS

A total of 3155 NUPT insertions were detected, and they represented approximated 0.06% of the nuclear genome. About 45% of the NUPTs were organized in clusters. These clusters were derived from various evolutionary events. The Y chromosome contained the highest number and largest proportion of NUPTs, suggesting more accumulation of NUPTs on sex chromosomes. NUPTs were distributed widely in all of the chromosomes, and some regions preferred these insertions. The highest density of NUPTs was found in a 47 kb region in the Y chromosome; more than 75% of this region was occupied by NUPTs. Further cytogenetic and sequence alignment analysis revealed that this region was likely the centromeric region of the sex chromosomes. On the other hand, the male-specific region of the Y chromosome (MSY) and the adjacent regions did not have NUPT insertions.

CONCLUSIONS

These results indicated that NUPTs were involved in shaping the genome of A. officinalis through complicated process. NUPTs may play important roles in the centromere shaping of the sex chromosomes of A. officinalis, but were not implicated in MSY formation.