An Excess of Risk-Increasing Low-Frequency Variants Can Be a Signal of Polygenic Inheritance in Complex Diseases

In most complex diseases, much of the heritability remains unaccounted for by common variants. It has been postulated that lower-frequency variants contribute to the remaining heritability. Here, we describe a method to test for polygenic inheritance from lower-frequency variants by using GWAS summary association statistics. We explored scenarios with many causal low-frequency variants and showed that there is more power to detect risk variants than to detect protective variants, resulting in an increase in the ratio of detected risk to protective variants (R/P ratio). Such an excess can also occur if risk variants are present and kept at lower frequencies because of negative selection. The R/P ratio can be falsely elevated because of reasons unrelated to polygenic inheritance, such as uneven sample sizes or asymmetric population stratification, so precautions to correct for these confounders are essential. We tested our method on published GWAS results and observed a strong signal in some diseases (schizophrenia and type 2 diabetes) but not others. We also explored the shared genetic component in overlapping phenotypes related to inflammatory bowel disease (Crohn disease [CD] and ulcerative colitis [UC]) and diabetic nephropathy (macroalbuminuria and end-stage renal disease [ESRD]). Although the signal was still present when both CD and UC were jointly analyzed, the signal was lost when macroalbuminuria and ESRD were jointly analyzed, suggesting that these phenotypes should best be studied separately. Thus, our method may also help guide the design of future genetic studies of various traits and diseases.

On the diversity of superluminous supernovae: ejected mass as the dominant factor

We assemble a sample of 24 hydrogen-poor superluminous supernovae(SLSNe). Parameterizing the light-curve shape through rise and declinetime-scales shows that the two are highly correlated. Magnetar-poweredmodels can reproduce the correlation, with the diversity in rise anddecline rates driven by the diffusion time-scale. Circumstellarinteraction models can exhibit a similar rise-decline relation, but onlyfor a narrow range of densities, which may be problematic for thesemodels. We find that SLSNe are approximately 3.5 mag brighter and havelight curves three times broader than SNe Ibc, but that the intrinsicshapes are similar. There are a number of SLSNe with particularly broadlight curves, possibly indicating two progenitor channels, butstatistical tests do not cleanly separate two populations. The generalspectral evolution is also presented. Velocities measured from Fe II aresimilar for SLSNe and SNe Ibc, suggesting that diffusion timedifferences are dominated by mass or opacity. Flat velocity evolution inmost SLSNe suggests a dense shell of ejecta. If opacities in SLSNe aresimilar to other SNe Ibc, the average ejected mass is higher by a factor2-3. Assuming κ = 0.1 cm2 g-1, we estimate amean (median) SLSN ejecta mass of 10 M⊙ (6M⊙), with a range of 3-30 M⊙. Doubling theassumed opacity brings the masses closer to normal SNe Ibc, but with ahigh-mass tail. The most probable mechanism for generating SLSNe seemsto be the core collapse of a very massive hydrogen-poor star, forming amillisecond magnetar.

The mitochondrial and autosomal mutation landscapes of prostate cancer

BACKGROUND: Prostate cancer (PCa) is the most common cancer in men. PCa is strongly age associated; low death rates in surveillance cohorts call into question the widespread use of surgery, which leads to overtreatment and a reduction in quality of life. There is a great need to increase the understanding of tumor characteristics in the context of disease progression.

OBJECTIVE: To perform the first multigenome investigation of PCa through analysis of both autosomal and mitochondrial DNA, and to integrate exome sequencing data, and RNA sequencing and copy-number alteration (CNA) data to investigate how various different tumor characteristics, commonly analyzed separately, are interconnected.

DESIGN, SETTING, AND PARTICIPANTS: Exome sequencing was applied to 64 tumor samples from 55 PCa patients with varying stage and grade. Integrated analysis was performed on a core set of 50 tumors from which exome sequencing, CNA, and RNA sequencing data were available.

OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS: Genes, mutated at a significantly higher rate relative to a genomic background, were identified. In addition, mitochondrial and autosomal mutation rates were correlated to CNAs and proliferation, assessed as a cell cycle gene expression signature.

