A sheep pangenome reveals the spectrum of structural variations and their effects on tail phenotypes.

Structural variations (SVs) are a major contributor to genetic diversity and phenotypic variations, but their prevalence and functions in domestic animals are largely unexplored. Here we generated high-quality genome assemblies for 15 individuals from genetically diverse sheep breeds using Pacific Biosciences (PacBio) high-fidelity sequencing, discovering 130.3 Mb nonreference sequences, from which 588 genes were annotated. A total of 149,158 biallelic insertions/deletions, 6531 divergent alleles, and 14,707 multiallelic variations with precise breakpoints were discovered. The SV spectrum is characterized by an excess of derived insertions compared to deletions (94,422 vs. 33,571), suggesting recent active LINE expansions in sheep. Nearly half of the SVs display low to moderate linkage disequilibrium with surrounding single-nucleotide polymorphisms (SNPs) and most SVs cannot be tagged by SNP probes from the widely used ovine 50K SNP chip. We identified 865 population-stratified SVs including 122 SVs possibly derived in the domestication process among 690 individuals from sheep breeds worldwide. A novel 168-bp insertion in the 5' untranslated region (5' UTR) of HOXB13 is found at high frequency in long-tailed sheep. Further genome-wide association study and gene expression analyses suggest that this mutation is causative for the long-tail trait. In summary, we have developed a panel of high-quality de novo assemblies and present a catalog of structural variations in sheep. Our data capture abundant candidate functional variations that were previously unexplored and provide a fundamental resource for understanding trait biology in sheep.

Epistasis and pleiotropy-induced variation for plant breeding.

Epistasis refers to nonallelic interaction between genes that cause bias in estimates of genetic parameters for a phenotype with interactions of two or more genes affecting the same trait. Partitioning of epistatic effects allows true estimation of the genetic parameters affecting phenotypes. Multigenic variation plays a central role in the evolution of complex characteristics, among which pleiotropy, where a single gene affects several phenotypic characters, has a large influence. While pleiotropic interactions provide functional specificity, they increase the challenge of gene discovery and functional analysis. Overcoming pleiotropy-based phenotypic trade-offs offers potential for assisting breeding for complex traits. Modelling higher order nonallelic epistatic interaction, pleiotropy and non-pleiotropy-induced variation, and genotype × environment interaction in genomic selection may provide new paths to increase the productivity and stress tolerance for next generation of crop cultivars. Advances in statistical models, software and algorithm developments, and genomic research have facilitated dissecting the nature and extent of pleiotropy and epistasis. We overview emerging approaches to exploit positive (and avoid negative) epistatic and pleiotropic interactions in a plant breeding context, including developing avenues of artificial intelligence, novel exploitation of large-scale genomics and phenomics data, and involvement of genes with minor effects to analyse epistatic interactions and pleiotropic quantitative trait loci, including missing heritability.

Phylogenomic analysis reveals five independently evolved African forage grass clades in the genus Urochloa.

