Analysis of gliomas DNA methylation: Assessment of pre-analytical variables.

Precision oncology is driven by biomarkers. For glioblastoma multiforme (GBM), the most common malignant adult primary brain tumor, O6-methylguanine-DNA methyltransferase (MGMT) gene promoter methylation is an important prognostic and treatment clinical biomarker. Time consuming pre-analytical steps such as biospecimen storage, fixation, sampling, and processing are sources of data irreproducibility, and all these pre-analytical variables are confounded by intratumor heterogeneity of MGMT promoter methylation. To assess the effect of pre-analytical variables on GBM DNA methylation, tissue storage/sampling (CryoGrid), sample preparation multi-sonicator (PIXUL), and 5-methylcytosine (5mC) DNA immunoprecipitation (Matrix MeDIP-qPCR/seq) platforms were used. MGMT promoter methylation status assayed by MeDIP-qPCR was validated with methylation specific PCR (MS-PCR). MGMT promoter methylation levels in frozen and formalin fixed paraffin embedded (FFPE) sample pairs were not statistically different, confirming reliability of FFPEs for MGMT promoter methylation analysis. Warm ex-vivo ischemia (up to 4hrs at 37oC) and 3 cycles of repeated sample thawing and freezing did not statistically impact 5mC at MGMT promoter, exon, and enhancer regions, indicating the resistance of DNA methylation to common variations in sample processing conditions that might be encountered in research and clinical settings. 26-34% of specimens exhibited intratumor heterogeneity in the MGMT DNA promoter methylation. These data demonstrate that variations in sample fixation, ischemia duration and temperature, and DNA methylation assay technique do not have a statistically significant impact on MGMT promoter methylation assessment. However, intratumor methylation heterogeneity underscores the value of multiple biopsies at different GBM geographic tumor sites in the evaluation of MGMT promoter methylation status. Matrix-MeDIP-seq analysis revealed that MGMT promoter methylation status clustered with other differentially methylated genomic loci (e.g. HOXA and lncRNAs) that are resilient to variation in the above pre-analytical conditions. These observations offer new opportunities to develop more granular data-based epigenetic GBM biomarkers. In this regard, the high throughput CryoGrid-PIXUL-Matrix toolbox could be useful.

A prognostic matrix gene expression signature defines functional glioblastoma phenotypes and niches.

Background: Interactions among tumor, immune, and vascular niches play major roles in driving glioblastoma (GBM) malignancy and treatment responses. The composition, heterogeneity, and localization of extracellular core matrix proteins (CMPs) that mediate such interactions, however, are not well understood. Methods: Here, through computational genomics and proteomics approaches, we analyzed the functional and clinical relevance of CMP expression in GBM at bulk, single cell, and spatial anatomical resolution. Results: We identified genes encoding CMPs whose expression levels categorize GBM tumors into CMP expression-high (M-H) and CMP expression-low (M-L) groups. CMP enrichment is associated with worse patient survival, specific driver oncogenic alterations, mesenchymal state, infiltration of pro-tumor immune cells, and immune checkpoint gene expression. Anatomical and single-cell transcriptome analyses indicate that matrisome gene expression is enriched in vascular and leading edge/infiltrative niches that are known to harbor glioma stem cells driving GBM progression. Finally, we identified a 17-gene CMP expression signature, termed Matrisome 17 (M17) signature that further refines the prognostic value of CMP genes. The M17 signature is a significantly stronger prognostic factor compared to MGMT promoter methylation status as well as canonical subtypes, and importantly, potentially predicts responses to PD1 blockade. Conclusion: The matrisome gene expression signature provides a robust stratification of GBM patients by survival and potential biomarkers of functionally relevant GBM niches that can mediate mesenchymal-immune cross talk. Patient stratification based on matrisome profiles can contribute to selection and optimization of treatment strategies.

Analysis of gliomas DNA methylation: Assessment of pre-analytical variables.

