The genome and population genomics of allopolyploid Coffea arabica reveal the diversification history of modern coffee cultivars.

Coffea arabica, an allotetraploid hybrid of Coffea eugenioides and Coffea canephora, is the source of approximately 60% of coffee products worldwide, and its cultivated accessions have undergone several population bottlenecks. We present chromosome-level assemblies of a di-haploid C. arabica accession and modern representatives of its diploid progenitors, C. eugenioides and C. canephora. The three species exhibit largely conserved genome structures between diploid parents and descendant subgenomes, with no obvious global subgenome dominance. We find evidence for a founding polyploidy event 350,000-610,000 years ago, followed by several pre-domestication bottlenecks, resulting in narrow genetic variation. A split between wild accessions and cultivar progenitors occurred ~30.5 thousand years ago, followed by a period of migration between the two populations. Analysis of modern varieties, including lines historically introgressed with C. canephora, highlights their breeding histories and loci that may contribute to pathogen resistance, laying the groundwork for future genomics-based breeding of C. arabica.

HERV-K (HML-2) insertion polymorphisms in the 8q24.13 region and their potential etiological associations with acute myeloid leukemia.

Human endogenous retroviruses (HERVs) are LTR retrotransposons that are present in the human genome. Among them, members of the HERV-K (HML-2) group are suspected to play a role in the development of different types of cancer, including lung, ovarian, and prostate cancer, as well as leukemia. Acute myeloid leukemia (AML) is an important disease that causes 1% of cancer deaths in the United States and has a survival rate of 28.7%. Here, we describe a method for assessing the statistical association between HERV-K (HML-2) transposable element insertion polymorphisms (or TIPs) and AML, using whole-genome sequencing and read mapping using TIP_finder software. Our results suggest that 101 polymorphisms involving HERV-K (HML-2) elements were correlated with AML, with a percentage between 44.4 to 56.6%, most of which (70) were located in the region from 8q24.13 to 8q24.21. Moreover, it was found that the TRIB1, LRATD2, POU5F1B, MYC, PCAT1, PVT1, and CCDC26 genes could be displaced or fragmented by TIPs. Furthermore, a general method was devised to facilitate analysis of the correlation between transposable element insertions and specific diseases. Finally, although the relationship between HERV-K (HML-2) TIPs and AML remains unclear, the data reported in this study indicate a statistical correlation, as supported by the χ2 test with p-values < 0.05.

A disease-associated mutation in thyroid hormone receptor α1 causes hearing loss and sensory hair cell patterning defects in mice.

Resistance to thyroid hormone due to mutations in THRA, which encodes the thyroid hormone receptor α (TRα1), shows variable clinical presentation. Mutations affecting TRß1 and TRß2 cause deafness in mice and have been associated with deafness in humans. To test whether TRα1 also affects hearing function, we used mice heterozygous for a frameshift mutation in Thra that is similar to human THRA mutations (ThraS1/+ mice) and reduces tissue sensitivity to thyroid hormone. Compared to wild-type littermates, ThraS1/+ mice showed moderate high-frequency sensorineural hearing loss as juveniles and increased age-related hearing loss. Ultrastructural examination revealed aberrant orientation of ~20% of sensory outer hair cells (OHCs), as well as increased numbers of mitochondria with fragmented morphology and autophagic vacuoles in both OHCs and auditory nerve fibers. Molecular dissection of the OHC lateral wall components revealed that the potassium ion channel Kcnq4 was aberrantly targeted to the cytoplasm of mutant OHCs. In addition, mutant cochleae showed increased oxidative stress, autophagy, and mitophagy associated with greater age-related cochlear cell damage, demonstrating that TRα1 is required for proper development of OHCs and for maintenance of OHC function. These findings suggest that patients with THRA mutations may present underdiagnosed, mild hearing loss and may be more susceptible to age-related hearing loss.

