Comparative Gene Expression Analysis Research Articles

A comparative analysis of gene expression induced by the embryo in the caprine endometrium.

Transcriptomics is an established powerful tool to identify potential mRNAs and ncRNAs (non‐coding RNAs) for endometrial receptivity. In this study, the goat endometrium at estrus day 5 (ED5) and estrus day 15 (ED15) were selected to systematically analyse the differential expressed genes (DEGs) what were induced by the embryo. There were 1,847 genes which were significantly differential expressed in endometrium induced by the embryo at ED5, and 1,346 at ED15 (p‐value < .05). Secreted phosphoprotein 1 (SPP) was the responsive genes for embryo in the goat endometrium during estrus cycle, neurotensis (NTS) and pleiotrophin (PTN) were the responsive genes for embryo in the goat endometrium at ED5, Testin (TES) and Phosphate and Tension Homology Deleted on Chromsome ten (PTEN) at ED15. Furthermore, Gene Ontology (GO) and Kyoto Encyclopedia of Genes (KEGG) analysis revealed cytoplasm and Endocytosis were indispensable for the endometrium development in dairy goat. In a word, this resulting view of the transcriptome greatly uncovered the global trends in mRNAs expression induced by the embryo in the endometrium of dairy goats.

Identification of common candidate genes and pathways for progression of ovarian, cervical and endometrial cancers

Ovarian, cervical and endometrial cancers are the most common gynecological cancers affecting women worldwide. Overall, the molecular mechanism of the disease is not well understood. These cancers have similar embryological origin which has led to the hypothesis that they might share common genes and pathways for disease development. In order to find out pathways and key candidate genes shared by these three lethal gynecological cancers, we performed comparative gene expression analysis for the three types of cancers using microarray data. Differentially expressed genes (DEGs) were identified by the combination of |log2fc| > 1 and adjusted p (padj) value<0.05. Our analysis revealed a set of 426 genes in ovarian cancer (OC), 128 common genes in cervical cancer (CC), and 288 genes in endometrial cancer (EC) were differentially expressed. Gene ontology (GO) and pathway enrichment analysis of the DEGs showed that cell cycle and p53 signaling pathways were conserved across the three cancer types. We constructed co-expression network of the DEGs and identified a subset of eight common genes namely ASPM, CEP55, ECT2, KIF4A, MELK, RRM2, TOP2A, and TTK that were differentially co-expressed in all the three cancer types formed a core circuitry called key candidate genes. Identified candidate genes were validated by using gene expression profiling interactive analysis (GEPIA) tool and found the candidate genes were highly expressed in OC, CC and EC compared to the respective normal tissues. Involvement of the candidate genes in the cancer progression were further supported by literature survey that might provide new insights on the prognosis and treatment strategies for the three female-specific cancers.

Comparative gene expression analysis in melanocytes driven by tumor cell-derived exosomes

Abundant with organelle-like membranous structures, the tumor microenvironment is composed of cancer cells that secrete exosomes. Studies have shown that these secreted exosomes transport RNA and active molecules to other cells to reshape the tumor microenvironment and promote tumor growth. In fact, we found that exosomes derived from melanoma cells drive pre-malignant transition in primary melanocytes. However, there is little available in the scientific literature on how exosomes modulate melanocytes in the microenvironment to optimize conditions for tumor progression and metastasis. We therefore focused this current study on identifying these conditions genetically. Through RNA sequencing, we analyzed gene expression levels of melanocytes driven by exosomes derived from melanoma and lung cancer cells compared with those without exosome controls. Significant differences were found in gene expression patterns of melanocytes driven by exosomes derived from melanoma and lung cancer cells. In the melanocytes responding to exosomes derived from melanoma cells, genes of lipopolysaccharide and regulation of leukocyte chemotaxis were predominant. In the melanocytes responding to exosomes derived from lung cancer cells, genes of DNA replication and mitotic nuclear division played an important role. These results provide further mechanistic understanding of tumor progression promoted by tumor-derived exosomes. This will also help identify potential therapeutic targets for melanoma progression.

Comparative Tumor RNA Sequencing Analysis for Difficult-to-Treat Pediatric and Young Adult Patients With Cancer

