Divergent sensory and immune gene evolution in sea turtles with contrasting demographic and life histories.

Sea turtles represent an ancient lineage of marine vertebrates that evolved from terrestrial ancestors over 100 Mya. The genomic basis of the unique physiological and ecological traits enabling these species to thrive in diverse marine habitats remains largely unknown. Additionally, many populations have drastically declined due to anthropogenic activities over the past two centuries, and their recovery is a high global conservation priority. We generated and analyzed high-quality reference genomes for the leatherback (Dermochelys coriacea) and green (Chelonia mydas) turtles, representing the two extant sea turtle families. These genomes are highly syntenic and homologous, but localized regions of noncollinearity were associated with higher copy numbers of immune, zinc-finger, and olfactory receptor (OR) genes in green turtles, with ORs related to waterborne odorants greatly expanded in green turtles. Our findings suggest that divergent evolution of these key gene families may underlie immunological and sensory adaptations assisting navigation, occupancy of neritic versus pelagic environments, and diet specialization. Reduced collinearity was especially prevalent in microchromosomes, with greater gene content, heterozygosity, and genetic distances between species, supporting their critical role in vertebrate evolutionary adaptation. Finally, diversity and demographic histories starkly contrasted between species, indicating that leatherback turtles have had a low yet stable effective population size, exhibit extremely low diversity compared with other reptiles, and harbor a higher genetic load compared with green turtles, reinforcing concern over their persistence under future climate scenarios. These genomes provide invaluable resources for advancing our understanding of evolution and conservation best practices in an imperiled vertebrate lineage.

Genome-wide association study meta-analysis of dizygotic twinning illuminates genetic regulation of female fecundity.

