Species co‐occurrence analysis: pairwise versus matrix‐level approaches

Veech (2013, Global Ecology and Biogeography, 22 , 252–260) introduced a formula to calculate the probability of two species co‐occurring in various sites under the assumption of statistical independence between the two distributional patterns. He presented his model as a new procedure, a ‘pairwise approach’, different from analyses of whole presence–absence matrices to examine patterns of co‐occurrence. Here I show that: (1) Veech's method is identical to Fisher's exact test, a standard procedure for measuring the statistical association between two discrete variables; (2) in a broad sense, the pairwise approach is very similar to early analyses of spatial association, such as the one advanced by Forbes in 1907; (3) implicit in Veech's formula is a sampling scheme that is indistinguishable from well‐known matrix‐level null models that randomize the distribution of species among equiprobable sites; (4) pairwise co‐occurrence patterns can be analysed using any matrix‐level null model, so pairwise comparisons are not limited to using Veech's formula. The methodological distinction that Veech proposed between pairwise and matrix‐level approaches does not in fact exist, although the conceptual distinction between the two approaches is still a debated topic.     

Phylogeography of the freshwater red alga B atrachospermum viride‐brasiliense (Rhodophyta,Batrachospermales)

Phylogeography of B atrachospermum viride‐brasiliense was investigated using two mitochondrial regions: the cox2‐3 spacer and the barcode region of cox1 gene. Eighty‐seven individuals were analyzed from nine stream segments throughout its distribution in Brazil. Ten cox2‐3 spacer and nine cox1 haplotypes were observed among the individuals studied (87 vs. 43, respectively). Divergences among haplotypes were relatively low (≤2.4% for cox2‐3 and ≤1.8% for cox1). Most locations have a single haplotype, whereas only two locations had two haplotypes for both markers. The haplotype network for cox2‐3 showed a phylogeographic trend from the south towards the southeast with haplotypes from the southeast more closely related. For cox1 a trend from the southeast spreading towards the south and north was revealed, with the southern haplotypes more closely associated. Results clearly indicated that B . viride‐brasiliense represents a single species and the phylogeographic pattern consisted of a closely connected group of haplotypes from southern and southeastern Brazil. Levels of intra‐ and inter‐population variation were similar for the two markers with slightly higher values for cox2‐3. The trend observed in this study is similar to that in other members of Batrachospermales with little variation within a stream segment (one or two haplotypes) and more distant haplotypes showing higher divergences. This pattern could be attributed to the fact that colonization of a site might be rare by a single event with subsequent proliferation of the population. The geographic distribution of B . viride‐brasiliense was interpreted according to the biogeographic models proposed for South America being limited to three morpho‐climatic domains or biogeographic provinces: tropical Atlantic rainforest, sub‐tropical rainforest and cerrado (Brazilian savannah).     

Differences in population structure estimated within maternally‐ and paternally‐inherited forms of mitochondria in Lampsilis siliquoidea (Bivalvia: Unionidae)

Mussels in several orders possess two separate mitochondrial lineages: a standard female‐inherited form and one inherited only through males. This system of doubly uniparental inheritance (DUI) for mitochondrial genes provides an opportunity to compare the population structure of gene‐lineages passed either mother‐to‐daughter or father‐to‐son. In the present study, we contrast variation in the male and female haplotype lineages of the American freshwater mussel species, Lampsilis siliquoidea (sometimes called Lampsilis radiata luteola), throughout the Lake Erie, Ohio River, and upper Mississippi River watersheds, and contrast variation with the sequences obtained for the related species/subspecies Lampsilis radiata radiata from Maine. The genetic markers were fragments of the cytochrome c oxidase subunit I gene (COI), which occurs in both mitochondrial types, F (female) and M (male). High haplotype diversity was found in the two independent lineages, although purifying selection against amino acid change appeared to be stronger in the female than the male lineage. Phylogeographical patterns also varied between mitochondria passing through females and males. The female lineage exhibited more population structure, with the occurrence of private or nearly‐private haplotypes within two streams, and three others showed restricted haplotype distributions. By contrast to the F‐haplotypes, complex phylogenetic structure occurred for M‐haplotypes, yet this phylogenetic variation coincided with almost no geographical pattern within haplotypes. Basically, F‐haplotypes showed isolation, especially above physical barriers, whereas M‐haplotypes did not. A few individuals in the eastern Lake Erie watershed even possessed M‐haplotypes of an Atlantic Slope (L. radiata radiata) origin, although their F‐haplotypes were typical of Midwestern L. siliquoidea. The finding that mussels package sperm as spermatozuegmata, which float downstream, may underlie greater gene mobility in male‐inherited mitochondria. © 2013 The Linnean Society of London, Biological Journal of the Linnean Society, 2013, 109 , 229–240.     