Erratum: In vitro assessment of pesticides capacity to act as agonists/antagonists of the thyroid hormone nuclear receptors.

[This corrects the article DOI: 10.1016/j.isci.2021.102957.].

TIP_finder: An HPC Software to Detect Transposable Element Insertion Polymorphisms in Large Genomic Datasets.

Transposable elements (TEs) are non-static genomic units capable of moving indistinctly from one chromosomal location to another. Their insertion polymorphisms may cause beneficial mutations, such as the creation of new gene function, or deleterious in eukaryotes, e.g., different types of cancer in humans. A particular type of TE called LTR-retrotransposons comprises almost 8% of the human genome. Among LTR retrotransposons, human endogenous retroviruses (HERVs) bear structural and functional similarities to retroviruses. Several tools allow the detection of transposon insertion polymorphisms (TIPs) but fail to efficiently analyze large genomes or large datasets. Here, we developed a computational tool, named TIP_finder, able to detect mobile element insertions in very large genomes, through high-performance computing (HPC) and parallel programming, using the inference of discordant read pair analysis. TIP_finder inputs are (i) short pair reads such as those obtained by Illumina, (ii) a chromosome-level reference genome sequence, and (iii) a database of consensus TE sequences. The HPC strategy we propose adds scalability and provides a useful tool to analyze huge genomic datasets in a decent running time. TIP_finder accelerates the detection of transposon insertion polymorphisms (TIPs) by up to 55 times in breast cancer datasets and 46 times in cancer-free datasets compared to the fastest available algorithms. TIP_finder applies a validated strategy to find TIPs, accelerates the process through HPC, and addresses the issues of runtime for large-scale analyses in the post-genomic era. TIP_finder version 1.0 is available at https://github.com/simonorozcoarias/TIP_finder.

Chloroplast genomes of Rubiaceae: Comparative genomics and molecular phylogeny in subfamily Ixoroideae.

In Rubiaceae phylogenetics, the number of markers often proved a limitation with authors failing to provide well-supported trees at tribal and generic levels. A robust phylogeny is a prerequisite to study the evolutionary patterns of traits at different taxonomic levels. Advances in next-generation sequencing technologies have revolutionized biology by providing, at reduced cost, huge amounts of data for an increased number of species. Due to their highly conserved structure, generally recombination-free, and mostly uniparental inheritance, chloroplast DNA sequences have long been used as choice markers for plant phylogeny reconstruction. The main objectives of this study are: 1) to gain insight in chloroplast genome evolution in the Rubiaceae (Ixoroideae) through efficient methodology for de novo assembly of plastid genomes; and, 2) to test the efficiency of mining SNPs in the nuclear genome of Ixoroideae based on the use of a coffee reference genome to produce well-supported nuclear trees. We assembled whole chloroplast genome sequences for 27 species of the Rubiaceae subfamily Ixoroideae using next-generation sequences. Analysis of the plastid genome structure reveals a relatively good conservation of gene content and order. Generally, low variation was observed between taxa in the boundary regions with the exception of the inverted repeat at both the large and short single copy junctions for some taxa. An average of 79% of the SNP determined in the Coffea genus are transferable to Ixoroideae, with variation ranging from 35% to 96%. In general, the plastid and the nuclear genome phylogenies are congruent with each other. They are well-resolved with well-supported branches. Generally, the tribes form well-identified clades but the tribe Sherbournieae is shown to be polyphyletic. The results are discussed relative to the methodology used and the chloroplast genome features in Rubiaceae and compared to previous Rubiaceae phylogenies.

Two Novel Cases of Resistance to Thyroid Hormone Due to THRA Mutation.

