The genome of a thorny species: comparative genomic analysis among South and North American Cactaceae.

MAIN CONCLUSION: The first South American cactus nuclear genome assembly associated with comparative genomic analyses provides insights into nuclear and plastid genomic features, such as size, transposable elements, and metabolic processes related to cactus development. Here, we assembled the partial genome, plastome, and transcriptome of Cereus fernambucensis (Cereeae, Cactaceae), a representative species of the South American core Cactoideae. We accessed other genomes and transcriptomes available for cactus species to compare the heterozygosity level, genome size, transposable elements, orthologous genes, and plastome structure. These estimates were obtained from the literature or using the same pipeline adopted for C. fermabucensis. In addition to the C. fernambucensis plastome, we also performed de novo plastome assembly of Pachycereus pringlei, Stenocereus thurberi, and Pereskia humboldtii based on the sequences available in public databases. We estimated a genome size of ~ 1.58 Gb for C. fernambucensis, the largest genome among the compared species. The genome heterozygosity was 0.88% in C. fernambucensis but ranged from 0.36 (Carnegiea gigantea) to 17.4% (Lophocereus schottii) in the other taxa. The genome lengths of the studied cacti are constituted by a high amount of transposable elements, ranging from ~ 57 to ~ 67%. Putative satellite DNAs are present in all species, excepting C. gigantea. The plastome of C. fernambucensis was ~ 104 kb, showing events of translocation, inversion, and gene loss. We observed a low number of shared unique orthologs, which may suggest gene duplication events and the simultaneous expression of paralogous genes. We recovered 37 genes that have undergone positive selection along the Cereus branch that are associated with different metabolic processes, such as improving photosynthesis during drought stress and nutrient absorption, which may be related to the adaptation to xeric areas of the Neotropics.

Early pregnancy-induced transcripts in peripheral blood immune cells in Bos indicus heifers.

Immune cells play a central role in early pregnancy establishment in cattle. We aimed to: (1) discover novel early-pregnancy-induced genes in peripheral blood mononuclear cells (PBMC); and (2) characterize the temporal pattern of early-pregnancy-induced transcription of select genes in PBMC and peripheral blood polymorphonuclear cells (PMN). Beef heifers were artificially inseminated on D0 and pregnancies were diagnosed on D28. On D10, 14, 16, 18, and 20, blood was collected for isolation of PBMC and PMN from heifers that were retrospectively classified as pregnant (P) or non-pregnant (NP). PBMC samples from D18 were submitted to RNAseq and 220 genes were differentially expressed between pregnant (P) and non-pregnant (NP) heifers. The temporal abundance of 20 transcripts was compared between P and NP, both in PBMC and PMN. In PBMC, pregnancy stimulated transcription of IFI6, RSAD2, IFI44, IFITM2, CLEC3B, OAS2, TNFSF13B, DMKN and LGALS3BP as early as D18. Expression of IFI44, RSAD2, OAS2, LGALS3BP, IFI6 and C1R in PMN was stimulated in the P group from D18. The novel early-pregnancy induced genes discovered in beef heifers will allow both the understanding of the role of immune cells during the pre-attachment period and the development of technologies to detect early pregnancies in beef cattle.

The potential of genome-wide RAD sequences for resolving rapid radiations: a case study in Cactaceae.

