Decrypting seasonal patterns of key pollen taxa in cool temperate Australia: A multi-barcode metabarcoding analysis

Pollen allergies pose a considerable global public health concern. Allergy risk can vary significantly within plant families, yet some key pollen allergens can only be identified to family level by current optical methods. Pollen information with greater taxonomic resolution is therefore required to best support allergy prevention and self-management. We used environmental DNA (eDNA) metabarcoding to deepen taxonomic insights into the seasonal composition of airborne pollen in cool temperate Australia, a region with high rates of allergic respiratory disease. In Hobart, Tasmania, we collected routine weekly air samples from December 2018 until October 2020 and sequenced the internal transcribed spacer 2 (ITS2) and chloroplastic tRNA-Leucine tRNA-Phenylalanine intergenic spacer (trnL-trnF) regions in order to address the following questions: a) What is the genus-level diversity of known and potential aeroallergens in Hobart, in particular, in the families Poaceae, Cupressaceae and Myrtaceae? b) How do the atmospheric concentrations of these taxa change over time, and c) Does trnL-trnF enhance resolution of biodiversity when used in addition to ITS2? Our results suggest that individuals in the region are exposed to temperate grasses including Poa and Bromus in the peak grass pollen season, however low levels of exposure to the subtropical grass Cynodon may occur in autumn and winter. Within Cupressaceae, both metabarcodes showed that exposure is predominantly to pollen from the introduced genera Cupressus and Juniperus. Only ITS2 detected the native genus, Callitris. Both metabarcodes detected Eucalyptus as the major Myrtaceae genus, with trnL-trnF exhibiting primer bias for this family. These findings help refine our understanding of allergy triggers in Tasmania and highlight the utility of multiple metabarcodes in aerobiome studies.

Identification of a novel recurrent EEF2 gene amplification in familial prostate tumors

Recurrent tumor copy number variations (CNVs) in prostate cancer (PrCa) have predominantly been discovered in sporadic tumor cohorts. Here, we examined familial prostate tumors for novel CNVs as prior studies suggest these harbor distinct CNVs. Array comparative genomic hybridization of 12 tumors from an Australian PrCa family, PcTas9, highlighted multiple recurrent CNVs, including amplification of EEF2 (19p13.3) in 100% of tumors. The EEF2 CNV was examined in a further 26 familial and seven sporadic tumors from the Australian cohort and in 494 tumors unselected for family history from The Cancer Genome Atlas (TCGA). EEF2 overexpression was observed in seven PcTas9 tumors, in addition to seven other predominantly familial tumors (n_total  = 34%). EEF2 amplification was only observed in 1.4% of TCGA tumors, however 7.5% harbored an EEF2 deletion. Analysis of genes co-expressed with EEF2 revealed significant upregulation of two genes, ZNF74 and ADSL, and downregulation of PLSCR1 in both EEF2 amplified familial tumors and EEF2 deleted TCGA tumors. Furthermore, in TCGA tumors, EEF2 amplification and deletion were significantly associated with a higher Gleason score. In summary, we identified a novel PrCa CNV that was predominantly amplified in familial tumors and deleted in unselected tumors. Our results provide further evidence that familial tumors harbor distinct CNVs, potentially due to an inherited predisposition, but also suggest that regardless of how EEF2 is dysregulated, a similar set of genes involved in key cancer pathways are impacted. Given the current lack of gene-based biomarkers and clinical targets in PrCa, further investigation of EEF2 is warranted.

From rivers to ocean basins: The role of ocean barriers and philopatry in the genetic structuring of a cosmopolitan coastal predator

