Orchestration of miRNA Patterns by Testosterone and Dietary Tomato Carotenoids during Early Prostate Carcinogenesis in TRAMP Mice.

BACKGROUND: The integrative effects of prostate cancer risk factors, such as diet and endocrine status, on cancer-associated miRNA expression are poorly defined. OBJECTIVES: This study aimed to define the influence of androgens and diet (tomato and lycopene) on prostatic miRNA expression during early carcinogenesis in the transgenic adenocarcinoma of the mouse prostate (TRAMP) model. METHODS: Wild type (WT) and TRAMP mice were fed control, tomato-containing, or lycopene-containing diets from 4 to 10 weeks of age. Mice underwent either sham (intact) or castration surgery at 8 wk, and half of the castrated mice received testosterone (2.5 mg/kg body weight/d) at 9 wk. Mice were killed at 10 wk, and dorsolateral prostate expression of 602 miRNAs was assessed. RESULTS: We detected expression of 88 miRNAs (15% of 602), all of which were present in the TRAMP, in comparison with 49 miRNAs being detectable (8%) in WT. Expression of 61 miRNAs differed by TRAMP genotype, with the majority upregulated in TRAMP. Of the 61 miRNAs, 42 were responsive to androgen status. Diet affected 41% of the miRNAs, which differed by genotype (25/61) and 48% of the androgen-sensitive miRNAs (20/42), indicating overlapping genetic and dietary influences on prostate miRNAs. Tomato and lycopene feeding influenced miRNAs previously associated with the regulation of androgen (miR-145 and let-7), MAPK (miR-106a, 204, 145/143, and 200b/c), and p53 signaling (miR-125 and miR-98) pathways. CONCLUSIONS: Expression of miRNAs in early prostate carcinogenesis is sensitive to genetic, endocrine, and diet drivers, suggesting novel mechanisms by which tomato and lycopene feeding modulate early prostate carcinogenesis.

Control Theory Forecasts of Optimal Training Dosage to Facilitate Children's Arithmetic Learning in a Digital Educational Application.

Education can be viewed as a control theory problem in which students seek ongoing exogenous input-either through traditional classroom teaching or other alternative training resources-to minimize the discrepancies between their actual and target (reference) performance levels. Using illustrative data from [Formula: see text] Dutch elementary school students as measured using the Math Garden, a web-based computer adaptive practice and monitoring system, we simulate and evaluate the outcomes of using off-line and finite memory linear quadratic controllers with constraintsto forecast students' optimal training durations. By integrating population standards with each student's own latent change information, we demonstrate that adoption of the control theory-guided, person- and time-specific training dosages could yield increased training benefits at reduced costs compared to students' actual observed training durations, and a fixed-duration training scheme. The control theory approach also outperforms a linear scheme that provides training recommendations based on observed scores under noisy and the presence of missing data. Design-related issues such as ways to determine the penalty cost of input administration and the size of the control horizon window are addressed through a series of illustrative and empirically (Math Garden) motivated simulations.

Evolutionary conservation and structural localizations suggest a physical trace of metabolism's progressive geochronological emergence.