The reconstruction of relationships within recently radiated groups is challenging even when massive amounts of sequencing data are available. The use of restriction site-associated DNA sequencing (RAD-Seq) to this end is promising. Here, we assessed the performance of RAD-Seq to infer the species-level phylogeny of the rapidly radiating genus Cereus (Cactaceae). To examine how the amount of genomic data affects resolution in this group, we used datasets and implemented different analyses. We sampled 52 individuals of Cereus, representing 18 of the 25 species currently recognized, plus members of the closely allied genera Cipocereus and Praecereus, and other 11 Cactaceae genera as outgroups. Three scenarios of permissiveness to missing data were carried out in iPyRAD, assembling datasets with 30% (333 loci), 45% (1440 loci), and 70% (6141 loci) of missing data. For each dataset, Maximum Likelihood (ML) trees were generated using two supermatrices, i.e., only SNPs and SNPs plus invariant sites. Accuracy and resolution were improved when the dataset with the highest number of loci was used (6141 loci), despite the high percentage of missing data included (70%). Coalescent trees estimated using SVDQuartets and ASTRAL are similar to those obtained by the ML reconstructions. Overall, we reconstruct a well-supported phylogeny of Cereus, which is resolved as monophyletic and composed of four main clades with high support in their internal relationships. Our findings also provide insights into the impact of missing data for phylogeny reconstruction using RAD loci.

Germline Mutation in MUS81 Resulting in Impaired Protein Stability is Associated with Familial Breast and Thyroid Cancer.

Multiple primary thyroid cancer (TC) and breast cancer (BC) are commonly diagnosed, and the lifetime risk for these cancers is increased in patients with a positive family history of both TC and BC. Although this phenotype is partially explained by TP53 or PTEN mutations, a significant number of patients are negative for these alterations. We judiciously recruited patients diagnosed with BC and/or TC having a family history of these tumors and assessed their whole-exome sequencing. After variant prioritization, we selected MUS81 c.1292G>A (p.R431H) for further investigation. This variant was genotyped in a healthy population and sporadic BC/TC tissues and investigated at the protein level and cellular models. MUS81 c.1292G>A was the most frequent variant (25%) and the strongest candidate due to its function of double-strand break repair. This variant was confirmed in four relatives from two families. MUS81 p.R431H protein exhibited lower expression levels in tumors from patients positive for the germline variant, compared with wild-type BC, and normal breast and thyroid tissues. Using cell line models, we showed that c.1292G>A induced protein instability and affected DNA damage response. We suggest that MUS81 is a novel candidate involved in familial BC/TC based on its low frequency in healthy individuals and proven effect in protein stability.

Beyond the limits of photoperception: constitutively active PHYTOCHROME B2 overexpression as a means of improving fruit nutritional quality in tomato.

Photoreceptor engineering has recently emerged as a means for improving agronomically beneficial traits in crop species. Despite the central role played by the red/far-red photoreceptor phytochromes (PHYs) in controlling fruit physiology, the applicability of PHY engineering for increasing fleshy fruit nutritional content remains poorly exploited. In this study, we demonstrated that the fruit-specific overexpression of a constitutively active GAF domain Tyr252 -to-His PHYB2 mutant version (PHYB2Y252H ) significantly enhances the accumulation of multiple health-promoting antioxidants in tomato fruits, without negative collateral consequences on vegetative development. Compared with the native PHYB2 overexpression, PHYB2Y252H -overexpressing lines exhibited more extensive increments in transcript abundance of genes associated with fruit plastid development, chlorophyll biosynthesis and metabolic pathways responsible for the accumulation of antioxidant compounds. Accordingly, PHYB2Y252H -overexpressing fruits developed more chloroplasts containing voluminous grana at the green stage and overaccumulated carotenoids, tocopherols, flavonoids and ascorbate in ripe fruits compared with both wild-type and PHYB2-overexpressing lines. The impacts of PHYB2 or PHYB2Y252H overexpression on fruit primary metabolism were limited to a slight promotion in lipid biosynthesis and reduction in sugar accumulation. Altogether, these findings indicate that mutation-based adjustments in PHY properties represent a valuable photobiotechnological tool for tomato biofortification, highlighting the potential of photoreceptor engineering for improving quality traits in fleshy fruits.

Transcriptional foliar profile of the C3-CAM bromeliad Guzmania monostachia.