Pediatric cancers are epigenetic diseases; therefore, considering tumor gene expression information is necessary for a complete understanding of the tumorigenic processes. To evaluate the feasibility and utility of incorporating comparative gene expression information into the precision medicine framework for difficult-to-treat pediatric and young adult patients with cancer. This cohort study was conducted as a consortium between the University of California, Santa Cruz (UCSC) Treehouse Childhood Cancer Initiative and clinical genomic trials. RNA sequencing (RNA-Seq) data were obtained from the following 4 clinical sites and analyzed at UCSC: British Columbia Children's Hospital (n = 31), Lucile Packard Children's Hospital at Stanford University (n = 80), CHOC Children's Hospital and Hyundai Cancer Institute (n = 46), and the Pacific Pediatric Neuro-Oncology Consortium (n = 24). The study dates were January 1, 2016, to March 22, 2017. Participants underwent tumor RNA-Seq profiling as part of 4 separate clinical trials at partner hospitals. The UCSC either downloaded RNA-Seq data from a partner institution for analysis in the cloud or provided a Docker pipeline that performed the same analysis at a partner institution. The UCSC then compared each participant's tumor RNA-Seq profile with more than 11 000 uniformly analyzed tumor profiles from pediatric and young adult patients with cancer, downloaded from public data repositories. These comparisons were used to identify genes and pathways that are significantly overexpressed in each patient's tumor. Results of the UCSC analysis were presented to clinical partners. Feasibility of a third-party institution (UCSC Treehouse Childhood Cancer Initiative) to obtain tumor RNA-Seq data from patients, conduct comparative analysis, and present analysis results to clinicians; and proportion of patients for whom comparative tumor gene expression analysis provided useful clinical and biological information. Among 144 samples from children and young adults (median age at diagnosis, 9 years; range, 0-26 years; 72 of 118 [61.0%] male [26 patients sex unknown]) with a relapsed, refractory, or rare cancer treated on precision medicine protocols, RNA-Seq-derived gene expression was potentially useful for 99 of 144 samples (68.8%) compared with DNA mutation information that was potentially useful for only 34 of 74 samples (45.9%). This study's findings suggest that tumor RNA-Seq comparisons may be feasible and highlight the potential clinical utility of incorporating such comparisons into the clinical genomic interpretation framework for difficult-to-treat pediatric and young adult patients with cancer. The study also highlights for the first time to date the potential clinical utility of harmonized publicly available genomic data sets.

Comparative proteomic analysis of cucumber powdery mildew resistance between a single-segment substitution line and its recurrent parent

Powdery mildew (PM) is considered a major cause of yield losses and reduced quality in cucumber worldwide, but the molecular basis of PM resistance remains poorly understood. A segment substitution line, namely, SSL508-28, was developed with dominant PM resistance in the genetic background of PM-susceptible cucumber inbred line D8. The substituted segment contains 860 genes. An iTRAQ-based comparative proteomic technology was used to map the proteomes of PM-inoculated and untreated (control) D8 and SSL508-28. The number of differentially regulated proteins (DRPs) in SSL508-28 was almost three times higher than that in D8. Fourteen DRPs were located in the substituted segment interval. Comparative gene expression analysis revealed that nodulin-related protein 1 (NRP1) may be a good candidate for PM resistance. Gene Ontology enrichment analysis showed that DRPs functioning in tetrapyrrole biosynthetic process, sulfur metabolic process and cell redox homeostasis were specifically enriched in the resistant line SSL508-28. DRPs categorized in the KEGG term photosynthesis increased in both lines upon PM infection, suggesting that the strategies used by cucumber may be different from those used by other crops to react to PM attacks at the initial stage. The measurement of hydrogen peroxide and superoxide anion production and net photosynthetic rate were consistent with the changes in protein abundance, suggesting that the proteomic results were reliable. There was a poor correlation between DRPs measured by iTRAQ and the corresponding gene expression changes measured by RNA-seq with the same experimental design. Taken together, these findings improve the understanding of the molecular mechanisms underlying the response of cucumber to PM infection.

Identification of a candidate gene underlying qKRN5b for kernel row number in Zea mays L.

A quantitative trait locus for kernel row number, qKRN5, was dissected into two tightly linked loci, qKRN5a and qKRN5b. Fine mapping, comparative analysis of nucleotide sequences and gene expression established the endonuclease/exonuclease/phosphatase family protein-encoding gene Zm00001d013603 as a causal gene of qKRN5b. Maize grain yield is determined by agronomically important traits that are controlled by interactions among and between genes and environmental factors. Considerable efforts have been made to identify major quantitative trait loci (QTLs) for yield-related traits; however, few were previously isolated and characterized in maize. In this study, we divided a QTL for kernel row number (KRN), qKRN5, into two tightly linked loci, qKRN5a and qKRN5b, using advanced backcross populations derived from near-isogenic lines. KRN was greater in individuals that were homozygous for the NX531 allele, which showed coupling-phase linkage. The major QTL qKRN5b had an additive effect of approximately one kernel row. Furthermore, fine mapping narrowed qKRN5b within a 147.2-kb region. The upstream sequence Zm00001d013603 and its expression in the ear inflorescence showed obvious differences between qKRN5b near-isogenic lines. In situ hybridization located Zm00001d013603 on the primordia of the spikelet pair meristems and spikelet meristems, but not in the inflorescence meristem, which indicates a role in regulating the initiation of reproductive axillary meristems of ear inflorescences. Expression analysis and nucleotide sequence alignment revealed that Zm00001d013603, which encodes an endonuclease/exonuclease/phosphatase family protein that hydrolyzes phosphatidyl inositol diphosphates, is the causal gene of qKRN5b. These results provide insight into the genetic basis of KRN and have potential value for enhancing maize grain yield.

