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.

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.

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.

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.

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.