Guzmania monostachia is an epiphytic tank bromeliad that displays the inducible CAM photosynthesis under stressful conditions and had the highest stomata density in the leaf apex, while the base portion has the highest density of trichomes, which are specialized structures used to acquire water and nutrients from the tank solution. In order to correlate the genetic factors behind these morpho-physiological characteristics along the leaf blade of G. monostachia, a comparative transcriptome analysis was performed to identify the functional enriched pathways and unigenes that could play a role in the apical, middle and basal leaf portions. A total of 653 million reads were used for de novo transcriptome assembly, resulting in 48,051 annotated unigenes. Analysis of differentially expressed genes (DEGs) among distinct leaf regions revealed that 806 DEGs were upregulated in the apex compared to the middle portion, while 9685 DEGs were upregulated in the apex and 9784 DEGs were upregulated in the middle portions compared to the base. Our outcomes correlated some DEGs and identified unigenes with their physiological functions, mainly suggesting that the leaf apex was related to the regulation of stomatal movement, production of chlorophyll, cellular response to stress, and H2O2 catabolic process. In contrast, the middle portion showed DEGs associated with the transport of amino acids. Furthermore, DEGs from the leaf base were mainly correlated with responses to nutrients and nitrogen compounds, regulation of potassium ion import, response to water deprivation, and trichome branching, indicating that, at least in part, this leaf portion can replace some of the root functions of terrestrial plants. Therefore, possibly candidate unigenes and enriched pathways presented here could be prospected in future experimental work, opening new possibilities to bioengineer non-inducible CAM plants and/or improve the fertilization use efficiency by increasing leaf nutrient acquisition of crop plants.

Analysis of the common bean (Phaseolus vulgaris L.) transcriptome regarding efficiency of phosphorus use.

Common bean is a highly important food in tropical regions, where most production occurs on small farms with limited use of technology and, consequently, greater vulnerability to abiotic stresses such as nutritional stress. Usually phosphorus (P) is the most limiting nutrient for crop growth in these regions. The aim of this study was to characterize the gene expression profiles of the genotypes of common bean IAC Imperador (P-responsive) and DOR 364 (P-unresponsive) under different P concentrations using RNA-seq transcriptome sequencing technology. Plants were grown hydroponically, with application of two P concentrations (4.00 mg L-1 restrictive level and 8.00 mg L-1 control level). Differential expression analyses, annotation, and functional classification were performed comparing genotypes within each P rate administered and comparing each genotype response to the different P levels. Considering differential expression analyses within genotypes, IAC Imperador exhibited 1538 up-regulated genes under P restriction and 1679 up-regulated genes in the control, while DOR 364 exhibited 13 up-regulated genes in the control and only 2 up-regulated genes under P restriction, strongly corroborating P-unresponsiveness of this genotype. Genes related to phosphorus restriction were identified among the differentially expressed genes, including transcription factors such as WRKY, ERF, and MYB families, phosphatase related genes such as pyrophosphatase, acid phosphatase, and purple acid phosphatase, and phosphate transporters. The enrichment test for the P restriction treatment showed 123 enriched gene ontologies (GO) for IAC Imperador, while DOR 364 enriched only 24. Also, the enriched GO correlated with P metabolism, compound metabolic processes containing phosphate, nucleoside phosphate binding, phosphorylation, and also response to stresses. Thus, this study proved to be informative to phosphorus limitation in common bean showing global changes at transcript level.

Gene Co-expression Analysis Indicates Potential Pathways and Regulators of Beef Tenderness in Nellore Cattle.