Swarming Inhibitory Potential of Cinnamtannin B1 from Cinnamomum tamala T. Nees and Eberm on Pseudomonas aeruginosa.

In a preliminary screening, the methanol extract of Cinnamomum tamala leaves was found to inhibit the swarming motility of Pseudomonas aeruginosa. Bioassay-guided fractionation by silica gel column chromatography led to the identification of cinnamtannin B1 (1) as one of the active components of the extract. It inhibited the swarming motility (at 12.5 μg/mL) and biofilm formation (at 25 μg/mL) ofP. aeruginosa. Comparative gene expression analysis revealed downregulation of rhlA and fliC genes upon treatment with the tannin. The tannin may be affecting rhamnolipid and flagellin production. Thus, cinnamtannin B1 is an active component of C. tamala responsible for inhibiting the swarming motility of P. aeruginosa.

Evolution of the visual sensory system in cichlid fishes from crater lake Barombi Mbo in Cameroon.

In deep-water animals, the visual sensory system is often challenged by the dim-light environment. Here, we focus on the molecular mechanisms involved in rapid deep-water adaptations. We examined visual system evolution in a small-scale yet phenotypically and ecologically diverse adaptive radiation, the species flock of cichlid fishes in deep crater lake Barombi Mbo in Cameroon, West Africa. We show that rapid adaptations of the visual system to the novel deep-water habitat primarily occurred at the level of gene expression changes rather than through nucleotide mutations, which is compatible with the young age of the radiation. Based on retinal bulk RNA sequencing of all eleven species, we found that the opsin gene expression pattern was substantially different for the deep-water species. The nine shallow-water species feature an opsin palette dominated by the red-sensitive (LWS) opsin, whereas the two unrelated deep-water species lack expression of LWS and the violet-sensitive (SWS2B) opsin, thereby shifting the cone sensitivity to the centre of the light spectrum. Deep-water species further predominantly express the green-sensitive RH2Aα over RH2Aβ. We identified one amino acid substitution in the RH2Aα opsin specific to the deep-water species. We finally performed a comparative gene expression analysis in retinal tissue of deep- vs. shallow-water species. We thus identified 46 differentially expressed genes, many of which are associated with functions in vision, hypoxia management or circadian clock regulation, with some of them being associated with human eye diseases.

Evaluation of suitable reference genes in Brassica juncea and its wild relative Camelina sativa for qRT-PCR analysis under various stress conditions

Quantitative real-time PCR (qRT-PCR) is an efficient method to estimate the gene expression levels but the accuracy of its result largely depends on the stability of the reference gene. Many studies have reported considerable variation in the expression of reference genes (RGs) in different tissue and conditions. Therefore, screening for appropriate RGs with stable expression is crucial for functional analysis of the target gene. Two closely related crucifers Brassica juncea (cultivated) and Camelina sativa (wild) respond differently towards various abiotic and biotic stress where C. sativa exhibits higher tolerance to various stress. Comparative gene expression analysis of the target genes between two such species is the key approach to understand the mechanism of a plant’s response to stress. However, using an unsuitable RG can lead to misinterpretation of expression levels of the target gene in such studies. In this investigation, the stability of seven candidate RGs including traditional housekeeping genes (HKGs) and novel candidate RGs were identified across diverse sample sets of B. juncea and C. sativa representing- hormone treated, wounded, Alternaria brassicae inoculated and combination treated samples (exogenous hormone treatment followed by A. brassicae inoculation). In this investigation, we identified stable RGs in both the species and the most suitable RGs to perform an unbiased comparative gene expression analysis between B. juncea and C. sativa. Results revealed that TIPS41 and PP2A were identified as the overall best performing RGs in both the species. However, the most suitable RG for each sample subset representing different condition must be individually selected. In Hormone treated and wounded samples TIPS41 expressed stably in both the species and in A. brassicae inoculated and combination treatment performance of PP2A was the best. In this study, for the first time, we have identified and validated stable reference gene in C. sativa for accurate normalization of gene expression data.

Comparative gene expression analysis of the engulfment and cell motility (ELMO) protein family in the mouse brain

