Evaluation of Algerian’s honey in terms of quality and authenticity based on the melissopalynology and physicochemical analysis and their antioxidant powers

BACKGROUND: Honey is a vegetable and animal product which comes from nectar and/or honeydew. It is used in different nutritional and therapeutic fields. OBJECTIVE: Melissopalynology and physicochemical analysis of Algerian honeys, determination of their phenolic compounds and authenticity parameters and the evaluation of their antioxidant properties. METHODS: Twenty Algerians honey were studied for their physicochemical parameters (moisture, pH, proteins, proline, hydroxymethylfurfural, ash, color, electrical conductivity, and optical rotation), floral origin and phenolic compounds contents. Antioxidant activities were tested too. RESULTS: Melissopalynologycal analyses revealed that the studied honeys were twelve multifloral, seven Fabaceae, and one Myrtaceae. All honeys were acidic (3.65≤pH≤4.35) and most of them were low in moisture content. The electrical conductivity varied between 0.29 mS/cm and 1.78 mS/cm. Ash, protein and proline contents results showed that the majority of honeys were in agreement with the legislation and were authentic. The color varied from mimosa yellow to dark brown. The specific rotation was levorotatory in most honey samples and the hydroxymethylfurfural values (from 1.5 mg/kg to 34.73 mg/kg) agreed with the international requirements. Honeys were rich in total phenolic compounds, 22.41 (Honey11) to 96.16 (Honey15)  mg gallic acid equivalents/100 g, and flavonoids, 8.90 (Honey11) to 80.02 (Honey02) mg quercetin equivalents/100 g. Honey samples 15, 03, 05, 01, and 06 exerted more than 50% reduction of 1,1-diphenyl-2-picrylhydrazyl and 2,2′-azino-bis(3-ethylbenzothiazoline-6-sulfonate) radicals and were able to reduce iron while honey samples 12, 18, 19, 14, and 11 chelate efficiently iron. High significant correlations between physicochemical parameters and antioxidant activities were found. CONCLUSION: The Algerian honeys analyzed were authentic and variations in their quality parameters and phenolics composition were directly associated with their demonstrated antioxidant properties.

Early diversification and permeable species boundaries in the Mediterranean firs

BACKGROUND AND AIMS

Inferring the evolutionary relationships of species and their boundaries is critical in order to understand patterns of diversification and their historical drivers. Despite Abies (Pinaceae) being the second most diverse group of conifers, the evolutionary history of Circum-Mediterranean firs (CMFs) remains under debate.

METHODS

We used restriction site-associated DNA sequencing (RAD-seq) on all proposed CMF taxa to investigate their phylogenetic relationships and taxonomic status.

KEY RESULTS

Based on thousands of genome-wide single nucleotide polymorphisms (SNPs), we present here the first formal test of species delimitation, and the first fully resolved, complete species tree for CMFs. We discovered that all previously recognized taxa in the Mediterranean should be treated as independent species, with the exception of Abies tazaotana and Abies marocana. An unexpectedly early pulse of speciation in the Oligocene-Miocene boundary is here documented for the group, pre-dating previous hypotheses by millions of years, revealing a complex evolutionary history encompassing both ancient and recent gene flow between distant lineages.

CONCLUSIONS

Our phylogenomic results contribute to shed light on conifers' diversification. Our efforts to resolve the CMF phylogenetic relationships help refine their taxonomy and our knowledge of their evolution.

Diversity of xerotolerant and xerophilic fungi in honey

Fungi can colonize most of the substrata on Earth. Honey, a sugary food produced by bees (and other insects) has been studied little in terms of its fungal diversity. We have surveyed and evaluated the presence of xerotolerant and xerophilic fungi in a set of honey bee samples collected from across Spain. From 84 samples, a total of 104 fungal strains were isolated, and morphologically and phylogenetically characterized. We identified 32 species distributed across 16 genera, most of them belonging to the ascomycetous genera Aspergillus, Bettsia, Candida, Eremascus, Monascus, Oidiodendron, Penicillium, Skoua, Talaromyces and Zygosaccharomyces. As a result of this survey, eight new taxa are proposed: i.e. the new family Helicoarthrosporaceae, two new genera, Helicoarthrosporum and Strongyloarthrosporum in Onygenales; three new species of Eurotiales, Talaromyces affinitatimellis, T. basipetosporus, and T. brunneosporus; and two new species of Myxotrichaceae, Oidiodendron mellicola, and Skoua asexualis.