Beef tenderness, a complex trait affected by many factors, is economically important to beef quality, industry, and consumer's palatability. In this study, RNA-Seq was used in network analysis to better understand the biological processes that lead to differences in beef tenderness. Skeletal muscle transcriptional profiles from 24 Nellore steers, selected by extreme estimated breeding values (EBVs) for shear force after 14 days of aging, were analyzed and 22 differentially expressed transcripts were identified. Among these were genes encoding ribosomal proteins, glutathione transporter ATP-binding cassette, sub-family C (CFTR/MRP), member 4 (ABCC4), and synaptotagmin IV (SYT4). Complementary co-expression analyses using Partial Correlation with Information Theory (PCIT), Phenotypic Impact Factor (PIF) and the Regulatory Impact Factor (RIF) methods identified candidate regulators and related pathways. The PCIT analysis identified ubiquitin specific peptidase 2 (USP2), growth factor receptor-bound protein 10 (GBR10), anoctamin 1 (ANO1), and transmembrane BAX inhibitor motif containing 4 (TMBIM4) as the most differentially hubbed (DH) transcripts. The transcripts that had a significant correlation with USP2, GBR10, ANO1, and TMBIM4 enriched for proteasome KEGG pathway. RIF analysis identified microRNAs as candidate regulators of variation in tenderness, including bta-mir-133a-2 and bta-mir-22. Both microRNAs have target genes present in the calcium signaling pathway and apoptosis. PIF analysis identified myoglobin (MB), enolase 3 (ENO3), and carbonic anhydrase 3 (CA3) as potentially having fundamental roles in tenderness. Pathways identified in our study impacted in beef tenderness included: calcium signaling, apoptosis, and proteolysis. These findings underscore some of the complex molecular mechanisms that control beef tenderness in Nellore cattle.

Impact of hormonal modulation at proestrus on ovarian responses and uterine gene expression of suckled anestrous beef cows.

BACKGROUND: This study evaluated the impact of hormonal modulation at the onset of proestrus on ovarian response and uterine gene expression of beef cows. METHODS: A total of 172 anestrous beef cows were assigned to one of four groups according to the treatment with estradiol cypionate (ECP) and/or equine chorionic gonadotropin (eCG) [CON (n = 43), ECP (n = 43), eCG (n = 44) and ECP + eCG (n = 42)]. RESULTS: ECP-treated cows (ECP and ECP + eCG groups) presented greater occurrence of estrus (44.6% vs. 65.4%; P = 0.01) and pregnancy per AI [47.1% vs. 33.3%; P = 0.07], but similar progesterone (P4) concentration at subsequent diestrus than cows not treated with ECP (CON and eCG groups). Nonetheless, eCG-treated cows (eCG and ECP + eCG groups) presented larger follicle at timed AI (12.6 ± 0.3 vs. 13.5 ± 0.3 mm; P = 0.03), greater ovulation rate (96.5% vs. 82.6%; P = 0.008) and greater P4 concentration at d 6 (3.9 ± 0.2 vs. 4.8 ± 0.2 ng/mL; P = 0.001) than cows not treated with eCG (CON and ECP groups). Next, cows with a new corpus luteum 6 d after TAI were submitted to uterine biopsy procedure. Uterine fragments [CON (n = 6), ECP (n = 6)] were analyzed by RNA-Seq and a total of 135 transcripts were differentially expressed between groups (73 genes up-regulated by ECP treatment). Subsequently, uterine samples were analyzed by qPCR (genes associated with cell proliferation). ECP treatment induced greater abundance of PTCH2 (P = 0.07) and COL4A1 (P = 0.02), whereas suppressed EGFR (P = 0.09) expression. Conversely, eCG treatment increased abundance of HB-EGF (P = 0.06), ESR2 (P = 0.09), and ITGB3 (P = 0.05), whereas it reduced transcription of ESR1 (P = 0.05). Collectively, supplementation with ECP or eCG at the onset of proestrous of anestrous beef cows influenced ovarian responses, global and specific endometrial gene expression. CONCLUSION: Proestrus estradiol regulate the endometrial transcriptome, particularly stimulating proliferative activity in the endometrium.

Oviductal transcriptional profiling of a bovine fertility model by next-generation sequencing.

