Gene Tissue Expression Research Articles

Komrower Memorial Lecture 2023. Molecular basis of phenotype expression in homocystinuria: Where are we 30 years later?

This review summarises progress in the research of homocystinuria (HCU) in the past three decades. HCU due to cystathionine β-synthase (CBS) was discovered in 1962, and Prof. Jan Peter Kraus summarised developments in the field in the first-ever Komrower lecture in 1993. In the past three decades, significant advancements have been achieved in the biology of CBS, including gene organisation, tissue expression, 3D structures, and regulatory mechanisms. Renewed interest in CBS arose in the late 1990s when this enzyme was implicated in biogenesis of H2S. Advancements in genetic and biochemical techniques enabled the identification of several hundreds of pathogenic CBS variants and the misfolding of missense mutations as a common mechanism. Several cellular, invertebrate and murine HCU models allowed us to gain insights into functional and metabolic pathophysiology of the disease. Establishing the E-HOD consortium and patient networks, HCU Network Australia and HCU Network America, offered new possibilities for acquiring clinical data in registries and data on patients' quality of life. A recent analysis of data from the E-HOD registry showed that the clinical variability of HCU is broad, extending from severe childhood disease to milder (late) adulthood forms, which typically respond to pyridoxine. Pyridoxine responsiveness appears to be the key factor determining the clinical course of HCU. Increased awareness about HCU played a role in developing novel therapies, such as gene therapy, correction of misfolding by chaperones, removal of methionine from the gut and enzyme therapies that decrease homocysteine or methionine in the circulation.

Co-expression analysis of transcriptomic data from cancer and healthy specimens reveals rewiring of proteasome genes and an interaction with the XPO1 gene across several tumour types

BackgroundThe 26S proteasome is a large intracellular multiprotein complex, that plays a homeostatic role by degrading proteins that have been tagged by ubiquitin. It is composed of 64 subunits assembled according to a well-defined structure and stoichiometry. Several proteasome subunits have been found to be overexpressed in tumours. However, comprehensive data are lacking on the relative abundance of each subunit and the impact on proteasome composition or stoichiometry. In cancer treatment, proteasome inhibitors and inhibitors of XPO1 (Exportin-1) have unexpectedly a similar range of activity, but the interaction between the two pathways has not been studied.MethodsWe performed gene co-expression analysis of 38 genes encoding proteasome subunits and 38 genes encoding proteins involved in nucleocytoplasmic transport in specimens from the Cancer Genome Atlas (33 tumour types) and from the Gene Tissue Expression database (32 healthy tissue types). We obtained 65 matrices, each containing Pearson correlation factors for 2964 gene pairs. We applied cluster analysis to the correlation matrices and compared the distribution of Pearson correlation coefficients of thirteen tumour types with their healthy tissue counterpart.ResultsStrong positive correlation (R Pearson correlation > 0.8) was observed for pairs of proteasome genes in the majority of healthy tissues, whereas the correlation for co-expression was significantly lower (R ≤ 0.50) for most gene pairs in the majority of cancer types. Cluster analysis based on gene co-expression allowed to distinguish cancers from healthy tissues in a clear-cut manner, and to identify the genes that contributed most to the separation. The crossed analysis between proteasome and nucleocytoplasmic transport genes showed that the expression of XPO1 and a subset of proteasome genes, including in particular PSMD14, is correlated in several cancer types and not in their healthy counterpart.ConclusionsThis analysis reveals that in cancer the co-expression of proteasome genes is significantly altered, highlighting the genes that are more often deregulated. In addition, it finds that XPO1 expression is often correlated with the expression of proteasome genes. From a therapeutic perspective, these findings support the investigation of novel targets within the proteasome and strategies of co-targeting of the proteasome and nucleocytoplasmic transport.

Genome-wide association study of traumatic brain injury in U.S. military veterans enrolled in the VA million veteran program.