BACKGROUND AND AIMS: The grass genus Urochloa (Brachiaria) sensu lato includes forage crops that are important for beef and dairy industries in tropical and sub-tropical Africa, South America and Oceania/Australia. Economically important species include U. brizantha, U. decumbens, U. humidicola, U. mutica, U. arrecta, U. trichopus, U. mosambicensis and Megathyrsus maximus, all native to the African continent. Perennial growth habits, large, fast growing palatable leaves, intra- and interspecific morphological variability, apomictic reproductive systems and frequent polyploidy are widely shared within the genus. The combination of these traits probably favoured the selection for forage domestication and weediness, but trait emergence across Urochloa cannot be modelled, as a robust phylogenetic assessment of the genus has not been conducted. We aim to produce a phylogeny for Urochloa that includes all important forage species, and identify their closest wild relatives (crop wild relatives). Finally, we will use our phylogeny and available trait data to infer the ancestral states of important forage traits across Urochloa s.l. and model the evolution of forage syndromes across the genus. METHODS: Using a target enrichment sequencing approach (Angiosperm 353), we inferred a species-level phylogeny for Urochloa s.l., encompassing 54 species (~40 % of the genus) and outgroups. Phylogenies were inferred using a multispecies coalescent model and maximum likelihood method. We determined the phylogenetic placement of agriculturally important species and identified their closest wild relatives, or crop wild relatives, based on well-supported monophyly. Further, we mapped key traits associated with Urochloa forage crops to the species tree and estimated ancestral states for forage traits along branch lengths for continuous traits and at ancestral nodes in discrete traits. KEY RESULTS: Agricultural species belong to five independent clades, including U. brizantha and U. decumbens lying in a previously defined species complex. Crop wild relatives were identified for these clades supporting previous sub-generic groupings in Urochloa based on morphology. Using ancestral trait estimation models, we find that five morphological traits that correlate with forage potential (perennial growth habits, culm height, leaf size, a winged rachis and large seeds) independently evolved in forage clades. CONCLUSIONS: Urochloa s.l. is a highly diverse genus that contains numerous species with agricultural potential, including crop wild relatives that are currently underexploited. All forage species and their crop wild relatives naturally occur on the African continent and their conservation across their native distributions is essential. Genomic and phenotypic diversity in forage clade species and their wild relatives need to be better assessed both to develop conservation strategies and to exploit the diversity in the genus for improved sustainability in Urochloa cultivar production.

Maintenance and expansion of genetic and trait variation following domestication in a clonal crop.

Clonal propagation enables favourable crop genotypes to be rapidly selected and multiplied. However, the absence of sexual propagation can lead to low genetic diversity and accumulation of deleterious mutations, which may eventually render crops less resilient to pathogens or environmental change. To better understand this trade-off, we characterize the domestication and contemporary genetic diversity of Enset (Ensete ventricosum), an indigenous African relative of bananas (Musa) and a principal starch staple for 20 million Ethiopians. Wild enset reproduction occurs strictly by sexual outcrossing, but for cultivation, it is propagated clonally and associated with diversification and specialization into hundreds of named landraces. We applied tGBS sequencing to generate genome-wide genotypes for 192 accessions from across enset's cultivated distribution, and surveyed 1340 farmers on enset agronomic traits. Overall, reduced heterozygosity in the domesticated lineage was consistent with a domestication bottleneck that retained 37% of wild diversity. However, an excess of putatively deleterious missense mutations at low frequency present as heterozygotes suggested an accumulation of mutational load in clonal domesticated lineages. Our evidence indicates that the major domesticated lineages initially arose through historic sexual recombination associated with a domestication bottleneck, followed by the amplification of favourable genotypes through an extended period of clonal propagation. Among domesticated lineages, we found a significant phylogenetic signal for multiple farmer-identified food, nutrition and disease resistance traits and little evidence of contemporary recombination. The development of future-climate adapted genotypes may require crop breeding, but outcrossing risks exposing deleterious alleles as homozygotes. This trade-off may partly explain the ubiquity and persistence of clonal propagation over recent centuries of comparative climate stability.

Anchoring alien chromosome segment substitutions bearing gene(s) for resistance to mustard aphid in Brassica juncea-B. fruticulosa introgression lines and their possible disruption through gamma irradiation.

KEY MESSAGE: Heavy doses of gamma irradiation can reduce linkage drag by disrupting large sized alien translocations and promoting exchanges between crop and wild genomes. Resistance to mustard aphid (Lipaphis erysimi) infestation was significantly improved in Brassica juncea through B. juncea-B. fruticulosa introgression. However, linkage drag caused by introgressed chromatin fragments has so far prevented the deployment of this resistance source in commercial cultivars. We investigated the patterns of donor chromatin segment substitutions in the introgression lines (ILs) through genomic in situ hybridization (GISH) coupled with B. juncea chromosome-specific oligonucleotide probes. These allowed identification of large chromosome translocations from B. fruticulosa in the terminal regions of chromosomes A05, B02, B03 and B04 in three founder ILs (AD-64, 101 and 104). Only AD-101 carried an additional translocation at the sub-terminal to intercalary position in both homologues of chromosome A01. We validated these translocations with a reciprocal blast hit analysis using shotgun sequencing of three ILs and species-specific contigs/scaffolds (kb sized) from a de novo assembly of B. fruticulosa. Alien segment substitution on chromosome A05 could not be validated. Current studies also endeavoured to break linkage drag by exposing seeds to a heavy dose (200kR) of gamma radiation. Reduction in the size of introgressed chromatin fragments was observed in many M3 plants. There was a complete loss of the alien chromosome fragment in one instance. A few M3 plants with novel patterns of chromosome segment substitutions displayed improved agronomic performance coupled with resistance to mustard aphid. SNPs in such genomic spaces should aid the development of markers to track introgressed DNA and allow application in plant breeding.

