Single cell tracing of Pomc neurons reveals recruitment of 'Ghost' subtypes with atypical identity in a mouse model of obesity

The hypothalamus contains a remarkable diversity of neurons that orchestrate behavioural and metabolic outputs in a highly plastic manner. Neuronal diversity is key to enabling hypothalamic functions and, according to the neuroscience dogma, it is predetermined during embryonic life. Here, by combining lineage tracing of hypothalamic pro-opiomelanocortin (Pomc) neurons with single-cell profiling approaches in adult male mice, we uncovered subpopulations of 'Ghost' neurons endowed with atypical molecular and functional identity. Compared to 'classical' Pomc neurons, Ghost neurons exhibit negligible Pomc expression and are 'invisible' to available neuroanatomical approaches and promoter-based reporter mice for studying Pomc biology. Ghost neuron numbers augment in diet-induced obese mice, independent of neurogenesis or cell death, but weight loss can reverse this shift. Our work challenges the notion of fixed, developmentally programmed neuronal identities in the mature hypothalamus and highlight the ability of specialised neurons to reversibly adapt their functional identity to adult-onset obesogenic stimuli.

Single-cell RNA sequencing identifies senescence as therapeutic target in rhabdomyolysis-induced acute kidney injury

BACKGROUND

The role of macrophages in the development of rhabdomyolysis-induced acute kidney injury (RM-AKI) has been established, but an in-depth understanding of the changes in the immune landscape could help to improve targeted strategies. Whereas senescence is usually associated with chronic kidney processes, we also wished to explore whether senescence could also occur in AKI and whether senolytics could act on immune cells.

METHODS

Single-cell RNA sequencing was used in the murine glycerol-induced RM-AKI model to dissect the transcriptomic characteristics of CD45+ live cells sorted from kidneys 2 days after injury. Public datasets from murine AKI models were reanalysed to explore cellular senescence signature in tubular epithelial cells (TECs). A combination of senolytics (dasatinib and quercetin, DQ) was administered to mice exposed or not to RM-AKI.

RESULTS

Unsupervised clustering of nearly 17 000 single-cell transcriptomes identified seven known immune cell clusters. Sub-clustering of the mononuclear phagocyte cells revealed nine distinct cell sub-populations differently modified with RM. One macrophage cluster was particularly interesting since it behaved as a critical node in a trajectory connecting one major histocompatibility complex class IIhigh (MHCIIhigh) cluster only present in Control to two MHCIIlow clusters only present in RM-AKI. This critical cluster expressed a senescence gene signature, that was very different from that of the TECs. Senolytic DQ treatment blocked the switch from a F4/80highCD11blow to F4/80lowCD11bhigh phenotype, which correlated with prolonged nephroprotection in RM-AKI.

CONCLUSIONS

Single-cell RNA sequencing unmasked novel transitional macrophage subpopulation associated with RM-AKI characterized by the activation of cellular senescence processes. This work provides a proof-of-concept that senolytics nephroprotective effects may rely, at least in part, on subtle immune modulation.

Eomes-dependent mitochondrial regulation promotes survival of pathogenic CD4+ T cells during inflammation

The mechanisms whereby Eomes controls tissue accumulation of T cells and strengthens inflammation remain ill-defined. Here, we show that Eomes deletion in antigen-specific CD4+ T cells is sufficient to protect against central nervous system (CNS) inflammation. While Eomes is dispensable for the initial priming of CD4+ T cells, it is required for long-term maintenance of CNS-infiltrating CD4+ T cells. We reveal that the impact of Eomes on effector CD4+ T cell longevity is associated with sustained expression of multiple genes involved in mitochondrial organization and functions. Accordingly, epigenetic studies demonstrate that Eomes supports mitochondrial function by direct binding to either metabolism-associated genes or mitochondrial transcriptional modulators. Besides, the significance of these findings was confirmed in CD4+ T cells from healthy donors and multiple sclerosis patients. Together, our data reveal a new mechanism by which Eomes promotes severity and chronicity of inflammation via the enhancement of CD4+ T cell mitochondrial functions and resistance to stress-induced cell death.

Macromolecular Complex Including MLL3, Carabin and Calcineurin Regulates Cardiac Remodeling

BACKGROUND