In cattle, the oviduct plays a fundamental role in the reproductive process. Oviductal functions are controlled by the ovarian sex steroids: estradiol and progesterone. Here, we tested the hypothesis that the exposure to contrasting sex steroid milieus differentially impacts the oviductal transcriptional profile. We manipulated growth of the pre-ovulatory follicle to obtain cows that ovulated a larger (LF group) or a smaller (SF group) follicle. The LF group presented greater proestrus/estrus concentrations of estradiol and metaestrus concentrations of progesterone (Gonella-Diaza et al. 2015 [1], Mesquita et al. 2014 [2]). Also, the LF group was associated with greater fertility in timed-artificial insemination programs (Pugliesi et al. 2016 [3]). Cows were slaughtered on day 4 of the estrous cycle and total RNA was extracted from ampulla and isthmus fragments and analyzed by RNAseq. The resulting reads were mapped to the bovine genome (Bos taurus UMD 3.1, NCBI). The differential expression analyses revealed that 325 and 367 genes in ampulla and 274 and 316 genes in the isthmus were up-regulated and down-regulated in LF samples, respectively. To validate the RNAseq results, transcript abundance of 23 genes was assessed by qPCR and expression patterns were consistent between the two techniques. A functional enrichment analysis was performed using Database for Annotation, Visualization and Integrated Discovery (DAVID) software. Processes enriched in the LF group included tissue morphology changes (extracellular matrix remodeling), cellular changes (proliferation), and secretion changes (growth factors, ions and metal transporters). An overview of the gene expression data was deposited in the NCBI's Gene Expression Omnibus (GEO) and is accessible through the accession number GSE65681. In conclusion, differences in the peri-ovulatory sex steroid milieu modify the oviductal gene expression profiles. Such differences may be associated with the greater fertility of the LF cows. This dataset is useful for further investigations of the oviductal biology and the impact of sex-steroid on the female reproductive tract.

Dynamic remodeling of endometrial extracellular matrix regulates embryo receptivity in cattle.

We aimed to evaluate in the bovine endometrium whether (1) key genes involved in endometrial extracellular matrix (ECM) remodeling are regulated by the endocrine peri-ovulatory milieu; and (2) specific endometrial ECM-related transcriptome can be linked to pregnancy outcome. In Experiment 1, pre-ovulatory follicle growth of cows was manipulated to obtain two groups with specific endocrine peri-ovulatory profiles: the Large Follicle-Large CL group (LF-LCL) served as a paradigm for greater receptivity and fertility and showed greater plasma pre-ovulatory estradiol and post-ovulatory progesterone concentrations when compared to the Small Follicle-Small CL group (SF-SCL). Endometrium was collected on days 4 and 7 of the estrous cycle. Histology revealed a greater abundance of total collagen content in SF-SCL on day 4 endometrium. In Experiment 2, cows were artificially inseminated and, six days later, endometrial biopsies were collected. Cows were retrospectively divided into pregnant and non-pregnant (P vs. NP) groups after diagnosis on day 30. In both experiments, expression of genes related to ECM remodeling in the endometrium was studied by RNAseq and qPCR. Gene ontology analysis showed an inhibition in the expression of ECM-related genes in the high receptivity groups (LF-LCL and P). Specifically, there was down-regulation of TGFB2, ADAMTS2, 5 and 14, TIMP3 and COL1A2, COL3A1, COL7A1 and COL3A3 in the LF-LCL and P groups. In summary, the overlapping set of genes differently expressed in both fertility models: (1) suggests that disregulation of ECM remodeling can impair receptivity and (2) can be used as markers to predict pregnancy outcome in cattle.

Integrative transcriptome analysis identifies deregulated microRNA-transcription factor networks in lung adenocarcinoma.