Isolation and characterization of microsatellite loci in Amazonian red‐handed howlers Alouatta belzebul (Primates,Plathyrrini)

Ten highly polymorphic microsatellite loci are described for Amazonian red‐handed howlers (Alouatta belzebul), an endemic Brazilian primate species subject to intense hunting pressure. The number of alleles observed in 30 individuals ranged from nine to 20, and observed heterozygosity varied from 0.35 to 0.93. No linkage associations were evident from pairwise comparisons of loci. These microsatellites offer a powerful tool for fine‐scale studies of genetic structure in both captive colonies and wild populations of red‐handed howlers.     

Within‐population polymorphism of sex‐determination systems in the common frog (Rana temporaria)

In sharp contrast with birds and mammals, the sex chromosomes of ectothermic vertebrates are often undifferentiated, for reasons that remain debated. A linkage map was recently published for Rana temporaria (Linnaeus, 1758) from Fennoscandia (Eastern European lineage), with a proposed sex‐determining role for linkage group 2 (LG2). We analysed linkage patterns in lowland and highland populations from Switzerland (Western European lineage), with special focus on LG2. Sibship analyses showed large differences from the Fennoscandian map in terms of recombination rates and loci order, pointing to large‐scale inversions or translocations. All linkage groups displayed extreme heterochiasmy (total map length was 12.2 cM in males, versus 869.8 cM in females). Sex determination was polymorphic within populations: a majority of families (with equal sex ratios) showed a strong correlation between offspring phenotypic sex and LG2 paternal haplotypes, whereas other families (some of which with female‐biased sex ratios) did not show any correlation. The factors determining sex in the latter could not be identified. This coexistence of several sex‐determination systems should induce frequent recombination of X and Y haplotypes, even in the absence of male recombination. Accordingly, we found no sex differences in allelic frequencies on LG2 markers among wild‐caught male and female adults, except in one high‐altitude population, where nonrecombinant Y haplotypes suggest sex to be entirely determined by LG2. Multifactorial sex determination certainly contributes to the lack of sex‐chromosome differentiation in amphibians.     

Association between polygenic liability for schizophrenia and substance involvement: A nationwide population‐based study in Taiwan

Schizophrenia and substance involvement frequently co‐occur in individuals, and a bidirectional relationship between the two has been proposed; shared underlying genetic factors could be an alternative explanation. This study investigated the genetic overlap between schizophrenia and substance involvement, including tobacco, alcohol and betel nut use. The study subjects were recruited from the Taiwan Biobank, and genome‐wide genotyping data was available for 18 327 participants without schizophrenia. We calculated the Psychiatric Genomics Consortium‐derived polygenic risk score (PRS) for schizophrenia in each participant. The significance of the schizophrenia PRS associated with substance involvement was evaluated using a regression model with adjustments for gender, age and population stratification components. The modified effect of gender or birth decade was also explored. The schizophrenia PRS was positively associated with lifetime tobacco smoking in women (OR in per SD increase in PRS = 1.12 with 95% CI 1.04‐1.20, P = .002), but not in men (OR = 0.99 with 95% CI 0.95‐1.04, P = .74), and the gender‐PRS interaction reached significance (P = .006). The OR between PRS and lifetime tobacco smoking increased with the birth decade (P of birth decade‐PRS interaction = .0002). In women, OR increased from 0.97 (P = .85) for subjects with a birth decade before 1950 to 1.21 (P = .04) for subjects with a birth decade after 1980; in men, the corresponding OR increased from 0.88 (P = .04) to 1.13 (P = .11). There was no association between schizophrenia PRS and alcohol/betel nut use phenotypes. This study provides evidence for the genetic overlap between schizophrenia and tobacco use in women, and this overlap was stronger in the younger population.     

Y‐specific microsatellites reveal an African subfamily in taurine (Bos taurus) cattle