Cardiac hypertrophy is an intermediate stage in the development of heart failure. The structural and functional processes occurring in cardiac hypertrophy include extensive gene reprogramming, which is dependent on epigenetic regulation and chromatin remodeling. However, the chromatin remodelers and their regulatory functions involved in the pathogenesis of cardiac hypertrophy are not well characterized.

METHODS

Protein interaction was determined by immunoprecipitation assay in primary cardiomyocytes and mouse cardiac samples subjected or not to transverse aortic constriction for 1 week. Chromatin immunoprecipitation and DNA sequencing (ChIP-seq) experiments were performed on the chromatin of adult mouse cardiomyocytes.

RESULTS

We report that the calcium-activated protein phosphatase CaN (calcineurin), its endogenous inhibitory protein carabin, the STK24 (STE20-like protein kinase 3), and the histone monomethyltransferase, MLL3 (mixed lineage leukemia 3) form altogether a macromolecular complex at the chromatin of cardiomyocytes. Under basal conditions, carabin prevents CaN activation while the serine/threonine kinase STK24 maintains MLL3 inactive via phosphorylation. After 1 week of transverse aortic constriction, both carabin and STK24 are released from the CaN-MLL3 complex leading to the activation of CaN, dephosphorylation of MLL3, and in turn, histone H3 lysine 4 monomethylation. Selective cardiac MLL3 knockdown mitigates hypertrophy, and chromatin immunoprecipitation and DNA sequencing analysis demonstrates that MLL3 is de novo recruited at the transcriptional start site of genes implicated in cardiomyopathy in stress conditions. We also show that CaN and MLL3 colocalize at chromatin and that CaN activates MLL3 histone methyl transferase activity at distal intergenic regions under hypertrophic conditions.

CONCLUSIONS

Our study reveals an unsuspected epigenetic mechanism of CaN that directly regulates MLL3 histone methyl transferase activity to promote cardiac remodeling.

The proteome and transcriptome of stress granules and P bodies during human T lymphocyte activation

Stress granules (SGs) and processing bodies (PBs) are membraneless cytoplasmic assemblies regulating mRNAs under environmental stress such as viral infections, neurological disorders, or cancer. Upon antigen stimulation, T lymphocytes mediate their immune functions under regulatory mechanisms involving SGs and PBs. However, the impact of T cell activation on such complexes in terms of formation, constitution, and relationship remains unknown. Here, by combining proteomic, transcriptomic, and immunofluorescence approaches, we simultaneously characterized the SGs and PBs from primary human T lymphocytes pre and post stimulation. The identification of the proteomes and transcriptomes of SGs and PBs indicate an unanticipated molecular and functional complementarity. Notwithstanding, these granules keep distinct spatial organizations and abilities to interact with mRNAs. This comprehensive characterization of the RNP granule proteomic and transcriptomic landscapes provides a unique resource for future investigations on SGs and PBs in T lymphocytes.

Tamoxifen Accelerates Endothelial Healing by Targeting ERα in Smooth Muscle Cells

RATIONALE

Tamoxifen prevents the recurrence of breast cancer and is also beneficial against bone demineralization and arterial diseases. It acts as an ER (estrogen receptor) α antagonist in ER-positive breast cancers, whereas it mimics the protective action of 17β-estradiol in other tissues such as arteries. However, the mechanisms of these tissue-specific actions remain unclear.

OBJECTIVE

Here, we tested whether tamoxifen is able to accelerate endothelial healing and analyzed the underlying mechanisms.

METHODS AND RESULTS

Using 3 complementary mouse models of carotid artery injury, we demonstrated that both tamoxifen and estradiol accelerated endothelial healing, but only tamoxifen required the presence of the underlying medial smooth muscle cells. Chronic treatment with 17β-estradiol and tamoxifen elicited differential gene expression profiles in the carotid artery. The use of transgenic mouse models targeting either whole ERα in a cell-specific manner or ERα subfunctions (membrane/extranuclear versus genomic/transcriptional) demonstrated that 17β-estradiol-induced acceleration of endothelial healing is mediated by membrane ERα in endothelial cells, while the effect of tamoxifen is mediated by the nuclear actions of ERα in smooth muscle cells.

CONCLUSIONS

Whereas tamoxifen acts as an antiestrogen and ERα antagonist in breast cancer but also on the membrane ERα of endothelial cells, it accelerates endothelial healing through activation of nuclear ERα in smooth muscle cells, inviting to revisit the mechanisms of action of selective modulation of ERα.