RESULTS AND LIMITATIONS: Genes not previously reported to be significantly mutated in PCa, such as cell division cycle 27 homolog (Saccharomyces cerevisiae) (CDC27), myeloid/lymphoid or mixed-lineage leukemia 3 (MLL3), lysine (K)-specific demethylase 6A (KDM6A), and kinesin family member 5A (KIF5A) were identified. The mutation rate in the mitochondrial genome was 55 times higher than that of the autosomes. Multilevel analysis demonstrated a tight correlation between high reactive-oxygen exposure, chromosomal damage, high proliferation, and in parallel, a transition from multiclonal indolent primary PCa to monoclonal aggressive disease. As we only performed targeted sequence analysis; copy-number neutral rearrangements recently described for PCa were not accounted for.

CONCLUSIONS: The mitochondrial genome displays an elevated mutation rate compared to the autosomal chromosomes. By integrated analysis, we demonstrated that different tumor characteristics are interconnected, providing an increased understanding of PCa etiology.

Adaptor protein-2 sigma subunit mutations causing familial hypocalciuric hypercalcaemia type 3 (FHH3) demonstrate genotype-phenotype correlations, codon bias and dominant-negative effects

The adaptor protein-2 sigma subunit (AP2sigma;2) is pivotal for clathrin-mediated endocytosis of plasma membrane constituents such as the calcium-sensing receptor (CaSR). Mutations of the AP2sigma;2 Arg15 residue result in familial hypocalciuric hypercalcaemia type 3 (FHH3), a disorder of extracellular calcium (Ca<inf>o</inf>²⁺) homeostasis. To elucidate the role of AP2sigma;2 in Ca<inf>o</inf>²⁺ regulation, we investigated 65 FHH probands, without other FHH-associated mutations, for AP2sigma;2 mutations, characterized their functional consequences and investigated the genetic mechanisms leading to FHH3. AP2sigma;2 mutations were identified in 17 probands, comprising 5 Arg15Cys, 4 Arg15His and 8 Arg15Leu mutations. A genotype-phenotype correlation was observed with the Arg15Leu mutation leading to marked hypercalcaemia. FHH3 probands harboured additional phenotypes such as cognitive dysfunction. All three FHH3-causing AP2sigma;2 mutations impaired CaSR signal transduction in a dominant-negative manner. Mutational bias was observed at the AP2sigma;2 Arg15 residue as other predicted missense substitutions (Arg15Gly, Arg15Pro and Arg15Ser), which also caused CaSR loss-of-function, were not detected in FHH probands, and these mutations were found to reduce the numbers of CaSR-expressing cells. FHH3 probands had significantly greater serum calcium (sCa) and magnesium (sMg) concentrations with reduced urinary calcium to creatinine clearance ratios (CCCR) in comparison with FHH1 probands with CaSR mutations, and a calculated index of sCa × sMg/100 × CCCR, which was ≥ 5.0, had a diagnostic sensitivity and specificity of 83 and 86%, respectively, for FHH3. Thus, our studies demonstrate AP2sigma;2 mutations to result in a more severe FHH phenotype with genotype-phenotype correlations, and a dominant-negative mechanism of action with mutational bias at the Arg15 residue.

Extracting statistically dominant modes of a steel truss bridge from vehicle-induced vibrations

This study is intended to investigate the validity of the stability diagram (SD) aided multivariate autoregressive (MAR) analysis for identifying modal parameters of a real truss bridge. The MAR models are adopted to fit the time series of the dynamic accelerations recorded from a number of observation points on the bridge; then the modal parameters are extracted from the MAR model coefficient matrix. The SD is adopted to determine statistically dominant modes. In plotting the SD, a number of stability criteria are further adopted for filtering out those modes with unstable modal parameters. By the present method, the first five modal frequencies and mode shapes are identified with very high precision, while the damping ratios are identified with high precision for the 1st mode but with poorer precision for higher modes. Moreover, the ability of the SD in selecting structural modes without getting involved in any model-order optimization problem is highlighted through a comparison study.

CD4+ T Cells Targeting Dominant and Cryptic Epitopes from Bacillus anthracis Lethal Factor