Resistance to thyroid hormone alpha (RTHα) is a rare and under-recognized genetic disease caused by mutations of THRA, the gene encoding thyroid hormone receptor α1 (TRα1). We report here two novel THRA missense mutations (M259T, T273A) in patients with RTHα. We combined biochemical and cellular assays with in silico modeling to assess the capacity of mutant TRα1 to bind triiodothyronine (T3), to heterodimerize with RXR, to interact with transcriptional coregulators, and to transduce a T3 transcriptional response. M259T, and to a lower extent T273A, reduces the affinity of TRα1 for T3. Their negative influence is only reverted by large excess of T3. The severity of the two novel RTHα cases originates from a reduction in the binding affinity of TRα1 mutants to T3 and thus correlates with the incapacity of corepressors to dissociate from TRα1 mutants in the presence of T3.

Evidence for hormonal control of heart regenerative capacity during endothermy acquisition.

Tissue regenerative potential displays striking divergence across phylogeny and ontogeny, but the underlying mechanisms remain enigmatic. Loss of mammalian cardiac regenerative potential correlates with cardiomyocyte cell-cycle arrest and polyploidization as well as the development of postnatal endothermy. We reveal that diploid cardiomyocyte abundance across 41 species conforms to Kleiber's law-the ¾-power law scaling of metabolism with bodyweight-and inversely correlates with standard metabolic rate, body temperature, and serum thyroxine level. Inactivation of thyroid hormone signaling reduces mouse cardiomyocyte polyploidization, delays cell-cycle exit, and retains cardiac regenerative potential in adults. Conversely, exogenous thyroid hormones inhibit zebrafish heart regeneration. Thus, our findings suggest that loss of heart regenerative capacity in adult mammals is triggered by increasing thyroid hormones and may be a trade-off for the acquisition of endothermy.

CRISPR/Cas9 Editing of the Mouse Thra Gene Produces Models with Variable Resistance to Thyroid Hormone.

BACKGROUND: Resistance to thyroid hormone due to THRA mutations (RTHα) is a recently discovered genetic disease, displaying important variability in its clinical presentation. The mutations alter the function of TRα1, one of the two nuclear receptors for thyroid hormone. METHODS: The aim of this study was to understand the relationship between specific THRA mutations and phenotype. CRISPR/Cas9 genome editing was used to generate five new mouse models of RTHα, with frameshift or missense mutations. RESULTS: Like human patients, mutant mice displayed a hypothyroid-like phenotype, with altered development. Phenotype severity varied between the different mouse models, mainly depending on the ability of the mutant receptor to interact with transcription corepressor in the presence of thyroid hormone. CONCLUSION: The present mutant mice represent highly relevant models for the human genetic disease which will be useful for future investigations.

Methylcytosine dioxygenase TET3 interacts with thyroid hormone nuclear receptors and stabilizes their association to chromatin.

Thyroid hormone receptors (TRs) are members of the nuclear hormone receptor superfamily that act as ligand-dependent transcription factors. Here we identified the ten-eleven translocation protein 3 (TET3) as a TR interacting protein increasing cell sensitivity to T3. The interaction between TET3 and TRs is independent of TET3 catalytic activity and specifically allows the stabilization of TRs on chromatin. We provide evidence that TET3 is required for TR stability, efficient binding of target genes, and transcriptional activation. Interestingly, the differential ability of different TRα1 mutants to interact with TET3 might explain their differential dominant activity in patients carrying TR germline mutations. So this study evidences a mode of action for TET3 as a nonclassical coregulator of TRs, modulating its stability and access to chromatin, rather than its intrinsic transcriptional activity. This regulatory function might be more general toward nuclear receptors. Indeed, TET3 interacts with different members of the superfamily and also enhances their association to chromatin.

Genome-wide analysis of LTR-retrotransposons in oil palm.

