@article{12193,
  abstract     = {DNA methylation regulates eukaryotic gene expression and is extensively reprogrammed during animal development. However, whether developmental methylation reprogramming during the sporophytic life cycle of flowering plants regulates genes is presently unknown. Here we report a distinctive gene-targeted RNA-directed DNA methylation (RdDM) activity in the Arabidopsis thaliana male sexual lineage that regulates gene expression in meiocytes. Loss of sexual-lineage-specific RdDM causes mis-splicing of the MPS1 gene (also known as PRD2), thereby disrupting meiosis. Our results establish a regulatory paradigm in which de novo methylation creates a cell-lineage-specific epigenetic signature that controls gene expression and contributes to cellular function in flowering plants.},
  author       = {Walker, James and Gao, Hongbo and Zhang, Jingyi and Aldridge, Billy and Vickers, Martin and Higgins, James D. and Feng, Xiaoqi},
  issn         = {1546-1718},
  journal      = {Nature Genetics},
  keywords     = {Genetics},
  number       = {1},
  pages        = {130--137},
  publisher    = {Nature Research},
  title        = {{Sexual-lineage-specific DNA methylation regulates meiosis in Arabidopsis}},
  doi          = {10.1038/s41588-017-0008-5},
  volume       = {50},
  year         = {2017},
}

@article{1228,
  abstract     = {Since 2006, reprogrammed cells have increasingly been used as a biomedical research technique in addition to neuro-psychiatric methods. These rapidly evolving techniques allow for the generation of neuronal sub-populations, and have sparked interest not only in monogenetic neuro-psychiatric diseases, but also in poly-genetic and poly-aetiological disorders such as schizophrenia (SCZ) and bipolar disorder (BPD). This review provides a summary of 19 publications on reprogrammed adult somatic cells derived from patients with SCZ, and five publications using this technique in patients with BPD. As both disorders are complex and heterogeneous, there is a plurality of hypotheses to be tested in vitro. In SCZ, data on alterations of dopaminergic transmission in vitro are sparse, despite the great explanatory power of the so-called DA hypothesis of SCZ. Some findings correspond to perturbations of cell energy metabolism, and observations in reprogrammed cells suggest neuro-developmental alterations. Some studies also report on the efficacy of medicinal compounds to revert alterations observed in cellular models. However, due to the paucity of replication studies, no comprehensive conclusions can be drawn from studies using reprogrammed cells at the present time. In the future, findings from cell culture methods need to be integrated with clinical, epidemiological, pharmacological and imaging data in order to generate a more comprehensive picture of SCZ and BPD.},
  author       = {Sauerzopf, Ulrich and Sacco, Roberto and Novarino, Gaia and Niello, Marco and Weidenauer, Ana and Praschak Rieder, Nicole and Sitte, Harald and Willeit, Matthaeus},
  journal      = {European Journal of Neuroscience},
  number       = {1},
  pages        = {45 -- 57},
  publisher    = {Wiley-Blackwell},
  title        = {{Are reprogrammed cells a useful tool for studying dopamine dysfunction in psychotic disorders? A review of the current evidence}},
  doi          = {10.1111/ejn.13418},
  volume       = {45},
  year         = {2017},
}

@article{123,
  abstract     = {The Leidenfrost effect occurs when an object near a hot surface vaporizes rapidly enough to lift itself up and hover. Although well understood for liquids and stiff sublimable solids, nothing is known about the effect with materials whose stiffness lies between these extremes. Here we introduce a new phenomenon that occurs with vaporizable soft solids - the elastic Leidenfrost effect. By dropping hydrogel spheres onto hot surfaces we find that, rather than hovering, they energetically bounce several times their diameter for minutes at a time. With high-speed video during a single impact, we uncover high-frequency microscopic gap dynamics at the sphere/substrate interface. We show how these otherwise-hidden agitations constitute work cycles that harvest mechanical energy from the vapour and sustain the bouncing. Our findings suggest a new strategy for injecting mechanical energy into a widely used class of soft materials, with potential relevance to fields such as active matter, soft robotics and microfluidics.},
  author       = {Waitukaitis, Scott R and Zuiderwijk, Antal and Souslov, Anton and Coulais, Corentin and Van Hecke, Martin},
  journal      = {Nature Physics},
  number       = {11},
  pages        = {1095 -- 1099},
  publisher    = {Nature Publishing Group},
  title        = {{Coupling the Leidenfrost effect and elastic deformations to power sustained bouncing}},
  doi          = {10.1038/nphys4194},
  volume       = {13},
  year         = {2017},
}

