In silico phylogenomics using complete genomes: a case study on the evolution of hominoids.

The increasing availability of complete genome data is facilitating the acquisition of phylogenomic data sets, but the process of obtaining orthologous sequences from other genomes and assembling multiple sequence alignments remains piecemeal and arduous. We designed software that performs these tasks and outputs anonymous loci (AL) or anchored enrichment/ultraconserved element loci (AE/UCE) data sets in ready-to-analyze formats. We demonstrate our program by applying it to the hominoids. Starting with human, chimpanzee, gorilla, and orangutan genomes, our software generated an exhaustive data set of 292 ALs (∼1 kb each) in ∼3 h. Not only did analyses of our AL data set validate the program by yielding a portrait of hominoid evolution in agreement with previous studies, but the accuracy and precision of our estimated ancestral effective population sizes and speciation times represent improvements. We also used our program with a published set of 512 vertebrate-wide AE "probe" sequences to generate data sets consisting of 171 and 242 independent loci (∼1 kb each) in 11 and 13 min, respectively. The former data set consisted of flanking sequences 500 bp from adjacent AEs, while the latter contained sequences bordering AEs. Although our AE data sets produced the expected hominoid species tree, coalescent-based estimates of ancestral population sizes and speciation times based on these data were considerably lower than estimates from our AL data set and previous studies. Accordingly, we suggest that loci subjected to direct or indirect selection may not be appropriate for coalescent-based methods. Complete in silico approaches, combined with the burgeoning genome databases, will accelerate the pace of phylogenomics.

Comparative mitogenomic analyses of Amazona parrots and Psittaciformes.

Amazon parrots are long-lived birds with highly developed cognitive skills, including vocal learning. Several parrot mitogenomes have been sequenced, but important aspects of their organization and evolution are not fully understood or have limited experimental support. The main aim of the present study was to describe the mitogenome of the blue-fronted Amazon, Amazona aestiva, and compare it to other mitogenomes from the genus Amazona and the order Psittaciformes. We observed that mitogenomes are highly conserved among Amazon parrots, and a detailed analysis of their duplicated control regions revealed conserved blocks. Population level analyses indicated that the specimen analyzed here seems to be close to A. aestiva individuals from Bahia state. Evolutionary relationships of 41 Psittaciformes species and three outgroups were inferred by BEAST. All relationships were retrieved with high support.

The complete mitochondrial genome of the threatened Neotropical catfish Lophiosilurus alexandri (Silurifomes: Pseudopimelodidae) and phylogenomic analysis indicate monophyly of Pimelodoidea.

Lophiosilurus alexandri is an endemic catfish from the São Francisco River Basin (Brazil) popularly known as pacamã, which has economic potential for aquaculture farming. The mitochondrial genome was sequenced for the threatened Neotropical catfish L. alexandri. Assembly into scaffolds using MIRA and MITObim software produced the whole, circularized mitochondrial genome, which comprises 16,445 bp and presents the typical gene arrangement of Teleostei mitochondria. A phylogenomic analysis was performed after the concatenation of all proteins obtained from whole mitogenomes of 20 Siluriformes and two outgroups. The results confirmed the monophyly of nine families of catfishes and also clustered L. alexandri as a sister group to the family Pimelodidae, thus confirming the monophyly of the superfamily Pimelodoidea. This is the first mitochondrial phylogenomics study for Pimelodoidea and the first mitogenome described for the Pseudopimelodidae family, representing an important resource for phylogeography, evolutionary biology, and conservation genetics studies in Neotropical fishes.

NRIP1 is activated by C-JUN/C-FOS and activates the expression of PGR, ESR1 and CCND1 in luminal A breast cancer.

Using chip array assays, we identified differentially expressed genes via a comparison between luminal A breast cancer subtype and normal mammary ductal cells from healthy donors. In silico analysis confirmed by western blot and immunohistochemistry revealed that C-JUN and C-FOS transcription factors are activated in luminal A patients as potential upstream regulators of these differentially expressed genes. Using a chip-on-chip assay, we identified potential C-JUN and C-FOS targets. Among these genes, the NRIP1 gene was revealed to be targeted by C-JUN and C-FOS. This was confirmed after identification and validation with transfection assays specific binding of C-JUN and C-FOS at consensus binding sites. NRIP1 is not only upregulated in luminal A patients and cell lines but also regulates breast cancer-related genes, including PR, ESR1 and CCND1. These results were confirmed by NRIP1 siRNA knockdown and chip array assays, thus highlighting the putative role of NRIP1 in PGR, ESR1 and CCND1 transcriptional regulation and suggesting that NRIP1 could play an important role in breast cancer ductal cell initiation.

Angiogenesis and evading immune destruction are the main related transcriptomic characteristics to the invasive process of oral tongue cancer.

Metastasis of head and neck tumors is responsible for a high mortality rate. Understanding its biochemistry may allow insights into tumorigenesis. To that end we carried out RNA-Seq analyses of 5 SCC9 derived oral cancer cell lines displaying increased invasive potential. Differentially expressed genes (DEGs) were annotated based on p-values and false discovery rate (q-values). All 292 KEGG pathways related to the human genome were compared in order to pinpoint the absolute and relative contributions to the invasive process considering the 8 hallmarks of cancer plus 2 new defined categories, as well as we made with our transcriptomic data. In terms of absolute contribution, the highest correlations were associated to the categories of evading immune destruction and energy metabolism and for relative contributions, angiogenesis and evading immune destruction. DEGs were distributed into each one of all possible modes of regulation, regarding up, down and continuum expression, along the 3 stages of metastatic progression. For p-values twenty-six genes were consistently present along the tumoral progression and 4 for q-values. Among the DEGs, we found 2 novel potentially informative metastatic markers: PIGG and SLC8B1. Furthermore, interactome analysis showed that MYH14, ANGPTL4, PPARD and ENPP1 are amenable to pharmacological interventions.