Visual analytics for large-scale bioinformatic data sets

Rapid advances in sequencing technologies (Next Generation Sequencing or NGS) have led to a vast increase in the quantity of bioinformatics data available, with this increasing scale presenting enormous challenges to researchers seeking to identify complex interactions. This paper is concerned with the domain of transcriptional regulation, and the use of visualisation to identify relationships between specific regulatory proteins (the transcription factors or TFs) and their associated target genes (TGs). We present preliminary work from an ongoing study which aims to determine the effectiveness of different visual representations and large scale displays in supporting discovery. Following an iterative process of implementation and evaluation, representations were tested by potential users in the bioinformatics domain to determine their efficacy, and to understand better the range of ad hoc practices among bioinformatics literate users. Results from two rounds of small scale user studies are considered with initial findings suggesting that bioinformaticians require richly detailed views of TF data, features to compare TF layouts between organisms quickly, and ways to keep track of interesting data points.

A simple method for estimating growth parameters from multiple length-frequency data in presence of continuous recruitment

The extended recruitment season for short-lived species such as prawns biases the estimation of growth parameters from length-frequency data when conventional methods are used. We propose a simple method for overcoming this bias given a time series of length-frequency data. The difficulties arising from extended recruitment are eliminated by predicting the growth of the succeeding samples and the length increments of the recruits in previous samples. This method requires that some maximum size at recruitment can be specified. The advantages of this multiple length-frequency method are: it is simple to use; it requires only three parameters; no specific distributions need to be assumed; and the actual seasonal recruitment pattern does not have to be specified. We illustrate the new method with length-frequency data on the tiger prawn Penaeus esculentus from the north-western Gulf of Carpentaria, Australia.

Maximum likelihood estimation of mortality and growth with individual variability from multiple length-frequency data

We consider estimation of mortality rates and growth parameters from length-frequency data of a fish stock and derive the underlying length distribution of the population and the catch when there is individual variability in the von Bertalanffy growth parameter L-infinity. The model is flexible enough to accommodate 1) any recruitment pattern as a function of both time and length, 2) length-specific selectivity, and 3) varying fishing effort over time. The maximum likelihood method gives consistent estimates, provided the underlying distribution for individual variation in growth is correctly specified. Simulation results indicate that our method is reasonably robust to violations in the assumptions. The method is applied to tiger prawn data (Penaeus semisulcatus) to obtain estimates of natural and fishing mortality.

A large-scale analysis of genetic variants within putative miRNA binding sites in prostate cancer

Prostate cancer is the second most common malignancy among men worldwide. Genome-wide association studies have identified 100 risk variants for prostate cancer, which can explain approximately 33% of the familial risk of the disease. We hypothesized that a comprehensive analysis of genetic variations found within the 3' untranslated region of genes predicted to affect miRNA binding (miRSNP) can identify additional prostate cancer risk variants. We investigated the association between 2,169 miRSNPs and prostate cancer risk in a large-scale analysis of 22,301 cases and 22,320 controls of European ancestry from 23 participating studies. Twenty-two miRSNPs were associated (P<2.3×10(-5)) with risk of prostate cancer, 10 of which were within 7 genes previously not mapped by GWAS studies. Further, using miRNA mimics and reporter gene assays, we showed that miR-3162-5p has specific affinity for the KLK3 rs1058205 miRSNP T-allele, whereas miR-370 has greater affinity for the VAMP8 rs1010 miRSNP A-allele, validating their functional role. SIGNIFICANCE Findings from this large association study suggest that a focus on miRSNPs, including functional evaluation, can identify candidate risk loci below currently accepted statistical levels of genome-wide significance. Studies of miRNAs and their interactions with SNPs could provide further insights into the mechanisms of prostate cancer risk.

Effect of individual variability on estimation of population parameters from length-frequency data

We consider estimation of mortality rates and growth parameters from length-frequency data of a fish stock when there is individual variability in the von Bertalanffy growth parameter L-infinity and investigate the possible bias in the estimates when the individual variability is ignored. Three methods are examined: (i) the regression method based on the Beverton and Holt's (1956, Rapp. P.V. Reun. Cons. Int. Explor. Mer, 140: 67-83) equation; (ii) the moment method of Powell (1979, Rapp. PV. Reun. Int. Explor. Mer, 175: 167-169); and (iii) a generalization of Powell's method that estimates the individual variability to be incorporated into the estimation. It is found that the biases in the estimates from the existing methods are, in general, substantial, even when individual variability in growth is small and recruitment is uniform, and the generalized method performs better in terms of bias but is subject to a larger variation. There is a need to develop robust and flexible methods to deal with individual variability in the analysis of length-frequency data.