Five cattle Y‐specific microsatellites, totalling six loci, were selected from a set of 44 markers and genotyped on 608 Bos taurus males belonging to 45 cattle populations from Europe and Africa. A total of 38 haplotypes were identified. Haplogroups (Y1 and Y2) previously defined using single nucleotide polymorphisms did not share haplotypes. Nine of the 27 Y2‐haplotypes were only present in African cattle. Network and correspondence analyses showed that this African‐specific subfamily clustered separately from the main Y2‐subfamily and the Y1 haplotypes. Within‐breed genetic variability was generally low, with most breeds (78%) showing haplotypes belonging to a single haplogroup. amova analysis showed that partitioning of genetic variation among breeds can be mainly explained by their geographical and haplogroup assignment. Between‐breed genetic variability summarized via Principal Component Analysis allowed the identification of three principal components explaining 94.2% of the available information. Projection of principal components on geographical maps illustrated that cattle populations located in mainland Europe, the three European Peninsulas and Mediterranean Africa presented similar genetic variation, whereas those breeds from Atlantic Europe and British Islands (mainly carrying Y1 haplotypes) and those from Sub‐Saharan Africa (belonging to Y2‐haplogroup) showed genetic variation of a different origin. Our study confirmed the existence of two large Y‐chromosome lineages (Y1 and Y2) in taurine cattle. However, Y‐specific microsatellites increased analytical resolution and allowed at least two different Y2‐haplotypic subfamilies to be distinguished, one of them restricted to the African continent.     

Population structure,population history and DNA barcoding of fruit fly Bactrocera latifrons (Hendel) (Diptera: Tephritidae)

The fruit fly Bactrocera latifrons (Hendel) is an important pest of commercially significant plants such as chili, tomato and eggplant. The species is native to South and Southeast Asia, but has now invaded Japan, Hawaii and Africa. In this study, mitochondrial DNA sequences were used to infer genetic structure and demographic history of B. latifrons. The efficiency of DNA barcodes for identification of B. latifrons was also tested. Ninety‐three specimens infesting four host‐plant species were obtained from 11 sampling locations in Thailand. The mitochondrial haplotype network revealed no major divergent lineage, which was consistent with a phylogenetic analysis that found strong support for the monophyly of B. latifrons. Population pairwise F_ST revealed that most (65%) comparisons were not significantly different, suggesting a high rate of gene flow. Analysis of molecular variance (amova ) found no significant genetic differentiation among populations from different host‐plant species. Sharing of several haplotypes among flies from different host‐plants indicates that the flies were moved freely across the plant species. Demographic history analysis revealed that the population has undergone recent expansion dating back to the end of the last glaciation. Thus, the results indicate that both ongoing and historical factors have played important roles in determining the genetic structure and diversity of B. latifrons. DNA barcoding analysis revealed that B. latifrons specimens were clearly differentiated from other species with 100% correct identification. Therefore, cytochrome oxidase I (COI) barcoding sequences could be effectively used to identify this important pest species, which could encourage monitoring and control efforts for this species.     

Microsatellite DNA markers for the study of Atlantic salmon (Salmo salar) kinship,population structure,and mixed‐fishery analyses

Eleven microsatellite DNA loci were identified and characterized for Atlantic salmon (Salmo salar) collected from the Penobscot River, Maine, USA and the River Nith, Scotland, UK. The markers revealed high levels of genetic diversity (seven to 48 alleles per locus), heterozygosity (to 100%), and allelic heterogeneity (all comparisons). Considerable differentiation was observed as the genetic distance (chord) between the two collections was 0.680 and the pairwise F_ST, 0.12, was highly significant. These findings are consistent with patterns of continental‐level differentiation observed previously using an alternate suite of microsatellite loci. Locus‐by‐locus analyses of molecular variance suggested that most markers were suitable for delineating kinships and population genetic structure.     

Nonparametric All‐Pairs Multiple Comparisons

Nonparametric all‐pairs multiple comparisons based on pairwise rankings can be performed in the one‐way design with the Steel‐Dwass procedure. To apply this test, Wilcoxon's rank sum statistic is calculated for all pairs of groups; the maximum of the rank sums is the test statistic. We provide exact calculations of the asymptotic critical values (and P‐values, respectively) even for unbalanced designs. We recommend this asymptotic method whenever large sample sizes are present. For small sample sizes we recommend the use of the new statistic according to Baumgartner , Weiss , and Schindler (1998, Biometrics 54 , 1129–1135) instead of Wilcoxon's rank sum for the multiple comparisons. We show that the resultant procedure can be less conservative and, according to simulation results, more powerful than the original Steel‐Dwass procedure. We illustrate the methods with a practical data set.     

Genetic structure and demographic history of Lymantria dispar (Linnaeus, 1758) (Lepidoptera: Erebidae) in its area of origin and adjacent areas