Engulfment and cell motility (ELMO) proteins bind to Dock180, a guanine nucleotide exchange factor (GEF) of the Rac family, and regulate GEF activity. The resultant ELMO/Dock180/Rac module regulates cytoskeletal reorganization responsible for the engulfment of apoptotic cells, cell migration, and neurite extension. The expression and function of Elmo family proteins in the nervous system, however, are not yet fully understood. Here, we characterize the comparative gene expression profiles of three Elmo family members (Elmo1, Elmo2, and Elmo3) in the brain of C57BL/6J mice, a widely used inbred strain, together with reeler mutant mice to understand gene expression in normal laminated brain areas compared with abnormal areas. Although all three Elmo genes showed widespread mRNA expression over various mouse tissues tested, Elmo1 and Elmo2 were the major types expressed in the brain, and three Elmo genes were up-regulated between the first postnatal week (infant stage) and the third postnatal week (juvenile, weaning stage). In addition, the mRNAs of Elmo genes showed distinct distribution patterns in various brain areas and cell-types; such as neurons including inhibitory interneurons as well as some non-neuronal cells. In the cerebral cortex, the three Elmo genes were widely expressed over many cortical regions, but the predominant areas of Elmo1 and Elmo2 expression tended to be distributed unevenly in the deep (a lower part of the VI) and superficial (II/III) layers, respectively, which also changed depending on the cortical areas and postnatal stages. In the dentate gyrus of the hippocampus, Elmo2 was expressed in dentate granule cells more in the mature stage rather than the immature-differentiating stage. In the thalamus, Elmo1 but not the other members was highly expressed in many nuclei. In the medial habenula, Elmo2 and Elmo3 were expressed at intermediate levels. In the cerebellar cortex, Elmo1 and Elmo2 were expressed in differentiating-mature granule cells and mature granule cells, respectively. In the Purkinje cell layer, Elmo1 and Elmo2 were expressed in Purkinje cells and Bergmann glia, respectively. Disturbed cellular distributions and laminar structures caused by the reeler mutation did not severely change expression in these cell types despite the disturbed cellular distributions and laminar structures, including those of the cerebrum, hippocampus, and cerebellum. Taken together, these results suggested that these three Elmo family members share their functional roles in various brain regions during prenatal-postnatal development.

ASSOCIATION OF EXPRESSION OF NFKB1, HIF1, VEGFA, VEGFВ, BAX, BCL2 GENES WITH BIOCHEMICAL REMEDIATION IN PATIENTS WITH LOCALIZED PROSTATE CANCER

Aim. To identify the association of NFKB1, HIF1, VEGFA, VEGFB, BAX, BCL2 gene expression in prostate adenocarcinoma cells with biochemical recurrence of localized prostate cancer. Patients and methods. Three groups of patients were formed in the study – the main one, the comparison group and the control group. In patients with prostate cancer (PC) in the main group (n = 56) with biochemical recurrence (BR) for two years after radical surgery, as well as in 60 patients without BR (experimental group) by real-time PCR in prostate cancer tissue the expression of genes NFKB1, HIF1, VEGFA, VEGFВ, BAX, BCL2 was determined. The control group consisted of 55 patients in whom, when performing diagnostic punctures for benign prostate tumors, biopsy specimens were taken in healthy tissues. The age of patients in the three groups ranged from 57 to 74 years (median 63 years). When quantifying expression of genes NFKB1, HIF1, VEGFA, VEGFВ, BAX, BCL2, the difference in the values of reaction threshold cycles (Ct) fixed for the studied and reference genes was determined. The relative level (Expr) was the ratio of Ct medians for each gene in two compared groups of the studied three ones: in the main group to the indicator in the control group, in the experimental group to the indicator in the control group, and also between the main group and the experimental group. Results. A comparative analysis of gene expression in prostate cancer tissue in the main group compared with the experimental group showed a statistically significant increase (p &lt; 0,05) in the relative index for the HIF1 gene (2,7 times), the VEGFA gene (2,4 times ) and the NFKB1 gene (2 times). Consequently, in patients with localized early recurrence prostate cancer, initially in the prostate tissue, a higher level of expression of the NFKB1, HIF1 and VEGFA genes was established. In the experimental group relative to the control group, the expression of the proapoptic gene BAX was 1,6 times higher (p &lt; 0,05), and for the antiapoptic gene BCL2 no changes were detected (p = 0,09). Thus, in patients with localized prostate cancer in the absence of BR, after radical prostatectomy, an initial increase in the expression of the BAX gene promoted the activation of apoptosis. In patients with localized prostate cancer, subsequent biochemical recurrence initially in the tissue of prostate adenocarcinoma inhibition of apoptosis due to increased expression of the BCL2 gene was observed. Conclusion. Enhancement of NFKB1, VEGFA, HIF1 and BCL2 gene expression in prostate tissue is associated with the development of BR in patients with localized prostate cancer.