BACKGROUND: The oil palm (Elaeis guineensis Jacq.) is a major cultivated crop and the world's largest source of edible vegetable oil. The genus Elaeis comprises two species E. guineensis, the commercial African oil palm and E. oleifera, which is used in oil palm genetic breeding. The recent publication of both the African oil palm genome assembly and the first draft sequence of its Latin American relative now allows us to tackle the challenge of understanding the genome composition, structure and evolution of these palm genomes through the annotation of their repeated sequences. METHODS: In this study, we identified, annotated and compared Transposable Elements (TE) from the African and Latin American oil palms. In a first step, Transposable Element databases were built through de novo detection in both genome sequences then the TE content of both genomes was estimated. Then putative full-length retrotransposons with Long Terminal Repeats (LTRs) were further identified in the E. guineensis genome for characterization of their structural diversity, copy number and chromosomal distribution. Finally, their relative expression in several tissues was determined through in silico analysis of publicly available transcriptome data. RESULTS: Our results reveal a congruence in the transpositional history of LTR retrotransposons between E. oleifera and E. guineensis, especially the Sto-4 family. Also, we have identified and described 583 full-length LTR-retrotransposons in the Elaeis guineensis genome. Our work shows that these elements are most likely no longer mobile and that no recent insertion event has occurred. Moreover, the analysis of chromosomal distribution suggests a preferential insertion of Copia elements in gene-rich regions, whereas Gypsy elements appear to be evenly distributed throughout the genome. CONCLUSIONS: Considering the high proportion of LTR retrotransposon in the oil palm genome, our work will contribute to a greater understanding of their impact on genome organization and evolution. Moreover, the knowledge gained from this study constitutes a valuable resource for both the improvement of genome annotation and the investigation of the evolutionary history of palms.

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.

A Novel Mutation in THRA Gene Associated With an Atypical Phenotype of Resistance to Thyroid Hormone.

CONTEXT: RTHα is a recently discovered resistance to thyroid hormone (RTH) due to mutation of THRA, the gene encoding TRα1, the thyroid hormone receptor. It has been described in a few patients with growth retardation, short stature, and a low free T4/free T3 (FT4/FT3) ratio. OBJECTIVE: A 27-year-old patient presenting with dwarfism and a low FT4/FT3 ratio was investigated. DESIGN: Clinical, biochemical, and radiological data were collected. Whole exome sequencing was performed in the patient and her relatives. RESULTS: The patient exhibited congenital macrocytic anemia and severe bone malformation with growth retardation, dwarfism, clavicular agenesis, and abnormalities of the fingers, toes, and elbow joints. In adulthood, she presented with active behavior, chronic motor diarrhea, and hypercalcemia. Treatment with T3 led to heart rate acceleration, worsening of diarrhea, and TSH suppression. Low resting energy expenditure normalized on T3. rT3, SHBG, and IGF-1 remained normal. A de novo monoallelic missense mutation in THRA was discovered, the N359Y amino acid substitution (c.1075A>T), which affected both the TRα1 and the non-receptor isoform TRα2. The mutant TRα1 had a decrease in transcriptional activity related to decreased T3 binding and a dominant-negative effect on the wild-type receptor. CONCLUSIONS: This patient presents a new phenotype including more significant bone abnormalities, lower TSH, and higher FT3 levels, without certainty of all her symptoms with the TRα1(N359Y) mutation. This case suggests that patients with a low FT4/FT3 ratio should be screened for THRA mutations, even if clinical and biological features differ from previous reported cases of RTHα.

The coffee genome provides insight into the convergent evolution of caffeine biosynthesis.

Coffee is a valuable beverage crop due to its characteristic flavor, aroma, and the stimulating effects of caffeine. We generated a high-quality draft genome of the species Coffea canephora, which displays a conserved chromosomal gene order among asterid angiosperms. Although it shows no sign of the whole-genome triplication identified in Solanaceae species such as tomato, the genome includes several species-specific gene family expansions, among them N-methyltransferases (NMTs) involved in caffeine production, defense-related genes, and alkaloid and flavonoid enzymes involved in secondary compound synthesis. Comparative analyses of caffeine NMTs demonstrate that these genes expanded through sequential tandem duplications independently of genes from cacao and tea, suggesting that caffeine in eudicots is of polyphyletic origin.