Workplace belongingness, psychological distress and resilience in emergency service workers

- Objective Ambulance personnel provide emergency medical services to the community, often attending to highly challenging and traumatic scenes in complex and chaotic circumstances. Currently the assessment of predictors of psychological well-being remains limited. The current study investigated whether workplace belongingness was significant in predicting psychological distress as well as the presence of resilience in ambulance personnel whilst controlling for more routinely examined factors. - Method Australian ambulance officers (N = 740) completed a survey battery including the Kessler 10 (Kessler & Mroczek, 1994), Brief Resilience Scale (Smith et al., 2008) and Psychological Sense of Organisational Membership (Cockshaw & Shochet, 2010) scale. - Results Controlling for more commonly examined factors such as severity of trauma exposure and length of service, hierarchical multiple regression analyses demonstrated that workplace belongingness was significantly associated with reduced distress levels and enhanced resilience levels. - Conclusions Results suggest that strategies to enhance a sense of workplace belongingness in emergency service organisations could promote the well-being of emergency workers despite routine exposure to potentially traumatic events.

A study of the use of the Edinburgh postnatal depression scale with parent groups outside the postpartum period

The Edinburgh Postnatal Depression Scale (EPDS) was sent by post to 206 mothers and 201 fathers of toddlers (aged between 19 and 22 months). At the same time these parents also completed subscales of the Crown—Crisp Experiential Index (CCEI). The responses were used to assess the feasibility of postal completion of the EPDS and its acceptability to parents outside the postpartum year, particularly fathers for whom there have been no previous reports of its use. On a small sub-group, the sensitivity, specificity and predictive values of the measures were assessed using the Present. State Examination. Answers to the depression subscale of the CCEI to the EPDS and to the Present State Examination were compared to assess validity. Completion of the postally-administered EPDS was satisfactory, though some difficulties were experienced in a second postal administration to a subsample. The scale was completed without obvious error or omission and this, combined with positive comments from parents, suggests the acceptability of the scale to both mothers and fathers. The mean scores were higher for mothers than for fathers, but the pattern of distribution was similar with a marked positive skew and a distinct decline in scores above 10. Because the subsample of parents interviewed was small the calculation of sensitivity and specificity has to be treated with caution. However, the results for mothers suggest that the EPDS has satisfactory validity for this group and one superior to the depression subscale of the CCEI. Among the fathers interviewed there were insufficient cases to enable calculation of sensitivity and specificity. Other results were encouraging, however, and suggest the merit of further studies of the application and validity of the EPDS with fathers.

Corneal confocal microscopy predicts 4-year incident peripheral neuropathy in Type 1 diabetes

OBJECTIVE This study determined if deficits in corneal nerve fiber length (CNFL) assessed using corneal confocal microscopy (CCM) can predict future onset of diabetic peripheral neuropathy (DPN). RESEARCH DESIGN AND METHODS CNFL and a range of other baseline measures were compared between 90 nonneuropathic patients with type 1 diabetes who did or did not develop DPN after 4 years. The receiver operator characteristic (ROC) curve was used to determine the capability of single and combined measures of neuropathy to predict DPN. RESULTS DPN developed in 16 participants (18%) after 4 years. Factors predictive of 4-year incident DPN were lower CNFL (P = 0.041); longer duration of diabetes (P = 0.002); higher triglycerides (P = 0.023); retinopathy (higher on the Early Treatment of Diabetic Retinopathy Study scale) (P = 0.008); nephropathy (higher albumin-to-creatinine ratio) (P = 0.001); higher neuropathy disability score (P = 0.037); lower cold sensation (P = 0.001) and cold pain (P = 0.027) thresholds; higher warm sensation (P = 0.008), warm pain (P = 0.024), and vibration (P = 0.003) thresholds; impaired monofilament response (P = 0.003); and slower peroneal (P = 0.013) and sural (P = 0.002) nerve conduction velocity. CCM could predict the 4-year incident DPN with 63% sensitivity and 74% specificity for a CNFL threshold cutoff of 14.1 mm/mm2 (area under ROC curve = 0.66, P = 0.041). Combining neuropathy measures did not improve predictive capability. CONCLUSIONS DPN can be predicted by various demographic, metabolic, and conventional neuropathy measures. The ability of CCM to predict DPN broadens the already impressive diagnostic capabilities of this novel ophthalmic marker.