Precision oncology is driven by molecular biomarkers. For glioblastoma multiforme (GBM), the most common malignant adult primary brain tumor, O6-methylguanine-DNA methyltransferase ( MGMT ) gene DNA promoter methylation is an important prognostic and treatment clinical biomarker. Time consuming pre-analytical steps such as biospecimen storage before fixing, sampling, and processing are major sources of errors and batch effects, that are further confounded by intra-tumor heterogeneity of MGMT promoter methylation. To assess the effect of pre-analytical variables on GBM DNA methylation, tissue storage/sampling (CryoGrid), sample preparation multi-sonicator (PIXUL) and 5-methylcytosine (5mC) DNA immunoprecipitation (Matrix MeDIP-qPCR/seq) platforms were used. MGMT promoter CpG methylation was examined in 173 surgical samples from 90 individuals, 50 of these were used for intra-tumor heterogeneity studies. MGMT promoter methylation levels in paired frozen and formalin fixed paraffin embedded (FFPE) samples were very close, confirming suitability of FFPE for MGMT promoter methylation analysis in clinical settings. Matrix MeDIP-qPCR yielded similar results to methylation specific PCR (MS-PCR). Warm ex-vivo ischemia (37°C up to 4hrs) and 3 cycles of repeated sample thawing and freezing did not alter 5mC levels at MGMT promoter, exon and upstream enhancer regions, demonstrating the resistance of DNA methylation to the most common variations in sample processing conditions that might be encountered in research and clinical settings. 20-30% of specimens exhibited intratumor heterogeneity in the MGMT DNA promoter methylation. Collectively these data demonstrate that variations in sample fixation, ischemia duration and temperature, and DNA methylation assay technique do not have significant impact on assessment of MGMT promoter methylation status. However, intratumor methylation heterogeneity underscores the need for histologic verification and value of multiple biopsies at different GBM geographic tumor sites in assessment of MGMT promoter methylation. Matrix-MeDIP-seq analysis revealed that MGMT promoter methylation status clustered with other differentially methylated genomic loci (e.g. HOXA and lncRNAs), that are likewise resilient to variation in above post-resection pre-analytical conditions. These MGMT -associated global DNA methylation patterns offer new opportunities to validate more granular data-based epigenetic GBM clinical biomarkers where the CryoGrid-PIXUL-Matrix toolbox could prove to be useful.

A prognostic matrix code defines functional glioblastoma phenotypes and niches.

Interactions among tumor, immune and vascular niches play major roles in driving glioblastoma (GBM) malignancy and treatment responses. The composition, heterogeneity, and localization of extracellular core matrix proteins (CMPs) that mediate such interactions, however, are not well understood. Here, we characterize functional and clinical relevance of genes encoding CMPs in GBM at bulk, single cell, and spatial anatomical resolution. We identify a "matrix code" for genes encoding CMPs whose expression levels categorize GBM tumors into matrisome-high and matrisome-low groups that correlate with worse and better patient survival, respectively. The matrisome enrichment is associated with specific driver oncogenic alterations, mesenchymal state, infiltration of pro-tumor immune cells and immune checkpoint gene expression. Anatomical and single cell transcriptome analyses indicate that matrisome gene expression is enriched in vascular and leading edge/infiltrative anatomic structures that are known to harbor glioma stem cells driving GBM progression. Finally, we identified a 17-gene matrisome signature that retains and further refines the prognostic value of genes encoding CMPs and, importantly, potentially predicts responses to PD1 blockade in clinical trials for GBM. The matrisome gene expression profiles provide potential biomarkers of functionally relevant GBM niches that contribute to mesenchymal-immune cross talk and patient stratification which could be applied to optimize treatment responses.

A prognostic matrix code defines functional glioblastoma phenotypes and niches.

Interactions among tumor, immune and vascular niches play major roles in driving glioblastoma (GBM) malignancy and treatment responses. The composition, heterogeneity, and localization of extracellular core matrix proteins (CMPs) that mediate such interactions, however, are not well understood. Here, we characterize functional and clinical relevance of genes encoding CMPs in GBM at bulk, single cell, and spatial anatomical resolution. We identify a "matrix code" for genes encoding CMPs whose expression levels categorize GBM tumors into matrisome-high and matrisome-low groups that correlate with worse and better survival, respectively, of patients. The matrisome enrichment is associated with specific driver oncogenic alterations, mesenchymal state, infiltration of pro-tumor immune cells and immune checkpoint gene expression. Anatomical and single cell transcriptome analyses indicate that matrisome gene expression is enriched in vascular and leading edge/infiltrative anatomic structures that are known to harbor glioma stem cells driving GBM progression. Finally, we identified a 17-gene matrisome signature that retains and further refines the prognostic value of genes encoding CMPs and, importantly, potentially predicts responses to PD1 blockade in clinical trials for GBM. The matrisome gene expression profiles may provide biomarkers of functionally relevant GBM niches that contribute to mesenchymal-immune cross talk and patient stratification to optimize treatment responses.