Blood CD9+ B cell, a biomarker of bronchiolitis obliterans syndrome after lung transplantation

Bronchiolitis obliterans syndrome is the main limitation for long-term survival after lung transplantation. Some specific B cell populations are associated with long-term graft acceptance. We aimed to monitor the B cell profile during early development of bronchiolitis obliterans syndrome after lung transplantation. The B cell longitudinal profile was analyzed in peripheral blood mononuclear cells from patients with bronchiolitis obliterans syndrome and patients who remained stable over 3 years of follow-up. CD24^hi CD38^hi transitional B cells were increased in stable patients only, and reached a peak 24 months after transplantation, whereas they remained unchanged in patients who developed a bronchiolitis obliterans syndrome. These CD24^hi CD38^hi transitional B cells specifically secrete IL-10 and express CD9. Thus, patients with a total CD9⁺ B cell frequency below 6.6% displayed significantly higher incidence of bronchiolitis obliterans syndrome (AUC = 0.836, PPV = 0.75, NPV = 1). These data are the first to associate IL-10-secreting CD24^hi CD38^hi transitional B cells expressing CD9 with better allograft outcome in lung transplant recipients. CD9-expressing B cells appear as a contributor to a favorable environment essential for the maintenance of long-term stable graft function and as a new predictive biomarker of bronchiolitis obliterans syndrome-free survival.

Draft genome and reference transcriptomic resources for the urticating pine defoliator Thaumetopoea pityocampa (Lepidoptera: Notodontidae)

