Capturing open ocean biodiversity: Comparing environmental DNA metabarcoding to the continuous plankton recorder.

Environmental DNA (eDNA) metabarcoding is emerging as a novel, objective tool for monitoring marine metazoan biodiversity. Zooplankton biodiversity in the vast open ocean is currently monitored through continuous plankton recorder (CPR) surveys, using ship-based bulk plankton sampling and morphological identification. We assessed whether eDNA metabarcoding (2 L filtered seawater) could capture similar Southern Ocean zooplankton biodiversity as conventional CPR bulk sampling (~1,500 L filtered seawater per CPR sample). We directly compared eDNA metabarcoding with (a) conventional morphological CPR sampling and (b) bulk DNA metabarcoding of CPR collected plankton (two transects for each comparison, 40 and 44 paired samples, respectively). A metazoan-targeted cytochrome c oxidase I (COI) marker was used to characterize species-level diversity. In the 2 L seawater eDNA samples, this marker amplified large amounts of non-metazoan picoplanktonic algae, but eDNA metabarcoding still detected up to 1.6 times more zooplankton species than morphologically analysed bulk CPR samples. COI metabarcoding of bulk DNA samples mostly avoided nonmetazoan amplifications and recovered more zooplankton species than eDNA metabarcoding. However, eDNA metabarcoding detected roughly two thirds of metazoan species and identified similar taxa contributing to community differentiation across the subtropical front separating transects. We observed a diurnal pattern in eDNA data for copepods which perform diel vertical migrations, indicating a surprisingly short temporal eDNA signal. Compared to COI, a eukaryote-targeted 18S ribosomal RNA marker detected a higher proportion, but lower diversity, of metazoans in eDNA. With refinement and standardization of methodology, eDNA metabarcoding could become an efficient tool for monitoring open ocean biodiversity.

Phylogenomic Resolution of the Cetacean Tree of Life Using Target Sequence Capture.

The evolution of cetaceans, from their early transition to an aquatic lifestyle to their subsequent diversification, has been the subject of numerous studies. However, although the higher-level relationships among cetacean families have been largely settled, several aspects of the systematics within these groups remain unresolved. Problematic clades include the oceanic dolphins (37 spp.), which have experienced a recent rapid radiation, and the beaked whales (22 spp.), which have not been investigated in detail using nuclear loci. The combined application of high-throughput sequencing with techniques that target specific genomic sequences provide a powerful means of rapidly generating large volumes of orthologous sequence data for use in phylogenomic studies. To elucidate the phylogenetic relationships within the Cetacea, we combined sequence capture with Illumina sequencing to generate data for $\sim $3200 protein-coding genes for 68 cetacean species and their close relatives including the pygmy hippopotamus. By combining data from $>$38,000 exons with existing sequences from 11 cetaceans and seven outgroup taxa, we produced the first comprehensive comparative genomic data set for cetaceans, spanning 6,527,596 aligned base pairs (bp) and 89 taxa. Phylogenetic trees reconstructed with maximum likelihood and Bayesian inference of concatenated loci, as well as with coalescence analyses of individual gene trees, produced mostly concordant and well-supported trees. Our results completely resolve the relationships among beaked whales as well as the contentious relationships among oceanic dolphins, especially the problematic subfamily Delphinidae. We carried out Bayesian estimation of species divergence times using MCMCTree and compared our complete data set to a subset of clocklike genes. Analyses using the complete data set consistently showed less variance in divergence times than the reduced data set. In addition, integration of new fossils (e.g., Mystacodon selenensis) indicates that the diversification of Crown Cetacea began before the Late Eocene and the divergence of Crown Delphinidae as early as the Middle Miocene. [Cetaceans; phylogenomics; Delphinidae; Ziphiidae; dolphins; whales.].

More than the eye can see: Genomic insights into the drivers of genetic differentiation in Royal/Macaroni penguins across the Southern Ocean.