Analysis of Multifloral Bee Pollen Pellets by Advanced Digital Imaging Applied to Functional Food Ingredients

Bee pollen is a hive product, resulting from floral pollen agglutination by worker bees and it is characterized by its excellent bioactive and nutritional composition. Currently, research is focused on bee pollen applications on food industry, because this product has been considered an excellent source of compounds for human nutrition. It is also important in some industries, where color and particle size are important characteristics for production. Due to the granular nature of bee pollen, conventional colorimetry does not allow describing color correctly; thus, digital image analysis is a better alternative. This technique could also allow classifying bee pollen according to its appearance beyond the color. Consequently, the aim of this work was to develop a novel methodology for image data processing to classify bee pollen as ingredient in food industry. Seven color groups in samples were established regarding harvest month and particle size. It was possible to calculate the percentage of each color group in all samples. This methodology also allowed selecting each fraction for different applications in food industry using colorimetry, granulometry and the relationship between both of them.

MICROBIOLOGICAL QUALITY OF MOROCCAN LABELED EUPHORBIA RESINIFERA HONEY

In the present work, microbiological profile of thirty-seven samples of labeled honey were collected in a Protected Geographical Indication “PGI” area of Tadla-Azilal region, which is an endemic zone of Euphorbia resinifera plant. A profile was assessed using conventional microbial methods, like enumeration, detection and/or germs identification, in accordance with ISO norms. This is the first study in which a honey with Moroccan “PGI” was tested, in order to assess its compliance with bacteriological recommendations. Coliforms (Total and fecal Coliforms), Salmonella spp., Shigella spp., Sporus of Bacillus cereus and Clostridium perfringens were not detected. The numbers of Standard Plate Count “SPC” were less than 102 CFU.g-1 for all samples. The molds and yeasts were found among samples and 32% and 40% of samples were positive, respectively. However, no samples showed a higher value than recommended limit [102 CFU.g-1]. We conclude that samples of labeled euphorbia honey of Tadla-Azilal analyzed present good commercial quality parameters (SPC, molds and yeasts “absence of unwanted fermentations”), a good sanitary quality (absence of coliforms and S. aureus) and are safe (Slam., Shig., Sporus of B. cereus and C. perf.). Standardization (regulation and specifications) and a rationalization of beekeeping techniques throughout Euphorbia “PGI” area studied may further sustainably improve the quality of this unique honey, and ensure it over the years.

Phylogeography above the species level for perennial species in a composite genus

In phylogeography, DNA sequence and fingerprint data at the population level are used to infer evolutionary histories of species. Phylogeography above the species level is concerned with the genealogical aspects of divergent lineages. Here, we present a phylogeographic study to examine the evolutionary history of a western Mediterranean composite, focusing on the perennial species of Helminthotheca (Asteraceae, Cichorieae). We used molecular markers (amplified fragment length polymorphism (AFLP), internal transcribed spacer and plastid DNA sequences) to infer relationships among populations throughout the distributional range of the group. Interpretation is aided by biogeographic and molecular clock analyses. Four coherent entities are revealed by Bayesian mixture clustering of AFLP data, which correspond to taxa previously recognized at the rank of subspecies. The origin of the group was in western North Africa, from where it expanded across the Strait of Gibraltar to the Iberian Peninsula and across the Strait of Sicily to Sicily. Pleistocene lineage divergence is inferred within western North Africa as well as within the western Iberian region. The existence of the four entities as discrete evolutionary lineages suggests that they should be elevated to the rank of species, yielding H. aculeata, H. comosa, H. maroccana and H. spinosa, whereby the latter two necessitate new combinations.

Big thistle eats the little thistle: does unidirectional introgressive hybridization endanger the conservation of Onopordum hinojense?