STUDY QUESTION: Which genetic factors regulate female propensity for giving birth to spontaneous dizygotic (DZ) twins? SUMMARY ANSWER: We identified four new loci, GNRH1, FSHR, ZFPM1, and IPO8, in addition to previously identified loci, FSHB and SMAD3. WHAT IS KNOWN ALREADY: The propensity to give birth to DZ twins runs in families. Earlier, we reported that FSHB and SMAD3 as associated with DZ twinning and female fertility measures. STUDY DESIGN, SIZE, DURATION: We conducted a genome-wide association meta-analysis (GWAMA) of mothers of spontaneous dizygotic (DZ) twins (8265 cases, 264 567 controls) and of independent DZ twin offspring (26 252 cases, 417 433 controls). PARTICIPANTS/MATERIALS, SETTING, METHODS: Over 700 000 mothers of DZ twins, twin individuals and singletons from large cohorts in Australia/New Zealand, Europe, and the USA were carefully screened to exclude twins born after use of ARTs. Genetic association analyses by cohort were followed by meta-analysis, phenome wide association studies (PheWAS), in silico and in vivo annotations, and Zebrafish functional validation. MAIN RESULTS AND THE ROLE OF CHANCE: This study enlarges the sample size considerably from previous efforts, finding four genome-wide significant loci, including two novel signals and a further two novel genes that are implicated by gene level enrichment analyses. The novel loci, GNRH1 and FSHR, have well-established roles in female reproduction whereas ZFPM1 and IPO8 have not previously been implicated in female fertility. We found significant genetic correlations with multiple aspects of female reproduction and body size as well as evidence for significant selection against DZ twinning during human evolution. The 26 top single nucleotide polymorphisms (SNPs) from our GWAMA in European-origin participants weakly predicted the crude twinning rates in 47 non-European populations (r = 0.23 between risk score and population prevalence, s.e. 0.11, 1-tail P = 0.058) indicating that genome-wide association studies (GWAS) are needed in African and Asian populations to explore the causes of their respectively high and low DZ twinning rates. In vivo functional tests in zebrafish for IPO8 validated its essential role in female, but not male, fertility. In most regions, risk SNPs linked to known expression quantitative trait loci (eQTLs). Top SNPs were associated with in vivo reproductive hormone levels with the top pathways including hormone ligand binding receptors and the ovulation cycle. LARGE SCALE DATA: The full DZT GWAS summary statistics will made available after publication through the GWAS catalog (https://www.ebi.ac.uk/gwas/). LIMITATIONS, REASONS FOR CAUTION: Our study only included European ancestry cohorts. Inclusion of data from Africa (with the highest twining rate) and Asia (with the lowest rate) would illuminate further the biology of twinning and female fertility. WIDER IMPLICATIONS OF THE FINDINGS: About one in 40 babies born in the world is a twin and there is much speculation on why twinning runs in families. We hope our results will inform investigations of ovarian response in new and existing ARTs and the causes of female infertility. STUDY FUNDING/COMPETING INTEREST(S): Support for the Netherlands Twin Register came from the Netherlands Organization for Scientific Research (NWO) and The Netherlands Organization for Health Research and Development (ZonMW) grants, 904-61-193, 480-04-004, 400-05-717, Addiction-31160008, 911-09-032, Biobanking and Biomolecular Resources Research Infrastructure (BBMRI.NL, 184.021.007), Royal Netherlands Academy of Science Professor Award (PAH/6635) to DIB, European Research Council (ERC-230374), Rutgers University Cell and DNA Repository (NIMH U24 MH068457-06), the Avera Institute, Sioux Falls, South Dakota (USA) and the National Institutes of Health (NIH R01 HD042157-01A1) and the Genetic Association Information Network (GAIN) of the Foundation for the National Institutes of Health and Grand Opportunity grants 1RC2 MH089951. The QIMR Berghofer Medical Research Institute (QIMR) study was supported by grants from the National Health and Medical Research Council (NHMRC) of Australia (241944, 339462, 389927, 389875, 389891, 389892, 389938, 443036, 442915, 442981, 496610, 496739, 552485, 552498, 1050208, 1075175). L.Y. is funded by Australian Research Council (Grant number DE200100425). The Minnesota Center for Twin and Family Research (MCTFR) was supported in part by USPHS Grants from the National Institute on Alcohol Abuse and Alcoholism (AA09367 and AA11886) and the National Institute on Drug Abuse (DA05147, DA13240, and DA024417). The Women's Genome Health Study (WGHS) was funded by the National Heart, Lung, and Blood Institute (HL043851 and HL080467) and the National Cancer Institute (CA047988 and UM1CA182913), with support for genotyping provided by Amgen. Data collection in the Finnish Twin Registry has been supported by the Wellcome Trust Sanger Institute, the Broad Institute, ENGAGE-European Network for Genetic and Genomic Epidemiology, FP7-HEALTH-F4-2007, grant agreement number 201413, National Institute of Alcohol Abuse and Alcoholism (grants AA-12502, AA-00145, AA-09203, AA15416, and K02AA018755) and the Academy of Finland (grants 100499, 205585, 118555, 141054, 264146, 308248, 312073 and 336823 to J. Kaprio). TwinsUK is funded by the Wellcome Trust, Medical Research Council, Versus Arthritis, European Union Horizon 2020, Chronic Disease Research Foundation (CDRF), Zoe Ltd and the National Institute for Health Research (NIHR) Clinical Research Network (CRN) and Biomedical Research Centre based at Guy's and St Thomas' NHS Foundation Trust in partnership with King's College London. For NESDA, funding was obtained from the Netherlands Organization for Scientific Research (Geestkracht program grant 10000-1002), the Center for Medical Systems Biology (CSMB, NVVO Genomics), Biobanking and Biomolecular Resources Research Infrastructure (BBMRI-NL), VU University's Institutes for Health and Care Research (EMGO+) and Neuroscience Campus Amsterdam, University Medical Center Groningen, Leiden University Medical Center, National Institutes of Health (NIH, ROI D0042157-01A, MH081802, Grand Opportunity grants 1 RC2 Ml-1089951 and IRC2 MH089995). Part of the genotyping and analyses were funded by the Genetic Association Information Network (GAIN) of the Foundation for the National Institutes of Health. Computing was supported by BiG Grid, the Dutch e-Science Grid, which is financially supported by NWO. Work in the Del Bene lab was supported by the Programme Investissements d'Avenir IHU FOReSIGHT (ANR-18-IAHU-01). C.R. was supported by an EU Horizon 2020 Marie Sklodowska-Curie Action fellowship (H2020-MSCA-IF-2014 #661527). H.S. and K.S. are employees of deCODE Genetics/Amgen. The other authors declare no competing financial interests. TRIAL REGISTRATION NUMBER: N/A.

Assessing the genetic risk of nodular melanoma using a candidate gene approach.

BACKGROUND: Nodular melanoma (NM) is a challenge to diagnose early due to its rapid growth and more atypical clinical presentation, making it the largest contributor to melanoma mortality. OBJECTIVES: Our study aim was to perform a rare-variant allele (RVA) analysis of whole-exome sequencing of patients with NM and non-NM (minor allele frequency ≤ 1% non-Finnish European) for a set of 500 candidate genes potentially implicated in melanoma. METHODS: This study recruited 131 participants with NM and 194 with non-NM from South-east Queensland and patients with NM from Victoria to perform a comparative analysis of possible genetic differences or similarities between the two melanoma cohorts. RESULTS: Phenotypic analysis revealed that a majority of patients diagnosed with NM were older males with a higher frequency of fair skin and red hair than is seen in the general population. The distribution of common melanoma polygenic risk scores was similar in patients with NM and non-NM, with over 28% in the highest quantile of scores. There was also a similar frequency of carriage of familial/high-penetrant melanoma gene and loss-of-function variants. We identified 39 genes by filtering 500 candidate genes based on the greatest frequency in NM compared with non-NM cases. The genes with RVAs of greatest frequency in NM included PTCH1, ARID2 and GHR. Rare variants in the SMO gene, which interacts with PTCH1 as ligand and receptor, were also identified, providing evidence that the Hedgehog pathway may contribute to NM risk. There was a cumulative effect in carrying multiple rare variants in the NM-associated genes. A 14.8-fold increased ratio for NM compared with non-NM was seen when two RVAs of the 39 genes were carried by a patient. CONCLUSIONS: This study highlights the importance of considering frequency of RVA to identify those at risk of NM in addition to known high penetrance genes.

Germline variants are associated with increased primary melanoma tumor thickness at diagnosis.

Germline genetic variants have been identified, which predispose individuals and families to develop melanoma. Tumor thickness is the strongest predictor of outcome for clinically localized primary melanoma patients. We sought to determine whether there is a heritable genetic contribution to variation in tumor thickness. If confirmed, this will justify the search for specific genetic variants influencing tumor thickness. To address this, we estimated the proportion of variation in tumor thickness attributable to genome-wide genetic variation (variant-based heritability) using unrelated patients with measured primary cutaneous melanoma thickness. As a secondary analysis, we conducted a genome-wide association study (GWAS) of tumor thickness. The analyses utilized 10 604 individuals with primary cutaneous melanoma drawn from nine GWAS datasets from eight cohorts recruited from the general population, primary care and melanoma treatment centers. Following quality control and filtering to unrelated individuals with study phenotypes, 8125 patients were used in the primary analysis to test whether tumor thickness is heritable. An expanded set of 8505 individuals (47.6% female) were analyzed for the secondary GWAS meta-analysis. Analyses were adjusted for participant age, sex, cohort and ancestry. We found that 26.6% (SE 11.9%, P = 0.0128) of variation in tumor thickness is attributable to genome-wide genetic variation. While requiring replication, a chromosome 11 locus was associated (P < 5 × 10-8) with tumor thickness. Our work indicates that sufficiently large datasets will enable the discovery of genetic variants associated with greater tumor thickness, and this will lead to the identification of host biological processes influencing melanoma growth and invasion.

Improving the Accuracy, Robustness, and Dynamic Range of Digital Bead Assays.

We report methods that improve the quantification of digital bead assays (DBA)─such as the digital enzyme-linked immunosorbent assay (ELISA)─that have found widespread use for high sensitivity measurement of proteins in clinical research and diagnostics. In digital ELISA, proteins are captured on beads, labeled with enzymes, individual beads are interrogated for activity from one or more enzymes, and the average number of enzymes per bead (AEB) is determined based on Poisson statistics. The widespread use of digital ELISA has revealed limitations to the original approaches to quantification that can lead to inaccurate AEB. Here, we have addressed the inaccuracy in AEB due to deviations from Poisson distribution in a digital ELISA for Aß-40 by changing the AEB calculation from a fixed threshold between digital counting and average normalized intensity to a smooth, continuous combination of digital counting and intensity. We addressed issues with determining the average product fluorescence intensity from single enzymes on beads by allowing outlier, high intensity arrays to be removed from average intensities, and by permitting the use of a wider range of arrays. These approaches improved the accuracy of a digital ELISA for tau protein that was affected by aggregated detection antibodies. We increased the dynamic range of a digital ELISA for IL-17A from AEB ∼25 to ∼130 by combining long and short exposure images at the product emission wavelength to create virtual images. The methods reported will significantly improve the accuracy and robustness of DBA based on imaging─such as single molecule arrays (Simoa)─and flow detection.

The MC1R r allele does not increase melanoma risk in MITF E318K carriers.

BACKGROUND: Population-wide screening for melanoma is not cost-effective, but genetic characterization could facilitate risk stratification and targeted screening. Common Melanocortin-1 receptor (MC1R) red hair colour (RHC) variants and Microphthalmia-associated transcription factor (MITF) E318K separately confer moderate melanoma susceptibility, but their interactive effects are relatively unexplored. OBJECTIVES: To evaluate whether MC1R genotypes differentially affect melanoma risk in MITF E318K+ vs. E318K- individuals. MATERIALS AND METHODS: Melanoma status (affected or unaffected) and genotype data (MC1R and MITF E318K) were collated from research cohorts (five Australian and two European). In addition, RHC genotypes from E318K+ individuals with and without melanoma were extracted from databases (The Cancer Genome Atlas and Medical Genome Research Bank, respectively). χ2 and logistic regression were used to evaluate RHC allele and genotype frequencies within E318K+/- cohorts depending on melanoma status. Replication analysis was conducted on 200 000 general-population exomes (UK Biobank). RESULTS: The cohort comprised 1165 MITF E318K- and 322 E318K+ individuals. In E318K- cases MC1R R and r alleles increased melanoma risk relative to wild type (wt), P < 0.001 for both. Similarly, each MC1R RHC genotype (R/R, R/r, R/wt, r/r and r/wt) increased melanoma risk relative to wt/wt (P < 0.001 for all). In E318K+ cases, R alleles increased melanoma risk relative to the wt allele [odds ratio (OR) 2.04 (95% confidence interval 1.67-2.49); P = 0.01], while the r allele risk was comparable with the wt allele [OR 0.78 (0.54-1.14) vs. 1.00, respectively]. E318K+ cases with the r/r genotype had a lower but not significant melanoma risk relative to wt/wt [OR 0.52 (0.20-1.38)]. Within the E318K+ cohort, R genotypes (R/R, R/r and R/wt) conferred a significantly higher risk compared with non-R genotypes (r/r, r/wt and wt/wt) (P < 0.001). UK Biobank data supported our findings that r did not increase melanoma risk in E318K+ individuals. CONCLUSIONS: RHC alleles/genotypes modify melanoma risk differently in MITF E318K- and E318K+ individuals. Specifically, although all RHC alleles increase risk relative to wt in E318K- individuals, only MC1R R increases melanoma risk in E318K+ individuals. Importantly, in the E318K+ cohort the MC1R r allele risk is comparable with wt. These findings could inform counselling and management for MITF E318K+ individuals.

GWAS of Dizygotic Twinning in an Enlarged Australian Sample of Mothers of DZ Twins.

Female fertility is a complex trait with age-specific changes in spontaneous dizygotic (DZ) twinning and fertility. To elucidate factors regulating female fertility and infertility, we conducted a genome-wide association study (GWAS) on mothers of spontaneous DZ twins (MoDZT) versus controls (3273 cases, 24,009 controls). This is a follow-up study to the Australia/New Zealand (ANZ) component of that previously reported (Mbarek et al., 2016), with a sample size almost twice that of the entire discovery sample meta-analysed in the previous article (and five times the ANZ contribution to that), resulting from newly available additional genotyping and representing a significant increase in power. We compare analyses with and without male controls and show unequivocally that it is better to include male controls who have been screened for recent family history, than to use only female controls. Results from the SNP based GWAS identified four genomewide significant signals, including one novel region, ZFPM1 (Zinc Finger Protein, FOG Family Member 1), on chromosome 16. Previous signals near FSHB (Follicle Stimulating Hormone beta subunit) and SMAD3 (SMAD Family Member 3) were also replicated (Mbarek et al., 2016). We also ran the GWAS with a dominance model that identified a further locus ADRB2 on chr 5. These results have been contributed to the International Twinning Genetics Consortium for inclusion in the next GWAS meta-analysis (Mbarek et al., in press).

Multiplex melanoma families are enriched for polygenic risk.

Cancers, including cutaneous melanoma, can cluster in families. In addition to environmental etiological factors such as ultraviolet radiation, cutaneous melanoma has a strong genetic component. Genetic risks for cutaneous melanoma range from rare, high-penetrance mutations to common, low-penetrance variants. Known high-penetrance mutations account for only about half of all densely affected cutaneous melanoma families, and the causes of familial clustering in the remainder are unknown. We hypothesize that some clustering is due to the cumulative effect of a large number of variants of individually small effect. Common, low-penetrance genetic risk variants can be combined into polygenic risk scores. We used a polygenic risk score for cutaneous melanoma to compare families without known high-penetrance mutations with unrelated melanoma cases and melanoma-free controls. Family members had significantly higher mean polygenic load for cutaneous melanoma than unrelated cases or melanoma-free healthy controls (Bonferroni-corrected t-test P = 1.5 × 10-5 and 6.3 × 10-45, respectively). Whole genome sequencing of germline DNA from 51 members of 21 families with low polygenic risk for melanoma identified a CDKN2A p.G101W mutation in a single family but no other candidate high-penetrance melanoma susceptibility genes. This work provides further evidence that melanoma, like many other common complex disorders, can arise from the joint action of multiple predisposing factors, including rare high-penetrance mutations, as well as via a combination of large numbers of alleles of small effect.

The Heritability of Twinning in Seven Large Historic Pedigrees.

It is widely recognized that dizygotic twinning (DZT) runs in families, but estimates of heritability from twin and family data are remarkably scarce and vary considerably. Here, we traced seven large, sometimes historical, multigeneration pedigrees from West Africans, fin de siècle French Jews, Canadians (two pedigrees), and the French royal family, in which twin births were recorded. We estimated heritability of twinning (of all types) as zygosity information was not available, diluting the true DZT heritability by a third or so. The estimates in the range 8-20% are remarkably consistent across time (8-19 generations) and ethnicities and also consistent with twin and family estimates.

Addressing Delicate and Variable Cancer Morphology in Spectral Histopathology Using Canine Visceral Hemangiosarcoma.

Spectral histopathology has shown promise for the classification and diagnosis of tumors with defined morphology, but application in tumors with variable or diffuse morphologies is yet to be investigated. To address this gap, we evaluated the application of Fourier transform infrared (FTIR) imaging as an accessory diagnostic tool for canine hemangiosarcoma (HSA), a vascular endothelial cell cancer that is difficult to diagnose. To preserve the delicate vascular tumor tissue structure, and potential classification of single endothelial cells, paraffin removal was not performed, and a partial least square discrimination analysis (PLSDA) and Random Forest (RF) models to classify different tissue types at individual pixel level were established using a calibration set (24 FTIR images from 13 spleen specimens). Next, the prediction capability of the PLSDA model was tested with an independent test set (n = 11), resulting in 74% correct classification of different tissue types at an individual pixel level. Finally, the performance of the FTIR spectropathology and chemometric algorithm for diagnosis of HSA was established in a blinded set of tissue samples (n = 24), with sensitivity and specificity of 80 and 81%, respectively. Taken together, these results show that FTIR imaging without paraffin removal can be applied to tumors with diffuse morphology, and this technique is a promising tool to assist in canine splenic HSA differential diagnosis.

Cell-type-specific eQTL of primary melanocytes facilitates identification of melanoma susceptibility genes.

Most expression quantitative trait locus (eQTL) studies to date have been performed in heterogeneous tissues as opposed to specific cell types. To better understand the cell-type-specific regulatory landscape of human melanocytes, which give rise to melanoma but account for <5% of typical human skin biopsies, we performed an eQTL analysis in primary melanocyte cultures from 106 newborn males. We identified 597,335 cis-eQTL SNPs prior to linkage disequilibrium (LD) pruning and 4997 eGenes (FDR < 0.05). Melanocyte eQTLs differed considerably from those identified in the 44 GTEx tissue types, including skin. Over a third of melanocyte eGenes, including key genes in melanin synthesis pathways, were unique to melanocytes compared to those of GTEx skin tissues or TCGA melanomas. The melanocyte data set also identified trans-eQTLs, including those connecting a pigmentation-associated functional SNP with four genes, likely through cis-regulation of IRF4 Melanocyte eQTLs are enriched in cis-regulatory signatures found in melanocytes as well as in melanoma-associated variants identified through genome-wide association studies. Melanocyte eQTLs also colocalized with melanoma GWAS variants in five known loci. Finally, a transcriptome-wide association study using melanocyte eQTLs uncovered four novel susceptibility loci, where imputed expression levels of five genes (ZFP90, HEBP1, MSC, CBWD1, and RP11-383H13.1) were associated with melanoma at genome-wide significant P-values. Our data highlight the utility of lineage-specific eQTL resources for annotating GWAS findings, and present a robust database for genomic research of melanoma risk and melanocyte biology.

'Essential Tremor' Phenotype in FMR1 Premutation/Gray Zone Sibling Series: Exploring Possible Genetic Modifiers.

Fragile X-associated tremor/ataxia syndrome (FXTAS) occurs in carriers of fragile X mental retardation 1 (FMR1) X-linked small CGG expansion (gray zone [GZ] and premutation [PM]) alleles, containing 41-200 repeats. Major features comprise kinetic tremor, gait ataxia, cognitive decline and cerebellar peduncular white matter lesions, but atypical/incomplete FXTAS may occur. We explored the possibility of polygenic effects modifying the FXTAS spectrum phenotypes. We used three motor scales and selected cognitive tests in a series of three males and three females from a single sibship carrying PM or GZ alleles (44 to 75 repeats). The molecular profiles from these siblings were determined by genomewide association study with single-nucleotide polymorphism (SNP) genotyping by Illumina Global Screening Array. Nonparametric linkage analysis was applied and Parkinson's disease (PD) polygenic risk scores (PRSs) were calculated for all the siblings, based on 107 known risk variants. All male and female siblings manifested similar kinetic tremor phenotypes. In contrast to FXTAS, they showed negligible gait ataxia, and few white matter lesions on MRI. Cognitive functioning was unaffected. Suggestive evidence of linkage to a broad region of the short arm of chromosome 10 was obtained, and median PD PRS for the sibship fell within the top 30% of a sample of over 500,000 UK and Australian controls. The genomewide study results are suggestive of modifying effects of genetic risk loci linked to PD, on the neurological phenotype of FMR1-CGG small expansion carriers, resulting in an oligosymptomatic kinetic tremor seen in FXTAS spectrum, but also consistent with essential tremor.

Toxin A-Predominant Pathogenic Clostridioides difficile: A Novel Clinical Phenotype.

BACKGROUND: Most Clostridioides difficile toxinogenic strains produce both toxins A and B (A+B+), but toxin A-negative, toxin B-positive (A-B+) variants also cause disease. We report the identification of a series of pathogenic clinical C. difficile isolates that produce high amounts of toxin A with low or nondetectable toxin B. METHODS: An ultrasensitive, quantitative immunoassay was used to measure toxins A and B in stool samples from 187 C. difficile infection (CDI) patients and 44 carriers. Isolates were cultured and assessed for in vitro toxin production and in vivo phenotypes (mouse CDI model). RESULTS: There were 7 CDI patients and 6 carriers who had stools with detectable toxin A (TcdA, range 23-17 422 pg/mL; 5.6% of samples overall) but toxin B (TcdB) below the clinical detection limit (<20 pg/mL; median TcdA:B ratio 17.93). Concentrations of toxin A far exceeded B in in vitro cultures of all 12 recovered isolates (median TcdA:B ratio 26). Of 8 toxin A>>B isolates tested in mice, 4 caused diarrhea, and 3 of those 4 caused lethal disease. Ribotyping demonstrated strain diversity. TcdA-predominant samples were also identified at 2 other centers, with similar frequencies (7.5% and 6.8%). CONCLUSIONS: We report the discovery of clinical pathogenic C. difficile strains that produce high levels of toxin A but minimal or no toxin B. This pattern of toxin production is not rare (>5% of isolates) and is consistently observed in vitro and in vivo in humans and mice. Our study highlights the significance of toxin A in human CDI pathogenesis and has important implications for CDI diagnosis, treatment, and vaccine development.

Customizable Multiplex Antibody Array Immunoassays with Attomolar Sensitivities.

We have developed a customizable contact printed multiplex immunoassay capable of simultaneously measuring up to five analytes with attomolar sensitivities. This enzyme-linked immunosorbent assay (ELISA) was based on spotting different antibodies in a circular pattern at the bottom of a microtiter plate well. Unlike traditional antibody printing for ELISA that prints a capture antibody specific to a target of interest, in this ELISA we printed unique "anchor" antibodies at the well surface, each having a high affinity for a specific peptide target. By coupling each peptide to a unique assay capture antibody, this array of anchor antibodies enabled a customizable contact printed multiplex immunoassay workflow. As a proof of concept, we developed a 5-plex assay measuring interleukin 5 (IL-5), interleukin 6 (IL-6), interleukin 10 (IL-10), interleukin 22 (IL-22), and tumor necrosis factor alpha (TNF-α). Measurements of these five analytes in serum and plasma correlated well between the method utilizing the anchor antibodies and peptides and the traditional capture antibody printing approach, with r2 values of 0.99, 0.93, 0.99, 0.96, and 0.75 for IL-5, IL-6, IL-10, IL-22, and TNFα, respectively. This approach makes customizable multiplex ultrasensitive ELISA available to laboratories without access to the precision printing instrumentation and will be useful for antibody screening, custom assay development, biomarker detection, and protein profiling for diagnostic applications.

Direct Detection of miR-122 in Hepatotoxicity Using Dynamic Chemical Labeling Overcomes Stability and isomiR Challenges.

Circulating microRNAs are biomarkers reported to be stable and translational across species. MicroRNA-122 (miR-122) is a hepatocyte-specific microRNA biomarker for drug-induced liver injury (DILI). We developed a single molecule, dynamic chemical labeling (DCL) assay to directly detect miR-122 in blood. The DCL assay specifically measured miR-122 directly from 10 µL of serum or plasma without any extraction steps, with a limit of detection of 1.32 pM that enabled the identification of DILI. Testing of 192 human serum samples showed that DCL accurately identified patients at risk of DILI after acetaminophen overdose (area under ROC curve 0.98 (95% CI; 0.96-1), P < 0.0001). The DCL assay also identified liver injury in rats and dogs. The use of specific captured beads had the additional benefit of stabilizing miR-122 after sample collection, with no signal loss after 14 days at room temperature, in contrast to PCR that showed significant loss of signal. RNA sequencing demonstrated the presence of multiple miR-122 isomiRs in the serum of patients with DILI that were at low concentration or not present in healthy individuals. Sample degradation over time produced more isomiRs, particularly rapidly with DILI. PCR was inaccurate when analyzing miR-122 isomiRs, whereas the DCL assay demonstrated accurate quantification. We conclude that the DCL assay can accurately measure miR-122 to diagnose liver injury in humans and other species and can overcome microRNA stability and isomiR challenges.

Gene Discovery Using Twins.

One of Nick's key early achievements at QIMR was to establish a twin study on melanoma risk factors. The Brisbane Twin Nevus Study (BTNS) had an initial focus on nevus (mole) count in adolescents but, reflecting Nick's broad interests, expanded in scope enormously over the decades. In the skin cancer arena, BTNS was essential to genetic discoveries in melanoma, eye color and pigmentation. Later studies amassed data on thousands of phenotypes, ranging from molecular phenotypes such as gene expression to studies where gene mapping findings in adolescents turned out to have translational potential in late-onset diseases. Nick's twin data have formed the basis for an enormous range of discoveries, with Nick and his colleagues continuing to capitalize on these data.

Excellent survival and outcomes with fixed-bearing medial UKA in young patients (≤ 60 years) at minimum 10-year follow-up.

PURPOSE: To evaluate whether long-term (10-year minimum) patient outcomes and survival of fixed-bearing medial unicompartmental knee arthroplasty (UKA) in patients aged ≤ 60 years were favorable despite non-conventional age criteria. METHODS: The authors reviewed the records of 91 consecutive medial UKAs performed in patients aged ≤ 60 by a single surgeon. All patients received the same fixed-bearing M/G Unicompartmental Knee System. Patients records were updated, noting complications or revisions, and Oxford Knee Scores and overall satisfaction collected. If deceased, the general practitioner or next of kin provided data. RESULTS: Of the initial 91 knees, 10 were revised, 6 were deceased, and 1 was lost to follow-up. The final cohort of 74 knees was aged 54.3 ± 4.3 years (range 41.8-60.6) at index surgery. Using revision of any component as endpoint, the present series had a KM survival of 92.9% (CI 84.8-96.7%) at 10 years, and 87.8% (CI 78.4-93.2%) at 15 years, and a single non-fatal DVT was reported. At final follow-up of 15 ± 1.3 years (range 11-18), OKS (available for all 74 knees) was 38.4 ± 8.4 (range 18-48). Overall patients were pleased or very pleased with 72 of the knees (97%). CONCLUSION: Fixed-bearing medial UKA yields favorable results in the treatment of single compartment osteoarthritis of the knee in patients ≤ 60 years. The present study demonstrates low complication rates, good-to-excellent long-term patient outcomes, and satisfactory implant survival for this age group considering the advantages of UKA. LEVEL OF EVIDENCE: Level IV, retrospective cohort study.

Variation at DENND1B and Asthma on the Island of Tristan da Cunha.

A high prevalence of asthma has been documented among the inhabitants of Tristan da Cunha, an isolated island in the South Atlantic. The population derives from just 28 founders. We performed lung function testing, including methacholine inhalation challenge, allergen skin prick testing, and collected DNA from essentially all of the current island population (269 individuals), and genotyped a panel of 43 single-nucleotide polymorphisms (SNPs) reported as associated with asthma and atopy. We carried out a mixed-model association analysis using the known pedigree. There were 96 individuals diagnosed as asthmatic (36%), and heritability estimates were similar to those from nonisolated population samples (multifactorial threshold model, h2 = 48%). The first component from a genetic principal components analysis using the entire SNP panel was nonlinearly associated with asthma, with the maximum risk to those intermediate to reference (Human Genome Diversity Project) European and African samples means. The single most strongly associated SNP was rs2786098 (p = 5.5 × 10-5), known to regulate the gene DENND1B. This explained approximately one-third of the trait heritability, with an allelic odds ratio for the A allele of 2.6. Among A/A carriers, 10 out of 12 individuals were asthmatic. The rs2786098*A variant was initially reported to decrease the risk of childhood (atopic) asthma in European but slightly increase the risk in African-descended populations, and does significantly alter Th2 cell function. Despite an absence of overall association with this variant in recent asthma genome wide association studies meta-analyses, an effect may exist on the particular genetic background of the Tristan da Cunha population.

Noncoding Variations in the Gene Encoding Ceramide Synthase 6 are Associated with Type 2 Diabetes in a Large Indigenous Australian Pedigree.

Type 2 diabetes (T2D) is a chronic disease that disproportionately affects Indigenous Australians. We have previously reported the localization of a novel T2D locus by linkage analysis to chromosome 2q24 in a large admixed Indigenous Australian pedigree (Busfield et al. (2002). American Journal of Human Genetics, 70, 349-357). Here we describe fine mapping of this region in this pedigree, with the identification of SNPs showing strong association with T2D: rs3845724 (diabetes p = 7 × 10-4), rs4668106 (diabetes p = 9 × 10-4) and rs529002 (plasma glucose p = 3 × 10-4). These associations were successfully replicated in an independent collection of Indigenous Australian T2D cases and controls. These SNPs all lie within the gene encoding ceramide synthase 6 (CERS6) and thus may regulate ceramide synthesis.

Problems, challenges and promises: perspectives on precision medicine.

The 'precision medicine (systems medicine)' concept promises to achieve a shift to future healthcare systems with a more proactive and predictive approach to medicine, where the emphasis is on disease prevention rather than the treatment of symptoms. The individualization of treatment for each patient will be at the centre of this approach, with all of a patient's medical data being computationally integrated and accessible. Precision medicine is being rapidly embraced by biomedical researchers, pioneering clinicians and scientific funding programmes in both the European Union (EU) and USA. Precision medicine is a key component of both Horizon 2020 (the EU Framework Programme for Research and Innovation) and the White House's Precision Medicine Initiative. Precision medicine promises to revolutionize patient care and treatment decisions. However, the participants in precision medicine are faced with a considerable central challenge. Greater volumes of data from a wider variety of sources are being generated and analysed than ever before; yet, this heterogeneous information must be integrated and incorporated into personalized predictive models, the output of which must be intelligible to non-computationally trained clinicians. Drawing primarily from the field of 'oncology', this article will introduce key concepts and challenges of precision medicine and some of the approaches currently being implemented to overcome these challenges. Finally, this article also covers the criticisms of precision medicine overpromising on its potential to transform patient care.