We analyzed the population genetic structure and demographic history of 20 Lymantria dispar populations from Far East Asia using microsatellite loci and mitochondrial genes. In the microsatellite analysis, the genetic distances based on pairwise F_ST values ranged from 0.0087 to 0.1171. A NeighborNet network based on pairwise F_ST genetic distances showed that the 20 regional populations were divided into five groups. Bayesian clustering analysis (K = 3) demonstrated the same groupings. The populations in the Korean Peninsula and adjacent regions, in particular, showed a mixed genetic pattern. In the mitochondrial genetic analysis based on 98 haplotypes, the median‐joining network exhibited a star shape that was focused on three high‐frequency haplotypes (Haplotype 1: central Korea and adjacent regions, Group 1; Haplotype 37: southern Korea, Group 2; and Haplotype 90: Hokkaido area, Group 3) connected by low‐frequency haplotypes. The mismatch distribution dividing the three groups was unimodal. In the neutral test, Tajima's D and Fu's FS tests were negative. We can thus infer that the Far East Asian populations of L. dispar underwent a sudden population expansion. Based on the age expansion parameter, the expansion time was inferred to be approximately 53,652 years before present (ybp) for Group 1, approximately 65,043 ybp for Group 2, and approximately 76,086 ybp for Group 3. We propose that the mixed genetic pattern of the inland populations of Far East Asia is due to these expansions and that the inland populations of the region should be treated as valid subspecies that are distinguishable from other subspecies by genetic traits.     

Species‐level phylogeographical history of the endemic species Calligonum roborovskii and its close relatives in Calligonum section Medusa (Polygonaceae) in arid north‐western China

Quaternary climatic oscillations appear to have influenced the genetic diversity and evolutionary history of arid‐adapted plants. To understand the processes involved and reveal evolutionary relationships, haplotypes were examined from Calligonum roborovskii, an endemic species occurring in the arid zones across the desert regions of north‐western China, and seven other species also from Calligonum section Medusa, including C. gobicum, C. mongolicum and the narrow endemic species C. ebi‐nuricum, C. pumilum, C. taklimakanense, C. trifarium and C. yengisaricum. Forty‐three haplotypes were identified in 422 individuals from 51 natural populations, from variation of two plastid DNA intergenic spacers (rpl32‐trnL and ycf6‐psbM). A high level of total genetic diversity was found across species for which more than two populations were examined, including C. gobicum, C. mongolicum, C. pumilum and C. roborovskii. A distinct isolation‐by‐distance pattern in each of these species was suggested by the Mantel test, indicating that restricted gene flow caused high genetic differentiation among populations. Three haplotypes were shared by two or three species each, but the other 40 haplotypes were species‐specific. The 43 haplotypes split into three major clades, but not species‐specific lineages; most of the Calligonum species were not reciprocally monophyletic, probably due to incomplete lineage sorting or introgression. The identified haplotypes were dated to 1.97 Mya (95% highest posterior density: 2.95–0.99 Mya) and diverged until the late Pleistocene, possibly linked to aridification and enlargement of deserts caused by climate changes. Variation of desert habitats during the Pleistocene might play a key role in causing the divergence.     

Single‐cell profiling screen identifies microtubule‐dependent reduction of variability in signaling

Populations of isogenic cells often respond coherently to signals, despite differences in protein abundance and cell state. Previously, we uncovered processes in the Saccharomyces cerevisiae pheromone response system (PRS) that reduced cell‐to‐cell variability in signal strength and cellular response. Here, we screened 1,141 non‐essential genes to identify 50 “variability genes”. Most had distinct, separable effects on strength and variability of the PRS, defining these quantities as genetically distinct “axes” of system behavior. Three genes affected cytoplasmic microtubule function: BIM1, GIM2, and GIM4. We used genetic and chemical perturbations to show that, without microtubules, PRS output is reduced but variability is unaffected, while, when microtubules are present but their function is perturbed, output is sometimes lowered, but its variability is always high. The increased variability caused by microtubule perturbations required the PRS MAP kinase Fus3 and a process at or upstream of Ste5, the membrane‐localized scaffold to which Fus3 must bind to be activated. Visualization of Ste5 localization dynamics demonstrated that perturbing microtubules destabilized Ste5 at the membrane signaling site. The fact that such microtubule perturbations cause aberrant fate and polarity decisions in mammals suggests that microtubule‐dependent signal stabilization might also operate throughout metazoans.     

