A common regulatory haplotype doubles lactoferrin concentration in milk.

BACKGROUND: Bovine lactoferrin (Lf) is an iron absorbing whey protein with antibacterial, antiviral, and antifungal activity. Lactoferrin is economically valuable and has an extremely variable concentration in milk, partly driven by environmental influences such as milking frequency, involution, or mastitis. A significant genetic influence has also been previously observed to regulate lactoferrin content in milk. Here, we conducted genetic mapping of lactoferrin protein concentration in conjunction with RNA-seq, ChIP-seq, and ATAC-seq data to pinpoint candidate causative variants that regulate lactoferrin concentrations in milk. RESULTS: We identified a highly-significant lactoferrin protein quantitative trait locus (pQTL), as well as a cis lactotransferrin (LTF) expression QTL (cis-eQTL) mapping to the LTF locus. Using ChIP-seq and ATAC-seq datasets representing lactating mammary tissue samples, we also report a number of regions where the openness of chromatin is under genetic influence. Several of these also show highly significant QTL with genetic signatures similar to those highlighted through pQTL and eQTL analysis. By performing correlation analysis between these QTL, we revealed an ATAC-seq peak in the putative promotor region of LTF, that highlights a set of 115 high-frequency variants that are potentially responsible for these effects. One of the 115 variants (rs110000337), which maps within the ATAC-seq peak, was predicted to alter binding sites of transcription factors known to be involved in lactation-related pathways. CONCLUSIONS: Here, we report a regulatory haplotype of 115 variants with conspicuously large impacts on milk lactoferrin concentration. These findings could enable the selection of animals for high-producing specialist herds.

Successful editing and maintenance of lactogenic gene expression in primary bovine mammary epithelial cells.

In vitro investigation of bovine lactation processes is limited by a lack of physiologically representative cell models. This deficiency is most evident through the minimal or absent expression of lactation-specific genes in cultured bovine mammary tissues. Primary bovine mammary epithelial cells (pbMECs) extracted from lactating mammary tissue and grown in culture initially express milk protein transcripts at relatively representative levels. However, expression drops dramatically after only three or four passages, which greatly reduces the utility of primary cells to model and further examine lactogenesis. To investigate the effects of alternate alleles in pbMECs including effects on transcription, we have developed methods to deliver CRISPR-Cas9 gene editing reagents to primary mammary cells, resulting in very high editing efficiencies. We have also found that culturing the cells on an imitation basement membrane composed of Matrigel, results in the restoration of a more representative lactogenic gene expression profile and the cells forming three-dimensional structures in vitro. Here, we present data from four pbMEC lines recovered from pregnant cows and detail the expression profile of five key milk synthesis genes in these MECs grown on Matrigel. Additionally, we describe an optimised method for preferentially selecting CRISPR-Cas9-edited cells conferring a knock-out of DGAT1, using fluorescence-activated cell sorting (FACS). The combination of these techniques facilitates the use of pbMECs as a model to investigate the effects of gene introgressions and genetic variation in lactating mammary tissue.

Simple combination of multiple somatic variant callers to increase accuracy.

Publications comparing variant caller algorithms present discordant results with contradictory rankings. Caller performances are inconsistent and wide ranging, and dependent upon input data, application, parameter settings, and evaluation metric. With no single variant caller emerging as a superior standard, combinations or ensembles of variant callers have appeared in the literature. In this study, a whole genome somatic reference standard was used to derive principles to guide strategies for combining variant calls. Then, manually annotated variants called from the whole exome sequencing of a tumor were used to corroborate these general principles. Finally, we examined the ability of these principles to reduce noise in targeted sequencing.

Predictors of success in establishing orthotopic patient-derived xenograft models of triple negative breast cancer.