The study of systematics in wide-ranging seabirds can be challenging due to the vast geographic scales involved, as well as the possible discordance between molecular, morphological and behavioral data. In the Southern Ocean, macaroni penguins (Eudyptes chrysolophus) are distributed over a circumpolar range including populations in Antarctic and sub-Antarctic areas. Macquarie Island, in its relative isolation, is home to a closely related endemic taxon - the royal penguin (Eudyptes schlegeli), which is distinguishable from E. chrysolophus mainly by facial coloration. Although these sister taxa are widely accepted as representing distinct species based on morphological grounds, the extent of their genome-wide differentiation remains uncertain. In this study, we use genome-wide Single Nucleotide Polymorphisms to test genetic differentiation between these geographically isolated taxa and evaluate the main drivers of population structure among breeding colonies of macaroni/royal penguins. Genetic similarity observed between macaroni and royal penguins suggests they constitute a single evolutionary unit. Nevertheless, royal penguins exhibited a tendency to cluster only with macaroni individuals from Kerguelen Island, suggesting that dispersal occurs mainly between these neighboring colonies. A stepping stone model of differentiation of macaroni/royal populations was further supported by a strong pattern of isolation by distance detected across its whole distribution range, possibly driven by large geographic distances between colonies as well as natal philopatry. However, we also detected intraspecific genomic differentiation between Antarctic and sub-Antarctic populations of macaroni penguins, highlighting the role of environmental factors together with geographic distance in the processes of genetic differentiation between Antarctic and sub-Antarctic waters.

Regulation of the ITGA2 gene by epigenetic mechanisms in prostate cancer.

BACKGROUND: Integrin alpha2 beta1 (α2 ß1 ) plays an integral role in tumour cell invasion, metastasis and angiogenesis, and altered expression of the receptor has been linked to tumour prognosis in several solid tumours. However, the relationship is complex, with both increased and decreased expression associated with different stages of tumour metastases in several tumour types. The ITGA2 gene, which codes for the α2 subunit, was examined to investigate whether a large CpG island associated with its promoter region is involved in the differential expression of ITGA2 observed in prostate cancer. METHODS: Bisulphite sequencing of the ITGA2 promoter was used to assess methylation in formalin-fixed paraffin-embedded (FFPE) prostate tumour specimens and prostate cancer cell lines, PC3, 22Rv1 and LNCaP. Changes in ITGA2 mRNA expression were measured using quantitative PCR. ITGA2 functionality was interrogated using cell migration scratch assays and siRNA knockdown experiments. RESULTS: Bisulphite sequencing revealed strikingly decreased methylation at key CpG sites within the promoter of tumour samples, when compared with normal prostate tissue. Altered methylation of this CpG island is also associated with differences in expression in the non-invasive LNCaP, and the highly metastatic PC3 and 22Rv1 prostate cancer cell lines. Further bisulphite sequencing confirmed that selected CpGs were highly methylated in LNCaP cells, whilst only low levels of methylation were observed in PC3 and 22Rv1 cells, correlating with ITGA2 transcript levels. Examination of the increased expression of ITGA2 was shown to influence migratory potential via scratch assay in PC3, 22Rv1 and LNCaP cells, and was confirmed by siRNA knockdown experiments. CONCLUSIONS: Taken together, our data supports the assertion that epigenetic modification of the ITGA2 promoter is a mechanism by which ITGA2 expression is regulated.

Molecular biomarkers for chronological age in animal ecology.

The chronological age of an individual animal predicts many of its biological characteristics, and these in turn influence population-level ecological processes. Animal age information can therefore be valuable in ecological research, but many species have no external features that allow age to be reliably determined. Molecular age biomarkers provide a potential solution to this problem. Research in this area of molecular ecology has so far focused on a limited range of age biomarkers. The most commonly tested molecular age biomarker is change in average telomere length, which predicts age well in a small number of species and tissues, but performs poorly in many other situations. Epigenetic regulation of gene expression has recently been shown to cause age-related modifications to DNA and to cause changes in abundance of several RNA types throughout animal lifespans. Age biomarkers based on these epigenetic changes, and other new DNA-based assays, have already been applied to model organisms, humans and a limited number of wild animals. There is clear potential to apply these marker types more widely in ecological studies. For many species, these new approaches will produce age estimates where this was previously impractical. They will also enable age information to be gathered in cross-sectional studies and expand the range of demographic characteristics that can be quantified with molecular methods. We describe the range of molecular age biomarkers that have been investigated to date and suggest approaches for developing the newer marker types as age assays in nonmodel animal species.

Australia's east coast humpback whales: Satellite tag-derived movements on breeding grounds, feeding grounds and along the northern and southern migration.