Large-scale genetic studies of traumatic brain injury (TBI) are lacking; thus, our understanding of the influence of genetic factors on TBI risk and recovery is incomplete. This study aimed to conduct a genome-wide association study (GWAS) of TBI in VA Million Veteran Program (MVP) enrollees. Participants included a multi-ancestry cohort (European, African, and Hispanic ancestries; N = 304,485; 111,494 TBI cases, 192,991 controls). TBI was assessed using MVP survey data and International Classification of Diseases (ICD) codes from the Veterans Health Administration's electronic health record. GWAS was performed using logistic regression in PLINK, and meta-analyzed in METAL. FUMA was used for post-GWAS analysis. Genomic structural equation modeling (gSEM) was conducted to investigate underlying genetic associations with TBI, and bivariate MiXeR was used to estimate phenotype specific and shared polygenicity. SNP-based heritability was 0.060 (SE = 0.004, p = 7.83×10-66). GWAS analysis identified 15 genome-wide significant (GWS) loci at p < 5×10-8. Gene-based analyses revealed 14 gene-wide significant genes; top genes included NCAM1, APOE, FTO, and FOXP2. Gene tissue expression analysis identified the brain as significantly enriched, particularly in the frontal cortex, anterior cingulate cortex, and nucleus accumbens. Genetic correlations with TBI were significant for risk-taking behaviors and psychiatric disorders, but generally not significant for the neurocognitive variables investigated. gSEM analysis revealed stronger associations with risk-taking traits than with psychiatric traits. Finally, the genetic architecture of TBI was similar to polygenic psychiatric disorders. Neurodegenerative disorders including Alzheimer's and Parkinson's disease showed much less polygenicity, however, the proportion of shared variance with TBI was high. This first well-powered GWAS of TBI identified 15 loci including genes relevant to TBI biology, and showed that TBI is a heritable trait with comparable genetic architecture and high genetic correlation with psychiatric traits. Our findings set the stage for future TBI GWASs that focus on injury severity and diversity and chronicity of symptom sequelae.

Screening of reliable reference genes for the normalization of RT-qPCR in chicken gastrointestinal tract

The application of reverse transcription quantitative real-time PCR technology for the production of gene tissue expression profiles is a widely employed approach in molecular biology research. It is imperative to ascertain internal reference genes that exhibit stable expression across diverse tissues to ensure the precision of tissue gene expression profiles. While there have been studies documenting the most suitable reference genes for various tissues in chickens, there is a dearth of research on the identification of reference genes in the gastrointestinal (GI) tract of chickens. This study utilized 4 different algorithms (Delta CT, BestKeeper, NormFinder, and Genorm) to assess the stability of 19 internal reference genes in various GI tract tissues, including individual GI tract tissues, the anterior and posterior GI tract, and the entire GI tissue. The RefFinder software was employed to comprehensively rank these genes. The research findings successfully identified the most appropriate internal reference genes for each type of GI tissue. Furthermore, TBP, DNAJC24, Polr2b, RPL13, andAp2m exhibited stable expression in the entire and posterior GI tract, whereas HMBS, TBP, Ap2m, GUSB, DNAJC24, and RPL13 demonstrated stable expression in the anterior GI tract. However, the internal reference genes commonly utilized, namely β-Actin, 18s RNA, and ALB, exhibit poor stability and are not advised for future investigations concerning gene expression in the GI region. Consequently, MUC2 and CDX1, 2 genes that specifically express in the gut, were chosen for examination to ascertain the stability of the aforementioned internal reference genes in this particular study. In summary, this study presents a relatively stable set of internal reference genes that can be employed to enhance the precision of quantifying mRNA expression levels in functional genes within the chicken GI tract.

A - 7 Genome-Wide Association Study of Traumatic Brain Injury in U.S. Military Veterans.