Emerging treatment strategies for breast cancer brain metastasis: from translational therapeutics to real-world experience.

Systemic therapies for primary breast cancer have made great progress over the past two decades. However, oncologists confront an insidious and particularly difficult problem: in those patients with metastatic breast cancer, up to 50% of human epidermal growth factor 2 (HER2)-positive and 25-40% of triple-negative subtypes, brain metastases (BM) kill most of them. Fortunately, standard- of-care treatments for BM have improved rapidly, with a decline in whole brain radiation therapy and use of fractionated stereotactic radiosurgery as well as targeted therapies and immunotherapies. Meanwhile, advances in fundamental understanding of the basic biological processes of breast cancer BM (BCBM) have led to many novel experimental therapeutic strategies. In this review, we describe the most recent clinical treatment options and emerging experimental therapeutic strategies that have the potential to combat BCBM.

The identification of a TNBC liver metastasis gene signature by sequential CTC-xenograft modeling.

Triple-negative breast cancer (TNBC) liver metastasis is associated with poor prognosis and low patient survival. It occurs when tumor cells disseminate from primary tumors, circulate in blood/lymph [circulating tumor cells (CTCs)], and acquire distinct characteristics during disease progression toward the metastatic phenotype. The purpose of this study was to decipher the genomic/transcriptomic properties of TNBC liver metastasis and its recurrence for potential therapeutic targeting. We employed a negative depletion strategy to isolate and interrogate CTCs from the blood of patients with TNBC, and to establish sequential generations of CTC-derived xenografts (CDXs) through injection of patient CTCs in immunodeficient mice. The isolation and validation of CDX-derived cell populations [analyses of CTCs were paired with bone marrow-resident cells (BMRTCs) and liver tissue cells obtained from the same animal] were performed by multiparametric flow cytometry, immune phenotyping, and genomic sequencing of putative CTCs. Comprehensive characterization of gene expression arrays from sequentially generated CDX-derived cell populations, online gene expression arrays, and TCGA databases were employed to discover a CTC-driven, liver metastasis-associated TNBC signature. We discovered a distinct transcriptomic signature of TNBC patient-isolated CTCs from primary TNBCs, which was consistent throughout sequential CDX modeling. We established a novel TNBC liver metastasis-specific CDX model that selectively recapitulates CTC biology for four sequential generations of mice. The evaluation of online databases and CDX-derived populations revealed 597 genes specific to the TNBC liver metastasis signatures. Further investigation of the TNBC liver metastasis signature predicted 16 hub genes, 6 biomarkers with clinically available drugs, and 22 survival genes. The sequential interrogation of CDX-CTCs is an innovative liquid biopsy-based approach for the discovery of organ metastasis-specific signatures of CTCs. This represents the first step for mechanistic and analytical validation in their application as prognostic indicators and therapeutic targets. Targeting CTC drug candidate biomarkers along with combination therapy can improve the clinical outcome of TNBC patients in general and recurrence of liver metastasis in particular.

PMN-MDSCs Enhance CTC Metastatic Properties through Reciprocal Interactions via ROS/Notch/Nodal Signaling.

