Animals modify acoustic communication signals in response to noise pollution, but consequences of these modifications are unknown. Vocalizations that transmit best in noise may not be those that best signal male quality, leading to potential conflict between selection pressures. For example, slow paced, narrow bandwidth songs transmit better in noise but are less effective in mate choice and competition than fast paced, wide bandwidth songs. We test the hypothesis that noise affects response to song pace and bandwidth in the context of competition using white-crowned sparrows (Zonotrichia leucophrys). We measure male response to song variation along a gradient of ambient noise levels in San Francisco, CA. We find that males discriminate between wide and narrow bandwidth songs but not between slow and fast paced songs. These findings are biologically relevant because songs in noisy areas tend to have narrow bandwidths. Therefore, this song phenotype potentially increases transmission distance in noise, but elicits weaker responses from competitors. Further, we find that males respond more strongly to stimuli in noisier conditions, supporting the 'urban anger' hypothesis. We suggest that noise affects male responsiveness to song, possibly leading to more territorial conflict in urban areas.

It is known that the allele ancestral to the variation at a polymorphic site cannot be assigned with certainty, and that the most frequently used method to assign the ancestral state - maximum parsimony - is prone to misinference. Estimates of counts of sites that have a certain number of copies of the derived allele in a sample (the unfolded site frequency spectrum, uSFS) made by parsimony are therefore also biased. We previously developed a maximum likelihood method to estimate the uSFS for a focal species using information from two outgroups while assuming simple models of nucleotide substitution. Here, we extend this approach to allow multiple outgroups (implemented for three outgroups), potentially any phylogenetic tree topology and more complex models of nucleotide substitution. We find, however, that two outgroups and the Kimura 2-parameter model are adequate for uSFS inference in most cases. We show that using parsimony to infer the ancestral state at a specific site seriously breaks down in two situations. The first is where the outgroups provide no information about the ancestral state of variation in the focal species. In this case, nucleotide variation will be under-estimated if such sites are excluded. The second is where the minor allele in the focal species agrees with the allelic state of the outgroups. In this situation, parsimony tends to over-estimate the probability of the major allele being derived, because it fails to account for the fact that sites with a high frequency of the derived allele tend to be rare. We present a method that corrects this deficiency, which is capable of providing nearly unbiased estimates of ancestral state probabilities on a site-by-site basis and the uSFS.

The process of speciation involves the accumulation of reproductive isolation (RI) between diverging lineages. Selection can favor increased RI via the process of reinforcement, whereby costs to hybridization impose selection for increased prezygotic RI. Reinforcement results in phenotypic divergence within at least one taxon, as a result of costly hybridization between sympatric taxa. The strength of selection driving reinforcement is determined by the cost of hybridization and the frequency of hybridization. We investigated the cost of hybridization by quantifying post-mating RI barriers among Phlox species that comprise one of the best-studied cases of reinforcement. We determined if the strength of RI differs among lineages that have and have not undergone reinforcement, how much variability there is within species in RI, and whether RI is associated with phylogenetic relatedness. We found high RI for the species that underwent phenotypic divergence due to reinforcement; however, RI was also high between other species pairs. We found extensive variability in RI among individuals within species, and no evidence that the strength of RI was associated with phylogenetic relatedness. We suggest that phenotypic divergence due to reinforcement is associated with the frequency of hybridization and introgression, and not the cost of hybridization in this clade. This article is protected by copyright. All rights reserved.

Clinal variation in body size and related life history traits is common and has stimulated the postulation of several eco-geographical rules. Whereas some clinal patterns are clearly adaptive, the causes of others remain obscure. We investigated intra-specific body size, development time and female fecundity (egg size and number) clines across 13 European populations of the dung fly Sepsis fulgens spanning 20° latitude from southern Italy to Estonia in a genetic common garden approach. Despite very short generation times (ca. 2 weeks at 24 °C), we found a converse Bergmann cline (smaller size at higher latitudes). As development time did not change with latitude (flat cline), integral growth rate thus likely declines towards the pole. At the same time, early fecundity, but not egg size, increased with latitude. Rather than being mediated by seasonal time constraints, the body size reduction in the northernmost flies from Estonia could suggest that these are marginal, edge populations, as when omitting them the body size cline became flat as well. Most of the other sepsid species investigated to date also show flat body size clines, a pattern that strikingly differs from Drosophila. We conclude that S. fulgens life history traits appear to be shaped by similar environmental pressures and selective mechanisms across Europe, be they adaptive or not. This reiterates the suggestion that body size clines can result as a secondary consequence of selection pressures shaping an entire life history syndrome, rendering them inconsistent and unpredictable in general.