Herein, we aimed at identifying global transcriptome microRNA (miRNA) changes and miRNA target genes in lung adenocarcinoma. Samples were selected as training (N = 24) and independent validation (N = 34) sets. Tissues were microdissected to obtain >90% tumor or normal lung cells, subjected to miRNA transcriptome sequencing and TaqMan quantitative PCR validation. We further integrated our data with published miRNA and mRNA expression datasets across 1,491 lung adenocarcinoma and 455 normal lung samples. We identified known and novel, significantly over- and under-expressed (p ≤ 0.01 and FDR≤0.1) miRNAs in lung adenocarcinoma compared to normal lung tissue: let-7a, miR-10a, miR-15b, miR-23b, miR-26a, miR-26b, miR-29a, miR-30e, miR-99a, miR-146b, miR-181b, miR-181c, miR-421, miR-181a, miR-574 and miR-1247. Validated miRNAs included let-7a-2, let-7a-3, miR-15b, miR-21, miR-155 and miR-200b; higher levels of miR-21 expression were associated with lower patient survival (p = 0.042). We identified a regulatory network including miR-15b and miR-155, and transcription factors with prognostic value in lung cancer. Our findings may contribute to the development of treatment strategies in lung adenocarcinoma.

Genome organization and host range of a Brazilian isolate of johnsongrass mosaic virus.

This work reports the complete genome sequence, production of a polyclonal antiserum, and host range of a Brazilian strain of johnsongrass mosaic virus (JGMV) found infecting Panicum maximum in the state of São Paulo, Brazil. The complete genome sequence of this potyvirus, comprising 9874 nucleotides, showed 82 % amino acid sequence identity in the polyprotein to that of an isolate of JGMV from Australia. The experimental host range of this virus included mainly fodder species. Cultivated species such as rice, oats, sugarcane, rye, corn and wheat were not infected, suggesting that current isolates of this potyvirus do not represent a threat to these crops in Brazil.

At least some meiofaunal species are not everywhere. Indication of geographic, ecological and geological barriers affecting the dispersion of species of Ototyphlonemertes (Nemertea, Hoplonemertea).

Most meiofaunal species are known to have a broad distribution with no apparent barriers to their dispersion. However, different morphological and/or molecular methods supported patterns of diversity and distribution that may be different among taxa while also conflicting within the same group. We accurately assessed the patterns of geographic distribution in actual genetic species of a marine meiofaunal animal model: Ototyphlonemertes. Specimens were collected from several sites around Europe, Northern and Central America, Southern America, Pacific Islands and Asia. We sequenced regions of two mitochondrial and two nuclear genes. Using single-gene, a concatenated data set, multilocus approaches and different DNA taxonomy methods, we disentangled the actual diversity and the spatial structures of haplotypes and tested the possible correlation between genetic diversity and geographic distance. The results show (i) the importance of using several genes to uncover both diversity and highlight phylogeographic relationships among species and that (ii) independent genetic evolutionary entities have a narrower distribution than morphological species. Moreover, (iii) a Mantel test supported a positive correlation between genetic and geographical distance. By sampling from the two sides of Isthmus of Panama, we were additionally able to identify lineage divergence times that are concordant with vicariance mechanisms caused by the geological closure of the seaway across the Isthmus. We therefore propose that in addition to distance, other geological and ecological conditions are also barriers to the dispersion of and gene flow in marine meiofaunal organisms.

Size of the Ovulatory Follicle Dictates Spatial Differences in the Oviductal Transcriptome in Cattle.

In cattle, molecular control of oviduct receptivity to the embryo is poorly understood. Here, we used a bovine model for receptivity based on size of the pre-ovulatory follicle to compare oviductal global and candidate gene transcript abundance on day 4 of the estrous cycle. Growth of the pre-ovulatory follicle (POF) of Nelore (Bos indicus) cows was manipulated to produce two groups: large POF large corpus luteum (CL) group (LF-LCL; greater receptivity) and small POF-small CL group (SF-SCL). Oviductal samples were collected four days after GnRH-induced ovulation. Ampulla and isthmus transcriptome was obtained by RNA-seq, regional gene expression was assessed by qPCR, and PGR and ERa protein distribution was evaluated by immunohistochemistry. There was a greater abundance of PGR and ERa in the oviduct of LF-LCL animals thus indicating a greater availability of receptors and possibly sex steroids stimulated signaling in both regions. Transcriptomic profiles indicated a series of genes associated with functional characteristics of the oviduct that are regulated by the periovulatory sex steroid milieu and that potentially affect oviductal receptivity and early embryo development. They include tissue morphology changes (extra cellular matrix remodeling), cellular changes (proliferation), and secretion changes (growth factors, ions and metal transporters), and were enriched for the genes with increased expression in the LF-LCL group. In conclusion, differences in the periovulatory sex steroid milieu lead to different oviductal gene expression profiles that could modify the oviductal environment to affect embryo survival and development.