Examination of a scale assessing attitudes towards individuals with intellectual disability in China

This study examined the applicability of the four-factor structure of the short form of the Community Living Attitudes Scale-Intellectual disability1 (CLAS-ID) in China, using a sample of 325 Chinese community members. Confirmatory factor analysis revealed that the original structure of the short form of the CLAS-ID did not adequately fit the data from the current sample. Most items of the Exclusion and Similarity subscales were retained while items on the Empowerment and Sheltering subscales were removed. Chinese community members held generally positive attitudes towards people with intellectual disability. However, a measurement tool originating from the Chinese context is needed to provide a better understanding of attitudes towards individuals with intellectual disability in mainland China.

The Connor-Davidson Resilience Scale (CD-RISC): Dimensionality and age-related measurement invariance with Australian cricketers

Objective This study explored the dimensionality and measurement invariance of the 25-item Connor-Davidson Resilience Scale (CD-RISC; Connor & Davidson, 2003) across samples of adult (n = 321; aged 20–36) and adolescent (n = 199; aged 12–18) Australian cricketers. Design Cross-sectional, self-report survey Methods An online, multi-section questionnaire. Results Confirmatory factor and item level analyses supported the psychometric superiority of a revised 10-item, unidimensional model of resilience over the original 25-item, five-factor measurement model. Positive and moderate correlations with hardiness as well as negative and moderate correlations with burnout components were evidenced thereby providing support for the convergent validity of the unidimensional model. Measurement invariance analyses of the unidimensional model across the two age-group samples supported configural (i.e., same factor structure across groups), metric (i.e., same pattern of factor loadings across the groups), and partial scalar invariance (i.e., mostly the same intercepts across the groups). Conclusion Evidence for a psychometrically sound measure of resilient qualities of the individual provides an important foundation upon which researchers can identify the antecedents to and outcomes of resilience in sport contexts.

Low birth weight in MZ twins discordant for birth weight is associated with shorter telomere length and lower IQ, but not anxiety/depression in later life

Shorter telomere length (TL) has found to be associated with lower birth weight and with lower cognitive ability and psychiatric disorders. However, the direction of causation of these associations and the extent to which they are genetically or environmentally mediated are unclear. Within-pair comparisons of monozygotic (MZ) and dizygotic (DZ) twins can throw light on these questions. We investigated correlations of within pair differences in telomere length, IQ, and anxiety/depression in an initial sample from Brisbane (242 MZ pairs, 245 DZ same sex (DZSS) pairs) and in replication samples from Amsterdam (514 MZ pairs, 233 DZSS pairs) and Melbourne (19 pairs selected for extreme high or low birth weight difference). Intra-pair differences of birth weight and telomere length were significantly correlated in MZ twins, but not in DZSS twins. Greater intra-pair differences of telomere length were observed in the 10% of MZ twins with the greatest difference in birth weight compared to the bottom 90% in both samples and also in the Melbourne sample. Intra-pair differences of telomere length and IQ, but not of TL and anxiety/depression, were correlated in MZ twins, and to a smaller extent in DZSS twins. Our findings suggest that the same prenatal effects that reduce birth weight also influence telomere length in MZ twins. The association between telomere length and IQ is partly driven by the same prenatal effects that decrease birth weight.

DCAF4, a novel gene associated with leucocyte telomere length