Patient-derived xenograft (PDX) models of breast cancer are an effective discovery platform and tool for preclinical pharmacologic testing and biomarker identification. We established orthotopic PDX models of triple negative breast cancer (TNBC) from the primary breast tumors of patients prior to and following neoadjuvant chemotherapy (NACT) while they were enrolled in the ARTEMIS trial (NCT02276443). Serial biopsies were obtained from patients prior to treatment (pre-NACT), from poorly responsive disease after four cycles of Adriamycin and cyclophosphamide (AC, mid-NACT), and in cases of AC-resistance, after a 3-month course of different experimental therapies and/or additional chemotherapy (post-NACT). Our study cohort includes a total of 269 fine needle aspirates (FNAs) from 217 women, generating a total of 62 PDX models (overall success-rate = 23%). Success of PDX engraftment was generally higher from those cancers that proved to be treatment-resistant, whether poorly responsive to AC as determined by ultrasound measurements mid-NACT (p = 0.063), RCB II/III status after NACT (p = 0.046), or metastatic relapse within 2 years of surgery (p = 0.008). TNBC molecular subtype determined from gene expression microarrays of pre-NACT tumors revealed no significant association with PDX engraftment rate (p = 0.877). Finally, we developed a statistical model predictive of PDX engraftment using percent Ki67 positive cells in the patient's diagnostic biopsy, positive lymph node status at diagnosis, and low volumetric reduction of the patient's tumor following AC treatment. This novel bank of 62 PDX models of TNBC provides a valuable resource for biomarker discovery and preclinical therapeutic trials aimed at improving neoadjuvant response rates for patients with TNBC.

Assessment of stained direct cytology smears of breast cancer for whole transcriptome and targeted messenger RNA sequencing.

BACKGROUND: Rather than surgical resection, cytologic specimens are often used as first-line clinical diagnostic procedures due to higher safety, speed, and cost-effectiveness. Archival diagnostic cytology slides containing cancer can be equivalent to tissue biopsies for DNA mutation testing, but the accuracy of transcriptomic profiling by RNA sequencing (RNA-seq) is less understood. METHODS: This study compares the results from whole transcriptome RNA-seq and a targeted RNA-seq assay of stained cytology smears (CS) versus matched tumor tissue samples preserved fresh-frozen (FF) and processed as formalin-fixed paraffin-embedded (FFPE) sections. Cellular cytology scrapes from all 11 breast cancers were fixed and stained using three common protocols: Carnoy's (CS_C) or 95% ethanol (CS_E) fixation and then Papanicolaou stain or air-dried then methanol fixation and DiffQuik stain (CS_DQ). Agreement between samples was assessed using Lin's concordance correlation coefficient. RESULTS: Library yield for CS_DQ was too low, therefore it was not sequenced. The distributions of concordance correlation coefficient of gene expression levels in comparison to FF were comparable between CS_C and CS_E, but expression of genes enriched in stroma was lower in cytosmear samples than in FF or FFPE. Six signatures showed similar concordance to FF for all methods and two were slightly worse in CS_C and CS_E. Genomic signatures were highly concordant using targeted RNA-seq. The allele fraction of selected mutations calculated on cytosmear specimens was highly correlated with FF tissues using both RNA-seq methods. CONCLUSION: RNA can be reliably extracted from cytology smears and is suitable for transcriptome profiling or mutation detection, except for signatures of tumor stroma.

Targeted RNAseq assay incorporating unique molecular identifiers for improved quantification of gene expression signatures and transcribed mutation fraction in fixed tumor samples.

BACKGROUND: Our objective was to assess whether modifications to a customized targeted RNA sequencing (RNAseq) assay to include unique molecular identifiers (UMIs) that collapse read counts to their source mRNA counts would improve quantification of transcripts from formalin-fixed paraffin-embedded (FFPE) tumor tissue samples. The assay (SET4) includes signatures that measure hormone receptor and PI3-kinase related transcriptional activity (SETER/PR and PI3Kges), and measures expression of selected activating point mutations and key breast cancer genes. METHODS: Modifications included steps to introduce eight nucleotides-long UMIs during reverse transcription (RT) in bulk solution, followed by polymerase chain reaction (PCR) of labeled cDNA in droplets, with optimization of the polymerase enzyme and reaction conditions. We used Lin's concordance correlation coefficient (CCC) to measure concordance, including precision (Rho) and accuracy (Bias), and nonparametric tests (Wilcoxon, Levene's) to compare the modified (NEW) SET4 assay to the original (OLD) SET4 assay and to whole transcriptome RNAseq using RNA from matched fresh frozen (FF) and FFPE samples from 12 primary breast cancers. RESULTS: The modified (NEW) SET4 assay measured single transcripts (p< 0.001) and SETER/PR (p=0.002) more reproducibly in technical replicates from FFPE samples. The modified SET4 assay was more precise for measuring single transcripts (Rho 0.966 vs 0.888, p< 0.01) but not multigene expression signatures SETER/PR (Rho 0.985 vs 0.968) or PI3Kges (Rho 0.985 vs 0.946) in FFPE, compared to FF samples. It was also more precise than wtRNAseq of FFPE for measuring transcripts (Rho 0.986 vs 0.934, p< 0.001) and SETER/PR (Rho 0.993 vs 0.915, p=0.004), but not PI3Kges (Rho 0.988 vs 0.945, p=0.051). Accuracy (Bias) was comparable between protocols. Two samples carried a PIK3CA mutation, and measurements of transcribed mutant allele fraction was similar in FF and FFPE samples and appeared more precise with the modified SET4 assay. Amplification efficiency (reads per UMI) was consistent in FF and FFPE samples, and close to the theoretically expected value, when the library size exceeded 400,000 aligned reads. CONCLUSIONS: Modifications to the targeted RNAseq protocol for SET4 assay significantly increased the precision of UMI-based and reads-based measurements of individual transcripts, multi-gene signatures, and mutant transcript fraction, particularly with FFPE samples.