The EU chemicals strategy for sustainability: in support of the BfR position.

The EU chemicals strategy for sustainability (CSS) asserts that both human health and the environment are presently threatened and that further regulation is necessary. In a recent Guest Editorial, members of the German competent authority for risk assessment, the BfR, raised concerns about the scientific justification for this strategy. The complexity and interdependence of the networks of regulation of chemical substances have ensured that public health and wellbeing in the EU have continuously improved. A continuous process of improvement in consumer protection is clearly desirable but any initiative directed towards this objective must be based on scientific knowledge. It must not confound risk with other factors in determining policy. This conclusion is fully supported in the present Commentary including the request to improve both, data collection and the time-consuming and bureaucratic procedures that delay the publication of regulations.

Human exposure to synthetic endocrine disrupting chemicals (S-EDCs) is generally negligible as compared to natural compounds with higher or comparable endocrine activity: how to evaluate the risk of the S-EDCs?

Theoretically, both synthetic endocrine disrupting chemicals (S-EDCs) and natural (exogenous and endogenous) endocrine disrupting chemicals (N-EDCs) can interact with endocrine receptors and disturb hormonal balance. However, compared to endogenous hormones, S-EDCs are only weak partial agonists with receptor affinities several orders of magnitude lower. Thus, to elicit observable effects, S-EDCs require considerably higher concentrations to attain sufficient receptor occupancy or to displace natural hormones and other endogenous ligands. Significant exposures to exogenous N-EDCs may result from ingestion of foods such as soy-based diets, green tea and sweet mustard. While their potencies are lower as compared to natural endogenous hormones, they usually are considerably more potent than S-EDCs. Effects of exogenous N-EDCs on the endocrine system were observed at high dietary intakes. A causal relation between their mechanism of action and these effects is established and biologically plausible. In contrast, the assumption that the much lower human exposures to S-EDCs may induce observable endocrine effects is not plausible. Hence, it is not surprising that epidemiological studies searching for an association between S-EDC exposure and health effects have failed. Regarding testing for potential endocrine effects, a scientifically justified screen should use in vitro tests to compare potencies of S-EDCs with those of reference N-EDCs. When the potency of the S-EDC is similar or smaller than that of the N-EDC, further testing in laboratory animals and regulatory consequences are not warranted.

Human exposure to synthetic endocrine disrupting chemicals (S-EDCs) is generally negligible as compared to natural compounds with higher or comparable endocrine activity. How to evaluate the risk of the S-EDCs?

Theoretically, both synthetic endocrine-disrupting chemicals (S-EDCs) and natural (exogenous and endogenous) endocrine-disrupting chemicals (N-EDCs) can interact with endocrine receptors and disturb hormonal balance. However, compared to endogenous hormones, S-EDCs are only weak partial agonists with receptor affinities several orders of magnitude lower than S-EDCs. Thus, to elicit observable effects, S-EDCs require considerably higher concentrations to attain sufficient receptor occupancy or to displace natural hormones and other endogenous ligands. Significant exposures to exogenous N-EDCs may result from ingestion of foods such as soy-based diets, green tea, and sweet mustard. While their potencies are lower as compared to natural endogenous hormones, they usually are considerably more potent than S-EDCs. Effects of exogenous N-EDCs on the endocrine system were observed at high dietary intakes. A causal relation between their mechanism of action and these effects is established and biologically plausible. In contrast, the assumption that the much lower human exposures to S-EDCs may induce observable endocrine effects is not plausible. Hence, it is not surprising that epidemiological studies searching for an association between S-EDC exposure and health effects have failed. Regarding testing for potential endocrine effects, a scientifically justified screen should use in vitro tests to compare potencies of S-EDCs with those of reference N-EDCs. When the potency of the S-EDC is similar or smaller than that of the N-EDC, further testing in laboratory animals and regulatory consequences are not warranted.