Toxicogenomic analysis of the ability of brominated flame retardants TBBPA and BDE-209 to disrupt thyroid hormone signaling in neural cells.

Brominated flame retardants are suspected to act as disruptors of thyroid hormone signaling. This raises the concern that they might affect children's cognitive functions by influencing thyroid hormone signaling in the developing brain. We present here an in vitro analysis of the ability of the most common compounds, tetrabromobisphenol A (TBBPA) and BDE-209, to alter thyroid hormone response based on a model neural cell line and genome-wide analysis of gene expression.

Genome-wide search reveals the existence of a limited number of thyroid hormone receptor alpha target genes in cerebellar neurons.

Thyroid hormone (T3) has a major influence on cerebellum post-natal development. The major phenotypic landmark of exposure to low levels of T3 during development (hypothyroidism) in the cerebellum is the retarded inward migration of the most numerous cell type, granular neurons. In order to identify the direct genetic regulation exerted by T3 on cerebellar neurons and their precursors, we used microarray RNA hybridization to perform a time course analysis of T3 induced gene expression in primary cultures of cerebellar neuronal cell. These experiments suggest that we identified a small set of genes which are directly regulated, both in vivo and in vitro, during cerebellum post-natal development. These modest changes suggest that T3 does not acts directly on granular neurons and mainly indirectly influences the cellular interactions taking place during development.

Micro-collinearity and genome evolution in the vicinity of an ethylene receptor gene of cultivated diploid and allotetraploid coffee species (Coffea).

Arabica coffee (Coffea arabica L.) is a self-compatible perennial allotetraploid species (2n=4x=44), whereas Robusta coffee (C. canephora L.) is a self-incompatible perennial diploid species (2n=2x=22). C. arabica (C(a) C(a) E(a) E(a) ) is derived from a spontaneous hybridization between two closely related diploid coffee species, C. canephora (CC) and C. eugenioides (EE). To investigate the patterns and degree of DNA sequence divergence between the Arabica and Robusta coffee genomes, we identified orthologous bacterial artificial chromosomes (BACs) from C. arabica and C. canephora, and compared their sequences to trace their evolutionary history. Although a high level of sequence similarity was found between BACs from C. arabica and C. canephora, numerous chromosomal rearrangements were detected, including inversions, deletions and insertions. DNA sequence identity between C. arabica and C. canephora orthologous BACs ranged from 93.4% (between E(a) and C(a) ) to 94.6% (between C(a) and C). Analysis of eight orthologous gene pairs resulted in estimated ages of divergence between 0.046 and 0.665 million years, indicating a recent origin of the allotetraploid species C. arabica. Analysis of transposable elements revealed differential insertion events that contributed to the size increase in the C(a) sub-genome compared to its diploid relative. In particular, we showed that insertion of a Ty1-copia LTR retrotransposon occurred specifically in C. arabica, probably shortly after allopolyploid formation. The two sub-genomes of C. arabica, C(a) and E(a) , showed sufficient sequence differences, and a whole-genome shotgun approach could be suitable for sequencing the allotetraploid genome of C. arabica.

Life and death among plant lysophosphatidic acid acyltransferases.

The tetraploid Brassica napus possesses several seed-expressed microsomal lysophosphatidic acid acyltransferases (LPAAT ) including BAT1.5, which has been retained after genome duplication as a consequence of a subfunctionalisation of the gene encoding the ubiquitously expressed Kennedy pathway enzyme BAT1.13. Next, cDNA BAT1.3, encoding a LPAAT was subsequently isolated from an embryo library. The rapeseed LPAAT encoded by BAT1.3 is orthologous to the Arabidopsis thaliana At1g51260 gene product possibly associated with tapetum development and male fertility. However, BAT1.3 expression is predominant during the mid stages of embryo development in seeds of Brassica napus. Functional characterisation of BAT1.3 provides further support for a hypothesis of gene dosage sensitivity of LPAATs as does an analysis of the chromosomal localisation of LPAAT genes in Arabidopsis thaliana. The pattern of retention or loss of LPAAT genes after polyploidisation or segmental duplication is consistent with a model of balanced gene drive.

