Cross-Platform Microarray Data Normalisation for Regulatory Network Inference

Background Inferring Gene Regulatory Networks (GRNs) from time course microarray data suffers from the dimensionality problem created by the short length of available time series compared to the large number of genes in the network. To overcome this, data integration from diverse sources is mandatory. Microarray data from different sources and platforms are publicly available, but integration is not straightforward, due to platform and experimental differences. Methods We analyse here different normalisation approaches for microarray data integration, in the context of reverse engineering of GRN quantitative models. We introduce two preprocessing approaches based on existing normalisation techniques and provide a comprehensive comparison of normalised datasets. Conclusions Results identify a method based on a combination of Loess normalisation and iterative K-means as best for time series normalisation for this problem

Comparison of evolutionary algorithms in gene regulatory network model inference

Background: The evolution of high throughput technologies that measure gene expression levels has created a data base for inferring GRNs (a process also known as reverse engineering of GRNs). However, the nature of these data has made this process very di±cult. At the moment, several methods of discovering qualitative causal relationships between genes with high accuracy from microarray data exist, but large scale quantitative analysis on real biological datasets cannot be performed, to date, as existing approaches are not suitable for real microarray data which are noisy and insu±cient. Results: This paper performs an analysis of several existing evolutionary algorithms for quantitative gene regulatory network modelling. The aim is to present the techniques used and o®er a comprehensive comparison of approaches, under a common framework. Algorithms are applied to both synthetic and real gene expression data from DNA microarrays, and ability to reproduce biological behaviour, scalability and robustness to noise are assessed and compared. Conclusions: Presented is a comparison framework for assessment of evolutionary algorithms, used to infer gene regulatory networks. Promising methods are identi¯ed and a platform for development of appropriate model formalisms is established

Genetic Architecture of a Reinforced, Postmating, Reproductive Isolation Barrier between Neurospora Species Indicates Evolution via Natural Selection

A role for natural selection in reinforcing premating barriers is recognized, but selection for reinforcement of postmating barriers remains controversial. Organisms lacking evolvable premating barriers can theoretically reinforce postmating isolation, but only under restrictive conditions: parental investment in hybrid progeny must inhibit subsequent reproduction, and selected postmating barriers must restore parents' capacity to reproduce successfully. We show that reinforced postmating isolation markedly increases maternal fitness in the fungus Neurospora crassa, and we detect the evolutionary genetic signature of natural selection by quantitative trait locus (QTL) analysis of the reinforced barrier. Hybrid progeny of N. crassa and N. intermedia are highly inviable. Fertilization by local N. intermedia results in early abortion of hybrid fruitbodies, and we show that abortion is adaptive because only aborted maternal colonies remain fully receptive to future reproduction. In the first QTL analysis of postmating reinforcement in microbial eukaryotes, we identify 11 loci for abortive hybrid fruitbody development, including three major QTLs that together explain 30% of trait variance. One of the major QTLs and six QTLs of lesser effect are found on the mating-type determining chromosome of Neurospora. Several reinforcement QTLs are flanked by genetic markers showing either segregation distortion or non-random associations with alleles at other loci in a cross between N. crassa of different clades, suggesting that the loci also are associated with local effects on same-species reproduction. Statistical analysis of the allelic effects distribution for abortive hybrid fruitbody development indicates its evolution occurred under positive selection. Our results strongly support a role for natural selection in the evolution of reinforced postmating isolation in N. crassa

Persistence of single species of symbionts across multiple closelyrelated host species

Some symbiont species are highly host-specific, inhabiting only one or a very few host species, and typically have limited dispersal abilities. When they do occur on multiple host species, populations of such symbionts are expected to become genetically structured across these different host species, and this may eventually lead to new symbiont species over evolutionary timescales. However, a low number of dispersal events of symbionts between host species across time might be enough to prevent population structure and species divergence. Overall, processes of evolutionary divergence and the species status of most putative multi-host symbiont systems are yet to be investigated. Here, we used DNA metabarcoding data of 6,023 feather mites (a total of 2,225 OTU representative sequences) from 147 infracommunities (i.e., the assemblage consisting of all mites of different species collected from the same bird host individual) to investigate patterns of population genetic structure and species status of three different putative multi-host feather mite species Proctophyllodes macedo Vitzthum, 1922, Proctophyllodes motacillae Gaud, 1953, and Trouessartia jedliczkai (Zimmerman, 1894), each of which inhabits a variable number of different closely related wagtail host species (genus Motacilla). We show that mite populations from different host species represent a single species. This pattern was found in all the mite species, suggesting that each of these species is a multi-host species in which dispersal of mites among host species prevents species divergence. Also, we found evidence of limited evolutionary divergence manifested by a low but significant level of population genetic structure among symbiont populations inhabiting different host species. Our study agrees with previous studies showing a higher than expected colonization opportunities in host-specific symbionts. Indeed, our results support that these dispersal events would allow the persistence of multi-host species even in symbionts with limited dispersal capabilities, though additional factors such as the geographical structure of some bird populations may also play a role.This work was supported by the MINECO CGL2011-24466 to RJ and CGL2015-69650-P to RJ and DS

In Silico Prediction and Analysis of Caenorhabditis EF-hand Containing Proteins