@inproceedings{12571,
  abstract     = {We consider the problems of maintaining approximate maximum matching and minimum vertex cover in a dynamic graph. Starting with the seminal work of Onak and Rubinfeld [STOC 2010], this problem has received significant attention in recent years. Very recently, extending the framework of Baswana, Gupta and Sen [FOCS 2011], Solomon [FOCS 2016] gave a randomized 2-approximation dynamic algorithm for this problem that has amortized update time of O(1) with high probability. We consider the natural open question of derandomizing this result. We present a new deterministic fully dynamic algorithm that maintains a O(1)-approximate minimum vertex cover and maximum fractional matching, with an amortized update time of O(1). Previously, the best deterministic algorithm for this problem was due to Bhattacharya, Henzinger and Italiano [SODA 2015]; it had an approximation ratio of (2+ϵ) and an amortized update time of O(logn/ϵ2). Our result can be generalized to give a fully dynamic O(f3)-approximation algorithm with O(f2) amortized update time for the hypergraph vertex cover and fractional matching problems, where every hyperedge has at most f vertices.},
  author       = {Bhattacharya, Sayan and Chakrabarty, Deeparnab and Henzinger, Monika H},
  booktitle    = {19th International Conference on Integer Programming and Combinatorial Optimization},
  isbn         = {9783319592497},
  issn         = {0302-9743},
  location     = {Waterloo, ON, Canada},
  pages        = {86--98},
  publisher    = {Springer Nature},
  title        = {{Deterministic fully dynamic approximate vertex cover and fractional matching in O(1) amortized update time}},
  doi          = {10.1007/978-3-319-59250-3_8},
  volume       = {10328},
  year         = {2017},
}

@article{12608,
  abstract     = {The spatio-temporal distribution of air temperature over mountain glaciers can demonstrate complex patterns, yet it is often represented simplistically using linear vertical temperature gradients (VTGs) extrapolated from off-glacier locations. We analyse a network of centreline and lateral air temperature observations at Tsanteleina Glacier, Italy, during summer 2015. On average, VTGs are steep (&lt;−0.0065 °C m<jats:sup>−1</jats:sup>), but they are shallow under warm ambient conditions when the correlation between air temperature and elevation becomes weaker. Published along-flowline temperature distribution methods explain centreline observations well, including warming on the lower glacier tongue, but cannot estimate lateral temperature variability. Application of temperature distribution methods improves simulation of melt rates (RMSE) in an energy-balance model by up to 36% compared to the environmental lapse rate extrapolated from an off-glacier station. However, results suggest that model parameters are not easily transferable to glaciers with a small fetch without recalibration. Such methods have potential to improve estimates of temperature across a glacier, but their parameter transferability should be further linked to the glacier and atmospheric characteristics. Furthermore, ‘cold spots’, which can be &gt;2°C cooler than expected for their elevation, whose occurrence is not predicted by the temperature distribution models, are identified at one-quarter of the measurement sites.},
  author       = {SHAW, THOMAS E. and BROCK, BEN W. and AYALA, ÁLVARO and RUTTER, NICK and Pellicciotti, Francesca},
  issn         = {1727-5652},
  journal      = {Journal of Glaciology},
  keywords     = {Earth-Surface Processes},
  number       = {242},
  pages        = {973--988},
  publisher    = {Cambridge University Press},
  title        = {{Centreline and cross-glacier air temperature variability on an Alpine glacier: Assessing temperature distribution methods and their influence on melt model calculations}},
  doi          = {10.1017/jog.2017.65},
  volume       = {63},
  year         = {2017},
}