Association of dopamine‐ and serotonin‐related genes with canine aggression

Human‐directed canine aggression was studied using 50 aggressive and 81 non‐aggressive dogs. We examined 62 single nucleotide polymorphisms (SNPs) occurring in or in the close vicinity of 16 neurotransmitter‐related genes. Allelic associations with aggression were identified for DRD1, HTR1D, HTR2C and SLC6A1. Risk or protective haplotypes for aggressive behaviour based on 2–5 SNPs were identified. The frequency of aggressive dogs varied significantly between the haplotypes within loci and the odds ratios of aggression in dogs with risk haplotypes compared with protective haplotypes varied from 4.4 (HTR2C) to 9.0 (SLC6A1). A risk haplotype across the neurotransmitter receptor gene HTR1D harboured a non‐synonymous SNP with a potential effect on protein function. We identified no haplotypes in complete association with the recorded phenotypes, supporting a complex inheritance of aggression.     

Isolation and characterization of 12 tetranucleotide repeat microsatellite loci from the green‐backed tit (Parus monticolus)

From a partial genomic library enriched for GATA short tandem repeats, we developed 12 polymorphic microsatellite loci from the green‐backed tit (Parus monticolus). We characterized these loci by genotyping 30 adult individuals with unknown relationship. The number of alleles ranged from four to 17 per locus (mean = 9.3 alleles) and the observed heterozygosity for each locus ranged from 0.633 to 0.933 (mean = 0.789). All loci conformed to Hardy–Weinberg expectations. Four of 66 possible pairwise comparisons between loci showed significant gametic disequilibrium.     

Phylogeography of the brown hare (Lepus europaeus) in Europe: a legacy of south‐eastern Mediterranean refugia?

Aim We analysed the population genetics of the brown hare (Lepus europaeus) in order to test the hypothesis that this species migrated into central Europe from a number of late glacial refugia, including some in Asia Minor. Location Thirty‐three localities in Greece, Bulgaria, Italy, Croatia, Serbia, Poland, Switzerland, Austria, France, Germany, the Netherlands, Spain, the United Kingdom, Turkey and Israel. Methods In total, 926 brown hares were analysed for mitochondrial DNA (mtDNA) variation by restriction fragment length polymorphism (RFLP) performed on polymerase chain reaction‐amplified products spanning cytochrome b (cyt b)/control region (CR), cytochrome oxidase I (COI) and 12S–16S rRNA. In addition, sequence analysis of the mtDNA CR‐I region was performed on 69 individuals, and the data were compared with 137 mtDNA CR‐I sequences retrieved from GenBank. Results The 112 haplotypes detected were partitioned into five phylogeographically well‐defined major haplogroups, namely the ‘south‐eastern European type haplogroup’ (SEEh), ‘Anatolian/Middle Eastern type haplogroup’ (AMh), ‘European type haplogroup, subgroup A’ (EUh‐A), ‘European type haplogroup, subgroup B’ (EUh‐B) and ‘Intermediate haplogroup’ (INTERh). Sequence data retrieved from GenBank were consistent with the haplogroups determined in this study. In Bulgaria and north‐eastern Greece numerous haplotypes of all five haplogroups were present, forming a large overlap zone. Main conclusions The mtDNA results allow us to infer post‐glacial colonization of large parts of Europe from a late glacial/early Holocene source population in the central or south‐central Balkans. The presence of Anatolian/Middle Eastern haplotypes in the large overlap zone in Bulgaria and north‐eastern Greece reveals gene flow from Anatolia to Europe across the late Pleistocene Bosporus land‐bridge. Although various restocking operations could be partly responsible for the presence of unexpected haplotypes in certain areas, we nevertheless trace a strong phylogeographic signal throughout all regions under study. Throughout Europe, mtDNA results indicate that brown hares are not separated into discernable phyletic groups.     

Haplotype‐based genotyping‐by‐sequencing in oat genome research

