Reconsidering Association Testing Methods Using Single-Variant Test Statistics as Alternatives to Pooling Tests for Sequence Data with Rare Variants

Association tests that pool minor alleles into a measure of burden at a locus have been proposed for case-control studies using sequence data containing rare variants. However, such pooling tests are not robust to the inclusion of neutral and protective variants, which can mask the association signal from risk variants. Early studies proposing pooling tests dismissed methods for locus-wide inference using nonnegative single-variant test statistics based on unrealistic comparisons. However, such methods are robust to the inclusion of neutral and protective variants and therefore may be more useful than previously appreciated. In fact, some recently proposed methods derived within different frameworks are equivalent to performing inference on weighted sums of squared single-variant score statistics. In this study, we compared two existing methods for locus-wide inference using nonnegative single-variant test statistics to two widely cited pooling tests under more realistic conditions. We established analytic results for a simple model with one rare risk and one rare neutral variant, which demonstrated that pooling tests were less powerful than even Bonferroni-corrected single-variant tests in most realistic situations. We also performed simulations using variants with realistic minor allele frequency and linkage disequilibrium spectra, disease models with multiple rare risk variants and extensive neutral variation, and varying rates of missing genotypes. In all scenarios considered, existing methods using nonnegative single-variant test statistics had power comparable to or greater than two widely cited pooling tests. Moreover, in disease models with only rare risk variants, an existing method based on the maximum single-variant Cochran-Armitage trend chi-square statistic in the locus had power comparable to or greater than another existing method closely related to some recently proposed methods. We conclude that efficient locus-wide inference using single-variant test statistics should be reconsidered as a useful framework for devising powerful association tests in sequence data with rare variants

Genetic architecture and evolution of the S locus supergene in Primula vulgaris

Darwin’s studies on heterostyly in Primula described two floral morphs, pin and thrum, with reciprocal anther and stigma heights that promote insect-mediated cross-pollination. This key innovation evolved independently in several angiosperm families. Subsequent studies on heterostyly in Primula contributed to the foundation of modern genetic theory and the neo-Darwinian synthesis. The established genetic model for Primula heterostyly involves a diallelic S locus comprising several genes, with rare recombination events that result in self-fertile homostyle flowers with anthers and stigma at the same height. Here we reveal the S locus supergene as a tightly-linked cluster of thrum-specific genes that are absent in pins. We show that thrums are hemizygous not heterozygous for the S locus, which suggests that homostyles do not arise by recombination between S locus haplotypes as previously proposed. Duplication of a floral homeotic gene 51.7 MYA, followed by its neofunctionalisation, created the current S locus assemblage which led to floral heteromorphy in Primula. Our findings provide new insights into the structure, function and evolution of this archetypal supergene

Conditions for the Evolution of Gene Clusters in Bacterial Genomes

Genes encoding proteins in a common pathway are often found near each other along bacterial chromosomes. Several explanations have been proposed to account for the evolution of these structures. For instance, natural selection may directly favour gene clusters through a variety of mechanisms, such as increased efficiency of coregulation. An alternative and controversial hypothesis is the selfish operon model, which asserts that clustered arrangements of genes are more easily transferred to other species, thus improving the prospects for survival of the cluster. According to another hypothesis (the persistence model), genes that are in close proximity are less likely to be disrupted by deletions. Here we develop computational models to study the conditions under which gene clusters can evolve and persist. First, we examine the selfish operon model by re-implementing the simulation and running it under a wide range of conditions. Second, we introduce and study a Moran process in which there is natural selection for gene clustering and rearrangement occurs by genome inversion events. Finally, we develop and study a model that includes selection and inversion, which tracks the occurrence and fixation of rearrangements. Surprisingly, gene clusters fail to evolve under a wide range of conditions. Factors that promote the evolution of gene clusters include a low number of genes in the pathway, a high population size, and in the case of the selfish operon model, a high horizontal transfer rate. The computational analysis here has shown that the evolution of gene clusters can occur under both direct and indirect selection as long as certain conditions hold. Under these conditions the selfish operon model is still viable as an explanation for the evolution of gene clusters

A Novel Adaptive Method for the Analysis of Next-Generation Sequencing Data to Detect Complex Trait Associations with Rare Variants Due to Gene Main Effects and Interactions