Background Leucocyte telomere length (LTL), which is fashioned by multiple genes, has been linked to a host of human diseases, including sporadic melanoma. A number of genes associated with LTL have already been identified through genome-wide association studies. The main aim of this study was to establish whether DCAF4 (DDB1 and CUL4-associated factor 4) is associated with LTL. In addition, using ingenuity pathway analysis (IPA), we examined whether LTL-associated genes in the general population might partially explain the inherently longer LTL in patients with sporadic melanoma, the risk for which is increased with ultraviolet radiation (UVR). Results Genome-wide association (GWA) meta-analysis and de novo genotyping of 20 022 individuals revealed a novel association (p=6.4×10−10) between LTL and rs2535913, which lies within DCAF4. Notably, eQTL analysis showed that rs2535913 is associated with decline in DCAF4 expressions in both lymphoblastoid cells and sun-exposed skin (p=4.1×10−3 and 2×10−3, respectively). Moreover, IPA revealed that LTL-associated genes, derived from GWA meta-analysis (N=9190), are over-represented among genes engaged in melanoma pathways. Meeting increasingly stringent p value thresholds (p<0.05, <0.01, <0.005, <0.001) in the LTL-GWA meta-analysis, these genes were jointly over-represented for melanoma at p values ranging from 1.97×10−169 to 3.42×10−24. Conclusions We uncovered a new locus associated with LTL in the general population. We also provided preliminary findings that suggest a link of LTL through genetic mechanisms with UVR and melanoma in the general population.

Identification of seven loci affecting mean telomere length and their association with disease

Interindividual variation in mean leukocyte telomere length (LTL) is associated with cancer and several age-associated diseases. We report here a genome-wide meta-analysis of 37,684 individuals with replication of selected variants in an additional 10,739 individuals. We identified seven loci, including five new loci, associated with mean LTL (P < 5 x 10(-8)). Five of the loci contain candidate genes (TERC, TERT, NAF1, OBFC1 and RTEL1) that are known to be involved in telomere biology. Lead SNPs at two loci (TERC and TERT) associate with several cancers and other diseases, including idiopathic pulmonary fibrosis. Moreover, a genetic risk score analysis combining lead variants at all 7 loci in 22,233 coronary artery disease cases and 64,762 controls showed an association of the alleles associated with shorter LTL with increased risk of coronary artery disease (21% (95% confidence interval, 5-35%) per standard deviation in LTL, P = 0.014). Our findings support a causal role of telomere-length variation in some age-related diseases.

Meta-analysis of telomere length in 19 713 subjects reveals high heritability, stronger maternal inheritance and a paternal age effect

Telomere length (TL) has been associated with aging and mortality, but individual differences are also influenced by genetic factors, with previous studies reporting heritability estimates ranging from 34 to 82%. Here we investigate the heritability, mode of inheritance and the influence of parental age at birth on TL in six large, independent cohort studies with a total of 19 713 participants. The meta-analysis estimate of TL heritability was 0.70 (95% CI 0.64–0.76) and is based on a pattern of results that is highly similar for twins and other family members. We observed a stronger mother–offspring (r=0.42; P-value=3.60 × 10−61) than father–offspring correlation (r=0.33; P-value=7.01 × 10−5), and a significant positive association with paternal age at offspring birth (β=0.005; P-value=7.01 × 10−5). Interestingly, a significant and quite substantial correlation in TL between spouses (r=0.25; P-value=2.82 × 10−30) was seen, which appeared stronger in older spouse pairs (mean age ≥55 years; r=0.31; P-value=4.27 × 10−23) than in younger pairs (mean age<55 years; r=0.20; P-value=3.24 × 10−10). In summary, we find a high and very consistent heritability estimate for TL, evidence for a maternal inheritance component and a positive association with paternal age.

A variant in LIN28B is associated with 2D:4D finger-length ratio, a putative retrospective biomarker of prenatal testosterone exposure

The ratio of the lengths of an individual's second to fourth digit (2D:4D) is commonly used as a noninvasive retrospective biomarker for prenatal androgen exposure. In order to identify the genetic determinants of 2D:4D, we applied a genome-wide association approach to 1507 11-year-old children from the Avon Longitudinal Study of Parents and Children (ALSPAC) in whom 2D:4D ratio had been measured, as well as a sample of 1382 12- to 16-year-olds from the Brisbane Adolescent Twin Study. A meta-analysis of the two scans identified a single variant in the LIN28B gene that was strongly associated with 2D:4D (rs314277: p = 4.1 x 10(-8)) and was subsequently independently replicated in an additional 3659 children from the ALSPAC cohort (p = 1.53 x 10(-6)). The minor allele of the rs314277 variant has previously been linked to increased height and delayed age at menarche, but in our study it was associated with increased 2D:4D in the direction opposite to that of previous reports on the correlation between 2D:4D and age at menarche. Our findings call into question the validity of 2D:4D as a simplistic retrospective biomarker for prenatal testosterone exposure.