The Bull Shark (Carcharhinus leucas) faces varying levels of exploitation around the world due to its coastal distribution. Information regarding population connectivity is crucial to evaluate its conservation status and local fishing impacts. In this study, we sampled 922 putative Bull Sharks from 19 locations in the first global assessment of population structure of this cosmopolitan species. Using a recently developed DNA-capture approach (DArTcap), samples were genotyped for 3400 nuclear markers. Additionally, full mitochondrial genomes of 384 Indo-Pacific samples were sequenced. Reproductive isolation was found between and across ocean basins (eastern Pacific, western Atlantic, eastern Atlantic, Indo-West Pacific) with distinct island populations in Japan and Fiji. Bull Sharks appear to maintain gene flow using shallow coastal waters as dispersal corridors, whereas large oceanic distances and historical land-bridges act as barriers. Females tend to return to the same area for reproduction, making them more susceptible to local threats and an important focus for management actions. Given these behaviors, the exploitation of Bull Sharks from insular populations, such as Japan and Fiji, may instigate local decline that cannot readily be replenished by immigration, which can in turn affect ecosystem dynamics and functions. These data also supported the development of a genetic panel to ascertain the population of origin, which will be useful in monitoring the trade of fisheries products and assessing population-level impacts of this harvest.

Rapid assessment of adult abundance and demographic connectivity from juvenile kin pairs in a critically endangered species

The viability of spatially structured populations depends on the abundance and connectivity between subpopulations of breeding adults. Yet, for many species, both are extremely difficult to assess. The speartooth shark is a critically endangered elasmobranch inhabiting tropical rivers with only three adults ever recorded in Australia. Close-kin mark-recapture models, informed by sibling pairs among 226 juveniles, were developed to estimate adult abundance and connectivity in two Australian river systems. Sixty-eight sibling pairs were found, and adult abundance was estimated at 892 for the Adelaide River and 1128 for the Alligator Rivers. We found strong evidence for female philopatry, with most females returning to the same river to pup. Adelaide River males appear largely philopatric, whereas Alligator Rivers males are highly connected to the Adelaide River. From only 4 years of sampling, our results demonstrate that juvenile-only kin pairs can inform simultaneous estimates of abundance and connectivity in a rare and threatened species.

Analysis of a large prostate cancer family identifies novel and recurrent gene fusion events providing evidence for inherited predisposition

There is strong interest in the characterisation of gene fusions and their use to enhance clinical practices in prostate cancer (PrCa). Significantly, ~50% of prostate tumours harbour a gene fusion. Inherited factors are thought to predispose to these events but, to date, only one study has investigated gene fusions in a familial context. Here, we examined the prevalence and diversity of gene fusions in 14 tumours from a single large PrCa family, PcTas9, using the TruSight® RNA Fusion Panel and Sanger sequencing validation. These fusions were then explored in The Cancer Genome Atlas (TCGA) PrCa data set (n = 494). Overall, 64.3% of PcTas9 tumours harboured a gene fusion, including known erythroblast transformation-specific (ETS) fusions involving ERG and ETV1, and two novel gene fusions, C19orf48:ETV4 and RYBP:FOXP1. Although 3' ETS genes were overexpressed in PcTas9 and TCGA tumour samples, 3' fusion of FOXP1 did not appear to alter its expression. In addition, PcTas9 fusion carriers were more likely to have lower-grade disease than noncarriers (p = 0.02). Likewise, TCGA tumours with high-grade disease were less likely to harbour fusions (p = 0.03). Our study further implicates an inherited predisposition to PrCa gene fusion events, which are associated with less aggressive tumours. This knowledge could lead to clinical strategies to predict men at risk for fusion-positive PrCa and, thus, identify patients who are more or less at risk of aggressive disease and/or responsive to particular therapies.

A rare variant in EZH2 is associated with prostate cancer risk