Large-scale genetic studies of traumatic brain injury (TBI) are lacking; thus, our understanding of the influence of genetic factors on TBI risk and recovery is incomplete. This study aimed to conduct a genome-wide association study (GWAS) of TBI in VA Million Veteran Program (MVP) enrollees. Participants included a multi-ancestry cohort (European, African, and Hispanic ancestries; N = 304,485; 111,494 TBI cases, 192,991 controls). TBI was assessed using MVP survey data and ICD codes from electronic health records. GWAS was performed using logistic regression in PLINK, and FUMA was used for post-GWAS analysis. Genetic correlations between TBI and several health domains (risk-taking behaviors, psychiatric disorders, neurocognitive variables, and brain morphometrics) were also evaluated. SNP-based heritability was 0.060 (SE = 0.004, p = 7.83x10-66). GWAS analysis identified 15 genome-wide significant loci at p < 5x10-8. Gene-based analyses revealed 14 gene-wide significant genes, including NCAM1, APOE, FTO, and FOXP2. Gene tissue expression analysis identified the brain as significantly enriched, particularly in the frontal cortex, anterior cingulate cortex, and nucleus accumbens. Genetic correlations with TBI were significant for risk-taking behaviors and psychiatric disorders, but generally non-significant for neurocognitive and brain morphometric variables. The strongest genetic correlation was with PTSD (rg = 0.69, SE = 0.04, p = 2.5x10-75). This first well-powered GWAS of TBI showed that TBI is a heritable trait and identified 15 loci including genes relevant to TBI biology. Findings also revealed high genetic correlations with psychiatric and risk-taking traits, significantly advancing understanding of the genetic basis of TBI. Results set the stage for future TBI GWASs that focus on diversity and chronicity of symptom sequelae.

Identification of single nucleotide polymorphisms (SNPs) associated with chronic graft-versus-host disease in patients undergoing allogeneic hematopoietic cell transplantation

IntroductionChronic graft-versus-host disease (cGVHD) is a debilitating side effect of allogeneic hematopoietic cell transplantation (HCT), affecting the quality of life of patients. We used whole exome sequencing to identify candidate SNPs and complete a multi-marker gene-level analysis using a cohort of cGVHD( +) (N = 16) and cGVHD( −) (N = 66) HCT patients.MethodsSaliva samples were collected from HCT patients (N = 82) pre-conditioning in a multi-center study from March 2011 to May 2018. Exome sequencing was performed and FASTQ files were processed for sequence alignments. Significant SNPs were identified by logistic regression using PLINK2v3.7 and Fisher’s exact test. One cGVHD( −) patient sample was excluded from further analysis since no SNP was present in at least 10% of the sample population. The FUMA platform’s SNP2GENE was utilized to annotate SNPs and generate a MAGMA output. Chromatin state visualization of lead SNPs was completed using Epilogos tool. FUMA’s GENE2FUNC was used to obtain gene function and tissue expression from lead genomic loci.ResultsLogistic regression classified 986 SNPs associated with cGVHD( +). SNP2GENE returned three genomic risk loci, four lead SNPs, 48 candidate SNPs, seven candidate GWAS tagged SNPs, and four mapped genes. Fisher’s exact test identified significant homozygous genotypes of four lead SNPs (p < 0.05). GENE2FUNC analysis of multi-marker SNP sets identified one positional gene set including lead SNPs for KANK1 and KDM4C and two curated gene sets including lead SNPs for PTPRD, KDM4C, and/or KANK1.ConclusionsOur data suggest that SNPs in three genes located on chromosome 9 confer genetic susceptibility to cGVHD in HCT patients. These genes modulate STAT3 expression and phosphorylation in cancer pathogenesis. The findings may have implications in the modulation of pathways currently targeted by JAK inhibitors in cGVHD clinical trials.

Prognostic significance of the rho GTPase RHOV and its role in tumor immune cell infiltration: a comprehensive pan-cancer analysis.