The banana (Musa acuminata) genome and the evolution of monocotyledonous plants.

Bananas (Musa spp.), including dessert and cooking types, are giant perennial monocotyledonous herbs of the order Zingiberales, a sister group to the well-studied Poales, which include cereals. Bananas are vital for food security in many tropical and subtropical countries and the most popular fruit in industrialized countries. The Musa domestication process started some 7,000 years ago in Southeast Asia. It involved hybridizations between diverse species and subspecies, fostered by human migrations, and selection of diploid and triploid seedless, parthenocarpic hybrids thereafter widely dispersed by vegetative propagation. Half of the current production relies on somaclones derived from a single triploid genotype (Cavendish). Pests and diseases have gradually become adapted, representing an imminent danger for global banana production. Here we describe the draft sequence of the 523-megabase genome of a Musa acuminata doubled-haploid genotype, providing a crucial stepping-stone for genetic improvement of banana. We detected three rounds of whole-genome duplications in the Musa lineage, independently of those previously described in the Poales lineage and the one we detected in the Arecales lineage. This first monocotyledon high-continuity whole-genome sequence reported outside Poales represents an essential bridge for comparative genome analysis in plants. As such, it clarifies commelinid-monocotyledon phylogenetic relationships, reveals Poaceae-specific features and has led to the discovery of conserved non-coding sequences predating monocotyledon-eudicotyledon divergence.

Chromosomal distribution and evolution of abundant retrotransposons in plants: gypsy elements in diploid and polyploid Brachiaria forage grasses.

Like other eukaryotes, the nuclear genome of plants consists of DNA with a small proportion of low-copy DNA (genes and regulatory sequences) and very abundant DNA sequence motifs that are repeated thousands up to millions of times in the genomes including transposable elements (TEs) and satellite DNA. Retrotransposons, one class of TEs, are sequences that amplify via an RNA intermediate and reinsert into the genome, are often the major fraction of a genome. Here, we put research on retrotransposons into the larger context of plant repetitive DNA and genome behaviour, showing features of genome evolution in a grass genus, Brachiaria, in relation to other plant species. We show the contrasting amplification of different retroelement fractions across the genome with characteristics for various families and domains. The genus Brachiaria includes both diploid and polyploid species, with similar chromosome types and chromosome basic numbers x = 6, 7, 8 and 9. The polyploids reproduce asexually and are apomictic, but there are also sexual species. Cytogenetic studies and flow cytometry indicate a large variation in DNA content (C-value), chromosome sizes and genome organization. In order to evaluate the role of transposable elements in the genome and karyotype organization of species of Brachiaria, we searched for sequences similar to conserved regions of TEs in RNAseq reads library produced in Brachiaria decumbens. Of the 9649 TE-like contigs, 4454 corresponded to LTR-retrotransposons, and of these, 79.5 % were similar to members of the gypsy superfamily. Sequences of conserved protein domains of gypsy were used to design primers for producing the probes. The probes were used in FISH against chromosomes of accesses of B. decumbens, Brachiaria brizantha, Brachiaria ruziziensis and Brachiaria humidicola. Probes showed hybridization signals predominantly in proximal regions, especially those for retrotransposons of the clades CRM and Athila, while elements of Del and Tat exhibited dispersed signals, in addition to those proximal signals. These results show that the proximal region of Brachiaria chromosomes is a hotspot for retrotransposon insertion, particularly for the gypsy family. The combination of high-throughput sequencing and a chromosome-centric cytogenetic approach allows the abundance, organization and nature of transposable elements to be characterized in unprecedented detail. By their amplification and dispersal, retrotransposons can affect gene expression; they can lead to rapid diversification of chromosomes between species and, hence, are useful for studies of genome evolution and speciation in the Brachiaria genus. Centromeric regions can be identified and mapped, and retrotransposon markers can also assisting breeders in the developing and exploiting interspecific hybrids.