Prostate cancer (PrCa) is highly heritable, and although rare variants contribute significantly to PrCa risk, few have been identified to date. Herein, whole-genome sequencing was performed in a large PrCa family featuring multiple affected relatives spanning several generations. A rare, predicted splice site EZH2 variant, rs78589034 (G > A), was identified as segregating with disease in all but two individuals in the family, one of whom was affected with lymphoma and bowel cancer and a female relative. This variant was significantly associated with disease risk in combined familial and sporadic PrCa datasets (n = 1551; odds ratio [OR] = 3.55, P = 1.20 × 10^-5 ). Transcriptome analysis was performed on prostate tumour needle biopsies available for two rare variant carriers and two wild-type cases. Although no allele-dependent differences were detected in EZH2 transcripts, a distinct differential gene expression signature was observed when comparing prostate tissue from the rare variant carriers with the wild-type samples. The gene expression signature comprised known downstream targets of EZH2 and included the top-ranked genes, DUSP1, FOS, JUNB and EGR1, which were subsequently validated by qPCR. These data provide evidence that rs78589034 is associated with increased PrCa risk in Tasmanian men and further, that this variant may be associated with perturbed EZH2 function in prostate tissue. Disrupted EZH2 function is a driver of tumourigenesis in several cancers, including prostate, and is of significant interest as a therapeutic target.

Multiple endocrine neoplasia type 1: clinical correlates of MEN1 gene methylation

Multiple endocrine neoplasia type 1 (MEN 1) has marked severity variation between individuals with the same mutation. To investigate any relationship between promoter methylation and clinical features, blood and tissue samples were collected from 16 members of the Tasman 1 MEN 1 kindred carrying a common splice site mutation and 7 patients with sporadic MEN 1. Methylation at 39 CpGs in the MEN1 promoter were assessed in formalin fixed, paraffin embedded parathyroid tissue. Clinical disease severity markers included age at first parathyroid operation, parathyroid hormone level and corrected serum calcium levels. Six patients with sporadic hyperparathyroidism were used for comparison. Minimal methylation was observed in all patients across CpG sites 1-23. In contrast, hypermethylation was observed at CpG sites 24-31 in MEN 1 patients, a pattern not observed in patients with non-MEN 1 parathyroid disease. Mean methylation at sites 24-31 was significantly correlated with age at first parathyroid operation (r = 0.652, p = 0.041). A permutation test, utilising the mean correlation coefficient (r = -0.401) revealed a possible association between relative PHPT severity and methylation score for each significant CpG site (p < 0.103). This novel study reveals evidence supporting a possible association between altered MEN1 promoter methylation and clinical severity of disease.

Development of 15 nuclear EST microsatellite markers for the paleoendemic conifer Pherosphaera hookeriana (Podocarpaceae)

PREMISE OF THE STUDY

Nuclear microsatellite markers were developed for population genetic analysis of the threatened paleoendemic conifer Pherosphaera hookeriana (Podocarpaceae).

METHODS AND RESULTS

Fifteen variable loci were identified showing one to 13 alleles per population, with seven loci displaying at least four alleles in all populations, and the average number of alleles per locus ranging from 4.80 to 5.93 per population. Levels of observed heterozygosity per locus varied from 0.00 to 0.91, while average heterozygosity across all loci varied from 0.54 to 0.63 between populations. All loci also amplified in the endangered congener P. fitzgeraldii, but only five of the loci had more than one allele.

CONCLUSIONS

These 15 loci are the first microsatellite markers developed in the genus Pherosphaera. These loci will be useful for investigating the species' extant genetic diversity and structure, the impact of past environmental change, and the significance of asexual reproduction.

Inferring contemporary and historical genetic connectivity from juveniles

Measuring population connectivity is a critical task in conservation biology. While genetic markers can provide reliable long-term historical estimates of population connectivity, scientists are still limited in their ability to determine contemporary patterns of gene flow, the most practical time frame for management. Here, we tackled this issue by developing a new approach that only requires juvenile sampling at a single time period. To demonstrate the usefulness of our method, we used the Speartooth shark (Glyphis glyphis), a critically endangered species of river shark found only in tropical northern Australia and southern Papua New Guinea. Contemporary adult and juvenile shark movements, estimated with the spatial distribution of kin pairs across and within three river systems, was contrasted with historical long-term connectivity patterns, estimated from mitogenomes and genome-wide SNP data. We found strong support for river fidelity in juveniles with the within-cohort relationship analysis. Male breeding movements were highlighted with the cross-cohort relationship analysis, and female reproductive philopatry to the river systems was revealed by the mitogenomic analysis. We show that accounting for juvenile river fidelity and female philopatry is important in population structure analysis and that targeted sampling in nurseries and juvenile aggregations should be included in the genomic toolbox of threatened species management.