Ras homolog gene family member V (RHOV) is an atypical Rho GTPase that participates in various important cellular processes. Although RHOV has been identified to play an oncogenic role in lung cancer and triple-negative breast cancer, its role in other types of tumors remains unknown. In this study, we investigated the expression of RHOV in pan-cancer analysis using The Cancer Genome Atlas (TCGA) and Gene-Tissue Expression datasets. RHOV mRNA levels were dysregulated in several types of tumors. RHOV expression was identified as an independent prognostic factor in 7 of 33 types of tumors; however, the relationship varied according to tumor type. Higher RHOV expression was associated with a favorable prognosis in kidney renal cell carcinoma and prostate adenocarcinoma, for which RHOV expression was downregulated, whereas RHOV expression was associated with a poor prognosis for patients with adenoid cystic carcinoma, lung adenocarcinoma, pancreatic ductal adenocarcinoma, skin cutaneous melanoma, and uveal melanoma with upregulated RHOV expression. Furthermore, RHOV expression was associated with various clinicopathological parameters in these tumors. RHOV expression showed varied associations with different types of tumor-infiltrating immune cells and demonstrated a potential impact on the response to immunotherapy depending on the cancer type. Additionally, functional enrichment analysis of RHOV-related genes demonstrated a role in a wide range of developmental and immune-related processes. This study provides valuable insights into the role of RHOV in pan-cancer development, indicating its role as a tumor suppressor or oncogene according to the cancer type and tumor microenvironment.

Genome-wide identification and expression analyses of the pectate lyase (PL) gene family in Fragaria vesca

BackgroundPectate lyase (PL, EC 4.2.2.2), as an endo-acting depolymerizing enzyme, cleaves α-1,4-glycosidic linkages in esterified pectin and involves a broad range of cell wall modifications. However, the knowledge concerning the genome-wide analysis of the PL gene family in Fragaria vesca has not been thoroughly elucidated.ResultsIn this study, sixteen PLs members in F. vesca were identified based on a genome-wide investigation. Substantial divergences existed among FvePLs in gene duplication, cis-acting elements, and tissue expression patterns. Four clusters were classified according to phylogenetic analysis. FvePL6, 8 and 13 in cluster II significantly contributed to the significant expansions during evolution by comparing orthologous PL genes from Malus domestica, Solanum lycopersicum, Arabidopsis thaliana, and Fragaria×ananassa. The cis-acting elements implicated in the abscisic acid signaling pathway were abundant in the regions of FvePLs promoters. The RNA-seq data and in situ hybridization revealed that FvePL1, 4, and 7 exhibited maximum expression in fruits at twenty days after pollination, whereas FvePL8 and FvePL13 were preferentially and prominently expressed in mature anthers and pollens. Additionally, the co-expression networks displayed that FvePLs had tight correlations with transcription factors and genes implicated in plant development, abiotic/biotic stresses, ions/Ca2+, and hormones, suggesting the potential roles of FvePLs during strawberry development. Besides, histological observations suggested that FvePL1, 4 and 7 enhanced cell division and expansion of the cortex, thus negatively influencing fruit firmness. Finally, FvePL1-RNAi reduced leaf size, altered petal architectures, disrupted normal pollen development, and rendered partial male sterility.ConclusionThese results provide valuable information for characterizing the evolution, expansion, expression patterns and functional analysis, which help to understand the molecular mechanisms of the FvePLs in the development of strawberries.

Predicting molecular mechanisms of hereditary diseases by using their tissue‐selective manifestation