Molecular characterization of the human lens epithelium-derived cell line SRA01/04

Cataract-associated gene discovery in human and animal models have informed on key aspects of human lens development, homeostasis and pathology. Additionally, in vitro models such as the culture of permanent human lens epithelium-derived cell lines (LECs) have also been utilized to understand the molecular biology of lens cells. However, these resources remain uncharacterized, specifically regarding their global gene expression and suitability to model lens cell biology. Therefore, we sought to molecularly characterize gene expression in the human LEC, SRA01/04, which is commonly used in lens studies. We first performed short tandem repeat (STR) analysis and validated SRA01/04 LEC for its human origin, as recommended by the eye research community. Next, we used Illumina HumanHT-12 v3.0 Expression BeadChip arrays to gain insights into the global gene expression profile of SRA01/04. Comparative analysis of SRA01/04 microarray data was performed using other resources such as the lens expression database iSyTE (integrated Systems Tool for Eye gene discovery), the cataract gene database Cat-Map and the published lens literature. This analysis showed that SRA01/04 significantly expresses >40% of the top iSyTE lens-enriched genes (313 out of 749) across different developmental stages. Further, SRA01/04 also significantly expresses ~53% (168 out of 318) of cataract-associated genes in Cat-Map. We also performed comparative gene expression analysis between SRA01/04 cells and the previously validated mouse LEC 21EM15. To gain insight into whether SRA01/04 reflects epithelial or fiber cell characteristics, we compared its gene expression profile to previously reported differentially expressed genes in isolated mouse lens epithelial and fiber cells. This analysis suggests that SRA01/04 has reduced expression of several fiber cell-enriched genes. In agreement with these findings, cell culture analysis demonstrates that SRA01/04 has reduced potential to initiate spontaneous lentoid body formation compared to 21EM15 cells. Next, to independently validate SRA01/04 microarray gene expression, we subjected several candidate genes to RT-PCR and RT-qPCR assays. This analysis demonstrates that SRA01/04 supports expression of many key genes associated with lens development and cataract, including CRYAB, CRYBB2, CRYGS, DKK3, EPHA2, ETV5, GJA1, HSPB1, INPPL1, ITGB1, PAX6, PVRL3, SFRP1, SPARC, TDRD7, and VIM, among others, and therefore can be relevant for understanding the mechanistic basis of these factors. At the same time, SRA01/04 cells do not exhibit robust expression of several genes known to be important to lens biology and cataract such as ALDH1A1, COL4A6, CP, CRYBA4, FOXE3, HMX1, HSF4, MAF, MEIS1, PITX3, PRX, SIX3, and TRPM3, among many others. Therefore, the present study offers a rich transcript-level resource for case-by-case evaluation of the potential advantages and limitations of SRA01/04 cells prior to their use in downstream investigations. In sum, these data show that the human LEC, SRA01/04, exhibits lens epithelial cell-like character reflected in the expression of several lens-enriched and cataract-associated genes, and therefore can be considered as a useful in vitro resource when combined with in vivo studies to gain insight into specific aspects of human lens epithelial cells.

Comparative Analysis of Gene Expression in Virulent and Attenuated Strains of Infectious Bronchitis Virus at Subcodon Resolution.

Like all coronaviruses, avian infectious bronchitis virus (IBV) possesses a long, single-stranded, positive-sense RNA genome (∼27 kb) and has a complex replication strategy that includes the production of a nested set of subgenomic mRNAs (sgmRNAs). Here, we used whole-transcriptome sequencing (RNASeq) and ribosome profiling (RiboSeq) to delineate gene expression in the IBV M41-CK and Beau-R strains at subcodon resolution. RNASeq facilitated a comparative analysis of viral RNA synthesis and revealed two novel transcription junction sites in the attenuated Beau-R strain, one of which would generate a sgmRNA encoding a ribosomally occupied open reading frame (dORF) located downstream of the nucleocapsid coding region. RiboSeq permitted quantification of the translational efficiency of virus gene expression and identified, for the first time, sites of ribosomal pausing on the genome. Quantification of reads flanking the programmed ribosomal frameshifting (PRF) signal at the genomic RNA ORF1a/ORF1b junction revealed that PRF in IBV is highly efficient (33 to 40%). Triplet phasing of RiboSeq data allowed precise determination of reading frames and revealed the translation of two ORFs (ORF4b and ORF4c on sgmRNA IR), which are widely conserved across IBV isolates. Analysis of differential gene expression in infected primary chick kidney cells indicated that the host cell response to IBV occurs primarily at the level of transcription, with global upregulation of immune-related mRNA transcripts following infection and comparatively modest changes in the translation efficiencies of host genes. Cellular genes and gene networks differentially expressed during virus infection were also identified, giving insights into the host cell response to IBV infection.IMPORTANCE IBV is a major avian pathogen and presents a substantial economic burden to the poultry industry. Improved vaccination strategies are urgently needed to curb the global spread of this virus, and the development of suitable vaccine candidates will be aided by an improved understanding of IBV molecular biology. Our high-resolution data have enabled a precise study of transcription and translation in cells infected with both pathogenic and attenuated forms of IBV and expand our understanding of gammacoronaviral gene expression. We demonstrate that gene expression shows considerable intraspecies variation, with single nucleotide polymorphisms being associated with altered production of sgmRNA transcripts, and our RiboSeq data sets enabled us to uncover novel ribosomally occupied ORFs in both strains. The numerous cellular genes and gene networks found to be differentially expressed during virus infection provide insights into the host cell response to IBV infection.