Comparison of pre-processing methodologies for Illumina 450k methylation array data in familial analyses

Background

Human methylome mapping in health and disease states has largely relied on Illumina Human Methylation 450k array (450k array) technology. Accompanying this has been the necessary evolution of analysis pipelines to facilitate data processing. The majority of these pipelines, however, cater for experimental designs where matched ‘controls’ or ‘normal’ samples are available. Experimental designs where no appropriate ‘reference’ exists remain challenging. Herein, we use data generated from our study of the inheritance of methylome profiles in families to evaluate the performance of eight normalisation pre-processing methods. Fifty individual samples representing four families were interrogated on five 450k array BeadChips. Eight normalisation methods were tested using qualitative and quantitative metrics, to assess efficacy and suitability.

Results

Stratified quantile normalisation combined with ComBat were consistently found to be the most appropriate when assessed using density, MDS and cluster plots. This was supported quantitatively by ANOVA on the first principal component where the effect of batch dropped from p < 0.01 to p = 0.97 after stratified QN and ComBat. Median absolute differences between replicated samples were the lowest after stratified QN and ComBat as were the standard error measures on known imprinted regions. Biological information was preserved after normalisation as indicated by the maintenance of a significant association between a known mQTL and methylation (p = 1.05e-05).

Conclusions

A strategy combining stratified QN with ComBat is appropriate for use in the analyses when no reference sample is available but preservation of biological variation is paramount. There is great potential for use of 450k array data to further our understanding of the methylome in a variety of similar settings. Such advances will be reliant on the determination of appropriate methodologies for processing these data such as established here.

Electronic supplementary material

The online version of this article (doi:10.1186/s13148-016-0241-2) contains supplementary material, which is available to authorized users.

Transient hybridization, not homoploid hybrid speciation, between ancient and deeply divergent conifers

PREMISE OF THE STUDY

Homoploid hybrid speciation is receiving growing attention due the increasing recognition of its role in speciation. We investigate if individuals intermediate in morphology between the two species of the conifer genus Athrotaxis represent a homoploid hybrid species, A. laxifolia, or are spontaneous F1 hybrids.

METHODS

A total of 1055 individuals of Athrotaxis cupressoides and A. selaginoides, morphologically intermediate individuals, and two putative hybrid swarms were sampled across the range of the genus and genotyped with 13 microsatellites. We used simulations to test the power of our data to identify the pure species, F1s, F2s, and backcross generations.

KEY RESULTS

We found that Athrotaxis cupressoides and A. selaginoides are likely the most divergent congeneric conifers known, but the intermediates are F1 hybrids, sharing one allele each from A. cupressoides and A. selaginoides at six loci with completely species specific alleles. The hybrid swarms contain wide genetic variation with stronger affinities to the locally dominant species, A. selaginoides and A. selaginoides backcrosses outnumbering A. cupressoides backcrosses. In addition, we observed evidence for isolated advanced generation backcrosses within the range of the pure species.

CONCLUSIONS

We conclude that, even though they can be large and long-lived, Athrotaxis hybrid swarms are on a trajectory of decline and will eventually be reabsorbed by the parental species. However, this process may take millennia and fossil evidence suggests that such events have occurred repeatedly since the early Quaternary. Given this timeline, our study highlights the many obstacles to homoploid hybrid speciation.

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.