Integration of genetic and physical maps of the Primula vulgaris S locus and localization by chromosome in situ hybridization.

Heteromorphic flower development in Primula is controlled by the S locus. The S locus genes, which control anther position, pistil length and pollen size in pin and thrum flowers, have not yet been characterized. We have integrated S-linked genes, marker sequences and mutant phenotypes to create a map of the P. vulgaris S locus region that will facilitate the identification of key S locus genes. We have generated, sequenced and annotated BAC sequences spanning the S locus, and identified its chromosomal location. We have employed a combination of classical genetics and three-point crosses with molecular genetic analysis of recombinants to generate the map. We have characterized this region by Illumina sequencing and bioinformatic analysis, together with chromosome in situ hybridization. We present an integrated genetic and physical map across the P. vulgaris S locus flanked by phenotypic and DNA sequence markers. BAC contigs encompass a 1.5-Mb genomic region with 1 Mb of sequence containing 82 S-linked genes anchored to overlapping BACs. The S locus is located close to the centromere of the largest metacentric chromosome pair. These data will facilitate the identification of the genes that orchestrate heterostyly in Primula and enable evolutionary analyses of the S locus.

Genome-wide analysis of flavonoid biosynthetic genes in Musaceae (Ensete, Musella, and Musa species) reveals amplification of flavonoid 3',5'-hydroxylase.

Flavonoids in Musaceae are involved in pigmentation and stress responses, including cold resistance, and are a component of the healthy human diet. Identification and analysis of the sequence and copy number of flavonoid biosynthetic genes are valuable for understanding the nature and diversity of flavonoid evolution in Musaceae species. In this study, we identified 71-80 flavonoid biosynthetic genes in chromosome-scale genome sequence assemblies of Musaceae, including those of Ensete glaucum, Musella lasiocarpa, Musa beccarii, M. acuminata, M. balbisiana and M. schizocarpa, checking annotations with BLAST and determining the presence of conserved domains. The number of genes increased through segmental duplication and tandem duplication. Orthologues of both structural and regulatory genes in the flavonoid biosynthetic pathway are highly conserved across Musaceae. The flavonoid 3',5'-hydroxylase gene F3'5'H was amplified in Musaceae and ginger compared with grasses (rice, Brachypodium, Avena longiglumis, and sorghum). One group of genes from this gene family amplified near the centromere of chromosome 2 in the x = 11 Musaceae species. Flavonoid biosynthetic genes displayed few consistent responses in the yellow and red bracts of Musella lasiocarpa when subjected to low temperatures. The expression levels of MlDFR2/3 (dihydroflavonol reductase) increased while MlLAR (leucoanthocyanidin reductase) was reduced by half. Overall, the results establish the range of diversity in both sequence and copy number of flavonoid biosynthetic genes during evolution of Musaceae. The combination of allelic variants of genes, changes in their copy numbers, and variation in transcription factors with the modulation of expression under cold treatments and between genotypes with contrasting bract-colours suggests the variation may be exploited in plant breeding programmes, particularly for improvement of stress-resistance in the banana crop.

The repetitive component of the A genome of peanut (Arachis hypogaea) and its role in remodelling intergenic sequence space since its evolutionary divergence from the B genome.