Hybridization is known to have a creative role in plant evolution. However, it can also have negative effects on parental species. Onopordum is a large genus whose species frequently hybridize. In the Southwest Iberian Peninsula, the rare O. hinojense co-occurs with the widely distributed O. nervosum, and hybrids between these two taxa have been described as O. × onubense. In this study we determine the extinction risk in a hybrid zone, both for hybrids and parentals, using analyses of morphological and cytogenetic traits as well as genetic markers and demographic models. To investigate the introgression process we used amplified fragment length polymorphism (AFLP) markers, Bayesian analyses and genome scan methods. Morphology, genome size and molecular markers confirmed homoploid hybridization and also indicated unidirectional backcrossing of F₁ hybrids with O. nervosum, which is likely to swamp O. hinojense, the parental with lower pollen size and a very low fruit set (8%). Genome scan methods revealed several loci significantly deviating from neutrality. Finally, our demographic modeling indicated that the higher fitness of O. nervosum threats the survival of O. hinojense by demographic swamping. Our study provides strong new evidence for a scenario of rapid extinction by unidirectional introgression and demographic swamping. The multifaceted approach used here sheds new light on the role of introgression in plant extinctions.

Effects of tree architecture on pollen dispersal and mating patterns in Abies pinsapo Boiss. (Pinaceae)

Plant architecture is crucial to pollination and mating in wind-pollinated species. We investigated the effect of crown architecture on pollen dispersal, mating system and offspring quality, combining phenotypic and genotypic analyses in a low-density population of the endangered species Abies pinsapo. A total of 598 embryos from three relative crown height levels (bottom, middle and top) in five mother plants were genotyped using eleven nuclear microsatellite markers (nSSRs). Paternity analysis and mating system models were used to infer mating and pollen dispersal parameters. In addition, seeds were weighed (N = 16 110) and germinated (N = 736), and seedling vigour was measured to assess inbreeding depression. Overall, A. pinsapo shows a fat-tailed dispersal kernel, with an average pollen dispersal distance of 113-227 m, an immigration rate of 0.84-26.92%, and a number of effective pollen donors (Nep ) ranging between 3.5 and 11.9. We found an effect of tree height and relative crown height levels on mating parameters. A higher proportion of seeds with embryo (about 50%) and a higher rate of self-fertilization (about 60%) were found at the bottom level in comparison with the top level. Seed weight and seedling vigour are positively related. Nevertheless, no differences were found in seed weight or in seedling-related variables such as weight and length of aerial and subterranean parts among the different relative crown height levels, suggesting that seeds from the more strongly inbred bottom level are not affected by inbreeding depression. Our results point to vertical isotropy for outcross-pollen and they suggest that self-pollen may ensure fertilization when outcross-pollen is not available in low-density population.

Nuclear microsatellite primers for the endangered relict fir, Abies pinsapo (Pinaceae) and cross-amplification in related Mediterranean Species

Twelve nuclear microsatellite primers (nSSR) were developed for the endangered species Abies pinsapo Boiss. to enable the study of gene flow and genetic structure in the remaining distribution areas. Microsatellite primers were developed using next-generation sequencing (454) data from a single Abies pinsapo individual. Primers were applied to thirty individuals from the three extant localities. The number of alleles per locus ranged from one to four. Cross-amplification was tested for other Abies species from the Mediterranean Basin, and most of the loci showed higher polymorphisms in the Mediterranean species than in A. pinsapo. These microsatellite markers provide tools for conservation genetic studies in Abies pinsapo as well other Abies species from the Mediterranean Basin.

Molecular phylogeny and systematics of the highly polymorphic Rumex bucephalophorus complex (Polygonaceae)

Rumex bucephalophorus is a very polymorphic species that has been subjected to various taxonomic studies in which diverse infraspecific taxa have been recognised on the basis of diaspore traits. In this study we used molecular markers (ITS and AFLP) to explore this remarkable diversity, to test previous hypotheses of classification, and attempt to explain biogeographic patterns. Results show that R. bucephalophorus forms a monophyletic group in which diversification began around 4.2 Mya, at the end of Messinian Salinity Crisis. The two molecular markers clearly show a deep divergence separating subsp. bucephalophorus from all other subspecific taxa, among which subsp. canariensis also constitutes a separate and well distinguishable unit. In contrast, subspecies hispanicus and subsp. gallicus constitute a monophyletic group in which three subgroups can be recognised: subsp. hispanicus, subsp. gallicus var. gallicus and subsp. gallicus var. subaegeus. However, these three subgroups are not clearly distinguished genetically or morphologically, so that in formal classification it would be preferable to treat them at the varietal level.

The Strait of Gibraltar as a major biogeographic barrier in Mediterranean conifers: a comparative phylogeographic survey