Exploring the molecular basis of neuronal excitability in a vocal learner

BackgroundVocal learning, the ability to learn to produce vocalizations through imitation, relies on specialized brain circuitry known in songbirds as the song system. While the connectivity and various physiological properties of this system have been characterized, the molecular genetic basis of neuronal excitability in song nuclei remains understudied. We have focused our efforts on examining voltage-gated ion channels to gain insight into electrophysiological and functional features of vocal nuclei. A previous investigation of potassium channel genes in zebra finches (Taeniopygia guttata) revealed evolutionary modifications unique to songbirds, as well as transcriptional specializations in the song system [Lovell PV, Carleton JB, Mello CV. BMC Genomics 14:470 2013]. Here, we expand this approach to sodium, calcium, and chloride channels along with their modulatory subunits using comparative genomics and gene expression analysis encompassing microarrays and in situ hybridization.ResultsWe found 23 sodium, 38 calcium, and 33 chloride channel genes (HGNC-based classification) in the zebra finch genome, several of which were previously unannotated. We determined 15 genes are missing relative to mammals, including several genes (CLCAs, BEST2) linked to olfactory transduction. The majority of sodium and calcium but few chloride channels showed differential expression in the song system, among them SCN8A and CACNA1E in the direct motor pathway, and CACNG4 and RYR2 in the anterior forebrain pathway. In several cases, we noted a seemingly coordinated pattern across multiple nuclei (SCN1B, SCN3B, SCN4B, CACNB4) or sparse expression (SCN1A, CACNG5, CACNA1B).ConclusionThe gene families examined are highly conserved between avian and mammalian lineages. Several cases of differential expression likely support high-frequency and burst firing in specific song nuclei, whereas cases of sparse patterns of expression may contribute to the unique electrophysiological signatures of distinct cell populations. These observations lay the groundwork for manipulations to determine how ion channels contribute to the neuronal excitability properties of vocal learning systems.

Familial Left Ventricular Non-Compaction Is Associated With a Rare p.V407I Variant in Bone Morphogenetic Protein 10.

Left ventricular non-compaction (LVNC) is a heritable cardiomyopathy characterized by hypertrabeculation, inter-trabecular recesses and thin compact myocardium, but the genetic basis and mechanisms remain unclear. This study identified novel LVNC-associated mutations inNOTCH-dependent genes and investigated their mutational effects.Methods and Results:High-resolution melting screening was performed in 230 individuals with LVNC, followed by whole exome and Sanger sequencing of available family members. Dimerization of bone morphogenetic protein 10 (BMP10) and its binding to BMP receptors (BMPRs) were evaluated. Cellular differentiation, proliferation and tolerance to mechanical stretch were assessed in H9C2 cardiomyoblasts, expressing wild-type (WT) or mutant BMP10 delivered by adenoviral vectors. Rare variants, p.W143*-NRG1and p.V407I-BMP10, were identified in 2 unrelated probands and their affected family members. Although dimerization of mutant V407I-BMP10 was preserved like WT-BMP10, V407I-BMP10 pulled BMPR1a and BMPR2 receptors more weakly compared with WT-BMP10. On comparative gene expression and siRNA analysis, expressed BMPR1a and BMPR2 receptors were responsive to BMP10 treatment in H9C2 cardiomyoblasts. Expression of V407I-BMP10 resulted in a significantly lower rate of proliferation in H9C2 cells compared with WT-BMP10. Cyclic stretch resulted in destruction and death of V407I-BMP10 cells. The W143*-NRG1and V470I-BMP10variants are associated with LVNC. Impaired BMPR-binding ability, perturbed proliferation and differentiation processes and intolerance to stretch in V407I-BMP10 mutant cardiomyoblasts may underlie myocardial non-compaction.