BACKGROUND AND AIMS: Peanut (Arachis hypogaea) is an allotetraploid (AABB-type genome) of recent origin, with a genome of about 2·8 Gb and a high repetitive content. This study reports an analysis of the repetitive component of the peanut A genome using bacterial artificial chromosome (BAC) clones from A. duranensis, the most probable A genome donor, and the probable consequences of the activity of these elements since the divergence of the peanut A and B genomes. METHODS: The repetitive content of the A genome was analysed by using A. duranensis BAC clones as probes for fluorescence in situ hybridization (BAC-FISH), and by sequencing and characterization of 12 genomic regions. For the analysis of the evolutionary dynamics, two A genome regions are compared with their B genome homeologues. KEY RESULTS: BAC-FISH using 27 A. duranensis BAC clones as probes gave dispersed and repetitive DNA characteristic signals, predominantly in interstitial regions of the peanut A chromosomes. The sequences of 14 BAC clones showed complete and truncated copies of ten abundant long terminal repeat (LTR) retrotransposons, characterized here. Almost all dateable transposition events occurred <3·5 million years ago, the estimated date of the divergence of A and B genomes. The most abundant retrotransposon is Feral, apparently parasitic on the retrotransposon FIDEL, followed by Pipa, also non-autonomous and probably parasitic on a retrotransposon we named Pipoka. The comparison of the A and B genome homeologous regions showed conserved segments of high sequence identity, punctuated by predominantly indel regions without significant similarity. CONCLUSIONS: A substantial proportion of the highly repetitive component of the peanut A genome appears to be accounted for by relatively few LTR retrotransposons and their truncated copies or solo LTRs. The most abundant of the retrotransposons are non-autonomous. The activity of these retrotransposons has been a very significant driver of genome evolution since the evolutionary divergence of the A and B genomes.

Evolutionary dynamics and adaptive benefits of deleterious mutations in crop gene pools.

Mutations with deleterious consequences in nature may be conditionally deleterious in crop plants. That is, while some genetic variants may reduce fitness under wild conditions and be subject to purifying selection, they can be under positive selection in domesticates. Such deleterious alleles can be plant breeding targets, particularly for complex traits. The difficulty of distinguishing favorable from unfavorable variants reduces the power of selection, while favorable trait variation and heterosis may be attributable to deleterious alleles. Here, we review the roles of deleterious mutations in crop breeding and discuss how they can be used as a new avenue for crop improvement with emerging genomic tools, including HapMaps and pangenome analysis, aiding the identification, removal, or exploitation of deleterious mutations.

A design principle underlying the synchronization of oscillations in cellular systems.

Biological oscillations are found ubiquitously in cells and are widely variable, with periods varying from milliseconds to months, and scales involving subcellular components to large groups of organisms. Interestingly, independent oscillators from different cells often show synchronization that is not the consequence of an external regulator. What is the underlying design principle of such synchronized oscillations, and can modeling show that the complex consequences arise from simple molecular or other interactions between oscillators? When biological oscillators are coupled with each other, we found that synchronization is induced when they are connected together through a positive feedback loop. Increasing the coupling strength of two independent oscillators shows a threshold beyond which synchronization occurs within a few cycles, and a second threshold where oscillation stops. The positive feedback loop can be composed of either double-positive (PP) or double-negative (NN) interactions between a node of each of the two oscillating networks. The different coupling structures have contrasting characteristics. In particular, PP coupling is advantageous with respect to stability of period and amplitude, when local oscillators are coupled with a short time delay, whereas NN coupling is advantageous for a long time delay. In addition, PP coupling results in more robust synchronized oscillations with respect to amplitude excursions but not period, with applied noise disturbances compared to NN coupling. However, PP coupling can induce a large fluctuation in the amplitude and period of the resulting synchronized oscillation depending on the coupling strength, whereas NN coupling ensures almost constant amplitude and period irrespective of the coupling strength. Intriguingly, we have also observed that artificial evolution of random digital oscillator circuits also follows this design principle. We conclude that a different coupling strategy might have been selected according to different evolutionary requirements.

Evolutionary design principles and functional characteristics based on kingdom-specific network motifs.

BACKGROUND: Network motifs within biological networks show non-random abundances in systems at different scales. Large directed protein networks at the cellular level are now well defined in several diverse species. We aimed to compare the nature of significantly observed two- and three-node network motifs across three different kingdoms (Arabidopsis thaliana for multicellular plants, Saccharomyces cerevisiae for unicellular fungi and Homo sapiens for animals). RESULTS: 'Two-node feedback' is the most significant motif in all three species. By considering the sign of each two-node feedback interaction, we examined the enrichment of the three types of two-node feedbacks [positive-positive (PP), negative-negative (NN) and positive-negative (PN)]. We found that PN is enriched in the network of A.thaliana, NN in the network of S.cerevisiae and PP and NN in the network of H.sapiens. Each feedback type has characteristic features of robustness, multistability and homeostasis. CONCLUSIONS: We suggest that amplification of particular network motifs emerges from contrasting dynamical and topological properties of the motifs, reflects the evolutionary design principles selected by the characteristic behavior of each species and provides a signature pointing to their behavior and function.