A meta-analysis of genome-wide association studies to identify prostate cancer susceptibility loci associated with aggressive and non-aggressive disease

Genome-wide association studies (GWAS) have identified multiple common genetic variants associated with an increased risk of prostate cancer (PrCa), but these explain less than one-third of the heritability. To identify further susceptibility alleles, we conducted a meta-analysis of four GWAS including 5953 cases of aggressive PrCa and 11 463 controls (men without PrCa). We computed association tests for approximately 2.6 million SNPs and followed up the most significant SNPs by genotyping 49 121 samples in 29 studies through the international PRACTICAL and BPC3 consortia. We not only confirmed the association of a PrCa susceptibility locus, rs11672691 on chromosome 19, but also showed an association with aggressive PrCa [odds ratio = 1.12 (95% confidence interval 1.03-1.21), P = 1.4 × 10(-8)]. This report describes a genetic variant which is associated with aggressive PrCa, which is a type of PrCa associated with a poorer prognosis.

Emerging putative biomarkers: the role of alpha 2 and 6 integrins in susceptibility, treatment, and prognosis

The genetic architecture underpinning prostate cancer is complex, polygenic and despite recent significant advances many questions remain. Advances in genetic technologies have greatly improved our ability to identify genetic variants associated with complex disease including prostate cancer. Genome-wide association studies (GWASs) and microarray gene expression studies have identified genetic associations with prostate cancer susceptibility and tumour development. The integrins feature prominently in both studies examining the underlying genetic susceptibility and mechanisms driving prostate tumour development. Integrins are cell adhesion molecules involved in extracellular and intracellular signalling and are imperative for tumour development, migration, and angiogenesis. Although several integrins have been implicated in tumour development, the roles of integrin α(2) and integrin α(6) are the focus of this paper as evidence is now emerging that these integrins are implicit in prostate cancer susceptibility, cancer stem cell biology, angiogenesis, cell migration, and metastases to bone and represent potential biomarkers and therapeutic targets. There currently exists an urgent need to develop tools that differentiate indolent from aggressive prostate cancers and predict how patients will respond to treatment. This paper outlines the evidence supporting the use of α(2) and α(6) integrins in clinical applications for tailored patient treatment.

Low but structured chloroplast diversity in Atherosperma moschatum (Atherospermataceae) suggests bottlenecks in response to the Pleistocene glacials

BACKGROUND AND AIMS

The cool temperate rainforests of Australia were much reduced in range during the cold and dry glacial periods, although genetic evidence indicates that two key rainforest species, Nothofagus cunninghamii and Tasmannia lanceolata, survived within multiple locations and underwent only local range expansions at the end of the Last Glacial. To better understand the glacial response of a co-occurring but wind-dispersed and less cold-tolerant rainforest tree species, Atherosperma moschatum, a chloroplast phylogeographic study was undertaken.

METHODS

A total of 3294 bp of chloroplast DNA sequence was obtained for 155 samples collected from across the species' range.

KEY RESULTS

The distribution of six haplotypes observed in A. moschatum was geographically structured with an inferred ancestral haplotype restricted to Tasmania, while three non-overlapping and endemic haplotypes were found on the mainland of south-eastern Australia. Last glacial refugia for A. moschatum are likely to have occurred in at least one location in western Tasmania and in Victoria and within at least two locations in the Great Dividing Range of New South Wales. Nucleotide diversity of A. moschatum was lower (π = 0·00021) than either N. cunninghamii (0·00101) or T. lanceolata (0·00073), and was amongst the lowest recorded for any tree species.

CONCLUSIONS

This study provides evidence for past bottlenecks having impacted the chloroplast diversity of A. moschatum as a result of the species narrower climatic niche during glacials. This hypothesis is supported by the star-like haplotype network and similar estimated rates of chloroplast DNA substitution for A. moschatum and the two more cold tolerant and co-occurring species that have higher chloroplast diversity, N. cunninghamii and T. lanceolata.