Human bone marrow adipocytes display distinct immune regulatory properties

BackgroundThe bone marrow (BM) is a major reservoir of resting memory T cells and long-lived plasma cells, capable of providing protection against recurrent infections. Whether the age-related accumulation of adipose tissue in the BM affects the functionality and maintenance of memory cells is not well understood. MethodsFor the first time, we compare human femur marrow adipose tissue (fMAT) and subcutaneous white adipose tissue of the thigh (tsWAT) obtained from the same donors. Therefore, we used microarrays for comparative global gene expression analysis, and employed assays to analyse parameters of adipocyte biology, inflammation and oxidative stress. FindingsWe show that fMAT adipocytes differ significantly from tsWAT adipocytes regarding specific gene expression profiles including inflammatory responses and adipogenesis/adipocyte phenotype. Concomitant with considerably lower levels of CD36, a membrane-associated protein important for long-chain fatty acid uptake that is used as maturation marker, fMAT adipocytes are smaller and contain less triglycerides. fMAT adipocytes secrete similar levels of adiponectin and leptin as tsWAT adipocytes, and express increased levels of pro-inflammatory molecules concomitant with an elevated generation of reactive oxygen species (ROS) and impaired function of plasma cells in the BM. InterpretationOur findings suggest that fMAT is a unique type of adipose tissue containing small adipocytes with lower CD36 protein and triglyceride levels than tsWAT but high adipokine secretion. Moreover, fMAT adipocytes secrete high levels of pro-inflammatory cytokines, contributing to inflammation and impairment of plasma cell function in the BM, suggesting that fMAT has more immune regulatory functions than tsWAT.

Comparative proteomics and gene expression analyses revealed responsive proteins and mechanisms for salt tolerance in chickpea genotypes

BackgroundSalinity is a major abiotic stress that limits the growth, productivity, and geographical distribution of plants. A comparative proteomics and gene expression analysis was performed to better understand salinity tolerance mechanisms in chickpea.ResultsTen days of NaCl treatments resulted in the differential expression of 364 reproducible spots in seedlings of two contrasting chickpea genotypes, Flip 97-43c (salt tolerant, T1) and Flip 97-196c (salt susceptible, S1). Notably, after 3 days of salinity, 80% of the identified proteins in T1 were upregulated, while only 41% in S2 had higher expression than the controls. The proteins were classified into eight functional categories, and three groups of co-expression profile. The second co-expressed group of proteins had higher and/or stable expression in T1, relative to S2, suggesting coordinated regulation and the importance of some processes involved in salinity acclimation. This group was mainly enriched in proteins associated with photosynthesis (39%; viz. chlorophyll a-b binding protein, oxygen-evolving enhancer protein, ATP synthase, RuBisCO subunits, carbonic anhydrase, and fructose-bisphosphate aldolase), stress responsiveness (21%; viz. heat shock 70 kDa protein, 20 kDa chaperonin, LEA-2 and ascorbate peroxidase), and protein synthesis and degradation (14%; viz. zinc metalloprotease FTSH 2 and elongation factor Tu). Thus, the levels and/or early and late responses in the activation of targeted proteins explained the variation in salinity tolerance between genotypes. Furthermore, T1 recorded more correlations between the targeted transcripts and their corresponding protein expression profiles than S2.ConclusionsThis study provides insight into the proteomic basis of a salt-tolerance mechanism in chickpea, and offers unexpected and poorly understood molecular resources as reliable starting points for further dissection.

Reduced expression of apolipoprotein E and immunoglobulin heavy constant gamma 1 proteins in Fuchs endothelial corneal dystrophy.

Fuchs endothelial corneal dystrophy (FECD) is a progressive and potentially a sight threatening disease, and a common indication for corneal grafting in the elderly. Aberrant thickening of Descemet's membrane, formation of microscopic excrescences (guttae) and gradual loss of corneal endothelial cells are the hallmarks of the disease. The aim of this study was to identify differentially abundant proteins between FECD-affected and unaffected Descemet's membrane. Label-free quantitative proteomics using nanoscale ultra-performance liquid chromatography-mass spectrometry (nUPLC-MSE ) was employed on affected and unaffected Descemet's membrane extracts, and interesting findings were further investigated using quantitative reverse transcription-polymerase chain reaction and immunohistochemical techniques. Quantitative proteomics revealed significantly lower abundance of apolipoprotein E (APOE) and immunoglobulin heavy constant gamma 1 protein (IGHG1) in affected Descemet's membrane. The difference in the distribution of APOE between affected and unaffected Descemet's membrane and of IGHG1 detected by immunohistochemistry support their down-regulation in the disease. Comparative gene expression analysis showed significantly lower APOE mRNA levels in FECD-affected than unaffected corneal endothelium. IGHG1 gene is expressed at extremely low levels in the corneal endothelium, precluding relative expression analysis. This is the first study to report comparative proteomics of Descemet's membrane tissue, and implicates dysregulation of APOE and IGHG1 proteins in the pathogenesis of Fuchs endothelial corneal dystrophy.