Intratumoral infiltration of myeloid-derived suppressor cells (MDSCs) is known to promote neoplastic growth by inhibiting the tumoricidal activity of T cells. However, direct interactions between patient-derived MDSCs and circulating tumors cells (CTCs) within the microenvironment of blood remain unexplored. Dissecting interplays between CTCs and circulatory MDSCs by heterotypic CTC/MDSC clustering is critical as a key mechanism to promote CTC survival and sustain the metastatic process. We characterized CTCs and polymorphonuclear-MDSCs (PMN-MDSCs) isolated in parallel from peripheral blood of metastatic melanoma and breast cancer patients by multi-parametric flow cytometry. Transplantation of both cell populations in the systemic circulation of mice revealed significantly enhanced dissemination and metastasis in mice co-injected with CTCs and PMN-MDSCs compared to mice injected with CTCs or MDSCs alone. Notably, CTC/PMN-MDSC clusters were detected in vitro and in vivo either in patients' blood or by longitudinal monitoring of blood from animals. This was coupled with in vitro co-culturing of cell populations, demonstrating that CTCs formed physical clusters with PMN-MDSCs; and induced their pro-tumorigenic differentiation through paracrine Nodal signaling, augmenting the production of reactive oxygen species (ROS) by PMN-MDSCs. These findings were validated by detecting significantly higher Nodal and ROS levels in blood of cancer patients in the presence of naïve, heterotypic CTC/PMN-MDSC clusters. Augmented PMN-MDSC ROS upregulated Notch1 receptor expression in CTCs through the ROS-NRF2-ARE axis, thus priming CTCs to respond to ligand-mediated (Jagged1) Notch activation. Jagged1-expressing PMN-MDSCs contributed to enhanced Notch activation in CTCs by engagement of Notch1 receptor. The reciprocity of CTC/PMN-MDSC bi-directional paracrine interactions and signaling was functionally validated in inhibitor-based analyses, demonstrating that combined Nodal and ROS inhibition abrogated CTC/PMN-MDSC interactions and led to a reduction of CTC survival and proliferation. This study provides seminal evidence showing that PMN-MDSCs, additive to their immuno-suppressive roles, directly interact with CTCs and promote their dissemination and metastatic potency. Targeting CTC/PMN-MDSC heterotypic clusters and associated crosstalks can therefore represent a novel therapeutic avenue for limiting hematogenous spread of metastatic disease.

Targeting USP7 Identifies a Metastasis-Competent State within Bone Marrow-Resident Melanoma CTCs.

Systemic metastasis is the major cause of death from melanoma, the most lethal form of skin cancer. Although most patients with melanoma exhibit a substantial gap between onset of primary and metastatic tumors, signaling mechanisms implicated in the period of metastatic latency remain unclear. We hypothesized that melanoma circulating tumor cells (CTC) home to and reside in the bone marrow during the asymptomatic phase of disease progression. Using a strategy to deplete normal cell lineages (Lin-), we isolated CTC-enriched cell populations from the blood of patients with metastatic melanoma, verified by the presence of putative CTCs characterized by melanoma-specific biomarkers and upregulated gene transcripts involved in cell survival and prodevelopment functions. Implantation of Lin- population in NSG mice (CTC-derived xenografts, i.e., CDX), and subsequent transcriptomic analysis of ex vivo bone marrow-resident tumor cells (BMRTC) versus CTC identified protein ubiquitination as a significant regulatory pathway of BMRTC signaling. Selective inhibition of USP7, a key deubiquinating enzyme, arrested BMRTCs in bone marrow locales and decreased systemic micrometastasis. This study provides first-time evidence that the asymptomatic progression of metastatic melanoma can be recapitulated in vivo using patient-isolated CTCs. Furthermore, these results suggest that USP7 inhibitors warrant further investigation as a strategy to prevent progression to overt clinical metastasis.Significance: These findings provide insights into mechanism of melanoma recurrence and propose a novel approach to inhibit systematic metastatic disease by targeting bone marrow-resident tumor cells through pharmacological inhibition of USP7.Graphical Abstract: http://cancerres.aacrjournals.org/content/canres/78/18/5349/F1.large.jpg Cancer Res; 78(18); 5349-62. ©2018 AACR.

Molecular characterization of breast cancer CTCs associated with brain metastasis.