Parasitic chytrid fungi have emerged as a significant threat to amphibian species worldwide, necessitating the development of techniques to isolate these pathogens into culture for research purposes. However, early methods of isolating chytrids from their hosts relied on killing amphibians. We modified a pre-existing protocol for isolating chytrids from infected animals to use toe clips and biopsies from toe webbing rather than euthanizing hosts, and distributed the protocol to researchers as part of the BiodivERsA project RACE; here called the RML protocol. In tandem, we developed a lethal procedure for isolating chytrids from tadpole mouthparts. Reviewing a database of use a decade after their inception, we find that these methods have been applied across 5 continents, 23 countries and in 62 amphibian species. Isolation of chytrids by the non-lethal RML protocol occured in 18% of attempts with 207 fungal isolates and three species of chytrid being recovered. Isolation of chytrids from tadpoles occured in 43% of attempts with 334 fungal isolates of one species (Batrachochytrium dendrobatidis) being recovered. Together, these methods have resulted in a significant reduction and refinement of our use of threatened amphibian species and have improved our ability to work with this group of emerging pathogens.

Cnidaria is an animal phylum, whose members probably have the most ancestral musculature. We prepared and characterized, for the first time to our knowledge, native actomyosin from the striated myoepithelium of the adult moon jelly Aurelia sp. The actomyosin contained myosin, paramyosin-like protein, Ser/Thr-kinase, actin, and two isoforms of tropomyosin, but not troponin, which is known to activate contraction dependent on intracellular Ca2+ signaling in almost all striated muscles of bilaterians. Notably, the myosin comprised striated muscle-type heavy chain and smooth muscle-type regulatory light chains. In the presence of Ca2+, the Mg-ATPase activity of actomyosin was stimulated and Ser21 of the regulatory light chain was concomitantly phosphorylated by the addition of calmodulin and myosin light chain kinase prepared from chicken smooth muscle. Collectively, these results suggest that, similar to smooth muscle, the contraction of jellyfish striated muscle is regulated by Ca2+-dependent phosphorylation of the myosin light chain.

Population structure, spatial diffusion, and climatic conditions mediate the spatiotemporal spread of seasonal influenza in temperate regions. However, much of our knowledge of these dynamics stems from a few well-studied countries, such as the United States (US), and the extent to which this applies in different demographic and climatic environments is not fully understood. Using novel data from Norway, Sweden, and Denmark, we applied wavelet analysis and non-parametric spatial statistics to explore the spatiotemporal dynamics of influenza transmission at regional and international scales. We found the timing and amplitude of epidemics were highly synchronized both within and between countries, despite the geographical isolation of many areas in our study. Within Norway, this synchrony was most strongly modulated by population size, confirming previous findings that hierarchical spread between larger populations underlies seasonal influenza dynamics at regional levels. However, we found no such association when comparing across countries, suggesting that other factors become important at the international scale. Finally, to frame our results within a wider global context, we compared our findings from Norway to those from the US. After correcting for differences in geographic scale, we unexpectedly found higher levels of synchrony in Norway, despite its smaller population size. We hypothesize that this greater synchrony may be driven by more favorable and spatially uniform climatic conditions, although there are other likely factors we were unable to consider (such as reduced variation in school term times and differences in population movements). Overall, our results highlight the importance of comparing influenza spread at different spatial scales and across diverse geographic regions in order to better understand the complex mechanisms underlying disease dynamics.