There is solid evidence that rare variants contribute to complex disease etiology. Next-generation sequencing technologies make it possible to uncover rare variants within candidate genes, exomes, and genomes. Working in a novel framework, the kernel-based adaptive cluster (KBAC) was developed to perform powerful gene/locus based rare variant association testing. The KBAC combines variant classification and association testing in a coherent framework. Covariates can also be incorporated in the analysis to control for potential confounders including age, sex, and population substructure. To evaluate the power of KBAC: 1) variant data was simulated using rigorous population genetic models for both Europeans and Africans, with parameters estimated from sequence data, and 2) phenotypes were generated using models motivated by complex diseases including breast cancer and Hirschsprung's disease. It is demonstrated that the KBAC has superior power compared to other rare variant analysis methods, such as the combined multivariate and collapsing and weight sum statistic. In the presence of variant misclassification and gene interaction, association testing using KBAC is particularly advantageous. The KBAC method was also applied to test for associations, using sequence data from the Dallas Heart Study, between energy metabolism traits and rare variants in ANGPTL 3,4,5 and 6 genes. A number of novel associations were identified, including the associations of high density lipoprotein and very low density lipoprotein with ANGPTL4. The KBAC method is implemented in a user-friendly R package

Aggregating behaviour in invasive Caribbean lionfish is driven by habitat complexity

Caribbean lionfish (Pterois spp.) are considered the most heavily impacting invasive marine vertebrate ever recorded. However, current management is largely inadequate, relying on opportunistic culling by recreational SCUBA divers. Culling efficiency could be greatly improved by exploiting natural aggregations, but to date this behaviour has only been recorded anecdotally, and the drivers are unknown. We found aggregations to be common in situ, but detected no conspecific attraction through visual or olfactory cues in laboratory experiments. Aggregating individuals were on average larger, but showed no further differences in morphology or life history. However, using visual assessments and 3D modelling we show lionfish prefer broad-scale, but avoid fine-scale, habitat complexity. We therefore suggest that lionfish aggregations are coincidental based on individuals’ mutual attraction to similar reef structure to maximise hunting efficiency. Using this knowledge, artificial aggregation devices might be developed to concentrate lionfish densities and thus improve culling efficiency

Statin Treatment Increases Lifespan and Improves Cardiac Health in Drosophila by Decreasing Specific Protein Prenylation

Statins such as simvastatin are 3-hydroxy-3-methylglutaryl coenzyme A (HMG-CoA) reductase inhibitors and standard therapy for the prevention and treatment of cardiovascular diseases in mammals. Here we show that simvastatin significantly increased the mean and maximum lifespan of Drosophila melanogaster (Drosophila) and enhanced cardiac function in aging flies by significantly reducing heart arrhythmias and increasing the contraction proportion of the contraction/relaxation cycle. These results appeared independent of internal changes in ubiquinone or juvenile hormone levels. Rather, they appeared to involve decreased protein prenylation. Simvastatin decreased the membrane association (prenylation) of specific small Ras GTPases in mice. Both farnesyl (L744832) and type 1 geranylgeranyl transferase (GGTI-298) inhibitors increased Drosophila lifespan. These data are the most direct evidence to date that decreased protein prenylation can increase cardiac health and lifespan in any metazoan species, and may explain the pleiotropic (non-cholesterol related) health effects of statins

Primula vulgaris (primrose) genome assembly, annotation and gene expression, with comparative genomics on the heterostyly supergene

Primula vulgaris (primrose) exhibits heterostyly: plants produce self-incompatible pin- or thrum-form flowers, with anthers and stigma at reciprocal heights. Darwin concluded that this arrangement promotes insect-mediated cross-pollination; later studies revealed control by a cluster of genes, or supergene, known as the S (Style length) locus. The P. vulgaris S locus is absent from pin plants and hemizygous in thrum plants (thrum-specific); mutation of S locus genes produces self-fertile homostyle flowers with anthers and stigma at equal heights. Here, we present a 411 Mb P. vulgaris genome assembly of a homozygous inbred long homostyle, representing ~87% of the genome. We annotate over 24,000 P. vulgaris genes, and reveal more genes up-regulated in thrum than pin flowers. We show reduced genomic read coverage across the S locus in other Primula species, including P. veris, where we define the conserved structure and expression of the S locus genes in thrum. Further analysis reveals the S locus has elevated repeat content (64%) compared to the wider genome (37%). Our studies suggest conservation of S locus genetic architecture in Primula, and provide a platform for identification and evolutionary analysis of the S locus and downstream targets that regulate heterostyly in diverse heterostylous species