Calcium (Ca+2) is a ubiquitous messenger in eukaryotes including Caenorhabditis. Ca+2-mediated signalling processes are usually carried out through well characterized proteins like calmodulin (CaM) and other Ca+2 binding proteins (CaBP). These proteins interact with different targets and activate it by bringing conformational changes. Majority of the EF-hand proteins in Caenorhabditis contain Ca+2 binding motifs. Here, we have performed homology modelling of CaM-like proteins using the crystal structure of Drosophila melanogaster CaM as a template. Molecular docking was applied to explore the binding mechanism of CaM-like proteins and IQ1 motif which is a ∼25 residues and conform to the consensus sequence (I, L, V)QXXXRXXXX(R,K) to serve as a binding site for different EF hand proteins. We made an attempt to identify all the EF-hand (a helix-loop-helix structure characterized by a 12 residues loop sequence involved in metal coordination) containing proteins and their Ca+2 binding affinity in Caenorhabditis by analysing the complete genome sequence. Docking studies revealed that F165, F169, L29, E33, F44, L57, M61, M96, M97, M108, G65, V115, F93, N104, E144 of CaM-like protein is involved in the interaction with IQ1 motif. A maximum of 170 EF-hand proteins and 39 non-EF-hand proteins with Ca+2/metal binding motif were identified. Diverse proteins including enzyme, transcription, translation and large number of unknown proteins have one or more putative EF-hands. Phylogenetic analysis revealed seven major classes/groups that contain some families of proteins. Various domains that we identified in the EF-hand proteins (uncharacterized) would help in elucidating their functions. It is the first report of its kind where calcium binding loop sequences of EF-hand proteins were analyzed to decipher their calcium affinities. Variation in Ca+2-binding affinity of EF-hand CaBP could be further used to study the behaviour of these proteins. Our analyses postulated that Ca+2 is likely to be key player in Caenorhabditis cell signalling

Connectivity and resilience of coral reef metapopulations in marine protected areas : matching empirical efforts to predictive needs

© 2009 The Authors. This is an open-access article distributed under the terms of the Creative Commons Attribution Noncommercial License. The definitive version was published in Coral Reefs 28 (2009): 327-337, doi:10.1007/s00338-009-0466-z.Design and decision-making for marine protected areas (MPAs) on coral reefs require prediction of MPA effects with population models. Modeling of MPAs has shown how the persistence of metapopulations in systems of MPAs depends on the size and spacing of MPAs, and levels of fishing outside the MPAs. However, the pattern of demographic connectivity produced by larval dispersal is a key uncertainty in those modeling studies. The information required to assess population persistence is a dispersal matrix containing the fraction of larvae traveling to each location from each location, not just the current number of larvae exchanged among locations. Recent metapopulation modeling research with hypothetical dispersal matrices has shown how the spatial scale of dispersal, degree of advection versus diffusion, total larval output, and temporal and spatial variability in dispersal influence population persistence. Recent empirical studies using population genetics, parentage analysis, and geochemical and artificial marks in calcified structures have improved the understanding of dispersal. However, many such studies report current self-recruitment (locally produced settlement/settlement from elsewhere), which is not as directly useful as local retention (locally produced settlement/total locally released), which is a component of the dispersal matrix. Modeling of biophysical circulation with larval particle tracking can provide the required elements of dispersal matrices and assess their sensitivity to flows and larval behavior, but it requires more assumptions than direct empirical methods. To make rapid progress in understanding the scales and patterns of connectivity, greater communication between empiricists and population modelers will be needed. Empiricists need to focus more on identifying the characteristics of the dispersal matrix, while population modelers need to track and assimilate evolving empirical results.Work by CB Paris was supported by the National Science Foundation grant NSF-OCE 0550732. Work by M-A Coffroth and SR Thorrold was supported by the National Science Foundation grant NSF-OCE 0424688. Work by TL Shearer was supported by an International Cooperative Biodiversity Group grant R21 TW006662-01 from the Fogarty International Center at the National Institutes of Health

Fusion–fission experiments in Aphidius: evolutionary split without isolation in response to environmental bimodality

Studying host-based divergence naturally maintained by a balance between selection and gene flow can provide valuable insights into genetic underpinnings of host adaptation and ecological speciation in parasites. Selection-gene flow balance is often postulated in sympatric host races, but direct experimental evidence is scarce. In this study, we present such evidence obtained in host races of Aphidius ervi, an important hymenopteran agent of biological control of aphids in agriculture, using a novel fusion–fission method of gene flow perturbation. In our study, between-race genetic divergence was obliterated by means of advanced hybridisation, followed by a multi-generation exposure of the resulting genetically uniform hybrid swarm to a two-host environment. This fusion–fission procedure was implemented under two contrasting regimes of between-host gene flow in two replicated experiments involving different racial pairs. Host-based genetic fission in response to environmental bimodality occurred in both experiments in as little as six generations of divergent adaptation despite continuous gene flow. We demonstrate that fission recovery of host-based divergence evolved faster and hybridisation-induced linkage disequilibrium decayed slower under restricted (6.7%) compared with unrestricted gene flow, directly pointing at a balance between gene flow and divergent selection. We also show, in four separate tests, that random drift had no or little role in the observed genetic split. Rates and patterns of fission divergence differed between racial pairs. Comparative linkage analysis of these differences is currently under way to test for the role of genomic architecture of adaptation in ecology-driven divergent evolution