Anthrax is an endemic infection in many countries, particularly in the developing world. The causative agent, Bacillus anthracis, mediates disease through the secretion of binary exotoxins. Until recently, research into adaptive immunity targeting this bacterial pathogen has largely focused on the humoral response to these toxins. There is, however, growing recognition that cellular immune responses involving IFNγ producing CD4+ T cells also contribute significantly to a protective memory response. An established concept in adaptive immunity to infection is that during infection of host cells, new microbial epitopes may be revealed, leading to immune recognition of so called 'cryptic' or 'subdominant' epitopes. We analyzed the response to both cryptic and immunodominant T cell epitopes derived from the toxin component lethal factor and presented by a range of HLA-DR alleles. Using IFNγ-ELISpot assays we characterized epitopes that elicited a response following immunization with synthetic peptide and the whole protein and tested their capacities to bind purified HLA-DR molecules in vitro. We found that DR1 transgenics demonstrated T cell responses to a greater number of domain III cryptic epitopes than other HLA-DR transgenics, and that this pattern was repeated with the immunodominant epitopes, as a greater proportion of these epitopes induced a T cell response when presented within the context of the whole protein. Immunodominant epitopes LF457-476 and LF467-487 were found to induce a T cell response to the peptide, as well as to the whole native LF protein in DR1 and DR15, but not in DR4 transgenics. The analysis of Domain I revealed the presence of several unique cryptic epitopes all of which showed a strong to moderate relative binding affinity to HLA-DR4 molecules. However, none of the cryptic epitopes from either domain III or I displayed notably high binding affinities across all HLA-DR alleles assayed. These responses were influenced by the specific HLA alleles presenting the peptide, and imply that construction of future epitope string vaccines which are immunogenic across a wide range of HLA alleles could benefit from a combination of both cryptic and immunodominant anthrax epitopes.

Quantification of HER2 heterogeneity in breast cancer-implications for identification of sub-dominant clones for personalised treatment

Breast cancer is a heterogeneous disease, at both an inter- and intra-tumoural level. Appreciating heterogeneity through the application of biomarkers and molecular signatures adds complexity to tumour taxonomy but is key to personalising diagnosis, treatment and prognosis. The extent to which heterogeneity exists, and its interpretation remains a challenge to pathologists. Using HER2 as an exemplar, we have developed a simple reproducible heterogeneity index. Cell-to-cell HER2 heterogeneity was extensive in a proportion of both reported 'amplified' and 'non-amplified' cases. The highest levels of heterogeneity objectively identified occurred in borderline categories and higher ratio non-amplified cases. A case with particularly striking heterogeneity was analysed further with an array of biomarkers in order to assign a molecular diagnosis. Broad biological complexity was evident. In essence, interpretation, depending on the area of tumour sampled, could have been one of three distinct phenotypes, each of which would infer different therapeutic interventions. Therefore, we recommend that heterogeneity is assessed and taken into account when determining treatment options.

Identity of a Tilapia Pheromone Released by Dominant Males that Primes Females for Reproduction

Knowledge of the chemical identity and role of urinary pheromones in fish is scarce, yet it is necessary in order to understand the integration of multiple senses in adaptive responses and the evolution of chemical communication [1]. In nature, Mozambique tilapia (Oreochromis mossambicus) males form hierarchies, and females mate preferentially with dominant territorial males, which they visit in aggregations or leks [2]. Dominant males have thicker urinary bladder muscular walls than subordinates or females and store large volumes of urine, which they release at increased frequency in the presence of subordinate males or preovulatory, but not postspawned, females [3–5]. Females exposed to dominant-male urine augment their release of the oocyte maturation-inducing steroid 17α,20β-dihydroxypregn-4-en-3-one (17,20β-P) [6]. Here we isolate and identify a male Mozambique tilapia urinary sex pheromone as two epimeric (20α- and 20β-) pregnanetriol 3-glucuronates. We show that both males and females have high olfactory sensitivity to the two steroids, which cross-adapt upon stimulation. Females exposed to both steroids show a rapid, 10-fold increase in production of 17,20β-P. Thus, the identified urinary steroids prime the female endocrine system to accelerate oocyte maturation and possibly promote spawning synchrony. Tilapia are globally important as a food source but are also invasive species, with devastating impact on local freshwater ecosystems [7, 8]. Identifying the chemical cues that mediate reproduction may lead to the development of tools for population control [9–11].