The Strait of Gibraltar (SG) is reputed for being both a bridge and a geographic barrier to biological exchanges between Europe and Africa. Major genetic breaks associated with this strait have been identified in various taxa, but it is unknown whether these disjunctions have been produced simultaneously or by independent biogeographic processes. Here, the genetic structure of five conifers distributed on both sides of the SG was investigated using mitochondrial (nad1 b/c, nad5-1, nad5-4 and nad7-1) and chloroplast (Pt1254, Pt15169, Pt30204, Pt36480, Pt71936 and Pt87268) DNA markers. The distribution of genetic variation was partially congruent between types of markers within the same species. Across taxa, there was a significant overlapping between the SG and the genetic breaks detected, especially for the four Tertiary species surveyed (Abies pinsapo complex, Pinus nigra, Pinus pinaster and Taxus baccata). For most of these taxa, the divergence of populations across the SG could date back to long before the Pleistocene glaciations. However, their strongly different cpDNA G(ST) and R(ST) values point out that they have had dissimilar population histories, which might include contrasting amounts of pollen-driven gene flow since their initial establishment in the region. The fifth species, Pinus halepensis, was genetically depauperated and homogenous on both sides of the SG. A further analysis of nuclear DNA sequences with coalescent-based isolation with migration models suggests a Pleistocene divergence of P. halepensis populations across the SG, which is in sharp contrast with the pre-Pleistocene divergence dates obtained for P. pinaster. Altogether, these results indicate that the genetic breaks observed across this putative biogeographical barrier have been produced by independent evolutionary processes related to the biological history of each individual species instead of a common vicariant phenomenon.

Pleistocene refugia and polytopic replacement of diploids by tetraploids in the Patagonian and Subantarctic plant Hypochaeris incana (Asteraceae, Cichorieae)

We report the phylogeographic pattern of the Patagonian and Subantarctic plant Hypochaeris incana endemic to southeastern South America. We applied amplified fragment length polymorphism (AFLP) and chloroplast DNA (cpDNA) analysis to 28 and 32 populations, respectively, throughout its distributional range and assessed ploidy levels using flow cytometry. While cpDNA data suggest repeated or simultaneous parallel colonization of Patagonia and Tierra del Fuego by several haplotypes and/or hybridization, AFLPs reveal three clusters corresponding to geographic regions. The central and northern Patagonian clusters (approximately 38-51 degrees S), which are closer to the outgroup, contain mainly tetraploid, isolated and highly differentiated populations with low genetic diversity. To the contrary, the southern Patagonian and Fuegian cluster (approximately 51-55 degrees S) contains mainly diploid populations with high genetic diversity and connected by high levels of gene flow. The data suggest that H. incana originated at the diploid level in central or northern Patagonia, from where it migrated south. All three areas, northern, central and southern, have similar levels of rare and private AFLP bands, suggesting that all three served as refugia for H. incana during glacial times. In southern Patagonia and Tierra del Fuego, the species seems to have expanded its populational system in postglacial times, when the climate became warmer and more humid. In central and northern Patagonia, the populations seem to have become restricted to favourable sites with increasing temperature and decreasing moisture and there was a parallel replacement of diploids by tetraploids in local populations.

Phylogeography of North African Atlas cedar (Cedrus atlantica, Pinaceae): Combined molecular and fossil data reveal a complex Quaternary history

Northwest Africa is a major hotspot of plant biodiversity, but very little is known about the Quaternary range dynamics of plant species in this region. Here we investigate the range-wide population structure and phylogeography of Atlas cedar (Cedrus atlantica), an emblematic forest tree endemic to Morocco and Algeria. We genotyped 261 individuals from 11 populations using AFLP markers. Data were analyzed using both conventional F(ST)-based techniques and Bayesian clustering. Overall population differentiation was high (F(ST) = 0.25). Two major groups of populations were identified, one distributed through the Rif and Middle Atlas mountains in Morocco and the other through the Algerian Tell Atlas and Aurès mountains as well as the Middle Atlas. Combined molecular and fossil data indicate that C. atlantica survived the Last Glacial Maximum in at least three disjunct refugia along the coast of the Mediterranean Sea, whereas the Middle Atlas, today the core of the species range, has been colonized relatively recently (<10000 yr BP). The colonization history of individual populations has left clear imprints in their present-day diversity, which may vary greatly even between nearby stands. Our study illustrates how integrating different data sources and analytical approaches can help elucidate complex range dynamics that would otherwise remain undeciphered.