Genetic variation and population structuring in two brooding coral species (Siderastrea stellata and Siderastrea radians) from Brazil.

Siderastrea stellata and S. radians are scleractinian coral species that present a remarkable overlap of diagnostic characteristics and sympatric distribution. Moreover, both are viviparous with similar reproductive strategies and with a gregarious larval behavior. Samples of both species from the Brazilian coast were analyzed using 18 isozymic loci to quantify their genetic variability and populational structure. Results confirmed species identity, high intrapopulational variability and revealed moderate genetic structuring among all samples (S. stellata: F (ST) = 0.070; S. radians: F (ST) = 0.092). Based on genotypic diversity analysis, there was evidence that local recruitment may have a minor role in the populations (mean, G ( o ):G ( e )= 1.00 +/- 0.0003 SD for S. stellata and 0.99 +/- 0.0023 SD for S. radians). Deviations towards heterozygote deficiencies found in both Siderastrea species could be explained by the Wahlund effect, since there was evidence that populations might be composed of colonies of different ages. In S. radians it is also likely that there is some inbreeding occurring in the studied populations. Despite the brooding pattern and the gregarious larval behavior, our data suggest the occurrence of gene flow along the Brazilian coast. This is the first study on population genetics of Brazilian reef corals.

Allozyme diversity and morphometrics of Melocactus paucispinus (Cactaceae) and evidence for hybridization with M. concinnus in the Chapada Diamantina, North-eastern Brazil.

BACKGROUND AND AIMS: Melocacatus paucispinus (Cactaceae) is endemic to the state of Bahia, Brazil, and due to its rarity and desirability to collectors it has been considered threatened with extinction. This species is usually sympatric and inter-fertile with M. concinnus, and morphological evidence for hybridization between them is present in some populations. Levels of genetic and morphological variation and sub-structuring in populations of these species were assessed and an attempt was made to verify the occurrence of natural hybridization between them. METHODS: Genetic variability was surveyed using allozymes (12 loci) and morphological variability using multivariate morphometric analyses (17 vegetative characters) in ten populations of M. paucispinus and three of M. concinnus occurring in the Chapada Diamantina, Bahia. KEY RESULTS: Genetic variability was low in both species (P = 0.0-33.3, A = 1.0-1.6, H(e) = 0.000-0.123 in M. paucispinus; P = 0.0-25.0, A = 1.0-1.4, H(e) = 0.000-0.104 in M. concinnus). Deficit of heterozygotes within the populations was detected in both species, with high values of F(IS) (0.732 and 0.901 in M. paucispinus and M. concinnus, respectively). Evidence of hybridization was detected by the relative allele frequency in the two diaphorase loci. High levels of genetic (F(ST) = 0.504 in M. paucispinus and 0.349 in M. concinnus) and morphological (A = 0.20 in M. paucispinus and 0.17 in M. concinnus) structuring among populations were found. CONCLUSIONS: The Melocactus spp. displayed levels of genetic variability lower than the values reported for other cactus species. The evidence indicates the occurrence of introgression in both species at two sites. The high F(ST) values cannot be explained by geographical substructuring, but are consistent with hybridization. Conversely, morphological differentiation in M. paucispinus, but not in M. concinnus, is probably due to isolation by distance.