In a de novo genotyping‐by‐sequencing (GBS) analysis of short, 64‐base tag‐level haplotypes in 4657 accessions of cultivated oat, we discovered 164741 tag‐level (TL) genetic variants containing 241224 SNPs. From this, the marker density of an oat consensus map was increased by the addition of more than 70000 loci. The mapped TL genotypes of a 635‐line diversity panel were used to infer chromosome‐level (CL) haplotype maps. These maps revealed differences in the number and size of haplotype blocks, as well as differences in haplotype diversity between chromosomes and subsets of the diversity panel. We then explored potential benefits of SNP vs. TL vs. CL GBS variants for mapping, high‐resolution genome analysis and genomic selection in oats. A combined genome‐wide association study (GWAS) of heading date from multiple locations using both TL haplotypes and individual SNP markers identified 184 significant associations. A comparative GWAS using TL haplotypes, CL haplotype blocks and their combinations demonstrated the superiority of using TL haplotype markers. Using a principal component‐based genome‐wide scan, genomic regions containing signatures of selection were identified. These regions may contain genes that are responsible for the local adaptation of oats to Northern American conditions. Genomic selection for heading date using TL haplotypes or SNP markers gave comparable and promising prediction accuracies of up to r = 0.74. Genomic selection carried out in an independent calibration and test population for heading date gave promising prediction accuracies that ranged between r = 0.42 and 0.67. In conclusion, TL haplotype GBS‐derived markers facilitate genome analysis and genomic selection in oat.     

The prion‐related protein (testis‐specific) gene (PRNT) is highly polymorphic in Portuguese sheep

The objective of this study was to search for polymorphisms in the ovine prion‐related protein (testis‐specific) gene (PRNT). Sampling included 567 sheep from eight Portuguese breeds. The PRNT gene‐coding region was analyzed by single‐strand conformation polymorphism and sequencing, allowing the identification of the first ovine PRNT polymorphisms, in codons 6, 38, 43 and 48: c.17C>T (p.Ser6Phe, which disrupts a consensus arginine‐X‐X‐serine/threonine motif); c.112G>C (p.Gly38>Arg); c.129T>C and c.144A>G (synonymous) respectively. Polymorphisms in codons 6, 38 and 48 occur simultaneously in 50.6% of the animals, 38.8% presenting as heterozygous. To study the distribution of the polymorphism in codon 43, a restriction fragment length polymorphism analysis was performed. Polymorphic variant c.129C, identified in 89.8% of the animals with 32.8% presented as heterozygous, was considered the wild genotype in Portuguese sheep. Eight different haplotypes which have comparable distribution in all breeds were identified for the PRNT gene. In conclusion, the PRNT coding region is highly polymorphic in sheep, unlike the prion protein 2 dublet gene (PRND), in which we previously found only one synonymous substitution (c.78G>A), in codon 26. The absence or reduced number of PRND heterozygotes (c.78G>A) was significantly associated with three PRNT haplotypes (17C‐112G‐129T‐144A,17CT‐112GC‐129CT‐144AG and 17T‐112C‐129C‐144G), and the only three animals found homozygous at c.78A had the 17C‐112G‐129C‐144A PRNT haplotype. These results constitute evidence of an association between polymorphic variation in PRND and PRNT genes, as has already been observed for PRND and prion protein gene (PRNP).     

A Survey of Sample Size Formulas for Pairwise and Many‐one Multiple Comparisons in the Parametric,Nonparametric and Binomial Case

We present a survey of sample size formulas derived in other papers for pairwise comparisons of k treatments and for comparisons of k treatments with a control. We consider the calculation of sample sizes with preassigned per‐pair, any‐pair and all‐pairs power for tests that control either the comparisonwise or the experimentwise type I error rate. A comparison exhibits interesting similarities between the parametric, nonparametric and binomial case.     

Using conditional power of network meta‐analysis (NMA) to inform the design of future clinical trials

Clinical trials are typically designed with an aim to reach sufficient power to test a hypothesis about relative effectiveness of two or more interventions. Their role in informing evidence‐based decision‐making demands, however, that they are considered in the context of the existing evidence. Consequently, their planning can be informed by characteristics of relevant systematic reviews and meta‐analyses. In the presence of multiple competing interventions the evidence base has the form of a network of trials, which provides information not only about the required sample size but also about the interventions that should be compared in a future trial. In this paper we present a methodology to evaluate the impact of new studies, their information size, the comparisons involved, and the anticipated heterogeneity on the conditional power (CP) of the updated network meta‐analysis. The methods presented are an extension of the idea of CP initially suggested for a pairwise meta‐analysis and we show how to estimate the required sample size using various combinations of direct and indirect evidence in future trials. We apply the methods to two previously published networks and we show that CP for a treatment comparison is dependent on the magnitude of heterogeneity and the ratio of direct to indirect information in existing and future trials for that comparison. Our methodology can help investigators calculate the required sample size under different assumptions about heterogeneity and make decisions about the number and design of future studies (set of treatments compared).