Ustilaginomycotina is home to a broad array of fungi including important plant pathogens collectively called smut fungi. Smuts are biotrophs that produce characteristic perennating propagules called teliospores, one of which, Ustilago maydis, is a model genetic organism. Broad exploration of smut biology has been hampered by limited phylogenetic resolution of Ustilaginiomycotina as well as an overall lack of genomic data for members of this subphylum. In this study, we sequenced eight Ustilaginomycotina genomes from previously unrepresented lineages, deciphered ordinal-level phylogenetic relationships for the subphylum, and performed comparative analyses. Unlike other Basidiomycota subphyla, all sampled Ustilaginomycotina genomes are relatively small and compact. Ancestral state reconstruction analyses indicate that teliospore formation was present at the origin of the subphylum. Divergence time estimation dates the divergence of most extant smut fungi after that of grasses (Poaceae). However, we found limited conservation of well-characterized genes related to smut pathogenesis from U. maydis, indicating dissimilar pathogenic mechanisms exist across other smut lineages. The genomes of Malasseziomycetes are highly diverged from the other sampled Ustilaginomycotina, likely due to their unique history as mammal-associated lipophilic yeasts. Despite the extensive genomic data, the phylogenetic placement of this class remains ambiguous. Although the sampled Ustilaginomycotina members lack many core enzymes for plant cell wall decomposition and starch catabolism, we identified several novel carbohydrate active enzymes potentially related to pectin breakdown. Finally, about 50% of Ustilaginomycotina species-specific genes are present in previously undersampled and rare lineages, highlighting the importance of exploring fungal diversity as a resource for novel gene discovery.

Stomatal conductance (gs) impacts both photosynthesis and transpiration, and is therefore fundamental to the global carbon and water cycles, food production, and ecosystem services. Mathematical models provide the primary means of analysing this important leaf gas exchange parameter. A nearly universal assumption in such models is that the vapour pressure inside leaves (ei) remains saturated under all conditions. The validity of this assumption has not been well tested, because so far ei cannot be measured directly. Here, we test this assumption using a novel technique, based on coupled measurements of leaf gas exchange and the stable isotope compositions of CO2 and water vapour passing over the leaf. We applied this technique to mature individuals of two semiarid conifer species. In both species, ei routinely dropped below saturation when leaves were exposed to moderate to high air vapour pressure deficits. Typical values of relative humidity in the intercellular air spaces were as low 0.9 in Juniperus monosperma and 0.8 in Pinus edulis. These departures of ei from saturation caused significant biases in calculations of gs and the intercellular CO2 concentration. Our results refute the longstanding assumption of saturated vapour pressure in plant leaves under all conditions.

Feeding preference of fossil herbivorous mammals, concerning the coevolution of mammalian and floral ecosystems, has become of key research interest. In this paper, phytoliths in dental calculus from two gomphotheriid proboscideans of the middle Miocene Junggar Basin, Central Asia, have been identified, suggesting that Gomphotherium connexum was a mixed feeder, while the phytoliths from G. steinheimense indicates grazing preference. This is the earliest-known proboscidean with a predominantly grazing habit. These results are further confirmed by microwear and isotope analyses. Pollen record reveals an open steppic environment with few trees, indicating an early aridity phase in the Asian interior during the Mid-Miocene Climate Optimum, which might urge a diet remodeling of G. steinheimense. Morphological and cladistic analyses show that G. steinheimense comprises the sister taxon of tetralophodont gomphotheres, which were believed to be the general ancestral stock of derived "true elephantids"; whereas G. connexum represents a more conservative lineage in both feeding behavior and tooth morphology, which subsequently became completely extinct. Therefore, grazing by G. steinheimense may have acted as a behavior preadaptive for aridity, and allowing its lineage evolving new morphological features for surviving later in time. This study displays an interesting example of behavioral adaptation prior to morphological modification.

Organic matter and nutrient loading into aquatic ecosystems affects ecosystem structure and function and can result in eutrophication and hypoxia. Hypoxia is often attributed to anthropogenic pollution and is not common in unpolluted rivers. Here we show that organic matter loading from hippopotami causes the repeated occurrence of hypoxia in the Mara River, East Africa. We documented 49 high flow events over 3 years that caused dissolved oxygen decreases, including 13 events resulting in hypoxia, and 9 fish kills over 5 years. Evidence from experiments and modeling demonstrates a strong mechanistic link between the flushing of hippo pools and decreased dissolved oxygen in the river. This phenomenon may have been more widespread throughout Africa before hippopotamus populations were severely reduced. Frequent hypoxia may be a natural part of tropical river ecosystem function, particularly in rivers impacted by large wildlife.

BACKGROUND: In phylogenetic reconstruction the result is a tree where all taxa are leaves and internal nodes are hypothetical ancestors. In a live phylogeny, both ancestral and living taxa may coexist, leading to a tree where internal nodes may be living taxa. The well-known Neighbor-Joining heuristic is largely used for phylogenetic reconstruction. RESULTS: We present Live Neighbor-Joining, a heuristic for building a live phylogeny. We have investigated Live Neighbor-Joining on datasets of viral genomes, a plausible scenario for its application, which allowed the construction of alternative hypothesis for the relationships among virus that embrace both ancestral and descending taxa. We also applied Live Neighbor-Joining on a set of bacterial genomes and to sets of images and texts. Non-biological data may be better explored visually when their relationship in terms of content similarity is represented by means of a phylogeny. CONCLUSION: Our experiments have shown interesting alternative phylogenetic hypothesis for RNA virus genomes, bacterial genomes and alternative relationships among images and texts, illustrating a wide range of scenarios where Live Neighbor-Joining may be used.