We studied multiple sequence alignment (MSA) consensus amino acid distributional patterns in 2844 amino acid sequences of the eight enzymes of the Kreb's oxidative tricarboxylic acid pathway (oTCA) in Archaea, Bacteria and Eukarya and 5545 sequences of 33 bacteria as geochronologically separated enzymes with MSA consensus site modal identities. The 33 bacteria were 20 presumptive examples of early-oldest (Hadean-Archaean) ('Epoch I') or 13 late-newest (contemporary) ('Epoch III') appearing enzymes on Earth. The enzyme's MSA consensus sites were identified by their modal identity, % Occupancy in one of nine-graded evolutionary-conservation zones (CZs) and the Euclidean distance (Å) from each of their consensus MSA Cɑs to the same atom (Anchor-atom) in their reported functional center. These MSA consensus sites are tetrad-data points called recovered-amino acids (RAA). Across Domains, the % Occupancies of the eight-dominant RAAs of the Kreb's cycle and the 33 bacteria were found to be similarly ranked. Compared to Trifonov's 'putative ranked temporal order of the appearance of amino acids on Earth' (TOAE), the greatest statistical concordance with tetrad-RAAs across Domains were those characterized as within the most-evolutionary conserved conservation zone (CZ9), typically nearest (Å) their enzyme's catalytic/active center. The geochronologically characterized early-oldest Hadean-Archaean Bacteria 'Epoch I' enzymes, compared to late-newest Bacteria enzymes, had greater average numbers of amino acid residues/sequence and a statistically significant larger variability in their RAA compositional-Å3-volumes. The late-newest 'Epoch III' enzymes had statistically significant lower volumetric values, specifically, their native Å3-volume, void-volume and volume change on unfolding. Our enzyme data suggest a geochronological trace of 'metabolism's progressive emergence'.Communicated by Ramaswamy H. Sarma.

Probability Distributome: A Web Computational Infrastructure for Exploring the Properties, Interrelations, and Applications of Probability Distributions.