Background: Satellite tags were deployed on 50 east Australian humpback whales (breeding stock E1) between 2008 and 2010 on their southward migration, northward migration and feeding grounds in order to identify and describe migratory pathways, feeding grounds and possible calving areas. At the time, these movements were not well understood and calving grounds were not clearly identified. To the best of our knowledge, this dataset details all long-term, implantable tag deployments that have occurred to date on breeding stock E1. As such, these data provide researchers, regulators and industry with clear and valuable insights into the spatial and temporal nature of humpback whale movements along the eastern coastline of Australia and into the Southern Ocean. As this population of humpback whales navigates an increasingly complex habitat undergoing various development pressures and anthropogenic disturbances, in addition to climate-mediated changes in their marine environment, this dataset may also provide a valuable baseline. New information: At the time these tracks were generated, these were the first satellite tag deployments intended to deliver long-term, detailed movement information on east Australian (breeding stock E1) humpback whales. The tracking data revealed previously unknown migratory pathways into the Southern Ocean, with 11 individuals tracked to their Antarctic feeding grounds. Once assumed to head directly south on their southern migration, five individuals initially travelled west towards New Zealand. Six tracks detailed the coastal movement of humpback whales migrating south. One tag transmitted a partial southern migration, then ceased transmissions only to begin transmitting eight months later as the animal was migrating north. Northern migration to breeding grounds was detailed for 13 individuals, with four tracks including turning points and partial southern migrations. Another 14 humpback whales were tagged in Antarctica, providing detailed Antarctic feeding ground movements.Broadly speaking, the tracking data revealed a pattern of movement where whales were at their northern limit in July and their southern limit in March. Migration north was most rapid across the months of May and June, whilst migration south was most rapid between November and December. Tagged humpback whales were located on their Antarctic feeding grounds predominantly between January and May and approached their breeding grounds between July and August. Tracking distances ranged from 68 km to 8580 km and 1 to 286 days. To the best of our knowledge, this dataset compiles all of the long-term tag deployments that have occurred to date on breeding stock E1.

Variation in the tyrosinase gene associated with a white humpback whale (Megaptera novaeangliae).

Tyrosinase-negative oculocutaneous albinism (OCA1A) is characterized by lifelong white hair and skin, a phenotype that has been described in most mammalian species worldwide. Tyrosinase is the key enzyme in melanin biosynthesis, and mutations in the tyrosinase gene result in OCA1A. We examined sequence variation at exon 1 of the tyrosinase gene in 66 humpback whale samples collected from the east coast of Australia, including an anomalously white humpback whale known as "Migaloo." We identified 3 novel variants, including a cytosine deletion that results in a premature stop codon in exon 1. The deletion truncates the tyrosinase protein including the putative catalytic domains that are essential for tyrosinase enzymatic activity. Migaloo was homozygous for this deletion, suggesting that the albino phenotype is a consequence of inactive tyrosinase caused by the frameshift in the tyrosinase gene.

Environmental DNA metabarcoding for monitoring metazoan biodiversity in Antarctic nearshore ecosystems.

Antarctic benthic ecosystems support high biodiversity but their characterization is limited to a few well-studied areas, due to the extreme environment and remoteness making access and sampling difficult. Our aim was to compare water and sediment as sources of environmental DNA (eDNA) to better characterise Antarctic benthic communities and further develop practical approaches for DNA-based biodiversity assessment in remote environments. We used a cytochrome c oxidase subunit I (COI) metabarcoding approach to characterise metazoan communities in 26 nearshore sites across 12 locations in the Vestfold Hills (East Antarctica) based on DNA extracted from either sediment cores or filtered seawater. We detected a total of 99 metazoan species from 12 phyla across 26 sites, with similar numbers of species detected in sediment and water eDNA samples. However, significantly different communities were detected in the two sample types at sites where both were collected (i.e., where paired samples were available). For example, nematodes and echinoderms were more likely to be detected exclusively in sediment and water eDNA samples, respectively. eDNA from water and sediment core samples are complementary sample types, with epifauna more likely to be detected in water column samples and infauna in sediment. More reference DNA sequences are needed for infauna/meiofauna to increase the proportion of sequences and number of taxa that can be identified. Developing a better understanding of the temporal and spatial dynamics of eDNA at low temperatures would also aid interpretation of eDNA signals from polar environments. Our results provide a preliminary scan of benthic metazoan communities in the Vestfold Hills, with additional markers required to provide a comprehensive biodiversity survey. However, our study demonstrates the choice of sample type for eDNA studies of benthic ecosystems (sediment, water or both) needs to be carefully considered in light of the research or monitoring question of interest.