The patterning of the spine of a zebrafish is controlled by the notochord, a rod-like structure that supports and instructs the developing embryo.

Modern genomics has shed light on many entomopathogenic fungi and expanded our knowledge widely; however, little is known about the genomic features of the insect-commensal fungi. Harpellales are obligate commensals living in the digestive tracts of disease-bearing insects (black flies, midges, and mosquitoes). In this study, we produced and annotated whole-genome sequences of nine Harpellales taxa and conducted the first comparative analyses to infer the genomic diversity within the members of the Harpellales. The genomes of the insect gut fungi feature low (26% to 37%) GC content and large genome size variations (25 to 102 Mb). Further comparisons with insect-pathogenic fungi (from both Ascomycota and Zoopagomycota), as well as with free-living relatives (as negative controls), helped to identify a gene toolbox that is essential to the fungus-insect symbiosis. The results not only narrow the genomic scope of fungus-insect interactions from several thousands to eight core players but also distinguish host invasion strategies employed by insect pathogens and commensals. The genomic content suggests that insect commensal fungi rely mostly on adhesion protein anchors that target digestive system, while entomopathogenic fungi have higher numbers of transmembrane helices, signal peptides, and pathogen-host interaction (PHI) genes across the whole genome and enrich genes as well as functional domains to inactivate the host inflammation system and suppress the host defense. Phylogenomic analyses have revealed that genome sizes of Harpellales fungi vary among lineages with an integer-multiple pattern, which implies that ancient genome duplications may have occurred within the gut of insects.IMPORTANCE Insect guts harbor various microbes that are important for host digestion, immune response, and disease dispersal in certain cases. Bacteria, which are among the primary endosymbionts, have been studied extensively. However, fungi, which are also frequently encountered, are poorly known with respect to their biology within the insect guts. To understand the genomic features and related biology, we produced the whole-genome sequences of nine gut commensal fungi from disease-bearing insects (black flies, midges, and mosquitoes). The results show that insect gut fungi tend to have low GC content across their genomes. By comparing these commensals with entomopathogenic and free-living fungi that have available genome sequences, we found a universal core gene toolbox that is unique and thus potentially important for the insect-fungus symbiosis. This comparative work also uncovered different host invasion strategies employed by insect pathogens and commensals, as well as a model system to study ancient fungal genome duplication within the gut of insects.

The K50 (lysine at amino acid position 50) homeodomain (HD) protein Orthodenticle (Otd) is critical for anterior patterning and brain and eye development in most metazoans. In Drosophila melanogaster, another K50HD protein, Bicoid (Bcd), has evolved to replace Otd's ancestral function in embryo patterning. Bcd is distributed as a long-range maternal gradient and activates transcription of a large number of target genes, including otd Otd and Bcd bind similar DNA sequences in vitro, but how their transcriptional activities are integrated to pattern anterior regions of the embryo is unknown. Here we define three major classes of enhancers that are differentially sensitive to binding and transcriptional activation by Bcd and Otd. Class 1 enhancers are initially activated by Bcd, and activation is transferred to Otd via a feed-forward relay (FFR) that involves sequential binding of the two proteins to the same DNA motif. Class 2 enhancers are activated by Bcd and maintained by an Otd-independent mechanism. Class 3 enhancers are never bound by Bcd, but Otd binds and activates them in a second wave of zygotic transcription. The specific activities of enhancers in each class are mediated by DNA motif variants preferentially bound by Bcd or Otd and the presence or absence of sites for cofactors that interact with these proteins. Our results define specific patterning roles for Bcd and Otd and provide mechanisms for coordinating the precise timing of gene expression patterns during embryonic development.