The pine processionary moth Thaumetopoea pityocampa (Lepidoptera: Notodontidae) is the main pine defoliator in the Mediterranean region. Its urticating larvae cause severe human and animal health concerns in the invaded areas. This species shows a high phenotypic variability for various traits, such as phenology, fecundity and tolerance to extreme temperatures. This study presents the construction and analysis of extensive genomic and transcriptomic resources, which are an obligate prerequisite to understand their underlying genetic architecture. Using a well-studied population from Portugal with peculiar phenological characteristics, the karyotype was first determined and a first draft genome of 537 Mb total length was assembled into 68,292 scaffolds (N50 = 164 kb). From this genome assembly, 29,415 coding genes were predicted. To circumvent some limitations for fine-scale physical mapping of genomic regions of interest, a 3X coverage BAC library was also developed. In particular, 11 BACs from this library were individually sequenced to assess the assembly quality. Additionally, de novo transcriptomic resources were generated from various developmental stages sequenced with HiSeq and MiSeq Illumina technologies. The reads were de novo assembled into 62,376 and 63,175 transcripts, respectively. Then, a robust subset of the genome-predicted coding genes, the de novo transcriptome assemblies and previously published 454/Sanger data were clustered to obtain a high-quality and comprehensive reference transcriptome consisting of 29,701 bona fide unigenes. These sequences covered 99% of the cegma and 88% of the busco highly conserved eukaryotic genes and 84% of the busco arthropod gene set. Moreover, 90% of these transcripts could be localized on the draft genome. The described information is available via a genome annotation portal (http://bipaa.genouest.org/sp/thaumetopoea_pityocampa/).

PNPLA1 defects in patients with autosomal recessive congenital ichthyosis and KO mice sustain PNPLA1 irreplaceable function in epidermal omega-O-acylceramide synthesis and skin permeability barrier

Autosomal recessive congenital ichthyosis (ARCI) is a heterogeneous group of monogenic genodermatoses that encompasses non-syndromic disorders of keratinization. The pathophysiology of ARCI has been linked to a disturbance in epidermal lipid metabolism that impaired the stratum corneum function, leading to permeability barrier defects. Functional characterization of some genes involved in ARCI contributed to the identification of molecular actors involved in epidermal lipid synthesis, transport or processing. Recently, PNPLA1 has been identified as a gene causing ARCI. While other members of PNPLA family are key elements in lipid metabolism, the function of PNPLA1 remained unclear. We identified 5 novel PNPLA1 mutations in ARCI patients, mainly localized in the putative active enzymatic domain of PNPLA1. To investigate Pnpla1 biological role, we analysed Pnpla1-deficient mice. KO mice died soon after birth from severe epidermal permeability defects. Pnpla1-deficient skin presented an important impairment in the composition and organization of the epidermal lipids. Quantification of epidermal ceramide species highlighted a blockade in the production of ω-O-acylceramides with a concomitant accumulation of their precursors in the KO. The virtually loss of ω-O-acylceramides in the stratum corneum was linked to a defective lipid coverage of the resistant pericellular shell encapsulating corneocytes, the so-called cornified envelope, and most probably disorganized the extracellular lipid matrix. Finally, these defects in ω-O-acylceramides synthesis and cornified envelope formation were also evidenced in the stratum corneum from PNPLA1-mutated patients. Overall, our data support that PNPLA1/Pnpla1 is a key player in the formation of ω-O-acylceramide, a crucial process for the epidermal permeability barrier function.

Does training respiratory physicians in clinical respiratory physiology and interpretation of pulmonary function tests improve core knowledge?

Lung function tests have a major role in respiratory medicine. Training in lung function tests is variable within the European Union. In this study, we have shown that an internship in a lung function tests laboratory significantly improved the technical and diagnostic skills of French respiratory trainees.

Lowering relative humidity level increases epidermal protein deimination and drives human filaggrin breakdown

BACKGROUND

Deimination (also known as citrullination), the conversion of arginine in a protein to citrulline, is catalyzed by a family of enzymes called peptidylarginine deiminases (PADs). Three PADs are expressed in the epidermis, one of their targets being filaggrin. Filaggrin plays a central role in atopic dermatitis and is a key protein for the epidermal barrier. It aggregates keratins and is cross-linked to cornified envelopes. Following its deimination, it is totally degraded to release free amino acids, contributing to the natural moisturizing factor (NMF). The mechanisms controlling this multistep catabolism in human are unknown.

OBJECTIVE

To test whether external humidity plays a role, and investigate the molecular mechanisms involved.

METHODS

Specimens of reconstructed human epidermis (RHEs) produced in humid or dry conditions (>95% or 30-50% relative humidity) were compared.

RESULTS

RHEs produced in the dry condition presented structural changes, including a thicker stratum corneum and a larger amount of keratohyalin granules. The transepidermal water loss and the stratum corneum pH were decreased whereas the quantity of NMF was greater. This highly suggested that filaggrin proteolysis was up-regulated. The expression/activity of the proteases involved in filaggrin breakdown did not increase while PAD1 expression and the deimination rate of proteins, including filaggrin, were drastically enhanced. Partial inhibition of PADs with Cl-amidine reversed the effect of dryness on filaggrin breakdown.

CONCLUSION

These results demonstrate the importance of external humidity in the control of human filaggrin metabolism, and suggest that deimination plays a major role in this regulation.

A novel locus on chromosome 1 underlies the evolution of a melanic plumage polymorphism in a wild songbird

Understanding the mechanisms responsible for phenotypic diversification within and among species ultimately rests with linking naturally occurring mutations to functionally and ecologically significant traits. Colour polymorphisms are of great interest in this context because discrete colour patterns within a population are often controlled by just a few genes in a common environment. We investigated how and why phenotypic diversity arose and persists in the Zosterops borbonicus white-eye of Reunion (Mascarene archipelago), a colour polymorphic songbird in which all highland populations contain individuals belonging to either a brown or a grey plumage morph. Using extensive phenotypic and genomic data, we demonstrate that this melanin-based colour polymorphism is controlled by a single locus on chromosome 1 with two large-effect alleles, which was not previously described as affecting hair or feather colour. Differences between colour morphs appear to rely upon complex cis-regulatory variation that either prevents the synthesis of pheomelanin in grey feathers, or increases its production in brown ones. We used coalescent analyses to show that, from a 'brown' ancestral population, the dominant 'grey' allele spread quickly once it arose from a new mutation. Since colour morphs are always found in mixture, this implies that the selected allele does not go to fixation, but instead reaches an intermediate frequency, as would be expected under balancing selection.

SNP discovery and genetic mapping using genotyping by sequencing of whole genome genomic DNA from a pea RIL population

BACKGROUND

Progress in genetics and breeding in pea still suffers from the limited availability of molecular resources. SNP markers that can be identified through affordable sequencing processes, without the need for prior genome reduction or a reference genome to assemble sequencing data would allow the discovery and genetic mapping of thousands of molecular markers. Such an approach could significantly speed up genetic studies and marker assisted breeding for non-model species.

RESULTS

A total of 419,024 SNPs were discovered using HiSeq whole genome sequencing of four pea lines, followed by direct identification of SNP markers without assembly using the discoSnp tool. Subsequent filtering led to the identification of 131,850 highly designable SNPs, polymorphic between at least two of the four pea lines. A subset of 64,754 SNPs was called and genotyped by short read sequencing on a subpopulation of 48 RILs from the cross 'Baccara' x 'PI180693'. This data was used to construct a WGGBS-derived pea genetic map comprising 64,263 markers. This map is collinear with previous pea consensus maps and therefore with the Medicago truncatula genome. Sequencing of four additional pea lines showed that 33 % to 64 % of the mapped SNPs, depending on the pairs of lines considered, are polymorphic and can therefore be useful in other crosses. The subsequent genotyping of a subset of 1000 SNPs, chosen for their mapping positions using a KASP™ assay, showed that almost all generated SNPs are highly designable and that most (95 %) deliver highly qualitative genotyping results. Using rather low sequencing coverages in SNP discovery and in SNP inferring did not hinder the identification of hundreds of thousands of high quality SNPs.

CONCLUSIONS

The development and optimization of appropriate tools in SNP discovery and genetic mapping have allowed us to make available a massive new genomic resource in pea. It will be useful for both fine mapping within chosen QTL confidence intervals and marker assisted breeding for important traits in pea improvement.

From museums to genomics: old herbarium specimens shed light on a C3 to C4 transition

Collections of specimens held by natural history museums are invaluable material for biodiversity inventory and evolutionary studies, with specimens accumulated over 300 years readily available for sampling. Unfortunately, most museum specimens yield low-quality DNA. Recent advances in sequencing technologies, so called next-generation sequencing, are revolutionizing phylogenetic investigations at a deep level. Here, the Illumina technology (HiSeq) was used on herbarium specimens of Sartidia (subfamily Aristidoideae, Poaceae), a small African-Malagasy grass lineage (six species) characteristic of wooded savannas, which is the C3 sister group of Stipagrostis, an important C4 genus from Africa and SW Asia. Complete chloroplast and nuclear ribosomal sequences were assembled for two Sartidia species, one of which (S. perrieri) is only known from a single specimen collected in Madagascar 100 years ago. Partial sequences of a few single-copy genes encoding phosphoenolpyruvate carboxylases (ppc) and malic enzymes (nadpme) were also assembled. Based on these data, the phylogenetic position of Malagasy Sartidia in the subfamily Aristidoideae was investigated and the biogeographical history of this genus was analysed with full species sampling. The evolutionary history of two genes for C4 photosynthesis (ppc-aL1b and nadpme-IV) in the group was also investigated. The gene encoding the C4 phosphoenolpyruvate caroxylase of Stipagrostis is absent from S. dewinteri suggesting that it is not essential in C3 members of the group, which might have favoured its recruitment into a new metabolic pathway. Altogether, the inclusion of historical museum specimens in phylogenomic analyses of biodiversity opens new avenues for evolutionary studies.

Gonad transcriptome analysis of pearl oyster Pinctada margaritifera: identification of potential sex differentiation and sex determining genes

BACKGROUND

Black pearl farming is based on culture of the blacklip pearl oyster Pinctada margaritifera (Mollusca, lophotrochozoa), a protandrous hermaphrodite species. At first maturation, all individuals are males. The female sex appears progressively from two years old, which represents a limitation for broodstock conditioning for aquaculture production. In marine mollusks displaying hermaphroditic features, data on sexual determinism and differentiation, including the molecular sex determining cascade, are scarce. To increase genomic resources and identify the molecular mechanisms whereby gene expression may act in the sexual dimorphism of P. margaritifera, we performed gonad transcriptome analysis.

RESULTS

The gonad transcriptome of P. margaritifera was sequenced from several gonadic samples of males and females at different development stages, using a Next-Generation-Sequencing method and RNAseq technology. After Illumina sequencing, assembly and annotation, we obtained 70,147 contigs of which 62.2% shared homologies with existing protein sequences, and 9% showed functional annotation with Gene Ontology terms. Differential expression analysis identified 1,993 differentially expressed contigs between the different categories of gonads. Clustering methods of samples revealed that the sex explained most of the variation in gonad gene expression. K-means clustering of differentially expressed contigs showed 815 and 574 contigs were more expressed in male and female gonads, respectively. The analysis of these contigs revealed the presence of known specific genes coding for proteins involved in sex determinism and/or differentiation, such as dmrt and fem-1 like for males, or foxl2 and vitellogenin for females. The specific gene expression profiles of pmarg-fem1-like, pmarg-dmrt and pmarg-foxl2 in different reproductive stages (undetermined, sexual inversion and regression) suggest that these three genes are potentially involved in the sperm-oocyte switch in P. margaritifera.

CONCLUSIONS

The study provides a new transcriptomic tool to study reproduction in hermaphroditic marine mollusks. It identifies sex differentiation and potential sex determining genes in P. margaritifera, a protandrous hermaphrodite species.

A DNA metabarcoding study of a primate dietary diversity and plasticity across its entire fragmented range

In tropical regions, most primary ecosystems have been replaced by mosaic landscapes in which species must cope with a large shift in the distribution of their habitat and associated food resources. Primates are particularly vulnerable to habitat modifications. Most species persist in small fragments surrounded by complex human-mediated matrices whose structure and connectivity may strongly influence their dispersal and feeding behavior. Behavioral plasticity appears to be a crucial parameter governing the ability of organisms to exploit the resources offered by new matrix habitats and thus to persist in fragmented habitats. In this study, we were interested in the dietary plasticity of the golden-crowned sifaka (Propithecus tattersalli), an endangered species of lemur, found only in the Daraina region in north-eastern Madagascar. We used a DNA-based approach combining the barcoding concept and Illumina next-generation sequencing to (i) describe the species diet across its entire range and (ii) evaluate the influence of landscape heterogeneity on diet diversity and composition. Faeces from 96 individuals were sampled across the entire species range and their contents were analyzed using the trnL metabarcoding approach. In parallel, we built a large DNA reference database based on a checklist of the plant species of the Daraina region. Our results suggest that golden-crowned sifakas exhibit remarkable dietary diversity with at least 130 plant species belonging to 80 genera and 49 different families. We highlighted an influence of both habitat type and openness on diet composition suggesting a high flexibility of foraging strategies. Moreover, we observed the presence of numerous cultivated and naturalized plants in the faeces of groups living in forest edge areas. Overall, our findings support our initial expectation that P. tattersalli is able to cope with the current level of alteration of the landscape and confirm our previous results on the distribution and the dispersal ability of this species.

Semiochemical compounds of preen secretion reflect genetic make-up in a seabird species

Several vertebrates choose their mate according to genetic heterozygosity and relatedness, and use odour cues to assess their conspecifics' genetic make-up. In birds, although several species (including the black-legged kittiwake) exhibit non-random mating according to genetic traits, the cues used to assess genetic characteristics remain unknown. The importance of olfaction in birds' social behaviour is gaining attention among researchers, and it has been suggested that, as in other vertebrates, bird body scent may convey information about genetic traits. Here, we combined gas chromatography data and genetic analyses at microsatellite loci to test whether semiochemical messages in preen secretion of kittiwakes carried information about genetic heterozygosity and relatedness. Semiochemical profile was correlated with heterozygosity in males and females, while semiochemical distance was correlated with genetic distance only in male-male dyads. Our study is the first to demonstrate a link between odour and genetics in birds, which sets the stage for the existence of sophisticated odour-based mechanisms of mate choice also in birds.

Isolation and characterization of microsatellite loci in Rhododendron ferrugineum (Ericaceae) using pyrosequencing technology

PREMISE OF THE STUDY

Microsatellite primers were developed for Rhododendron ferrugineum (Ericaceae) to evaluate genetic diversity, population genetic structure, and mating system of this self-compatible species.

METHODS AND RESULTS

The new-generation 454 FLX Titanium pyrosequencing technology was used to isolate 102 novel microsatellite loci. Two multiplex PCR sets were optimized to genotype nine polymorphic microsatellite loci. The level of genetic diversity was assessed in two populations from the Pyrenees (France). The mean number of alleles per locus ranged from 1.5 to 7. The mean observed and expected heterozygosities ranged from 0 to 0.76 and from 0.03 to 0.66, respectively. Cross-species amplification was successful for 13 other Rhododendron species and two additional genera of Ericaceae, with an average of seven pairs of primers amplifying per species.

CONCLUSIONS

These markers will facilitate further studies on the evolutionary history of the large Ericaceae family and Rhododendron in particular.