The transforming acidic coiled coil (TACC1) protein modulates the transcriptional activity of the nuclear receptors TR and RAR.

BACKGROUND: The transcriptional activity of Nuclear hormone Receptors (NRs) is regulated by interaction with coactivator or corepressor proteins. Many of these cofactors have been shown to have a misregulated expression or to show a subcellular mislocalization in cancer cell lines or primary tumors. Therefore they can be factors involved in the process of oncogenesis. RESULTS: We describe a novel NR coregulator, TACC1, which belongs to the Transforming Acidic Coiled Coil (TACC) family. The interaction of TACC1 with Thyroid Hormone Receptors (TR) and several other NRs has been shown in a yeast two-hybrid screen and confirmed by GST pulldown, colocalization and co-immunoprecipitation experiments. TACC1 interacts preferentially with unliganded NRs. In F9 cells, endogenous TACC1 localized in the chromatin-enriched fraction of the nucleus and interacted with Retinoid Acid Receptors (RARalpha) in the nucleus. TACC1 depletion in the cell led to decreased RARalpha and TRalpha ligand-dependent transcriptional activity and to delocalization of TR from the nucleus to the cytoplasm. CONCLUSIONS: From these experimental studies we propose that TACC1 might be a scaffold protein building up a transcriptional complex around the NRs we studied. This function of TACC1 might account for its involvement in several forms of tumour development.

Microcollinearity in an ethylene receptor coding gene region of the Coffea canephora genome is extensively conserved with Vitis vinifera and other distant dicotyledonous sequenced genomes.

BACKGROUND: Coffea canephora, also called Robusta, belongs to the Rubiaceae, the fourth largest angiosperm family. This diploid species (2x = 2n = 22) has a fairly small genome size of approximately 690 Mb and despite its extreme economic importance, particularly for developing countries, knowledge on the genome composition, structure and evolution remain very limited. Here, we report the 160 kb of the first C. canephora Bacterial Artificial Chromosome (BAC) clone ever sequenced and its fine analysis. RESULTS: This clone contains the CcEIN4 gene, encoding an ethylene receptor, and twenty other predicted genes showing a high gene density of one gene per 7.8 kb. Most of them display perfect matches with C. canephora expressed sequence tags or show transcriptional activities through PCR amplifications on cDNA libraries. Twenty-three transposable elements, mainly Class II transposon derivatives, were identified at this locus. Most of these Class II elements are Miniature Inverted-repeat Transposable Elements (MITE) known to be closely associated with plant genes. This BAC composition gives a pattern similar to those found in gene rich regions of Solanum lycopersicum and Medicago truncatula genomes indicating that the CcEIN4 regions may belong to a gene rich region in the C. canephora genome. Comparative sequence analysis indicated an extensive conservation between C. canephora and most of the reference dicotyledonous genomes studied in this work, such as tomato (S. lycopersicum), grapevine (V. vinifera), barrel medic M. truncatula, black cottonwood (Populus trichocarpa) and Arabidopsis thaliana. The higher degree of microcollinearity was found between C. canephora and V. vinifera, which belong respectively to the Asterids and Rosids, two clades that diverged more than 114 million years ago. CONCLUSION: This study provides a first glimpse of C. canephora genome composition and evolution. Our data revealed a remarkable conservation of the microcollinearity between C. canephora and V. vinifera and a high conservation with other distant dicotyledonous reference genomes. Altogether, these results provide valuable information to identify candidate genes in C. canephora genome and serve as a foundation to establish strategies for whole genome sequencing. Future large-scale sequence comparison between C. canephora and reference sequenced genomes will help in understanding the evolutionary history of dicotyledonous plants.