Transcriptome Analysis of Chrysanthemum lavandulifolium Response to Salt Stress and Overexpression a K+ Transport ClAKT Gene-enhanced Salt Tolerance in Transgenic Arabidopsis

Plant growth and development are significantly affected by salt stress. Chrysanthemum lavandulifolium is a halophyte species and one of the ancestors of chrysanthemum (C. ×morifolium). Understanding how this species tolerates salt stress could provide vital insight for clarifying the salt response systems of higher plants, and chrysanthemum-breeding programs could be improved. In this study, salt tolerance was compared among C. lavandulifolium and three chrysanthemum cultivars by physiological experiments, among which C. lavandulifolium and Jinba displayed better tolerance to salt stress than the other two cultivars, whereas Xueshan was a salt-sensitive cultivar. Using the transcriptome database of C. lavandulifolium as a reference, we used digital gene expression technology to analyze the global gene expression changes in C. lavandulifolium seedlings treated with 200 mm NaCl for 12 hours compared with seedlings cultured in normal conditions. In total, 2254 differentially expressed genes (DEGs), including 1418 up-regulated and 836 down-regulated genes, were identified. These DEGs were significantly enriched in 35 gene ontology terms and 29 Kyoto Encyclopedia of Genes and Genomes pathways. Genes related to signal transduction, ion transport, proline biosynthesis, reactive oxygen species scavenging systems, and flavonoid biosynthesis pathways were relevant to the salt tolerance of C. lavandulifolium. Furthermore, comparative gene expression analysis was conducted using reverse transcription polymerase chain reaction to compare the transcriptional levels of significantly up-regulated DEGs in C. lavandulifolium and the salt-sensitive cultivar Xueshan, and species-specific differences were observed. The analysis of one of the DEGs, ClAKT, an important K+ transport gene, was found to enable transgenic Arabidopsis thaliana to absorb K+ and efflux Na+ under salt stress and to absorb K+ under drought stress. The present study investigated potential genes and pathways involved in salt tolerance in C. lavandulifolium and provided a hereditary resource for the confinement of genes and pathways responsible for salt tolerance in this species. This study provided a valuable source of reference genes for chrysanthemum cultivar transgenesis breeding.

Abstract 4689: Subclone-specific evolution of tumor phenotypes – A framework to study subclone-specific gene expression from a combination of bulk DNA and single cell RNA sequencing data

Abstract Several approaches are now available for subclonal reconstruction of heterogeneous tumor biopsies from somatic variant allele frequencies measured in bulk DNA sequencing datasets, including our own SubcloneSeeker program. However, to understand how chemoresistance, relapse, or metastasis evolves at a subclonal level, we also need to study the molecular phenotype corresponding to each subclone. Single-cell RNAseq technologies now allow one to study the transcriptomic characteristics and phenotypic behaviors of individual cells, but new algorithms are needed to link these cells to the genomically defined tumor subclones they represent. Here we present a computational approach to make such cell assignments, using a combination of bulk DNA sequencing data (WGS or WES), and single cell RNA sequencing (scRNAseq) collected from the same tissues. Subclones constructed from bulk DNA sequencing data are defined by specific combinations of somatic mutations (SNVs and CNVs). Our algorithm uses the scRNA-seq reads to assess which of these subclone-defining mutations are present in each individual cell, then utilizes a Bayesian probabilistic framework for making the cell-to-subclone assignment. This framework overcomes the challenges of sparse scRNA-seq read coverage, and maximizes the accuracy and efficiency of cell assignment. We have successfully applied this method to multiple longitudinal and metastatic cancer patient datasets, representing both WES and WGS bulk DNA sequencing, as well as datasets collected using the 10X Chromium and Fluidigm C1 technologies. Our algorithm uses both somatically acquired CNV and SNV events in the tumor for cell assignment. Using our approach, as many of 80% of tumor cells could be assigned to specific subclones, enabling comparative gene expression and pathway analysis across subclones. Citation Format: Yi Qiao, Xiaomeng Huang, Samuel Brady, Andrea Bild, David Bowtell, William Johnson, Gabor Marth. Subclone-specific evolution of tumor phenotypes – A framework to study subclone-specific gene expression from a combination of bulk DNA and single cell RNA sequencing data [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2019; 2019 Mar 29-Apr 3; Atlanta, GA. Philadelphia (PA): AACR; Cancer Res 2019;79(13 Suppl):Abstract nr 4689.