The enumeration of EpCAM-positive circulating tumor cells (CTCs) has allowed estimation of overall metastatic burden in breast cancer patients. However, a thorough understanding of CTCs associated with breast cancer brain metastasis (BCBM) is necessary for early identification and evaluation of treatment response to BCBM. Here we report that BCBM CTCs is enriched in a distinct sub-population of cells identifiable by their biomarker expression and mutational content. Deriving from a comprehensive analysis of CTC transcriptomes, we discovered a unique "circulating tumor cell gene signature" that is distinct from primary breast cancer tissues. Further dissection of the circulating tumor cell gene signature identified signaling pathways associated with BCBM CTCs that may have roles in potentiating BCBM. This study proposes CTC biomarkers and signaling pathways implicated in BCBM that may be used either as a screening tool for brain micro-metastasis detection or for making rational treatment decisions and monitoring therapeutic response in patients with BCBM.Characterization of CTCs derived from breast cancer patients with brain metastasis (BCBM) may allow for early diagnosis of brain metastasis and/or help for treatment choice and its efficacy. In this study, the authors identify a unique signature, based on patient-derived CTCs transcriptomes, for BCBM- CTCs that is different from primary tumors.

The isolation and characterization of CTC subsets related to breast cancer dormancy.

Uncovering CTCs phenotypes offer the promise to dissect their heterogeneity related to metastatic competence. CTC survival rates are highly variable and this can lead to many questions as yet unexplored properties of CTCs responsible for invasion and metastasis vs dormancy. We isolated CTC subsets from peripheral blood of patients diagnosed with or without breast cancer brain metastasis. CTC subsets were selected for EpCAM negativity but positivity for CD44(+)/CD24(-) stem cell signature; along with combinatorial expression of uPAR and int ß1, two markers directly implicated in breast cancer dormancy mechanisms. CTC subsets were cultured in vitro generating 3D CTC tumorspheres which were interrogated for biomarker profiling and biological characteristics. We identified proliferative and invasive properties of 3D CTC tumorspheres distinctive upon uPAR/int ß1 combinatorial expression. The molecular characterization of uPAR/int ß1 CTC subsets may enhance abilities to prospectively identify patients who may be at high risk of developing BCBM.

In vivo functional characterization of the transmembrane histidine kinase KinC in Bacillus subtilis.

In response to starvation, Bacillus subtilis cells differentiate into different subsets, undergoing cannibalism, biofilm formation or sporulation. These processes require a multiple component phosphorelay, wherein the master regulator Spo0A is activated upon phosphorylation by one or a combination of five histidine kinases (KinA-KinE) via two intermediate phosphotransferases, Spo0F and Spo0B. In this study, we focused on KinC, which was originally identified as a sporulation kinase and was later shown to regulate cannibalism and biofilm formation. First, genetic experiments using both the domesticated and undomesticated (biofilm forming) strains revealed that KinC activity and the membrane localization are independent of both the lipid raft marker proteins FloTA and cytoplasmic potassium concentration, which were previously shown to be required for the kinase activity. Next, we demonstrated that KinC controls cannibalism and biofilm formation in a manner dependent on phosphorelay. For further detailed characterization of KinC, we established an IPTG-inducible expression system in the domesticated strain, in which biofilm formation is defective, for simplicity of study. Using this system, we found that the N-terminal transmembrane domain is dispensable but the PAS domain is needed for the kinase activity. An in vivo chemical cross-linking experiment demonstrated that the soluble and functional KinC (KinC(ΔTM1+2)) forms a tetramer. Based on these results, we propose a revised model in which KinC becomes active by forming a homotetramer via the N-terminal PAS domain, but its activity is independent of both the lipid raft and the potassium leakage, which was previously suggested to be induced by surfactin.

Circulating tumor cells from a 4-dimensional lung cancer model are resistant to cisplatin.

OBJECTIVE: To determine the effect of cisplatin on circulatory tumor cells (CTC) and tumor nodules in a four-dimensional (4D) lung cancer model. METHODS: CTCs from the 4D model seeded with H1299, A549, or H460 and respective cells that were grown under two-dimensional conditions in a Petri dish were treated with 50 µM cisplatin for 24 and 48 hours and cell viability was determined. The lung nodules in the 4D model were then treated with different continuous or intermittent doses of cisplatin and the nodule size, the number of CTCs, and the level of matrix metalloproteinase (MMP) were determined. RESULTS: Cisplatin led to a significant decrease in the viability of tumor cells grown under 2D conditions (P < .01) but not in CTCs from the 4D model after both 24 hours and 48 hours. Cisplatin led to regression of tumor nodules with both the continuous and intermittent treatments. Moreover, there was a significantly higher number of CTCs per tumor area (P < .05) and MMP-2 production per tumor area (P = .007) for all human lung cancer cell lines grown in the 4D model when treated with cisplatin. CONCLUSIONS: The 4D lung cancer model allows for the isolation of CTCs that are resistant to cisplatin treatment. The model may allow us to better understand the biology of cisplatin resistance.