Phylogeography of the invasive weed Hypochaeris radicata (Asteraceae): from Moroccan origin to worldwide introduced populations

In an attempt to delineate the area of origin and migratory expansion of the highly successful invasive weedy species Hypochaeris radicata, we analysed amplified fragment length polymorphisms from samples taken from 44 populations. Population sampling focused on the central and western Mediterranean area, but also included sites from Northern Spain, Western and Central Europe, Southeast Asia and South America. The six primer combinations applied to 213 individuals generated a total of 517 fragments of which 513 (99.2%) were polymorphic. The neighbour-joining tree presented five clusters and these divisions were supported by the results of Bayesian analyses: plants in the Moroccan, Betic Sierras (Southern Spain), and central Mediterranean clusters are all heterocarpic. The north and central Spanish, southwestern Sierra Morena, and Central European, Asian and South American cluster contain both heterocarpic (southwestern Sierra Morena) and homocarpic populations (all other populations). The Doñana cluster includes two homocarpic populations. Analyses of fragment parameters indicate that the oldest populations of H. radicata are located in Morocco and that the species expanded from this area in the Late Quaternary via at least three migratory routes, the earliest of which seems to have been to the southwestern Iberian Peninsula, with subsequent colonizations to the central Mediterranean area and the Betic Sierras. Homocarpic populations originated in the southwestern Iberian Peninsula and subsequently spread across north and central Spain, Central Europe and worldwide, where they became a highly successful weed.

Range-wide phylogeography of Juniperus thurifera L., a presumptive keystone species of western Mediterranean vegetation during cold stages of the Pleistocene

We investigate the range-wide population structure and phylogeography of thuriferous juniper (Juniperus thurifera L.), a species with a highly disjunct distribution in the western Mediterranean. We genotyped a total of 327 individuals from 20 populations using amplified fragment length polymorphisms (AFLP). Different analyses such as principal co-ordinate analysis (PCoA), nonmetric multidimensional scaling of F(ST) distances among populations, unweighted pair group method with arithmetic mean (UPGMA), and Bayesian clustering revealed that the Strait of Gibraltar acted as an efficient barrier against gene flow between the Moroccan and European populations for a very long time, and consequently support that the Moroccan populations should be recognised as a distinct subspecies (J. thurifera L. subsp. africana (Maire) Romo and Boratyńsky). The Algerian population was genetically more closely related to the European than to the Moroccan ones, probably due to dispersal events from Europe to Algeria. With respect to the mainland European populations, our data are not conclusive to reject any of the two following hypotheses: (1) the Iberian Peninsula was subdivided into different gene pools, and was the source for the colonisation of the Pyrenees and the Alps; and (2) the pattern we see today is partly the result of immigration into the Iberian Peninsula, e.g. from the Alps. Finally, the Corsican population was closely related genetically to two northern Iberian populations most probably due to relatively recent long-distance dispersal.

Genetic diversity and population structure in natural populations of Moroccan Atlas cedar (Cedrus atlantica; Pinaceae) determined with cpSSR markers

Atlas cedar (Cedrus atlantica) is an ecologically and economically important forest tree species of northern Africa and is considered one of the endangered conifer species in the region. Chloroplast microsatellites (cpSSR) were used to study genetic variation within and among populations and geographical structure in natural populations of C. atlantica throughout its entire distribution range in Morocco. A total of 25 chloroplast haplotypes and 66 cpSSR alleles were found among 162 individuals. The cpSSRs indicate that C. atlantica appears to maintain a high level of genetic diversity (mean H(e) = 0.95), as observed in most coniferous species. Values of mean pairwise distance within a population (D(2)(SH)) were related to the size and location of the populations. AMOVA analysis showed that most of the variation in C. atlantica occurs within populations and confirmed the general tendency of gymnosperms to display lower values of population differentiation than angiosperms. The distance-based clustering method (PCoA and neighbor-joining analysis) and the geographical structure revealed a poor structure among the six populations of Cedrus atlantica. Also, a Mantel test indicated a weak correlation between geographic and genetic distances (P = 0.106, r = 0.363). These results are also interpreted in the context of postglacial history of the region plus human impacts.