Age estimation in a long-lived seabird (Ardenna tenuirostris) using DNA methylation-based biomarkers.

Age structure is a fundamental aspect of animal population biology. Age is strongly related to individual physiological condition, reproductive potential and mortality rate. Currently, there are no robust molecular methods for age estimation in birds. Instead, individuals must be ringed as chicks to establish known-age populations, which is a labour-intensive and expensive process. The estimation of chronological age using DNA methylation (DNAm) is emerging as a robust approach in mammals including humans, mice and some non-model species. Here, we quantified DNAm in whole blood samples from a total of 71 known-age Short-tailed shearwaters (Ardenna tenuirostris) using digital restriction enzyme analysis of methylation (DREAM). The DREAM method measures DNAm levels at thousands of CpG dinucleotides throughout the genome. We identified seven CpG sites with DNAm levels that correlated with age. A model based on these relationships estimated age with a mean difference of 2.8 years to known age, based on validation estimates from models created by repeated sampling of training and validation data subsets. Longitudinal observation of individuals re-sampled over 1 or 2 years generally showed an increase in estimated age (6/7 cases). For the first time, we have shown that epigenetic changes with age can be detected in a wild bird. This approach should be of broad interest to researchers studying age biomarkers in non-model species and will allow identification of markers that can be assessed using targeted techniques for accurate age estimation in large population studies.

Sequence variants of alpha-methylacyl-CoA racemase are associated with prostate cancer risk: a replication study in an ethnically homogeneous population.

BACKGROUND: Examination of variants of the alpha-methylacyl-CoA racemase (AMACR) gene, as genetic contributors to prostate cancer risk, has been of considerable interest given the gene's recently established role as a diagnostic biomarker for prostate cancer. METHODS: The AMACR gene variants, M9V and D175G, were genotyped in a familial dataset comprising 127 cases and in a second sporadic prostate cancer dataset comprising 414 cases and 319 controls. Genotype-disease associations were examined employing the M(QLS) test and unconditional logistic regression. Differences in allele frequencies were examined using the Fisher's exact test. Association between the AMACR haplotypes and prostate cancer risk was also investigated using haplo.score. RESULTS: Significant evidence for association with prostate cancer risk for both the M9V and D175G variants was observed in the Tasmanian prostate cancer dataset. Whilst this association remained significant, it was diminished when relatedness amongst the familial prostate cancer cases was considered. CONCLUSION: This study, performed in a relatively genetically homogenous Tasmanian population, provides further evidence for a significant association between variants within the AMACR gene and prostate cancer risk. Risk was found to be more significantly associated with AMACR gene variants in sporadic compared to familial prostate cancer cases. These findings again highlight that genetic heterogeneity in the study population should be considered when examining genetic risk factors in prostate cancer.

Genetic monitoring of open ocean biodiversity: An evaluation of DNA metabarcoding for processing continuous plankton recorder samples.

DNA metabarcoding is an efficient method for measuring biodiversity, but the process of initiating long-term DNA-based monitoring programmes, or integrating with conventional programs, is only starting. In marine ecosystems, plankton surveys using the continuous plankton recorder (CPR) have characterized biodiversity along transects covering millions of kilometres with time-series spanning decades. We investigated the potential for use of metabarcoding in CPR surveys. Samples (n = 53) were collected in two Southern Ocean transects and metazoans identified using standard microscopic methods and by high-throughput sequencing of a cytochrome c oxidase subunit I marker. DNA increased the number of metazoan species identified and provided high-resolution taxonomy of groups problematic in conventional surveys (e.g., larval echinoderms and hydrozoans). Metabarcoding also generally produced more detections than microscopy, but this sensitivity may make cross-contamination during sampling a problem. In some samples, the prevalence of DNA from large plankton such as krill masked the presence of smaller species. We investigated adding a fixed amount of exogenous DNA to samples as an internal control to allow determination of relative plankton biomass. Overall, the metabarcoding data represent a substantial shift in perspective, making direct integration into current long-term time-series challenging. We discuss a number of hurdles that exist for progressing DNA metabarcoding from the current snapshot studies to the requirements of a long-term monitoring programme. Given the power and continually increasing efficiency of metabarcoding, it is almost certain this approach will play an important role in future plankton monitoring.