@article{12609,
  abstract     = {Previous estimates of melt and surface sublimation on glaciers of the subtropical semiarid Andes (29–34°S) have been obtained at few specific locations, but it is not clear how ablation components vary across the entire extent of a glacier in this dry environment. Here, we simulate the distributed energy and mass balance of Juncal Norte Glacier (33°S) during a 2-month summer period. Forcing fields of near-surface air temperature and wind speed are generated using two methods accounting for the main physical processes that shape their spatial variations. Simulated meteorological variables and ablation agree well with observations on the glacier tongue and reveal complex patterns of energy and mass fluxes. Ablation decreases from 70 mm w.e. d<jats:sup>−1</jats:sup> at the low-albedo glacier terminus (~3000 m), where almost 100% of total ablation corresponds to melt, to &lt;5 mm w.e. d<jats:sup>−1</jats:sup> at wind-exposed, strong-radiated sites above 5500 m, where surface sublimation represents &gt;75% of total ablation. Our simulations provide the first glacier-scale estimates of ablation components on a glacier in the study region and better reproduce the observed and expected spatial variations of melt and surface sublimation, in comparison with more simple assumptions, such as linear gradients and uniform wind speeds.},
  author       = {AYALA, A. and Pellicciotti, Francesca and PELEG, N. and BURLANDO, P.},
  issn         = {1727-5652},
  journal      = {Journal of Glaciology},
  keywords     = {Earth-Surface Processes},
  number       = {241},
  pages        = {803--822},
  publisher    = {Cambridge University Press},
  title        = {{Melt and surface sublimation across a glacier in a dry environment: distributed energy-balance modelling of Juncal Norte Glacier, Chile}},
  doi          = {10.1017/jog.2017.46},
  volume       = {63},
  year         = {2017},
}