Genome assembly of Musa beccarii shows extensive chromosomal rearrangements and genome expansion during evolution of Musaceae genomes.

BACKGROUND: Musa beccarii (Musaceae) is a banana species native to Borneo, sometimes grown as an ornamental plant. The basic chromosome number of Musa species is x = 7, 10, or 11; however, M. beccarii has a basic chromosome number of x = 9 (2n = 2x = 18), which is the same basic chromosome number of species in the sister genera Ensete and Musella. Musa beccarii is in the section Callimusa, which is sister to the section Musa. We generated a high-quality chromosome-scale genome assembly of M. beccarii to better understand the evolution and diversity of genomes within the family Musaceae. FINDINGS: The M. beccarii genome was assembled by long-read and Hi-C sequencing, and genes were annotated using both long Iso-seq and short RNA-seq reads. The size of M. beccarii was the largest among all known Musaceae assemblies (∼570 Mbp) due to the expansion of transposable elements and increased 45S ribosomal DNA sites. By synteny analysis, we detected extensive genome-wide chromosome fusions and fissions between M. beccarii and the other Musa and Ensete species, far beyond those expected from differences in chromosome number. Within Musaceae, M. beccarii showed a reduced number of terpenoid synthase genes, which are related to chemical defense, and enrichment in lipid metabolism genes linked to the physical defense of the cell wall. Furthermore, type III polyketide synthase was the most abundant biosynthetic gene cluster (BGC) in M. beccarii. BGCs were not conserved in Musaceae genomes. CONCLUSIONS: The genome assembly of M. beccarii is the first chromosome-scale genome assembly in the Callimusa section in Musa, which provides an important genetic resource that aids our understanding of the evolution of Musaceae genomes and enhances our knowledge of the pangenome.

A chromosome-level reference genome of Ensete glaucum gives insight into diversity and chromosomal and repetitive sequence evolution in the Musaceae.

BACKGROUND: Ensete glaucum (2n = 2x = 18) is a giant herbaceous monocotyledonous plant in the small Musaceae family along with banana (Musa). A high-quality reference genome sequence assembly of E. glaucum is a resource for functional and evolutionary studies of Ensete, Musaceae, and the Zingiberales. FINDINGS: Using Oxford Nanopore Technologies, chromosome conformation capture (Hi-C), Illumina and RNA survey sequence, supported by molecular cytogenetics, we report a high-quality 481.5 Mb genome assembly with 9 pseudo-chromosomes and 36,836 genes. A total of 55% of the genome is composed of repetitive sequences with predominantly LTR-retroelements (37%) and DNA transposons (7%). The single 5S ribosomal DNA locus had an exceptionally long monomer length of 1,056 bp, more than twice that of the monomers at multiple loci in Musa. A tandemly repeated satellite (1.1% of the genome, with no similar sequence in Musa) was present around all centromeres, together with a few copies of a long interspersed nuclear element (LINE) retroelement. The assembly enabled us to characterize in detail the chromosomal rearrangements occurring between E. glaucum and the x = 11 species of Musa. One E. glaucum chromosome has the same gene content as Musa acuminata, while others show multiple, complex, but clearly defined evolutionary rearrangements in the change between x= 9 and 11. CONCLUSIONS: The advance towards a Musaceae pangenome including E. glaucum, tolerant of extreme environments, makes a complete set of gene alleles, copy number variation, and a reference for structural variation available for crop breeding and understanding environmental responses. The chromosome-scale genome assembly shows the nature of chromosomal fusion and translocation events during speciation, and features of rapid repetitive DNA change in terms of copy number, sequence, and genomic location, critical to understanding its role in diversity and evolution.