Probability distributions are useful for modeling, simulation, analysis, and inference on varieties of natural processes and physical phenomena. There are uncountably many probability distributions. However, a few dozen families of distributions are commonly defined and are frequently used in practice for problem solving, experimental applications, and theoretical studies. In this paper, we present a new computational and graphical infrastructure, the Distributome, which facilitates the discovery, exploration and application of diverse spectra of probability distributions. The extensible Distributome infrastructure provides interfaces for (human and machine) traversal, search, and navigation of all common probability distributions. It also enables distribution modeling, applications, investigation of inter-distribution relations, as well as their analytical representations and computational utilization. The entire Distributome framework is designed and implemented as an open-source, community-built, and Internet-accessible infrastructure. It is portable, extensible and compatible with HTML5 and Web2.0 standards (http://Distributome.org). We demonstrate two types of applications of the probability Distributome resources: computational research and science education. The Distributome tools may be employed to address five complementary computational modeling applications (simulation, data-analysis and inference, model-fitting, examination of the analytical, mathematical and computational properties of specific probability distributions, and exploration of the inter-distributional relations). Many high school and college science, technology, engineering and mathematics (STEM) courses may be enriched by the use of modern pedagogical approaches and technology-enhanced methods. The Distributome resources provide enhancements for blended STEM education by improving student motivation, augmenting the classical curriculum with interactive webapps, and overhauling the learning assessment protocols.

Examining Gender Differences in Written Assessment Tasks in Biology: A Case Study of Evolutionary Explanations.

Understanding sources of performance bias in science assessment provides important insights into whether science curricula and/or assessments are valid representations of student abilities. Research investigating assessment bias due to factors such as instrument structure, participant characteristics, and item types are well documented across a variety of disciplines. However, the relationships among these factors are unclear for tasks evaluating understanding through performance on scientific practices, such as explanation. Using item-response theory (Rasch analysis), we evaluated differences in performance by gender on a constructed-response (CR) assessment about natural selection (ACORNS). Three isomorphic item strands of the instrument were administered to a sample of undergraduate biology majors and nonmajors (Group 1: n = 662 [female = 51.6%]; G2: n = 184 [female = 55.9%]; G3: n = 642 [female = 55.1%]). Overall, our results identify relationships between item features and performance by gender; however, the effect is small in the majority of cases, suggesting that males and females tend to incorporate similar concepts into their CR explanations. These results highlight the importance of examining gender effects on performance in written assessment tasks in biology.

Suppression of Proinflammatory and Prosurvival Biomarkers in Oral Cancer Patients Consuming a Black Raspberry Phytochemical-Rich Troche.

Black raspberries (BRB) demonstrate potent inhibition of aerodigestive tract carcinogenesis in animal models. However, translational clinical trials evaluating the ability of BRB phytochemicals to impact molecular biomarkers in the oral mucosa remain limited. The present phase 0 study addresses a fundamental question for oral cancer food-based prevention: Do BRB phytochemicals successfully reach the targeted oral tissues and reduce proinflammatory and antiapoptotic gene expression profiles? Patients with biopsy-confirmed oral squamous cell carcinomas (OSCC) administered oral troches containing freeze-dried BRB powder from the time of enrollment to the date of curative intent surgery (13.9 ± 1.27 days). Transcriptional biomarkers were evaluated in patient-matched OSCCs and noninvolved high at-risk mucosa (HARM) for BRB-associated changes. Significant expression differences between baseline OSCC and HARM tissues were confirmed using a panel of genes commonly deregulated during oral carcinogenesis. Following BRB troche administration, the expression of prosurvival genes (AURKA, BIRC5, EGFR) and proinflammatory genes (NFKB1, PTGS2) were significantly reduced. There were no BRB-associated grade 3-4 toxicities or adverse events, and 79.2% (N = 30) of patients successfully completed the study with high levels of compliance (97.2%). The BRB phytochemicals cyanidin-3-rutinoside and cyanidin-3-xylosylrutinoside were detected in all OSCC tissues analyzed, demonstrating that bioactive components were successfully reaching targeted OSCC tissues. We confirmed that hallmark antiapoptotic and proinflammatory molecular biomarkers were overexpressed in OSCCs and that their gene expression was significantly reduced following BRB troche administration. As these molecular biomarkers are fundamental to oral carcinogenesis and are modifiable, they may represent emerging biomarkers of molecular efficacy for BRB-mediated oral cancer chemoprevention.

Dietary tomato and lycopene impact androgen signaling- and carcinogenesis-related gene expression during early TRAMP prostate carcinogenesis.

Consumption of tomato products containing the carotenoid lycopene is associated with a reduced risk of prostate cancer. To identify gene expression patterns associated with early testosterone-driven prostate carcinogenesis, which are impacted by dietary tomato and lycopene, wild-type (WT) and transgenic adenocarcinoma of the mouse prostate (TRAMP) mice were fed control or tomato- or lycopene-containing diets from 4 to 10 weeks of age. Eight-week-old mice underwent sham surgery, castration, or castration followed by testosterone repletion (2.5 mg/kg/d initiated 1 week after castration). Ten-week-old intact TRAMP mice exhibit early multifocal prostatic intraepithelial neoplasia. Of the 200 prostate cancer-related genes measured by quantitative NanoString, 189 are detectable, 164 significantly differ by genotype, 179 by testosterone status, and 30 by diet type (P < 0.05). In TRAMP, expression of Birc5, Mki67, Aurkb, Ccnb2, Foxm1, and Ccne2 is greater compared with WT and is decreased by castration. In parallel, castration reduces Ki67-positive staining (P < 0.0001) compared with intact and testosterone-repleted TRAMP mice. Expression of genes involved in androgen metabolism/signaling pathways is reduced by lycopene feeding (Srd5a1) and by tomato feeding (Srd5a2, Pxn, and Srebf1). In addition, tomato feeding significantly reduced expression of genes associated with stem cell features, Aldh1a and Ly6a, whereas lycopene feeding significantly reduced expression of neuroendocrine differentiation-related genes, Ngfr and Syp. Collectively, these studies demonstrate a profile of testosterone-regulated genes associated with early prostate carcinogenesis that are potential mechanistic targets of dietary tomato components. Future studies on androgen signaling/metabolism, stem cell features, and neuroendocrine differentiation pathways may elucidate the mechanisms by which dietary tomato and lycopene impact prostate cancer risk.

The selective Aurora-A kinase inhibitor MLN8237 (alisertib) potently inhibits proliferation of glioblastoma neurosphere tumor stem-like cells and potentiates the effects of temozolomide and ionizing radiation.

The selective Aurora-A kinase inhibitor MLN8237 is in clinical trials for hematologic malignancies, ovarian cancer and other solid tumors. We previously showed that MLN8237 is potently antiproliferative toward standard monolayer-cultured glioblastoma cells. We have now investigated the effect of MLN8237 with and without temozolomide or ionizing radiation on the proliferation of glioblastoma tumor stem-like cells (neurospheres) using soft agar colony formation assays and normal human astrocytes by MTT assay. Western blotting was utilized to compare MLN8237 IC50s to cellular Aurora-A and phosphoThr(288)Aurora-A levels. MLN8237 was more potently antiproliferative to neurosphere cells than to standard monolayer glioma cells, and was non-toxic to normal human astrocytes. Western blot analysis revealed that MLN8237 treatment inhibits phosphoThr(288)Aurora-A levels providing proof of drug target-hit in glioblastoma cells. Furthermore, phosphoThr(288)Aurora-A levels partially predicted the antiproliferative efficacy of MLN8237. We also found that Aurora-A inhibition by MLN8237 was synergistic with temozolomide and potentiated the effects of ionizing radiation on colony formation in neurosphere glioblastoma tumor stem-like cells. These results further support the potential of Aurora-A inhibitors as primary chemotherapy agents or biologic response modifiers in glioblastoma patients.

ß-Carotene-9',10'-oxygenase status modulates the impact of dietary tomato and lycopene on hepatic nuclear receptor-, stress-, and metabolism-related gene expression in mice.

Tomato and lycopene (ψ,ψ-carotene) consumption is hypothesized to protect against nonalcoholic steatohepatitis and hepatocarcinogenesis, processes that may depend upon diet and gene interactions. To investigate the interaction of tomato or lycopene feeding with ß-carotene-9',10'-monooxygenase (Bco2) on hepatic metabolic and signaling pathways, male wild-type (WT) and Bco2(-/-) mice (3-wk-old; n = 36) were fed semi-purified control, 10% tomato powder-containing, or 0.25% lycopene beadlet-containing diets for 3 wk. Serum lycopene concentrations were higher in lycopene- and tomato-fed Bco2(-/-) mice compared with WT (P = 0.03). Tomato- and lycopene-fed mice had detectable hepatic apolipoprotein (apo)-6'-, apo-8'-, and apo-12'-lycopenal concentrations. Hepatic expression of ß-carotene-15,15'-monooxygenase was increased in Bco2(-/-) mice compared with WT (P = 0.02), but not affected by diet. Evaluation of hepatic gene expression by focused quantitative reverse transcriptase-polymerase chain reaction arrays for nuclear receptors and coregulators (84 genes) and stress and metabolism (82 genes) genes indicates that tomato feeding affected 31 genes (≥1.5-fold, P < 0.05) and lycopene feeding affected 19 genes, 16 of which were affected by both diets. Lycopene down-regulation of 7 nuclear receptors and coregulators, estrogen-related receptor-α, histone deacetylase 3, nuclear receptor coactivator 4, RevErbA-ß, glucocorticoid receptor, peroxisome proliferator-activated receptor (PPAR)-α, and PPAR-γ, coactivator 1 ß was dependent upon interaction with Bco2 status. Lycopene and tomato feeding induced gene expression patterns consistent with decreased lipid uptake, decreased cell proliferation and mitosis, down-regulated aryl hydrocarbon receptor signaling, and decreased expression of genes involved in retinoid X receptor heterodimer activation. Tomato feeding also caused expression changes consistent with down-regulation of DNA synthesis and terpenoid metabolism. These data suggest tomato components, particularly lycopene, affect hepatic gene expression, potentially affecting hepatic responses to metabolic, infectious, or chemical stress.

The choice of reference gene affects statistical efficiency in quantitative PCR data analysis.

Quantitative polymerase chain reaction (qPCR), a highly sensitive method of measuring gene expression, is widely used in biomedical research. To produce reliable results, it is essential to use stably expressed reference genes (RGs) for data normalization so that sample-to-sample variation can be controlled. In this study, we examine the effect of different RGs on statistical efficiency by analyzing a qPCR data set that contains 12 target genes and 3 RGs. Our results show that choosing the most stably expressed RG for data normalization does not guarantee reduced variance or improved statistical efficiency. We also provide a formula for determining when data normalization will improve statistical efficiency and hence increase the power of statistical tests in data analysis.

Uniquely localized intra-molecular amino acid concentrations at the glycolytic enzyme catalytic/active centers of Archaea, Bacteria and Eukaryota are associated with their proposed temporal appearances on earth.

The distributions of amino acids at most-conserved sites nearest catalytic/active centers (C/AC) in 4,645 sequences of ten enzymes of the glycolytic Embden-Meyerhof-Parnas pathway in Archaea, Bacteria and Eukaryota are similar to the proposed temporal order of their appearance on Earth. Glycine, isoleucine, leucine, valine, glutamic acid and possibly lysine often described as prebiotic, i.e., existing or occurring before the emergence of life, were localized in positional and conservational defined aggregations in all enzymes of all Domains. The distributions of all 20 biologic amino acids in most-conserved sites nearest their C/ACs were quite different either from distributions in sites less-conserved and further from their C/ACs or from all amino acids regardless of their position or conservation. The major concentrations of glycine, e.g., perhaps the earliest prebiotic amino acid, occupies ≈ 16 % of all the most-conserved sites within a volume of ≈ 7-8 Å radius from their C/ACs and decreases linearly towards the molecule's peripheries. Spatially localized major concentrations of isoleucine, leucine and valine are in the mid-conserved and mid-distant sites from their C/ACs in protein interiors. Lysine and glutamic acid comprise ≈ 25-30 % of all amino acids within an irregular volume bounded by ≈ 24-28 Å radii from their C/ACs at the most-distant least-conserved sites. The unreported characteristics of these amino acids: their spatially and conservationally identified concentrations in Archaea, Bacteria and Eukaryota, suggest some common structural organization of glycolytic enzymes that may be relevant to their evolution and that of other proteins. We discuss our data in relation to enzyme evolution, their reported prebiotic putative temporal appearances on Earth, abundances, biological "cost", neighbor-sequence preferences or "ordering" and some thermodynamic parameters.

miR-181b is a biomarker of disease progression in chronic lymphocytic leukemia.

MicroRNAs play a crucial role in chronic lymphocytic leukemia. We investigated whether microRNAs can discriminate patients with a progressive disease from patients with a stable disease. We analyzed microRNA expression on leukemic cells isolated from 358 sequential samples of 114 patients with either stable or progressive disease. We found that during the course of the disease the expression values of miR-181b, the most dysregulated microRNA, decreased in samples of patients with a progressive (P < .001, training and validation sets) but not in samples of patients with a stable disease (P = .3, training set; P = .2, validation set) over time. A drop of ≥ 50% between sequential samples and/or a miR-181b value ≤ 0.005 at the starting time point were significant to differentiate progressive from stable disease (P = .004, training set; P < .001, validation set). These parameters were associated with high risk of requiring treatment (risk ratio, 5.8; 95% confidence interval, 2.5-14.9). We also observed that miR-181b targets Mcl-1 protein and that the decrease of its expression inversely correlated with increased protein levels of MCL1 and BCL2 target genes. We conclude that parameters defined on the basis of the miR-181b expression values specify disease progression in chronic lymphocytic leukemia and are associated with clinical outcome.

Newborn and carrier screening for spinal muscular atrophy.

Spinal muscular atrophy (SMA) is a common autosomal recessive neuromuscular disorder caused by mutations in the survival motor neuron (SMN1) gene, affecting approximately 1 in 10,000 live births. The homozygous absence of SMN1 exon 7 has been observed in the majority of patients and is being utilized as a reliable and sensitive SMA diagnostic test. Treatment and prevention of SMA are complementary responses to the challenges presented by SMA. Even though a specific therapy for SMA is not currently available, a newborn screening test may allow the child to be enrolled in a clinical trial before irreversible neuronal loss occurs and enable patients to obtain more proactive treatments. Until an effective treatment is found to cure or arrest the progression of the disease, prevention of new cases through accurate diagnosis and carrier and prenatal diagnosis is of the utmost importance. The goal of population-based SMA carrier screening is to identify couples at risk for having a child with SMA, thus allowing carriers to make informed reproductive choices. During this study we performed two pilot projects addressing the clinical applicability of testing in the newborn period and carrier screening in the general population. We have demonstrated that an effective technology does exist for newborn screening of SMA. We also provide an estimate of the carrier frequency among individuals who accepted carrier screening, and report on patient's knowledge and attitudes toward SMA testing.

Concentration of specific amino acids at the catalytic/active centers of highly-conserved "housekeeping" enzymes of central metabolism in archaea, bacteria and Eukaryota: is there a widely conserved chemical signal of prebiotic assembly?

In alignments of 1969 protein sequences the amino acid glycine and others were found concentrated at most-conserved sites within approximately 15 A of catalytic/active centers (C/AC) of highly conserved kinases, dehydrogenases or lyases of Archaea, Bacteria and Eukaryota. Lysine and glutamic acid were concentrated at least-conserved sites furthest from their C/ACs. Logistic-regression analyses corroborated the "movement" of glycine towards and lysine away from their C/ACs: the odds of a glycine occupying a site were decreased by 19%, while the odds for a lysine were increased by 53%, for every 10 A moving away from the C/AC. Average conservation of MSA consensus sites was highest surrounding the C/AC and directly decreased in transition toward model's peripheries. Findings held with statistical confidence using sequences restricted to individual Domains or enzyme classes or to both. Our data describe variability in the rate of mutation and likelihoods for phylogenetic trees based on protein sequence data and endorse the extension of substitution models by incorporating data on conservation and distance to C/ACs rather than only using cumulative levels. The data support the view that in the most-conserved environment immediately surrounding the C/AC of taxonomically distant and highly conserved essential enzymes of central metabolism there are amino acids whose identity and degree of occupancy is similar to a proposed amino acid set and frequency associated with prebiotic evolution.

Maximum tree: a consistent estimator of the species tree.

We propose a model based approach to use multiple gene trees to estimate the species tree. The coalescent process requires that gene divergences occur earlier than species divergences when there is any polymorphism in the ancestral species. Under this scenario, speciation times are restricted to be smaller than the corresponding gene split times. The maximum tree (MT) is the tree with the largest possible speciation times in the space of species trees restricted by available gene trees. If all populations have the same population size, the MT is the maximum likelihood estimate of the species tree. It can be shown the MT is a consistent estimator of the species tree even when the MT is built upon the estimates of the true gene trees if the gene tree estimates are statistically consistent. The MT converges in probability to the true species tree at an exponential rate.

Coalescent methods for estimating phylogenetic trees.

We review recent models to estimate phylogenetic trees under the multispecies coalescent. Although the distinction between gene trees and species trees has come to the fore of phylogenetics, only recently have methods been developed that explicitly estimate species trees. Of the several factors that can cause gene tree heterogeneity and discordance with the species tree, deep coalescence due to random genetic drift in branches of the species tree has been modeled most thoroughly. Bayesian approaches to estimating species trees utilizes two likelihood functions, one of which has been widely used in traditional phylogenetics and involves the model of nucleotide substitution, and the second of which is less familiar to phylogeneticists and involves the probability distribution of gene trees given a species tree. Other recent parametric and nonparametric methods for estimating species trees involve parsimony criteria, summary statistics, supertree and consensus methods. Species tree approaches are an appropriate goal for systematics, appear to work well in some cases where concatenation can be misleading, and suggest that sampling many independent loci will be paramount. Such methods can also be challenging to implement because of the complexity of the models and computational time. In addition, further elaboration of the simplest of coalescent models will be required to incorporate commonly known issues such as deviation from the molecular clock, gene flow and other genetic forces.

Sphingosine-1-phosphate regulates glioblastoma cell invasiveness through the urokinase plasminogen activator system and CCN1/Cyr61.

Glioblastoma multiforme (GBM) is an aggressively invasive brain neoplasm with poor patient prognosis. We have previously shown that the bioactive lipid sphingosine-1-phosphate (S1P) stimulates in vitro invasiveness of GBM cells and that high expression levels of the enzyme that forms S1P, sphingosine kinase-1 (SphK1), correlate with shorter survival time of GBM patients. We also recently showed that S1P induces expression of CCN1 (also known as Cyr61), a matricellular protein known to correlate with poor patient prognosis, in GBM cells. In this study, we further explored the role of CCN1 as well as the urokinase plasminogen activator (uPA), a protein known to stimulate GBM cell invasiveness, in S1P-induced invasion using a spheroid invasion assay. We also investigated the roles of various S1P receptors in stimulating invasiveness through these pathways. S1P induced expression of uPA and its receptor, uPAR, in GBM cells. Whereas S1P(1), S1P(2), and S1P(3) receptors all contribute, at least partially, S1P(1) overexpression led to the most dramatic induction of the uPA system and of spheroid invasion, even in the absence of added S1P. Furthermore, neutralizing antibodies directed against uPA or CCN1 significantly decreased both basal and S1P-stimulated GBM cell invasiveness. Inhibition of SphK blocked basal expression of uPA and uPAR, as well as glioma cell invasion; however, overexpression of SphK did not augment S1P receptor-mediated enhancement of uPA activity or invasion. Thus, SphK is necessary for basal activity of the uPA system and glioma cell invasion, whereas S1P receptor signaling enhances invasion, partially through uPA and CCN1.

Estimating species phylogenies using coalescence times among sequences.

The estimation of species trees (phylogenies) is one of the most important problems in evolutionary biology, and recently, there has been greater appreciation of the need to estimate species trees directly rather than using gene trees as a surrogate. A Bayesian method constructed under the multispecies coalescent model can consistently estimate species trees but involves intensive computation, which can hinder its application to the phylogenetic analysis of large-scale genomic data. Many summary statistics-based approaches, such as shallowest coalescences (SC) and Global LAteSt Split (GLASS), have been developed to infer species phylogenies for multilocus data sets. In this paper, we propose 2 methods, species tree estimation using average ranks of coalescences (STAR) and species tree estimation using average coalescence times (STEAC), based on the summary statistics of coalescence times. It can be shown that the 2 methods are statistically consistent under the multispecies coalescent model. STAR uses the ranks of coalescences and is thus resistant to variable substitution rates along the branches in gene trees. A simulation study suggests that STAR consistently outperforms STEAC, SC, and GLASS when the substitution rates among lineages are highly variable. Two real genomic data sets were analyzed by the 2 methods and produced species trees that are consistent with previous results.

Estimating species trees using multiple-allele DNA sequence data.

Several techniques, such as concatenation and consensus methods, are available for combining data from multiple loci to produce a single statement of phylogenetic relationships. However, when multiple alleles are sampled from individual species, it becomes more challenging to estimate relationships at the level of species, either because concatenation becomes inappropriate due to conflicts among individual gene trees, or because the species from which multiple alleles have been sampled may not form monophyletic groups in the estimated tree. We propose a Bayesian hierarchical model to reconstruct species trees from multiple-allele, multilocus sequence data, building on a recently proposed method for estimating species trees from single allele multilocus data. A two-step Markov Chain Monte Carlo (MCMC) algorithm is adopted to estimate the posterior distribution of the species tree. The model is applied to estimate the posterior distribution of species trees for two multiple-allele datasets--yeast (Saccharomyces) and birds (Manacus-manakins). The estimates of the species trees using our method are consistent with those inferred from other methods and genetic markers, but in contrast to other species tree methods, it provides credible regions for the species tree. The Bayesian approach described here provides a powerful framework for statistical testing and integration of population genetics and phylogenetics.