We present AAI-profiler, a web server for exploratory analysis and quality control in comparative genomics. AAI-profiler summarizes proteome-wide sequence search results to identify novel species, assess the need for taxonomic reclassification and detect multi-isolate and contaminated samples. AAI-profiler visualises results using a scatterplot that shows the Average Amino-acid Identity (AAI) from the query proteome to all similar species in the sequence database. Taxonomic groups are indicated by colour and marker styles, making outliers easy to spot. AAI-profiler uses SANSparallel to perform high-performance homology searches, making proteome-wide analysis possible. We demonstrate the efficacy of AAI-profiler in the discovery of a close relationship between two bacterial symbionts of an omnivorous pirate bug (Orius) and a thrip (Frankliniella occidentalis), an important pest in agriculture. The symbionts represent novel species within the genus Rosenbergiella so far described only in floral nectar. AAI-profiler is easy to use, the analysis presented only required two mouse clicks and was completed in a few minutes. AAI-profiler is available at http://ekhidna2.biocenter.helsinki.fi/AAI.

Inbreeding depression plays a major role in shaping mating systems: in particular, inbreeding avoidance is often proposed as a mechanism explaining extra-pair reproduction in socially-monogamous species. This suggestion relies on assumptions which are rarely comprehensively tested: that inbreeding depression is present, that higher kinship between social partners increases infidelity, and that infidelity reduces the frequency of inbreeding. Here, we test these assumptions using 26 years of data for a cooperatively-breeding, socially-monogamous bird with high female infidelity, the superb fairy-wren (Malurus cyaneus). Although inbred individuals were rare (∼6% of offspring), we found evidence of inbreeding depression in nestling mass (but not in fledgling survival). Mother-son social pairings resulted in 100% infidelity, but kinship between a social pair did not otherwise predict female infidelity. Nevertheless, extra-pair offspring were less likely to be inbred than within-pair offspring. Finally, the social environment (the number of helpers in a group) did not affect offspring inbreeding coefficients or inbreeding depression levels. In conclusion, despite some agreement with the assumptions that are necessary for inbreeding avoidance to drive infidelity, the apparent scarcity of inbreeding events and the observed levels of inbreeding depression seem insufficient to explain the ubiquitous infidelity in this system, beyond the mother-son mating avoidance. This article is protected by copyright. All rights reserved.

The extent to which diet and environment influence gut community membership (presence or absence of taxa) and structure (individual taxon abundance) is the subject of growing interest in microbiome research. Here, we examined the gut bacterial communities of three cricket groups: (1) wild caught field crickets, (2) laboratory-reared crickets fed cat chow, and (3) laboratory-reared crickets fed chemically defined diets. We found that both environment and diet greatly altered the structure of the gut bacterial community. Wild crickets had greater gut microbial diversity and higher Firmicutes to Bacteroidetes ratios, in contrast to laboratory-reared crickets. Predictive metagenomes revealed that laboratory-reared crickets were significantly enriched in amino acid degradation pathways, while wild crickets had a higher relative abundance of peptidases that would aid in amino acid release. Although wild and laboratory animals differ greatly in their bacterial communities, we show that the community proportional membership remains stable from Phylum to Family taxonomic levels regardless of differences in environment and diet, suggesting that endogenous factors, such as host genetics, have greater control in shaping gut community membership.

There have been over 25 independent unicellular to multicellular evolutionary transitions, which have been transformational in the complexity of life. All of these transitions likely occurred in communities numerically dominated by unicellular organisms, mostly bacteria. Hence, it is reasonable to expect that bacteria were involved in generating the ecological conditions that promoted the stability and proliferation of the first multicellular forms as protective units. In this study, we addressed this problem by analyzing the occurrence of multicellularity in an experimental phylogeny of yeasts (Sacharomyces cerevisiae) a model organism that is unicellular but can generate multicellular clusters under some conditions. We exposed a single ancestral population to periodic divergences, coevolving with a cocktail of environmental bacteria that were inoculated to the environment of the ancestor, and compared to a control (no bacteria). We quantified culturable microorganisms to the level of genera, finding up to 20 taxa (all bacteria) that competed with the yeasts during diversification. After 600 generations of coevolution, the yeasts produced two types of multicellular clusters: clonal and aggregative. Whereas clonal clusters were present in both treatments, aggregative clusters were only present under the bacteria treatment and showed significant phylogenetic signal. However, clonal clusters showed different properties if bacteria were present as follows: They were more abundant and significantly smaller than in the control. These results indicate that bacteria are important modulators of the occurrence of multicellularity, providing support to the idea that they generated the ecological conditions-promoting multicellularity.