Triggering sporulation in Bacillus subtilis with artificial two-component systems reveals the importance of proper Spo0A activation dynamics.

Sporulation initiation in Bacillus subtilis is controlled by the phosphorylated form of the master regulator Spo0A which controls transcription of a multitude of sporulation genes. In this study, we investigated the importance of temporal dynamics of phosphorylated Spo0A (Spo0A∼P) accumulation by rewiring the network controlling its phosphorylation. We showed that simultaneous induction of KinC, a kinase that can directly phosphorylate Spo0A, and Spo0A itself from separately controlled inducible promoters can efficiently trigger sporulation even under nutrient rich conditions. However, the sporulation efficiency in this artificial two-component system was significantly impaired when KinC and/or Spo0A induction was too high. Using mathematical modelling, we showed that gradual accumulation of Spo0A∼P is essential for the proper temporal order of the Spo0A regulon expression, and that reduction in sporulation efficiency results from the reversal of that order. These insights led us to identify premature repression of DivIVA as one possible explanation for the adverse effects of accelerated accumulation of Spo0A∼P on sporulation. Moreover, we found that positive feedback resulting from autoregulation of the native spo0A promoter leads to robust control of Spo0A∼P accumulation kinetics. Thus we propose that a major function of the conserved architecture of the sporulation network is controlling Spo0A activation dynamics.

Stallion sperm transcriptome comprises functionally coherent coding and regulatory RNAs as revealed by microarray analysis and RNA-seq.

Mature mammalian sperm contain a complex population of RNAs some of which might regulate spermatogenesis while others probably play a role in fertilization and early development. Due to this limited knowledge, the biological functions of sperm RNAs remain enigmatic. Here we report the first characterization of the global transcriptome of the sperm of fertile stallions. The findings improved understanding of the biological significance of sperm RNAs which in turn will allow the discovery of sperm-based biomarkers for stallion fertility. The stallion sperm transcriptome was interrogated by analyzing sperm and testes RNA on a 21,000-element equine whole-genome oligoarray and by RNA-seq. Microarray analysis revealed 6,761 transcripts in the sperm, of which 165 were sperm-enriched, and 155 were differentially expressed between the sperm and testes. Next, 70 million raw reads were generated by RNA-seq of which 50% could be aligned with the horse reference genome. A total of 19,257 sequence tags were mapped to all horse chromosomes and the mitochondrial genome. The highest density of mapped transcripts was in gene-rich ECA11, 12 and 13, and the lowest in gene-poor ECA9 and X; 7 gene transcripts originated from ECAY. Structural annotation aligned sperm transcripts with 4,504 known horse and/or human genes, rRNAs and 82 miRNAs, whereas 13,354 sequence tags remained anonymous. The data were aligned with selected equine gene models to identify additional exons and splice variants. Gene Ontology annotations showed that sperm transcripts were associated with molecular processes (chemoattractant-activated signal transduction, ion transport) and cellular components (membranes and vesicles) related to known sperm functions at fertilization, while some messenger and micro RNAs might be critical for early development. The findings suggest that the rich repertoire of coding and non-coding RNAs in stallion sperm is not a random remnant from spermatogenesis in testes but a selectively retained and functionally coherent collection of RNAs.

Genome-wide association study implicates testis-sperm specific FKBP6 as a susceptibility locus for impaired acrosome reaction in stallions.