The impact of RNA extraction method on accurate RNA sequencing from formalin-fixed paraffin-embedded tissues.

BACKGROUND: Utilization of RNA sequencing methods to measure gene expression from archival formalin-fixed paraffin-embedded (FFPE) tumor samples in translational research and clinical trials requires reliable interpretation of the impact of pre-analytical variables on the data obtained, particularly the methods used to preserve samples and to purify RNA. METHODS: Matched tissue samples from 12 breast cancers were fresh frozen (FF) and preserved in RNAlater or fixed in formalin and processed as FFPE tissue. Total RNA was extracted and purified from FF samples using the Qiagen RNeasy kit, and in duplicate from FFPE tissue sections using three different kits (Norgen, Qiagen and Roche). All RNA samples underwent whole transcriptome RNA sequencing (wtRNAseq) and targeted RNA sequencing for 31 transcripts included in a signature of sensitivity to endocrine therapy. We assessed the effect of RNA extraction kit on the reliability of gene expression levels using linear mixed-effects model analysis, concordance correlation coefficient (CCC) and differential analysis. All protein-coding genes in the wtRNAseq and three gene expression signatures for breast cancer were assessed for concordance. RESULTS: Despite variable quality of the RNA extracted from FFPE samples by different kits, all had similar concordance of overall gene expression from wtRNAseq between matched FF and FFPE samples (median CCC 0.63-0.66) and between technical replicates (median expression difference 0.13-0.22). More than half of genes were differentially expressed between FF and FFPE, but with low fold change (median |LFC| 0.31-0.34). Two out of three breast cancer signatures studied were highly robust in all samples using any kit, whereas the third signature was similarly discordant irrespective of the kit used. The targeted RNAseq assay was concordant between FFPE and FF samples using any of the kits (CCC 0.91-0.96). CONCLUSIONS: The selection of kit to purify RNA from FFPE did not influence the overall quality of results from wtRNAseq, thus variable reproducibility of gene signatures probably relates to the reliability of individual gene selected and possibly to the algorithm. Targeted RNAseq showed promising performance for clinical deployment of quantitative assays in breast cancer from FFPE samples, although numerical scores were not identical to those from wtRNAseq and would require calibration.

MelanomaDB: A Web Tool for Integrative Analysis of Melanoma Genomic Information to Identify Disease-Associated Molecular Pathways.

Despite on-going research, metastatic melanoma survival rates remain low and treatment options are limited. Researchers can now access a rapidly growing amount of molecular and clinical information about melanoma. This information is becoming difficult to assemble and interpret due to its dispersed nature, yet as it grows it becomes increasingly valuable for understanding melanoma. Integration of this information into a comprehensive resource to aid rational experimental design and patient stratification is needed. As an initial step in this direction, we have assembled a web-accessible melanoma database, MelanomaDB, which incorporates clinical and molecular data from publically available sources, which will be regularly updated as new information becomes available. This database allows complex links to be drawn between many different aspects of melanoma biology: genetic changes (e.g., mutations) in individual melanomas revealed by DNA sequencing, associations between gene expression and patient survival, data concerning drug targets, biomarkers, druggability, and clinical trials, as well as our own statistical analysis of relationships between molecular pathways and clinical parameters that have been produced using these data sets. The database is freely available at http://genesetdb.auckland.ac.nz/melanomadb/about.html. A subset of the information in the database can also be accessed through a freely available web application in the Illumina genomic cloud computing platform BaseSpace at http://www.biomatters.com/apps/melanoma-profiler-for-research. The MelanomaDB database illustrates dysregulation of specific signaling pathways across 310 exome-sequenced melanomas and in individual tumors and identifies the distribution of somatic variants in melanoma. We suggest that MelanomaDB can provide a context in which to interpret the tumor molecular profiles of individual melanoma patients relative to biological information and available drug therapies.