Polyandry, where multiple mating by females results in the temporal and spatial overlap of ejaculates from two or more males, is taxonomically widespread and occurs in varying frequencies within and among species. In decapods (crabs, lobsters, crayfish, and prawns), rates of polyandry are likely to be variable, but the extent to which patterns of multiple paternity reflect multiple mating, and thus are shaped by postmating processes that bias fertilization toward one or a subset of mated males, is unclear. Here, we use microsatellite markers to examine the frequency of multiple mating (the presence of spermatophores from two or more males) and patterns of paternity in wild populations of western rock lobster (Panulirus cygnus). Our data confirm that >45% of females had attached spermatophores arising from at least two males (i.e., confirming polyandry), but we found very limited evidence for multiple paternity; among 24 clutches sampled in this study, only two arose from fertilizations by two or more males. Single inferred paternal genotypes accounted for all remaining progeny genotypes in each clutch, including several instances when the mother had been shown to mate with two or more males. These findings highlight the need for further work to understand whether polyandry is adaptive and to uncover the mechanisms underlying postmating paternity biases in this system.

Background: Glycogen synthase kinase 3 (GSK3) regulates many cell fate decisions in animal development. In multicellular structures of the group 4 dictyostelid Dictyostelium discoideum, GSK3 promotes spore over stalk-like differentiation. We investigated whether, similar to other sporulation-inducing genes such as cAMP-dependent protein kinase (PKA), this role of GSK3 is derived from an ancestral role in encystation of unicellular amoebas. Results: We deleted GSK3 in Polysphondylium pallidum, a group 2 dictyostelid which has retained encystation as an alternative survival strategy. Loss of GSK3 inhibited cytokinesis of cells in suspension, as also occurs in D. discoideum, but did not affect spore or stalk differentiation in P. pallidum. However, gsk3- amoebas entered into encystation under conditions that in wild type favour aggregation and fruiting body formation. The gsk3- cells were hypersensitive to osmolytes, which are known to promote encystation, and to cyst-inducing factors that are secreted during starvation. GSK3 was not itself regulated by these factors, but inhibited their effects. Conclusions: Our data show that GSK3 has a deeply conserved role in controlling cytokinesis, but not spore differentiation in Dictyostelia. Instead, in P. pallidum, one of many Dictyostelia that like their solitary ancestors can still encyst to survive starvation, GSK3 promotes multicellular development into fruiting bodies over unicellular encystment.

The genetic code used in nuclear genes is almost universal, but here we report that it changed three times in parallel during the evolution of budding yeasts. All three changes were reassignments of the codon CUG, which is translated as serine (in 2 yeast clades), alanine (1 clade), or the 'universal' leucine (2 clades). The newly discovered Ser2 clade is in the final stages of a genetic code transition. Most species in this clade have genes for both a novel tRNASer(CAG) and an ancestral tRNALeu(CAG) to read CUG, but only tRNASer(CAG) is used in standard growth conditions. The coexistence of these alloacceptor tRNA genes indicates that the genetic code transition occurred via an ambiguous translation phase. We propose that the three parallel reassignments of CUG were not driven by natural selection in favor of their effects on the proteome, but by selection to eliminate the ancestral tRNALeu(CAG).

Modularity is a widespread property in biological systems. It implies that interactions occur mainly within groups of system elements. A modular arrangement facilitates adjustment of one module without perturbing the rest of the system. Therefore, modularity of developmental mechanisms is a major factor for evolvability, the potential to produce beneficial variation from random genetic change. Understanding how modularity evolves in gene regulatory networks, that create the distinct gene activity patterns that characterize different parts of an organism, is key to developmental and evolutionary biology. One hypothesis for the evolution of modules suggests that interactions between some sets of genes become maladaptive when selection favours additional gene activity patterns. The removal of such interactions by selection would result in the formation of modules. A second hypothesis suggests that modularity evolves in response to sparseness, the scarcity of interactions within a system. Here I simulate the evolution of gene regulatory networks and analyse diverse experimentally sustained networks to study the relationship between sparseness and modularity. My results suggest that sparseness alone is neither sufficient nor necessary to explain modularity in gene regulatory networks. However, sparseness amplifies the effects of forms of selection that, like selection for additional gene activity patterns, already produce an increase in modularity. That evolution of new gene activity patterns is frequent across evolution also supports that it is a major factor in the evolution of modularity. That sparseness is widespread across gene regulatory networks indicates that it may have facilitated the evolution of modules in a wide variety of cases.