@article{12610,
  abstract     = {The hydrological systems of heavily-downwasted debris-covered glaciers differ from those of clean-ice glaciers due to the hummocky surface and debris mantle of such glaciers, leading to a relatively limited understanding of drainage pathways. Supraglacial ponds represent sinks within the discontinuous supraglacial drainage system, and occasionally drain englacially. To assess pond dynamics, we made pond water level measurements on Lirung Glacier, Nepal, during May and October of 2013 and 2014. Simultaneously, aerial, satellite, and terrestrial orthoimages and digital elevation models were obtained, providing snapshots of the ponds and their surroundings. We performed a DEM-based analysis of the glacier's closed surface catchments to identify surface drainage pathways and englacial drainage points, and compared this to field observations of surface and near-surface water flow. The total ponded area was higher in the pre-monsoon than post-monsoon, with individual ponds filling and draining seasonally associated with the surface exposure of englacial conduit segments. We recorded four pond drainage events, all of which occurred gradually (duration of weeks), observed diurnal fluctuations indicative of varying water supply and outflow discharge, and we documented instances of interaction between distant ponds. The DEM drainage analysis identified numerous sinks >3 m in depth across the glacier surface, few of which exhibited ponds (23%), while the field survey highlighted instances of surface water only explicable via englacial routes. Taken together, our observations provide evidence for widespread supraglacial-englacial connectivity of meltwater drainage paths. Results suggest that successive englacial conduit collapse events, themselves likely driven by supraglacial pond drainage, cause the glacier surface drainage system to evolve into a configuration following relict englacial conduit systems. Within this system, ponds form in depressions of reduced drainage efficiency and link the supraglacial and englacial drainage networks.},
  author       = {Miles, Evan S. and Steiner, Jakob and Willis, Ian and Buri, Pascal and Immerzeel, Walter W. and Chesnokova, Anna and Pellicciotti, Francesca},
  issn         = {2296-6463},
  journal      = {Frontiers in Earth Science},
  keywords     = {General Earth and Planetary Sciences},
  publisher    = {Frontiers Media},
  title        = {{Pond dynamics and supraglacial-englacial connectivity on debris-covered Lirung Glacier, Nepal}},
  doi          = {10.3389/feart.2017.00069},
  volume       = {5},
  year         = {2017},
}

@article{12612,
  abstract     = {Supraglacial ponds play a key role in absorbing atmospheric energy and directing it to the ice of debris-covered glaciers, but the spatial and temporal distribution of these features is not well documented. We analyse 172 Landsat TM/ETM+ scenes for the period 1999–2013 to identify thawed supraglacial ponds for the debris-covered tongues of five glaciers in the Langtang Valley of Nepal. We apply an advanced atmospheric correction routine (Landcor/6S) and use band ratio and image morphological techniques to identify ponds and validate our results with 2.5 m Cartosat-1 observations. We then characterize the spatial, seasonal and interannual patterns of ponds. We find high variability in pond incidence between glaciers (May–October means of 0.08–1.69% of debris area), with ponds most frequent in zones of low surface gradient and velocity. The ponds show pronounced seasonality, appearing in the pre-monsoon as snow melts, peaking at the monsoon onset at 2% of debris-covered area, then declining in the post-monsoon as ponds drain or freeze. Ponds are highly recurrent and persistent, with 40.5% of pond locations occurring for multiple years. Rather than a trend in pond cover over the study period, we find high interannual variability for each glacier after controlling for seasonality.},
  author       = {MILES, EVAN S. and WILLIS, IAN C. and ARNOLD, NEIL S. and STEINER, JAKOB and Pellicciotti, Francesca},
  issn         = {1727-5652},
  journal      = {Journal of Glaciology},
  keywords     = {Earth-Surface Processes},
  number       = {237},
  pages        = {88--105},
  publisher    = {Cambridge University Press},
  title        = {{Spatial, seasonal and interannual variability of supraglacial ponds in the Langtang Valley of Nepal, 1999–2013}},
  doi          = {10.1017/jog.2016.120},
  volume       = {63},
  year         = {2017},
}

@inproceedings{12905,
  author       = {Schlögl, Alois and Kiss, Janos},
  booktitle    = {AHPC17 – Austrian HPC Meeting 2017},
  location     = {Grundlsee, Austria},
  pages        = {28},
  publisher    = {FSP Scientific Computing},
  title        = {{Scientific Computing at IST Austria}},
  year         = {2017},
}

@article{1294,
  abstract     = {We study controller synthesis problems for finite-state Markov decision processes, where the objective is to optimize the expected mean-payoff performance and stability (also known as variability in the literature). We argue that the basic notion of expressing the stability using the statistical variance of the mean payoff is sometimes insufficient, and propose an alternative definition. We show that a strategy ensuring both the expected mean payoff and the variance below given bounds requires randomization and memory, under both the above definitions. We then show that the problem of finding such a strategy can be expressed as a set of constraints.},
  author       = {Brázdil, Tomáš and Chatterjee, Krishnendu and Forejt, Vojtěch and Kučera, Antonín},
  journal      = {Journal of Computer and System Sciences},
  pages        = {144 -- 170},
  publisher    = {Elsevier},
  title        = {{Trading performance for stability in Markov decision processes}},
  doi          = {10.1016/j.jcss.2016.09.009},
  volume       = {84},
  year         = {2017},
}

@misc{9707,
  abstract     = {Branching morphogenesis of the epithelial ureteric bud forms the renal collecting duct system and is critical for normal nephron number, while low nephron number is implicated in hypertension and renal disease. Ureteric bud growth and branching requires GDNF signaling from the surrounding mesenchyme to cells at the ureteric bud tips, via the Ret receptor tyrosine kinase and coreceptor Gfrα1; Ret signaling up-regulates transcription factors Etv4 and Etv5, which are also critical for branching. Despite extensive knowledge of the genetic control of these events, it is not understood, at the cellular level, how renal branching morphogenesis is achieved or how Ret signaling influences epithelial cell behaviors to promote this process. Analysis of chimeric embryos previously suggested a role for Ret signaling in promoting cell rearrangements in the nephric duct, but this method was unsuited to study individual cell behaviors during ureteric bud branching. Here, we use Mosaic Analysis with Double Markers (MADM), combined with organ culture and time-lapse imaging, to trace the movements and divisions of individual ureteric bud tip cells. We first examine wild-type clones and then Ret or Etv4 mutant/wild-type clones in which the mutant and wild-type sister cells are differentially and heritably marked by green and red fluorescent proteins. We find that, in normal kidneys, most individual tip cells behave as self-renewing progenitors, some of whose progeny remain at the tips while others populate the growing UB trunks. In Ret or Etv4 MADM clones, the wild-type cells generated at a UB tip are much more likely to remain at, or move to, the new tips during branching and elongation, while their Ret−/− or Etv4−/− sister cells tend to lag behind and contribute only to the trunks. By tracking successive mitoses in a cell lineage, we find that Ret signaling has little effect on proliferation, in contrast to its effects on cell movement. Our results show that Ret/Etv4 signaling promotes directed cell movements in the ureteric bud tips, and suggest a model in which these cell movements mediate branching morphogenesis.},
  author       = {Riccio, Paul and Cebrián, Christina and Zong, Hui and Hippenmeyer, Simon and Costantini, Frank},
  publisher    = {Dryad},
  title        = {{Data from: Ret and Etv4 promote directed movements of progenitor cells during renal branching morphogenesis}},
  doi          = {10.5061/dryad.pk16b},
  year         = {2017},
}

@misc{9709,
  abstract     = {Across the nervous system, certain population spiking patterns are observed far more frequently than others. A hypothesis about this structure is that these collective activity patterns function as population codewords–collective modes–carrying information distinct from that of any single cell. We investigate this phenomenon in recordings of ∼150 retinal ganglion cells, the retina’s output. We develop a novel statistical model that decomposes the population response into modes; it predicts the distribution of spiking activity in the ganglion cell population with high accuracy. We found that the modes represent localized features of the visual stimulus that are distinct from the features represented by single neurons. Modes form clusters of activity states that are readily discriminated from one another. When we repeated the same visual stimulus, we found that the same mode was robustly elicited. These results suggest that retinal ganglion cells’ collective signaling is endowed with a form of error-correcting code–a principle that may hold in brain areas beyond retina.},
  author       = {Prentice, Jason and Marre, Olivier and Ioffe, Mark and Loback, Adrianna and Tkačik, Gašper and Berry, Michael},
  publisher    = {Dryad},
  title        = {{Data from: Error-robust modes of the retinal population code}},
  doi          = {10.5061/dryad.1f1rc},
  year         = {2017},
}

@misc{9842,
  abstract     = {Mathematica notebooks used to generate figures.},
  author       = {Etheridge, Alison and Barton, Nicholas H},
  publisher    = {Mendeley Data},
  title        = {{Data for: Establishment in a new habitat by polygenic adaptation}},
  doi          = {10.17632/nw68fxzjpm.1},
  year         = {2017},
}

@misc{9847,
  abstract     = {information on culture conditions, phage mutagenesis, verification and lysate preparation; Raw data},
  author       = {Pleska, Maros and Guet, Calin C},
  publisher    = {The Royal Society},
  title        = {{Supplementary materials and methods; Full data set from effects of mutations in phage restriction sites during escape from restriction–modification}},
  doi          = {10.6084/m9.figshare.5633917.v1},
  year         = {2017},
}

@misc{9853,
  abstract     = {Egg laying rates and infection loads of C. obscurior queens},
  author       = {Giehr, Julia and Grasse, Anna V and Cremer, Sylvia and Heinze, Jürgen and Schrempf, Alexandra},
  publisher    = {The Royal Society},
  title        = {{Raw data from ant queens increase their reproductive efforts after pathogen infection}},
  doi          = {10.6084/m9.figshare.5117788.v1},
  year         = {2017},
}

@misc{9859,
  abstract     = {Lists of all differentially expressed genes in the different priming-challenge treatments (compared to the fully naïve control; xlsx file). Relevant columns include the following: sample_1 and sample_2 – treatment groups being compared; Normalised FPKM sample_1 and sample_2 – FPKM of samples being compared; log2(fold_change) – log2(FPKM sample 2/FPKM sample 1), i.e. negative means sample 1 upregulated compared with sample 2, positive means sample 2 upregulated compared with sample 1; cuffdiff test_statistic – test statistic of differential expression test; p_value – p-value of differential expression test; q_value (FDR correction) – adjusted P-value of differential expression test. (XLSX 598 kb)},
  author       = {Greenwood, Jenny and Milutinovic, Barbara and Peuß, Robert and Behrens, Sarah and Essar, Daniela and Rosenstiel, Philip and Schulenburg, Hinrich and Kurtz, Joachim},
  publisher    = {Springer Nature},
  title        = {{Additional file 1: Table S1. of Oral immune priming with Bacillus thuringiensis induces a shift in the gene expression of Tribolium castaneum larvae}},
  doi          = {10.6084/m9.figshare.c.3756974_d1.v1},
  year         = {2017},
}

@misc{9860,
  author       = {Greenwood, Jenny and Milutinovic, Barbara and Peuß, Robert and Behrens, Sarah and Essar, Daniela and Rosenstiel, Philip and Schulenburg, Hinrich and Kurtz, Joachim},
  publisher    = {Springer Nature},
  title        = {{Additional file 5: Table S3. of Oral immune priming with Bacillus thuringiensis induces a shift in the gene expression of Tribolium castaneum larvae}},
  doi          = {10.6084/m9.figshare.c.3756974_d5.v1},
  year         = {2017},
}

@misc{9861,
  abstract     = {As a consequence of its difference in copy number between males and females, the X chromosome is subject to unique evolutionary forces and gene regulatory mechanisms. Previous studies of Drosophila melanogaster have shown that the expression of X-linked, testis-specific reporter genes is suppressed in the male germline. However, it is not known whether this phenomenon is restricted to testis-expressed genes or if it is a more general property of genes with tissue-specific expression, which are also underrepresented on the X chromosome. To test this, we compared the expression of three tissue-specific reporter genes (ovary, accessory gland and Malpighian tubule) inserted at various autosomal and X-chromosomal locations. In contrast to testis-specific reporter genes, we found no reduction of X-linked expression in any of the other tissues. In accessory gland and Malpighian tubule, we detected higher expression of the X-linked reporter genes, which suggests that they are at least partially dosage compensated. We found no difference in the tissue-specificity of X-linked and autosomal reporter genes. These findings indicate that, in general, the X chromosome is not a detrimental environment for tissue-specific gene expression and that the suppression of X-linked expression is limited to the male germline.},
  author       = {Argyridou, Eliza and Huylmans, Ann K and Königer, Annabella and Parsch, John},
  publisher    = {Dryad},
  title        = {{Data from: X-linkage is not a general inhibitor of tissue-specific gene expression in Drosophila melanogaster}},
  doi          = {10.5061/dryad.02f6r},
  year         = {2017},
}

@article{988,
  abstract     = {The current-phase relation (CPR) of a Josephson junction (JJ) determines how the supercurrent evolves with the superconducting phase difference across the junction. Knowledge of the CPR is essential in order to understand the response of a JJ to various external parameters. Despite the rising interest in ultraclean encapsulated graphene JJs, the CPR of such junctions remains unknown. Here, we use a fully gate-tunable graphene superconducting quantum intereference device (SQUID) to determine the CPR of ballistic graphene JJs. Each of the two JJs in the SQUID is made with graphene encapsulated in hexagonal boron nitride. By independently controlling the critical current of the JJs, we can operate the SQUID either in a symmetric or asymmetric configuration. The highly asymmetric SQUID allows us to phase-bias one of the JJs and thereby directly obtain its CPR. The CPR is found to be skewed, deviating significantly from a sinusoidal form. The skewness can be tuned with the gate voltage and oscillates in antiphase with Fabry-Pérot resistance oscillations of the ballistic graphene cavity. We compare our experiments with tight-binding calculations that include realistic graphene-superconductor interfaces and find a good qualitative agreement.},
  author       = {Nanda, Gaurav and Aguilera Servin, Juan L and Rakyta, Péter and Kormányos, Andor and Kleiner, Reinhold and Koelle, Dieter and Watanabe, Kazuo and Taniguchi, Takashi and Vandersypen, Lieven and Goswami, Srijit},
  issn         = {15306984},
  journal      = {Nano Letters},
  number       = {6},
  pages        = {3396 -- 3401},
  publisher    = {American Chemical Society},
  title        = {{Current-phase relation of ballistic graphene Josephson junctions}},
  doi          = {10.1021/acs.nanolett.7b00097},
  volume       = {17},
  year         = {2017},
}

@article{990,
  abstract     = {Assortative mating is an important driver of speciation in populations with gene flow and is predicted to evolve under certain conditions in few-locus models. However, the evolution of assortment is less understood for mating based on quantitative traits, which are often characterized by high genetic variability and extensive linkage disequilibrium between trait loci. We explore this scenario for a two-deme model with migration, by considering a single polygenic trait subject to divergent viability selection across demes, as well as assortative mating and sexual selection within demes, and investigate how trait divergence is shaped by various evolutionary forces. Our analysis reveals the existence of sharp thresholds of assortment strength, at which divergence increases dramatically. We also study the evolution of assortment via invasion of modifiers of mate discrimination and show that the ES assortment strength has an intermediate value under a range of migration-selection parameters, even in diverged populations, due to subtle effects which depend sensitively on the extent of phenotypic variation within these populations. The evolutionary dynamics of the polygenic trait is studied using the hypergeometric and infinitesimal models. We further investigate the sensitivity of our results to the assumptions of the hypergeometric model, using individual-based simulations.},
  author       = {Sachdeva, Himani and Barton, Nicholas H},
  issn         = {00143820},
  journal      = {Evolution; International Journal of Organic Evolution},
  number       = {6},
  pages        = {1478 -- 1493 },
  publisher    = {Wiley-Blackwell},
  title        = {{Divergence and evolution of assortative mating in a polygenic trait model of speciation with gene flow}},
  doi          = {10.1111/evo.13252},
  volume       = {71},
  year         = {2017},
}