Impaired acrosomal reaction (IAR) of sperm causes male subfertility in humans and animals. Despite compelling evidence about the genetic control over acrosome biogenesis and function, the genomics of IAR is as yet poorly understood, providing no molecular tools for diagnostics. Here we conducted Equine SNP50 Beadchip genotyping and GWAS using 7 IAR-affected and 37 control Thoroughbred stallions. A significant (P<6.75E-08) genotype-phenotype association was found in horse chromosome 13 in FK506 binding protein 6 (FKBP6). The gene belongs to the immunophilins FKBP family known to be involved in meiosis, calcium homeostasis, clathrin-coated vesicles, and membrane fusions. Direct sequencing of FKBP6 exons in cases and controls identified SNPs g.11040315G>A and g.11040379C>A (p.166H>N) in exon 4 that were significantly associated with the IAR phenotype both in the GWAS cohort (n = 44) and in a large multi-breed cohort of 265 horses. All IAR stallions were homozygous for the A-alleles, while this genotype was found only in 2% of controls. The equine FKBP6 was exclusively expressed in testis and sperm and had 5 different transcripts, of which 4 were novel. The expression of this gene in AC/AG heterozygous controls was monoallelic, and we observed a tendency for FKBP6 up-regulation in IAR stallions compared to controls. Because exon 4 SNPs had no effect on the protein structure, it is likely that FKBP6 relates to the IAR phenotype via regulatory or modifying functions. In conclusion, FKBP6 was considered a susceptibility gene of incomplete penetrance for IAR in stallions and a candidate gene for male subfertility in mammals. FKBP6 genotyping is recommended for the detection of IAR-susceptible individuals among potential breeding stallions. Successful use of sperm as a source of DNA and RNA propagates non-invasive sample procurement for fertility genomics in animals and humans.

Expression level of a chimeric kinase governs entry into sporulation in Bacillus subtilis.

Upon starvation, Bacillus subtilis cells switch from growth to sporulation. It is believed that the N-terminal sensor domain of the cytoplasmic histidine kinase KinA is responsible for detection of the sporulation-specific signal(s) that appears to be produced only under starvation conditions. Following the sensing of the signal, KinA triggers autophosphorylation of the catalytic histidine residue in the C-terminal domain to transmit the phosphate moiety, via phosphorelay, to the master regulator for sporulation, Spo0A. However, there is no direct evidence to support the function of the sensor domain, because the specific signal(s) has never been found. To investigate the role of the N-terminal sensor domain, we replaced the endogenous three-PAS repeat in the N-terminal domain of KinA with a two-PAS repeat derived from Escherichia coli and examined the function of the resulting chimeric protein. Despite the introduction of a foreign domain, we found that the resulting chimeric protein, in a concentration-dependent manner, triggered sporulation by activating Spo0A through phosphorelay, irrespective of nutrient availability. Further, by using chemical cross-linking, we showed that the chimeric protein exists predominantly as a tetramer, mediated by the N-terminal domain, as was found for KinA. These results suggest that tetramer formation mediated by the N-terminal domain, regardless of the origin of the protein, is important and sufficient for the kinase activity catalyzed by the C-terminal domain. Taken together with our previous observations, we propose that the primary role of the N-terminal domain of KinA is to form a functional tetramer, but not for sensing an unknown signal.

Total RNA isolation from stallion sperm and testis biopsies.

Sperm mRNA transcriptional profiles can be used to evaluate male fertility, yet differences between species in sperm attributes and packaging require adjustments in sperm RNA isolation protocols. The objective was to optimize RNA isolation methodology for fresh, frozen, and extended ejaculates, and epididymal sperm of stallions. Additionally, a protocol for RNA isolation from testis biopsies was established. Separation of sperm from somatic cells was critical for assuring the isolation of sperm-specific RNA. The highest purity was obtained by centrifuging ejaculates and epididymal sperm at 200 x g for 30 min through a 40% Equipure silanized silica particle solution. Sperm RNA isolation was more efficient with TRIzol reagent than with a spin-column based method; it resulted in 2 microg of total RNA per 100 x 10(6) sperm. To evaluate RNA quantity and quality, we used a NanoDrop spectrophotometer and Agilent Bioanalyzer. A protocol for reverse transcriptase PCR with equine primers for PRM2 and PTPRC genes was developed to determine sperm RNA contamination with genomic DNA or RNA from somatic cells. By these methods, hybridization- and sequencing-quality RNA was isolated from 11 samples of stallion sperm. Stallion testis biopsy with a 14 gauge 22 mm deep biopsy needle yielded approximately 12 microg of good quality total RNA, and could serve as an alternative to excision surgery for sample procurement. Compared to RNA isolation from testis, the sperm required advanced processing and RNA quality control. The described methodologies provided a foundation to establish functional genomic studies of stallion fertility.