DNA methylation levels in candidate genes associated with chronological age in mammals are not conserved in a long-lived seabird.

Most seabirds do not have any outward identifiers of their chronological age, so estimation of seabird population age structure generally requires expensive, long-term banding studies. We investigated the potential to use a molecular age biomarker to estimate age in short-tailed shearwaters (Ardenna tenuirostris). We quantified DNA methylation in several A. tenuirostris genes that have shown age-related methylation changes in mammals. In birds ranging from chicks to 21 years of age, bisulphite treated blood and feather DNA was sequenced and methylation levels analysed in 67 CpG sites in 13 target gene regions. From blood samples, five of the top relationships with age were identified in KCNC3 loci (CpG66: R2 = 0.325, p = 0.019). In feather samples ELOVL2 (CpG42: R2 = 0.285, p = 0.00048) and EDARADD (CpG46: R2 = 0.168, p = 0.0067) were also weakly correlated with age. However, the majority of markers had no clear association with age (of 131 comparisons only 12 had a p-value < 0.05) and statistical analysis using a penalised lasso approach did not produce an accurate ageing model. Our data indicate that some age-related signatures identified in orthologous mammalian genes are not conserved in the long-lived short tailed shearwater. Alternative molecular approaches will be required to identify a reliable biomarker of chronological age in these seabirds.

A phylogenomic analysis of the role and timing of molecular adaptation in the aquatic transition of cetartiodactyl mammals.

Recent studies have reported multiple cases of molecular adaptation in cetaceans related to their aquatic abilities. However, none of these has included the hippopotamus, precluding an understanding of whether molecular adaptations in cetaceans occurred before or after they split from their semi-aquatic sister taxa. Here, we obtained new transcriptomes from the hippopotamus and humpback whale, and analysed these together with available data from eight other cetaceans. We identified more than 11 000 orthologous genes and compiled a genome-wide dataset of 6845 coding DNA sequences among 23 mammals, to our knowledge the largest phylogenomic dataset to date for cetaceans. We found positive selection in nine genes on the branch leading to the common ancestor of hippopotamus and whales, and 461 genes in cetaceans compared to 64 in hippopotamus. Functional annotation revealed adaptations in diverse processes, including lipid metabolism, hypoxia, muscle and brain function. By combining these findings with data on protein-protein interactions, we found evidence suggesting clustering among gene products relating to nervous and muscular systems in cetaceans. We found little support for shared ancestral adaptations in the two taxa; most molecular adaptations in extant cetaceans occurred after their split with hippopotamids.

Epigenetic estimation of age in humpback whales.

Age is a fundamental aspect of animal ecology, but is difficult to determine in many species. Humpback whales exemplify this as they have a lifespan comparable to humans, mature sexually as early as 4 years and have no reliable visual age indicators after their first year. Current methods for estimating humpback age cannot be applied to all individuals and populations. Assays for human age have recently been developed based on age-induced changes in DNA methylation of specific genes. We used information on age-associated DNA methylation in human and mouse genes to identify homologous gene regions in humpbacks. Humpback skin samples were obtained from individuals with a known year of birth and employed to calibrate relationships between cytosine methylation and age. Seven of 37 cytosines assayed for methylation level in humpback skin had significant age-related profiles. The three most age-informative cytosine markers were selected for a humpback epigenetic age assay. The assay has an R(2) of 0.787 (P = 3.04e-16) and predicts age from skin samples with a standard deviation of 2.991 years. The epigenetic method correctly determined which of parent-offspring pairs is the parent in more than 93% of cases. To demonstrate the potential of this technique, we constructed the first modern age profile of humpback whales off eastern Australia and compared the results to population structure 5 decades earlier. This is the first epigenetic age estimation method for a wild animal species and the approach we took for developing it can be applied to many other nonmodel organisms.

Adélie penguin population diet monitoring by analysis of food DNA in scats.