How do aberrations in widely expressed genes lead to tissue‐selective hereditary diseases? Previous attempts to answer this question were limited to testing a few candidate mechanisms. To answer this question at a larger scale, we developed “Tissue Risk Assessment of Causality by Expression” (TRACE), a machine learning approach to predict genes that underlie tissue‐selective diseases and selectivity‐related features. TRACE utilized 4,744 biologically interpretable tissue‐specific gene features that were inferred from heterogeneous omics datasets. Application of TRACE to 1,031 disease genes uncovered known and novel selectivity‐related features, the most common of which was previously overlooked. Next, we created a catalog of tissue‐associated risks for 18,927 protein‐coding genes (https://netbio.bgu.ac.il/trace/). As proof‐of‐concept, we prioritized candidate disease genes identified in 48 rare‐disease patients. TRACE ranked the verified disease gene among the patient's candidate genes significantly better than gene prioritization methods that rank by gene constraint or tissue expression. Thus, tissue selectivity combined with machine learning enhances genetic and clinical understanding of hereditary diseases.

Genome-Wide Identification and Characterization of Copper Chaperone for Superoxide Dismutase (CCS) Gene Family in Response to Abiotic Stress in Soybean

Copper Chaperone For Superoxide Dismutase (CCS) genes encode copper chaperone for Superoxide dismutase (SOD) and dramatically affect the activity of SOD through regulating copper delivery from target to SOD. SOD is the effective component of the antioxidant defense system in plant cells to reduce oxidative damage by eliminating Reactive oxygen species (ROS), which are produced during abiotic stress. CCS might play an important role in abiotic stress to eliminate the damage caused by ROS, however, little is known about CCS in soybean in abiotic stress regulation. In this study, 31 GmCCS gene family members were identified from soybean genome. These genes were classified into 4 subfamilies in the phylogenetic tree. Characteristics of 31 GmCCS genes including gene structure, chromosomal location, collinearity, conserved domain, protein motif, cis-elements, and tissue expression profiling were systematically analyzed. RT-qPCR was used to analyze the expression of 31 GmCCS under abiotic stress, and the results showed that 5 GmCCS genes(GmCCS5, GmCCS7, GmCCS8, GmCCS11 and GmCCS24) were significantly induced by some kind of abiotic stress. The functions of these GmCCS genes in abiotic stress were tested using yeast expression system and soybean hairy roots. The results showed that GmCCS7/GmCCS24 participated in drought stress regulation. Soybean hairy roots expressing GmCCS7/GmCCS24 showed improved drought stress tolerance, with increased SOD and other antioxidant enzyme activities. The results of this study provide reference value in-depth study CCS gene family, and important gene resources for the genetic improvement of soybean drought stress tolerance.

102P Co-expression analysis of genes encoding proteasome subunits and XPO1-related proteins in the Cancer Genome Atlas (TCGA) and the Gene Tissue Expression (GTEx) databases as a tool to devise new treatment strategies

102P Co-expression analysis of genes encoding proteasome subunits and XPO1-related proteins in the Cancer Genome Atlas (TCGA) and the Gene Tissue Expression (GTEx) databases as a tool to devise new treatment strategies

Mitochondrial TXNRD2 and ME3 Genetic Risk Scores Are Associated with Specific Primary Open-Angle Glaucoma Phenotypes

Genetic variants in regions that include the mitochondrial genes thioredoxin reductase 2 (TXNRD2) and malic enzyme 3 (ME3) are associated with primary open-angle glaucoma (POAG) in genome-wide association studies (GWASs). To assess their clinical impact, we investigated whether TXNRD2 and ME3 genetic risk scores (GRSs) are associated with specific glaucoma phenotypes. Cross-sectional study. A total of 2617 patients with POAG and 2634 control participants from the National Eye Institute Glaucoma Human Genetics Collaboration Hereditable Overall Operational Database (NEIGHBORHOOD) consortium. All POAG-associated single nucleotide polymorphisms (SNPs) in the TXNRD2 and ME3 loci were identified using GWAS data (P < 0.05). Of these, 20 TXNRD2 and 24 ME3 SNPs were selected after adjusting for linkage disequilibrium. The correlation between SNP effect size and gene expression levels was investigated using the Gene-Tissue Expression database. Genetic risk scores were constructed for each individual using the unweighted sum of TXNRD2, ME3, and TXNRD2 + ME3 combined risk alleles. Age- and sex-adjusted odds ratios (ORs) for POAG diagnosis were calculated per decile for each GRS. Additionally, the clinical features of patients with POAG in the top 1%, 5%, and 10% of each GRS were compared with those in the bottom 1%, 5%, and 10%, respectively. Primary open-angle glaucoma OR per GRS decile, maximum treated intraocular pressure (IOP), and prevalence of paracentral visual field loss among patients with POAG with high versus low GRSs. A larger SNP effect size strongly correlated with higher TXNRD2 and lower ME3 expression levels (r= 0.95 and r= -0.97, respectively; P < 0.05 for both). Individuals in decile 10 of the TXNRD2 + ME3 GRS had the highest odds of POAG diagnosis (OR, 1.79 compared with decile 1; 95% confidence interval, 1.39-2.30; P<0.001). Patients with POAG in the top 1% of the TXNRD2 GRS showed higher mean maximum treated IOP compared with the bottom 1% (19.9 mmHg vs. 15.6 mmHg; adjusted P= 0.03). Patients with POAG in the top 1% of the ME3 and TXNRD2 + ME3 GRS showed a higher prevalence of paracentral field loss than the bottom 1% (72.7% vs. 14.3% for ME3 GRS and 88.9% vs. 33.3% for TXNRD2+ME3 GRS; adjusted P= 0.03 for both). Patients with POAG with higher TXNRD2 and ME3 GRSs showed higher treated IOP and a greater prevalence of paracentral field loss. Functional studies exploring how these variants impact mitochondrial function in patients with glaucoma are warranted. Proprietary or commercial disclosure may be found after the references.

Proteolysis of CD44 at the cell surface controls a downstream protease network.

The cell surface receptor cluster of differentiation 44 (CD44) is the main hyaluronan receptor of the human body. At the cell surface, it can be proteolytically processed by different proteases and was shown to interact with different matrix metalloproteinases. Upon proteolytic processing of CD44 and generation of a C-terminal fragment (CTF), an intracellular domain (ICD) is released after intramembranous cleavage by the γ-secretase complex. This intracellular domain then translocates to the nucleus and induces transcriptional activation of target genes. In the past CD44 was identified as a risk gene for different tumor entities and a switch in CD44 isoform expression towards isoform CD44s associates with epithelial to mesenchymal transition (EMT) and cancer cell invasion. Here, we introduce meprin β as a new sheddase of CD44 and use a CRISPR/Cas9 approach to deplete CD44 and its sheddases ADAM10 and MMP14 in HeLa cells. We here identify a regulatory loop at the transcriptional level between ADAM10, CD44, MMP14 and MMP2. We show that this interplay is not only present in our cell model, but also across different human tissues as deduced from GTEx (Gene Tissue Expression) data. Furthermore, we identify a close relation between CD44 and MMP14 that is also reflected in functional assays for cell proliferation, spheroid formation, migration and adhesion.

Systemic interindividual epigenetic variation in humans is associated with transposable elements and under strong genetic control

BackgroundGenetic variants can modulate phenotypic outcomes via epigenetic intermediates, for example at methylation quantitative trait loci (mQTL). We present the first large-scale assessment of mQTL at human genomic regions selected for interindividual variation in CpG methylation, which we call correlated regions of systemic interindividual variation (CoRSIVs). These can be assayed in blood DNA and do not reflect interindividual variation in cellular composition.ResultsWe use target-capture bisulfite sequencing to assess DNA methylation at 4086 CoRSIVs in multiple tissues from each of 188 donors in the NIH Gene-Tissue Expression (GTEx) program. At CoRSIVs, DNA methylation in peripheral blood correlates with methylation and gene expression in internal organs. We also discover unprecedented mQTL at these regions. Genetic influences on CoRSIV methylation are extremely strong (median R2=0.76), cumulatively comprising over 70-fold more human mQTL than detected in the most powerful previous study. Moreover, mQTL beta coefficients at CoRSIVs are highly skewed (i.e., the major allele predicts higher methylation). Both surprising findings are independently validated in a cohort of 47 non-GTEx individuals. Genomic regions flanking CoRSIVs show long-range enrichments for LINE-1 and LTR transposable elements; the skewed beta coefficients may therefore reflect evolutionary selection of genetic variants that promote their methylation and silencing. Analyses of GWAS summary statistics show that mQTL polymorphisms at CoRSIVs are associated with metabolic and other classes of disease.ConclusionsA focus on systemic interindividual epigenetic variants, clearly enhanced in mQTL content, should likewise benefit studies attempting to link human epigenetic variation to the risk of disease.

Gene identification and tissue expression analysis inform the floral organization and color in the basal angiosperm Magnolia polytepala (Magnoliaceae).

In Magnolia polytepala, the formation of floral organization and color was attributed to tissue-dependent differential expression levels of MADS-box genes and anthocyanin biosynthetic genes. In angiosperms, the diversity of floral morphology and organization suggests its value in exploring plant evolution. Magnolia polytepala, an endemic basal angiosperm species in China, possesses three green sepal-like tepals in the outermost whorl and pink petal-like tepals in the inner three whorls, forming unique floral morphology and organization. However, we know little about its underlying molecular regulatory mechanism. Here, we first reported the full-length transcriptome of M. polytepala using PacBio sequencing. A total of 16 MADS-box transcripts were obtained from the transcriptome data, including floral homeotic genes (e.g., MpAPETALA3) and other non-floral homeotic genes (MpAGL6, etc.). Phylogenetic analysis and spatial expression pattern reflected their putative biological function as their homologues in Arabidopsis. In addition, nine structural genes involved in anthocyanin biosynthesis pathway had been screened out, and tepal color difference was significantly associated with their tissue-dependent differential expression levels. This study provides a relatively comprehensive investigation of the MADS-box family and anthocyanin biosynthetic genes in M. polytepala, and will facilitate our understanding of the regulatory mechanism underlying floral organization and color in basal angiosperms.

Characterization of CRH-Binding Protein (CRHBP) in Chickens: Molecular Cloning, Tissue Distribution and Investigation of Its Role as a Negative Feedback Regulator within the Hypothalamus-Pituitary-Adrenal Axis.

Corticotropin (ACTH) is a pituitary hormone playing important roles in stress response within the hypothalamus–pituitary–adrenal (HPA) axis. The biosynthesis and secretion of ACTH are controlled by multiple factors, including corticotropin-releasing hormone (CRH). As a key hypothalamus-derived regulator, CRH binds to corticotropin-releasing hormone receptor 1 (CRHR1) in the anterior pituitary gland to regulate ACTH synthesis and release. Thus, CRH-binding protein (CRHBP), which binds CRH with high affinity to inhibit CRH-induced ACTH secretion from pituitary cells, draws wide attention. In contrast to the extensive investigation of CRHBP in mammals and other lower vertebrates, the gene structure, tissue expression and physiological functions of CRHBP in birds remain largely unknown. In the present study, using chicken (c-) as our animal model, we examined the gene structure, tissue expression and functionality of CRHBP. Our results showed that: (1) cCRHBP cDNA encodes a 345 amino acid precursor, which shares high sequence identity with that of mammals, reptiles, frogs and fish; (2) cCRHBP is abundantly expressed in the brain (cerebrum and hypothalamus), pituitary and ovary; (3) cCRHBP inhibits the signaling of cCRHRs induced by cCRH, thus reducing the cCRH-induced ACTH secretion from cultured chick pituitary cells; (4) stress mediators (e.g., glucocorticoids) and stress significantly upregulate CRHBP mRNA expression in chickens, supporting its role as a negative feedback regulator in the HPA axis. The present study enriches our understanding of the conserved roles of CRHBP across vertebrates. In addition, chicken is an important poultry animal with multiple economic traits which are tightly controlled by the HPA axis. The characterization of the chicken CRHBP gene helps to reveal the molecular basis of the chicken HPA axis and is thus beneficial to the poultry industry.

Calmodulin and calmodulin-like gene family in barley: Identification, characterization and expression analyses.

Calmodulin (CaM) and calmodulin-like (CML) proteins are Ca2+ relays and play diverse and multiple roles in plant growth, development and stress responses. However, CaM/CML gene family has not been identified in barley (Hordeum vulgare). In the present study, 5 HvCaMs and 80 HvCMLs were identified through a genome-wide analysis. All HvCaM proteins possessed 4 EF-hand motifs, whereas HvCMLs contained 1 to 4 EF-hand motifs. HvCaM2, HvCaM3 and HvCaM5 coded the same polypeptide although they differed in nucleotide sequence, which was identical to the polypeptides coded by OsCaM1-1, OsCaM1-2 and OsCaM1-3. HvCaMs/CMLs were unevenly distributed over barley 7 chromosomes, and could be phylogenetically classified into 8 groups. HvCaMs/CMLs differed in gene structure, cis-acting elements and tissue expression patterns. Segmental and tandem duplication were observed among HvCaMs/CMLs during evolution. HvCML16, HvCML18, HvCML50 and HvCML78 were dispensable genes and the others were core genes in barley pan-genome. In addition, 14 HvCaM/CML genes were selected to examine their responses to salt, osmotic and low potassium stresses by qRT-PCR, and their expression were stress-and time-dependent. These results facilitate our understanding and further functional identification of HvCaMs/CMLs.

Balancing the transcriptome: leveraging sample similarity to improve measures of gene specificity.

The spatial and temporal domain of a gene’s expression can range from ubiquitous to highly specific. Quantifying the degree to which this expression is unique to a specific tissue or developmental timepoint can provide insight into the etiology of genetic diseases. However, quantifying specificity remains challenging as measures of specificity are sensitive to similarity between samples in the sample set. For example, in the Gene-Tissue Expression project (GTEx), brain subregions are overrepresented at 13 of 54 (24%) unique tissues sampled. In this dataset, existing specificity measures have a decreased ability to identify genes specific to the brain relative to other organs. To solve this problem, we leverage sample similarity information to weight samples such that overrepresented tissues do not have an outsized effect on specificity estimates. We test this reweighting procedure on 4 measures of specificity, Z-score, Tau, Tsi and Gini, in the GTEx data and in single cell datasets for zebrafish and mouse. For all of these measures, incorporating sample similarity information to weight samples results in greater stability of sets of genes called as specific and decreases the overall variance in the change of specificity estimates as sample sets become more unbalanced. Furthermore, the genes with the largest improvement in their specificity estimate’s stability are those with functions related to the overrepresented sample types. Our results demonstrate that incorporating similarity information improves specificity estimates’ stability to the choice of the sample set used to define the transcriptome, providing more robust and reproducible measures of specificity for downstream analyses.

Concurrent outcomes from multiple approaches of epistasis analysis for human body mass index associated loci provide insights into obesity biology

Genome wide association studies (GWAS) have focused on elucidating the genetic architecture of complex traits by assessing single variant effects in additive genetic models, albeit explaining a fraction of the trait heritability. Epistasis has recently emerged as one of the intrinsic mechanisms that could explain part of this missing heritability. We conducted epistasis analysis for genome-wide body mass index (BMI) associated SNPs in Alzheimer’s Disease Neuroimaging Initiative (ADNI) and followed up top significant interacting SNPs for replication in the UK Biobank imputed genotype dataset. We report two pairwise epistatic interactions, between rs2177596 (RHBDD1) and rs17759796 (MAPK1), rs1121980 (FTO) and rs6567160 (MC4R), obtained from a consensus of nine different epistatic approaches. Gene interaction maps and tissue expression profiles constructed for these interacting loci highlights co-expression, co-localisation, physical interaction, genetic interaction, and shared pathways emphasising the neuronal influence in obesity and implicating concerted expression of associated genes in liver, pancreas, and adipose tissues insinuating to metabolic abnormalities characterized by obesity. Detecting epistasis could thus be a promising approach to understand the effect of simultaneously interacting multiple genetic loci in disease aetiology, beyond single locus effects.

Association mapping and domestication analysis to dissect genetic improvement process of upland cotton yield-related traits in China

BackgroundCotton fiber yield is a complex trait, which can be influenced by multiple agronomic traits. Unravelling the genetic basis of cotton fiber yield-related traits contributes to genetic improvement of cotton.ResultsIn this study, 503 upland cotton varieties covering the four breeding stages (BS1–BS4, 1911–2011) in China were used for association mapping and domestication analysis. One hundred and forty SSR markers significantly associated with ten fiber yield-related traits were identified, among which, 29 markers showed an increasing trend contribution to cotton yield-related traits from BS1 to BS4, and 26 markers showed decreased trend effect. Four favorable alleles of 9 major loci (R2 ≥ 3) were strongly selected during the breeding stages, and the candidate genes of the four strongly selected alleles were predicated according to the gene function annotation and tissue expression data.ConclusionsThe study not only uncovers the genetic basis of 10 cotton yield-related traits but also provides genetic evidence for cotton improvement during the cotton breeding process in China.