Measuring local depletion of terrestrial game vertebrates by central-place hunters in rural Amazonia

The degree to which terrestrial vertebrate populations are depleted in tropical forests occupied by human communities has been the subject of an intense polarising debate that has important conservation implications. Conservation ecologists and practitioners are divided over the extent to which community-based subsistence offtake is compatible with ecologically functional populations of tropical forest game species. To quantify depletion envelopes of forest vertebrates around human communities, we deployed a total of 383 camera trap stations and 78 quantitative interviews to survey the peri-community areas controlled by 60 semi-subsistence communities over a combined area of over 3.2 million hectares in the Médio Juruá and Uatumã regions of Central-Western Brazilian Amazonia. Our results largely conform with prior evidence that hunting large-bodied vertebrates reduces wildlife populations near settlements, such that they are only found at a distance to settlements where they are hunted less frequently. Camera trap data suggest that a select few harvest-sensitive species, including lowland tapir, are either repelled or depleted by human communities. Nocturnal and cathemeral species were detected relatively more frequently in disturbed areas close to communities, but individual species did not necessarily shift their activity patterns. Group biomass of all species was depressed in the wider neighbourhood of urban areas rather than communities. Interview data suggest that species traits, especially group size and body mass, mediate these relationships. Large-bodied, large-group-living species are detected farther from communities as reported by experienced informants. Long-established communities in our study regions have not “emptied” the surrounding forest. Low human population density and low hunting offtake due to abundant sources of alternative aquatic protein, suggest that these communities represent a best-case scenario for sustainable hunting of wildlife for food, thereby providing a conservative assessment of game depletion. Given this ‘best-case’ camera trap and interview-based evidence for hunting depletion, regions with higher human population densities, external trade in wildlife and limited access to alternative protein will likely exhibit more severe depletion

Making sense of policy choices: understanding the roles of value predispositions, mass media, and cognitive processing in public attitudes toward nanotechnology

Using a nationally representative telephone survey of 1,015 adults in the United States, this study examines how value predispositions, communication variables, and perceptions of risks and benefits are associated with public support for federal funding of nanotechnology. Our findings show that highly religious individuals were less supportive of funding of nanotech than less religious individuals, whereas individuals who held a high deference for scientific authority were more supportive of funding of the emerging technology than those low in deference. Mass media use and elaborative processing of scientific news were positively associated with public support for funding, whereas factual scientific knowledge had no significant association with policy choices. The findings suggest that thinking about and reflecting upon scientific news promote better understanding of the scientific world and may provide a more sophisticated cognitive structure for the public to form opinions about nanotech than factual scientific knowledge. Finally, heuristic cues including trust in scientists and perceived risks and benefits of nanotech were found to be associated with public support for nanotech funding. We conclude with policy implications that will be useful for policymakers and science communication practitioners

Complex SUMO-1 Regulation of Cardiac Transcription Factor Nkx2-5

Reversible post-translational protein modifications such as SUMOylation add complexity to cardiac transcriptional regulation. The homeodomain transcription factor Nkx2-5/Csx is essential for heart specification and morphogenesis. It has been previously suggested that SUMOylation of lysine 51 (K51) of Nkx2-5 is essential for its DNA binding and transcriptional activation. Here, we confirm that SUMOylation strongly enhances Nkx2-5 transcriptional activity and that residue K51 of Nkx2-5 is a SUMOylation target. However, in a range of cultured cell lines we find that a point mutation of K51 to arginine (K51R) does not affect Nkx2-5 activity or DNA binding, suggesting the existence of additional Nkx2-5 SUMOylated residues. Using biochemical assays, we demonstrate that Nkx2-5 is SUMOylated on at least one additional site, and this is the predominant site in cardiac cells. The second site is either non-canonical or a “shifting” site, as mutation of predicted consensus sites and indeed every individual lysine in the context of the K51R mutation failed to impair Nkx2-5 transcriptional synergism with SUMO, or its nuclear localization and DNA binding. We also observe SUMOylation of Nkx2-5 cofactors, which may be critical to Nkx2-5 regulation. Our data reveal highly complex regulatory mechanisms driven by SUMOylation to modulate Nkx2-5 activity