The Adélie penguin is the most important animal currently used for ecosystem monitoring in the Southern Ocean. The diet of this species is generally studied by visual analysis of stomach contents; or ratios of isotopes of carbon and nitrogen incorporated into the penguin from its food. There are significant limitations to the information that can be gained from these methods. We evaluated population diet assessment by analysis of food DNA in scats as an alternative method for ecosystem monitoring with Adélie penguins as an indicator species. Scats were collected at four locations, three phases of the breeding cycle, and in four different years. A novel molecular diet assay and bioinformatics pipeline based on nuclear small subunit ribosomal RNA gene (SSU rDNA) sequencing was used to identify prey DNA in 389 scats. Analysis of the twelve population sample sets identified spatial and temporal dietary change in Adélie penguin population diet. Prey diversity was found to be greater than previously thought. Krill, fish, copepods and amphipods were the most important food groups, in general agreement with other Adélie penguin dietary studies based on hard part or stable isotope analysis. However, our DNA analysis estimated that a substantial portion of the diet was gelatinous groups such as jellyfish and comb jellies. A range of other prey not previously identified in the diet of this species were also discovered. The diverse prey identified by this DNA-based scat analysis confirms that the generalist feeding of Adélie penguins makes them a useful indicator species for prey community composition in the coastal zone of the Southern Ocean. Scat collection is a simple and non-invasive field sampling method that allows DNA-based estimation of prey community differences at many temporal and spatial scales and provides significant advantages over alternative diet analysis approaches.

The effect of input DNA copy number on genotype call and characterising SNP markers in the humpback whale genome using a nanofluidic array.

Recent advances in nanofluidic technologies have enabled the use of Integrated Fluidic Circuits (IFCs) for high-throughput Single Nucleotide Polymorphism (SNP) genotyping (GT). In this study, we implemented and validated a relatively low cost nanofluidic system for SNP-GT with and without Specific Target Amplification (STA). As proof of principle, we first validated the effect of input DNA copy number on genotype call rate using well characterised, digital PCR (dPCR) quantified human genomic DNA samples and then implemented the validated method to genotype 45 SNPs in the humpback whale, Megaptera novaeangliae, nuclear genome. When STA was not incorporated, for a homozygous human DNA sample, reaction chambers containing, on average 9 to 97 copies, showed 100% call rate and accuracy. Below 9 copies, the call rate decreased, and at one copy it was 40%. For a heterozygous human DNA sample, the call rate decreased from 100% to 21% when predicted copies per reaction chamber decreased from 38 copies to one copy. The tightness of genotype clusters on a scatter plot also decreased. In contrast, when the same samples were subjected to STA prior to genotyping a call rate and a call accuracy of 100% were achieved. Our results demonstrate that low input DNA copy number affects the quality of data generated, in particular for a heterozygous sample. Similar to human genomic DNA, a call rate and a call accuracy of 100% was achieved with whale genomic DNA samples following multiplex STA using either 15 or 45 SNP-GT assays. These calls were 100% concordant with their true genotypes determined by an independent method, suggesting that the nanofluidic system is a reliable platform for executing call rates with high accuracy and concordance in genomic sequences derived from biological tissue.

Identification of a prostate cancer susceptibility gene on chromosome 5p13q12 associated with risk of both familial and sporadic disease.

Genetic heterogeneity is a difficulty frequently encountered in the search for genes conferring susceptibility to prostate cancer. To circumvent this issue, we selected a large prostate cancer pedigree for genome-wide linkage analysis from a population that is genetically homogeneous. Selected cases and first-degree relatives were genotyped with Affymetrix 10K SNP arrays, identifying a 14 Mb haplotype on chromosome 5 (5p13-q12) inherited identical-by-descent (IBD) by multiple cases. Microsatellite genotyping of additional deceased case samples confirmed that a total of eight cases inherited the common haplotype (P=0.0017). Re-sequencing of eight prioritised candidate genes in the region in six selected individuals identified 15 SNPs segregating with the IBD haplotype, located within the ITGA2 gene. Three of these polymorphisms were selected for genotyping in an independent Tasmanian data set comprising 127 cases with familial prostate cancer, 412 sporadic cases and 319 unaffected controls. Two were associated with prostate cancer risk: rs3212649 (OR=1.67 (1.07-2.6), P=0.0009) and rs1126643 (OR=1.52 (1.01-2.28), P=0.0088). Significant association was observed in both familial and sporadic prostate cancer. Although the functional SNP remains to be identified, considerable circumstantial evidence, provided by in vivo and in vitro studies, supports a role for ITGA2 in tumour development.