Multivariate correlation between color and mineral composition of honeys and by their botanical origin

The mineral content and color characteristics of 77 honey samples were analyzed. Eighteen minerals were quantified for each honey. Multiple linear regression (MLR) was used to establish equations relating the colorimetric CIELAB coordinates to the mineral data. The results obtained shown that lightness (L) was significantly correlated with S, Ca, Fe, As, Pb, and Cd for the dark honey types (avocado, heather, chestnut, and honeydew). For the light and brown honey types (citrus, rosemary, lavender, eucalyptus, and thyme), C(ab) and b showed the lower correlation with the mineral content of the honeys; their regression functions involve a few independent variables (Mg and Al for b and only Al for C(ab)). Furthermore, by means of application of linear discriminant analysis to the mineral content, it was possible to obtain a model that classifies the honeys by their lightness. The prediction ability of the built model, determined with the test set method, was 85%.

Inductively coupled plasma optical emission spectrometric determination of minerals in thyme honeys and their contribution to geographical discrimination

Twenty-four Spanish thyme honey samples were analyzed using inductively coupled plasma optical emission spectrometry (ICP-OES). Twenty-four minerals were quantified for each honey. The elements Al, As, Ba, Ca, Cr, Cu, Fe, K, Li, Mg, Mn, Na, P, Pb, S, Se, Si, Sr, and Zn were detected in all samples; seven elements are very abundant (Ca, K, Mg, Na, P, S, and Si), and six are not abundant (Al, Cu, Fe, Li, Mn, and Zn). Eleven of them are trace elements (As, Ba, Cd, Co, Cr, Ni, Mo, Pb, Se, Sr, and V) at <1 mg kg(-)(1). Classification of thyme honeys according to their origin (coast, mountains) was achieved by pattern recognition techniques on the mineral data. By means of principal component analysis, a good separation by geographical origin is obtained when scores for the two first principal components are plotted. Classification functions of 11 metals (Al, As, Cr, Cu, K, Li, Mg, Na, P, S, and V) were obtained using stepwise discriminant analysis and applied to classify correctly approximately 100% of the honey samples.

Análisis polínico de mieles en las regiones de Ottazzane y Costa Atlántica (Noroeste de Marruecos)

RESUMEN. Análisis polínico de mieles en las regiones de Ottazzane y Costa Atlántica (Noroeste de Marruecos). Se ha realizado el análisis microscópico de 13 muestras de miel de las regiones de Ouazzane y Costa Atlántica. Las muestras fueron proporcionadas directamente por los apicultores, en su mayoría aficionados. Los resultados reflejan que el néctar de las flores es la principal fuente de miel en el territorio y que siete de las muestras son pobres polinicamente, con 4.600-47.800 GP (Clase I y II de Maurizio), dos muestras presentan una riqueza media, con 189.000-209.700 GP (Clase III), y cuatro son ricas o muy ricas, con 872.000-2.950.000 GP (Clases IV y V). Se han identificado 63 taxones por el análisis microscópico, resultando seis de las mieles monoflorales: dos de Eucalyptus sp., dos de Lythrum sp., una de Leucojum sp. y una de Citrus sp.Palabras clave. Marruecos, Ouazzane, Costa Atlántica, polen, miel, monofloral, Eucalyptus sp., Lythrum sp., Leucojunz sp.,Citrus sp.ABSTRACT. Pollen analysis of honey from Ouazzane and Athlantic coast of Morocco. Thirteen honey samples from Ouazzane and Atlantic Coast are analysed. The samples were directly provided by amateur beekeepers. A total of 63 plants have been identified by microscopic analysis. Results show that nectar is the main source for honey in this territory and that seven samples have pollen gathering poor, with 4.600- 47.800 PG (Maurizio's Classes 1 and II), two samples have pollen gathering medium, with 189.000-209.000 PG (Class III), and four have pollen gathering high, with 872.000-2.950.000 PG (Classes IV and V). By microscopic analysis 63 taxa were identified and six of the honey samples are monofloral: two of Eucalyptus sp., two of Lythrum sp., one of Leucojum sp. and one of Citrus sp.Key words. Morocco, Ouazzane, Atlantic Coast, pollen, honey, unifloral, Eucalyptus sp., Lythrum sp., Leucojunz sp.,Citrus sp.