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 (ntotal  = 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.

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.

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.

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.

Runs of homozygosity and a cluster of vulvar cancer in young Australian Aboriginal women.

OBJECTIVE: A cluster of vulvar cancer exists in young Aboriginal women living in remote communities in Arnhem Land, Australia. A genetic case-control study was undertaken involving 30 cases of invasive vulvar cancer and its precursor lesion, high-grade vulvar intraepithelial neoplasia (VIN), and 61 controls, matched for age and community of residence. It was hypothesized that this small, isolated population may exhibit increased autozygosity, implicating recessive effects as a possible mechanism for increased susceptibility to vulvar cancer. METHODS: Genotyping data from saliva samples were used to identify runs of homozygosity (ROH) in order to calculate estimates of genome-wide homozygosity. RESULTS: No evidence of an effect of genome-wide homozygosity on vulvar cancer and VIN in East Arnhem women was found, nor was any individual ROH found to be significantly associated with case status. This study found further evidence supporting an association between previous diagnosis of CIN and diagnosis of vulvar cancer or VIN, but found no association with any other medical history variable. CONCLUSIONS: These findings do not eliminate the possibility of genetic risk factors being involved in this cancer cluster, but rather suggest that alternative analytical strategies and genetic models should be explored.

Genetic and epigenetic variation in vulvar cancer: current research and future clinical practice.

Vulvar cancer is a relatively rare gynaecological malignancy, the treatment of which is associated with significant patient morbidity. With reports that the incidence of vulvar cancer is increasing, there is a rising need for improved preventive, diagnostic and therapeutic tools. Recent advances within genetics and epigenetics present possible approaches for addressing this need, by contributing to the clarification of the aetiology of this disease, identifying screening and drug targets and introducing the potential for personalised treatments. This paper reviews the genetic and epigenetic research undertaken to date within vulvar cancer, evaluates its potential for clinical application and identifies directions for future research.

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.

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.

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.

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.

Impact of the G84E variant on HOXB13 gene and protein expression in formalin-fixed, paraffin-embedded prostate tumours.

The HOXB13 G84E variant is associated with risk of prostate cancer (PCa), however the role this variant plays in PCa development is unknown. This study examined 751 cases, 450 relatives and 355 controls to determine the contribution of this variant to PCa risk in Tasmania and investigated HOXB13 gene and protein expression in tumours from nine G84E heterozygote variant and 13 wild-type carriers. Quantitative PCR and immunohistochemistry showed that HOXB13 gene and protein expression did not differ between tumour samples from variant and wild-type carriers. Allele-specific transcription revealed that two of seven G84E carriers transcribed both the variant and wild-type allele, while five carriers transcribed the wild-type allele. Methylation of surrounding CpG sites was lower in the variant compared to the wild-type allele, however overall methylation across the region was very low. Notably, tumour characteristics were less aggressive in the two variant carriers that transcribed the variant allele compared to the five that did not. This study has shown that HOXB13 expression does not differ between tumour tissue of G84E variant carriers and non-carriers. Intriguingly, the G84E variant allele was rarely transcribed in carriers, suggesting that HOXB13 expression may be driven by the wild-type allele in the majority of carriers.

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.

A retrospective examination of mean relative telomere length in the Tasmanian Familial Hematological Malignancies Study.

Telomere length has a biological link to cancer, with excessive telomere shortening leading to genetic instability and resultant malignant transformation. Telomere length is heritable and genetic variants determining telomere length have been identified. Telomere biology has been implicated in the development of hematological malignancies (HMs), therefore, closer examination of telomere length in HMs may provide further insight into genetic etiology of disease development and support for telomere length as a prognostic factor in HMs. We retrospectively examined mean relative telomere length in the Tasmanian Familial Hematological Malignancies Study using a quantitative PCR method on genomic DNA from peripheral blood samples. Fifty-five familial HM cases, 191 unaffected relatives of familial HM cases and 75 non-familial HM cases were compared with 758 population controls. Variance components modeling was employed to identify factors influencing variation in telomere length. Overall, HM cases had shorter mean relative telomere length (p=2.9×10-6) and this was observed across both familial and non-familial HM cases (p=2.2x10-4 and 2.2x10-5, respectively) as well as additional subgroupings of HM cases according to broad subtypes. Mean relative telomere length was also significantly heritable (62.6%; p=4.7x10-5) in the HM families in the present study. We present new evidence of significantly shorter mean relative telomere length in both familial and non-familial HM cases from the same population adding further support to the potential use of telomere length as a prognostic factor in HMs. Whether telomere shortening is the cause of or the result of HMs is yet to be determined, but as telomere length was found to be highly heritable in our HM families this suggests that genetics driving the variation in telomere length is related to HM disease risk.