Selection pressure on the hepatitis B virus pre-S/S and P open reading frames in Tongan subjects with a chronic hepatitis B virus infection.

Identification of the full repertoire of hepatitis B virus (HBV) peptides that are presented to CD8+ T cells by common HLA class I alleles will be useful for designing immunotherapies for chronic hepatitis B. One hundred and seventy five cloned sequences containing the pre-S/S and P open reading frames (ORF) of the HBV were obtained from serum HBV-DNA of HBeAg-positive (n=4) and HBeAg-negative (inactive healthy carriers (IHC), n=16) Tongan subjects with an inactive chronic HBV infection. In addition, 34 and 32 sequences were obtained 5.2±1.4 (mean±SD) years apart from eight subjects. PAML was used to identify codons in the pre-S/S and P ORFs that were under positive selection pressure (ω>1). The number of non-synonymous substitutions in these codons was compared in IHC who were homozygous for either HLA-B∗4001 (n=9) or HLA-B*5602 (n=7), and who were either positive (n=6) or negative (n=10) for HLA-A*02. 34 codons in the pre-S/S and 11 codons in the P ORFs were under positive selection pressure. There was a higher number of non-synonymous substitutions in these codons in HBeAg-negative versus HBeAg-positive subjects in the P (p=0.02) but not the pre-S/S (p=0.64) ORF. There was no association between any HLA class I allele and non-synonymous substitutions in these codons. There was no increase in positive selection pressure on the pre-S/S and P ORFs with time. In conclusion, we could not find HLA class I-restricted selection pressure on any pre-S/S or P ORF amino acid; raising the possibility that peptide-based immunotherapies for chronic hepatitis B may not require peptides from these ORFs.

Involvement of TRP-like channels in the acute ischemic response of hippocampal CA1 neurons in brain slices.

During a period of acute ischemia in vivo or oxygen-glucose deprivation (OGD) in vitro, CA1 neurons depolarize, swell and become overloaded with calcium. Our aim was to test the hypothesis that the initial responses to OGD are at least partly due to transient receptor potential (TRP) channel activation. As some TRP channels are temperature-sensitive, we also compared the effects of pharmacological blockade of the channels with the effects of reducing temperature. Acute hippocampal slices (350 mum) obtained from Wistar rats were submerged in ACSF at 36 degrees C. CA1 neurons were monitored electrophysiologically using extracellular, intracellular or whole-cell patch-clamp recordings. Cell swelling was assessed by recording changes in relative tissue resistance, and changes in intracellular calcium were measured after loading neurons with fura-2 dextran. Blockers of TRP channels (ruthenium red, La3+, Gd3+, 2-APB) or lowering temperature by 3 degrees C reduced responses to OGD. This included: (a) an increased delay to negative shifts of extracellular DC potential; (b) reduction in rate of the initial slow membrane depolarization, slower development of OGD-induced increase in cell input resistance and slower development of whole-cell inward current; (c) reduced tissue swelling; and (d) a smaller rise in intracellular calcium. Mild hypothermia (33 degrees C) and La3+ or Gd3+ (100 microM) showed an occlusion effect when delay to extracellular DC shifts was measured. Expression of TRPM2/TRPM7 (oxidative stress-sensitive) and TRPV3/TRPV4 (temperature-sensitive) channels was demonstrated in the CA1 subfield with RT-PCR. These results indicate that TRP or TRP-like channels are activated by cellular stress and contribute to ischemia-induced membrane depolarization, intracellular calcium accumulation and cell swelling. We also hypothesize that closing of some TRP channels (TRPV3 and/or TRPV4) by lowering temperature may be partly responsible for the neuroprotective effect of hypothermia.