Tissue-specific Expression Patterns Research Articles

Genome-wide identification and characterization of polycomb repressive complex 2 core components in upland cotton (Gossypium hirsutum L.)

BackgroundThe evolutionarily conserved Polycomb Repressive Complex 2 (PRC2) plays a vital role in epigenetic gene repression by depositing tri-methylation on lysine residue K27 of histone H3 (H3K27me3) at the target loci, thus participating in diverse biological processes. However, few reports about PRC2 are available in plant species with large and complicated genomes, like cotton.ResultsHere, we performed a genome-wide identification and comprehensive analysis of cotton PRC2 core components, especially in upland cotton (Gossypium hirsutum). Firstly, a total of 8 and 16 PRC2 core components were identified in diploid and tetraploid cotton species, respectively. These components were classified into four groups, E(z), Su(z)12, ESC and p55, and the members in the same group displayed good collinearity, similar gene structure and domain organization. Next, we cloned G. hirsutum PRC2 (GhPRC2) core components, and found that most of GhPRC2 proteins were localized in the nucleus, and interacted with each other to form multi-subunit complexes. Moreover, we analyzed the expression profile of GhPRC2 genes. The transcriptome data and quantitative real-time PCR (qRT-PCR) assays indicated that GhPRC2 genes were ubiquitously but differentially expressed in various tissues, with high expression levels in reproductive organs like petals, stamens and pistils. And the expressions of several GhPRC2 genes, especially E(z) group genes, were responsive to various abiotic and biotic stresses, including drought, salinity, extreme temperature, and Verticillium dahliae (Vd) infection.ConclusionWe identified PRC2 core components in upland cotton, and systematically investigated their classifications, phylogenetic and synteny relationships, gene structures, domain organizations, subcellular localizations, protein interactions, tissue-specific and stresses-responsive expression patterns. Our results will provide insights into the evolution and composition of cotton PRC2, and lay the foundation for further investigation of their biological functions and regulatory mechanisms.

Genome-Wide Identification and Expression Analysis of the bHLH Transcription Factor Family and Its Response to Abiotic Stress in Mongolian Oak (Quercus mongolica)

The basic helix-loop-helix (bHLH) family, one of the largest families of transcription factors in plants, is extensively involved in the growth, development, and stress response of several woody plants. However, no systematic analysis of the bHLH gene family in Quercus mongolica has been reported. We characterize QmbHLH genes and identify the functions of QmbHLH proteins in Q. mongolica. We used bioinformatics approaches, qRT-PCR analysis, and RNA sequencing data to examine chromosomal distributions, gene structures, and conserved patterns, and identified 89 QmbHLH genes, which were divided into 21 subgroups based on the phylogenetic analysis of bHLH genes in Arabidopsis thaliana. Segmental replication played a more prominent role than tandem duplication in the expansion of the QmbHLH gene family. Based on patterns of tissue-specific expression, protein interactions, and cis-element analysis, QmbHLH genes may be extensively involved in the growth and development of Q. mongolica. In leaves, stems, and roots, 12 selected QmbHLH genes exhibited responsiveness to abiotic stresses (salt, cold, weak light, and drought). Our study facilitates follow-up functional investigations of the bHLH gene family in Q. mongolica and provides novel insights into bHLH superfamilies in woody plants.

Genome-Wide Analysis and Expression of Cyclic Nucleotide-Gated Ion Channel (CNGC) Family Genes under Cold Stress in Mango (Mangifera indica).

The 'king of fruits' mango (Mangifera indica) is widely cultivated in tropical areas and has been threatened by frequent extreme cold weather. Cyclic nucleotide-gated ion channel (CNGC) genes have an important function in the calcium-mediated development and cold response of plants. However, few CNGC-related studies are reported in mango, regardless of the mango cold stress response. In this study, we identified 43 CNGC genes in mango showing tissue-specific expression patterns. Five MiCNGCs display more than 3-fold gene expression induction in the fruit peel and leaf under cold stress. Among these, MiCNGC9 and MiCNGC13 are significantly upregulated below 6 °C, suggesting their candidate functions under cold stress. Furthermore, cell membrane integrity was damaged at 2 °C in the mango leaf, as shown by the content of malondialdehyde (MDA), and eight MiCNGCs are positively correlated with MDA contents. The high correlation between MiCNGCs and MDA implies MiCNGCs might regulate cell membrane integrity by regulating MDA content. Together, these findings provide a valuable guideline for the functional characterization of CNGC genes and will benefit future studies related to cold stress and calcium transport in mango.

Characterization of the soybean KRP gene family reveals a key role for GmKRP2a in root development.

Kip-related proteins (KRPs), as inhibitory proteins of cyclin-dependent kinases, are involved in the growth and development of plants by regulating the activity of the CYC-CDK complex to control cell cycle progression. The KRP gene family has been identified in several plants, and several KRP proteins from Arabidopsis thaliana have been functionally characterized. However, there is little research on KRP genes in soybean, which is an economically important crop. In this study, we identified nine GmKRP genes in the Glycine max genome using HMM modeling and BLASTP searches. Protein subcellular localization and conserved motif analysis showed soybean KRP proteins located in the nucleus, and the C-terminal protein sequence was highly conserved. By investigating the expression patterns in various tissues, we found that all GmKRPs exhibited transcript abundance, while several showed tissue-specific expression patterns. By analyzing the promoter region, we found that light, low temperature, an anaerobic environment, and hormones-related cis-elements were abundant. In addition, we performed a co-expression analysis of the GmKRP gene family, followed by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) set enrichment analysis. The co-expressing genes were mainly involved in RNA synthesis and modification and energy metabolism. Furthermore, the GmKRP2a gene, a member of the soybean KRP family, was cloned for further functional analysis. GmKRP2a is located in the nucleus and participates in root development by regulating cell cycle progression. RNA-seq results indicated that GmKRP2a is involved in cell cycle regulation through ribosome regulation, cell expansion, hormone response, stress response, and plant pathogen response pathways. To our knowledge, this is the first study to identify and characterize the KRP gene family in soybean.

Optimization and application of non-native Phragmites australis transcriptome assemblies.

Phragmites australis (common reed) has a cosmopolitan distribution and has been suggested as a model organism for the study of invasive plant species. In North America, the non-native subspecies (ssp. australis) is widely distributed across the contiguous 48 states in the United States and large parts of Canada. Even though millions of dollars are spent annually on Phragmites management, insufficient knowledge of P. australis impeded the efficiency of management. To solve this problem, transcriptomic information generated from multiple types of tissue could be a valuable resource for future studies. Here, we constructed forty-nine P. australis transcriptomes assemblies via different assembly tools and multiple parameter settings. The optimal transcriptome assembly for functional annotation and downstream analyses was selected among these transcriptome assemblies by comprehensive assessments. For a total of 422,589 transcripts assembled in this transcriptome assembly, 319,046 transcripts (75.5%) have at least one functional annotation. Within the transcriptome assembly, we further identified 1,495 transcripts showing tissue-specific expression pattern, 10,828 putative transcription factors, and 72,165 candidates for simple sequence repeats markers. The identification and analyses of predicted transcripts related to herbicide- and salinity-resistant genes were shown as two applications of the transcriptomic information to facilitate further research on P. australis. Transcriptome assembly and selection would be important for the transcriptome annotation. With this optimal transcriptome assembly and all relative information from downstream analyses, we have helped to establish foundations for future studies on the mechanisms underlying the invasiveness of non-native P. australis subspecies.

Identification of the Genome-Wide Expression Patterns of Non-Coding RNAs Associated with Tanshinones Synthesis Pathway in Salvia miltiorrhiza

The red root of Salvia miltiorrhiza Bunge, a famous traditional Chinese medicine (TCM), was caused by tanshinone in epidermis cells. In order to study the biological function of ncRNAs in the tanshinone synthesis, the expression patterns of mRNA and ncRNAs were comprehensively analyzed in red (high tanshinone content) and white root (low tanshinone content) tissues derived from the same plant. A total of 731 differentially expressed genes (DEGs) were mainly enriched in primary metabolic pathways such as galactose and nitrogen, and some secondary metabolic pathways such as phenylpropanoid and terpenoids. A total of 70 miRNAs, 48 lncRNAs, and 26 circRNAs were identified as differentially expressed (DE) ones. The enrichment pathway of the targets of DE-lncRNA were mainly in ribosome, carbon metabolism, plant hormone signal transduction, and glycerophospholipid metabolism. The function of the targets genes of 59 miRNAs combined with DE-circRNAs was mainly involved in plant–pathogen interaction, endocytosis, phenylpropanoid biosynthesis, and sesquiterpenoid and triterpenoid biosynthesis pathways. Most genes of the tanshinone synthesis pathway had a higher expression. Some ncRNAs were predicted to regulate several key enzyme genes of the tanshinone synthesis pathway, such as SmDXS2, SmGGPPS1, SmKSL. Furthermore, most target genes were related to the resistance of pathogens. The present study exhibited the tissue-specific expression patterns of ncRNAs, which would provide a basis for further research into the regulation mechanism of ncRNAs in the tanshinone synthesis process.

Investigating nicotine pathway-related long non-coding RNAs in tobacco.

Long non-coding RNAs (lncRNAs) are transcripts longer than 200bp with low or no protein-coding ability, which play essential roles in various biological processes in plants. Tobacco is an ideal model plant for studying nicotine biosynthesis and metabolism, and there is little research on lncRNAs in this field. Therefore, how to take advantage of the mature tobacco system to profoundly investigate the lncRNAs involved in the nicotine pathway is intriguing. By exploiting 549 public RNA-Seq datasets of tobacco, 30,212 lncRNA candidates were identified, including 24,084 large intervening non-coding RNAs (lincRNAs), 5,778 natural antisense transcripts (NATs) and 350 intronic non-coding RNAs (incRNAs). Compared with protein-coding genes, lncRNAs have distinct properties in terms of exon number, sequence length, A/U content, and tissue-specific expression pattern. lincRNAs showed an asymmetric evolutionary pattern, with a higher proportion (68.71%) expressed from the Nicotiana sylvestris (S) subgenome. We predicted the potential cis/trans-regulatory effects on protein-coding genes. One hundred four lncRNAs were detected as precursors of 30 known microRNA (miRNA) family members, and 110 lncRNAs were expected to be the potential endogenous target mimics for 39 miRNAs. By combining the results of weighted gene co-expression network analysis with the differentially expressed gene analysis of topping RNA-seq data, we constructed a sub-network containing eight lncRNAs and 25 nicotine-related coding genes. We confirmed that the expression of seven lncRNAs could be affected by MeJA treatment and may be controlled by the transcription factor NtMYC2 using a quantitative PCR assay and gene editing. The results suggested that lncRNAs are involved in the nicotine pathway. Our findings further deepened the understanding of the features and functions of lncRNAs and provided new candidates for regulating nicotine biosynthesis in tobacco.

Insulin-like peptide 8 (Ilp8) regulates female fecundity in flies.

Introduction: Insulin-like peptides (Ilps) play crucial roles in nearly all life stages of insects. Ilp8 is involved in developmental stability, stress resistance and female fecundity in several insect species, but the underlying mechanisms are not fully understood. Here we report the functional characterization of Ilp8s in three fly species, including Bactrocera dorsalis, Drosophila mercatorum and Drosophila melanogaster. Methods: Phylogenetic analyses were performed to identify and characterize insect Ilp8s. The amino acid sequences of fly Ilp8s were aligned and the three-dimensional structures of fly Ilp8s were constructed and compared. The tissue specific expression pattern of fly Ilp8s were examined by qRT-PCR. In Bactrocera dorsalis and Drosophila mercatorum, dsRNAs were injected into virgin females to inhibit the expression of Ilp8 and the impacts on female fecundity were examined. In Drosophila melanogaster, the female fecundity of Ilp8 loss-of-function mutant was compared with wild type control flies. The mutant fruit fly strain was also used for sexual behavioral analysis and transcriptomic analysis. Results: Orthologs of Ilp8s are found in major groups of insects except for the lepidopterans and coleopterans, and Ilp8s are found to be well separated from other Ilps in three fly species. The key motif and the predicted three-dimensional structure of fly Ilp8s are well conserved. Ilp8 are specifically expressed in the ovary and are essential for female fecundity in three fly species. Behavior analysis demonstrates that Ilp8 mutation impairs female sexual attractiveness in fruit fly, which results in decreased mating success and is likely the cause of fecundity reduction. Further transcriptomic analysis indicates that Ilp8 might influence metabolism, immune activity, oocyte development as well as hormone homeostasis to collectively regulate female fecundity in the fruit fly. Discussion: Our findings support a universal role of insect Ilp8 in female fecundity, and also provide novel clues for understanding the modes of action of Ilp8.

The chromosome-scale genome assembly of the yellowtail clownfish Amphiprion clarkii provides insights into the melanic pigmentation of anemonefish.

Anemonefish are an emerging group of model organisms for studying genetic, ecological, evolutionary, and developmental traits of coral reef fish. The yellowtail clownfish Amphiprion clarkii possesses species-specific characteristics such as inter-species co-habitation, high intra-species color variation, no anemone specificity, and a broad geographic distribution, that can increase our understanding of anemonefish evolutionary history, behavioral strategies, fish-anemone symbiosis, and color pattern evolution. Despite its position as an emerging model species, the genome of A. clarkii is yet to be published. Using PacBio long-read sequencing and Hi-C chromatin capture technology, we generated a high-quality chromosome-scale genome assembly initially comprised of 1,840 contigs with an N50 of 1,203,211 bp. These contigs were successfully anchored into 24 chromosomes of 843,582,782 bp and annotated with 25,050 protein-coding genes encompassing 97.0% of conserved actinopterygian genes, making the quality and completeness of this genome the highest among all published anemonefish genomes to date. Transcriptomic analysis identified tissue-specific gene expression patterns, with the brain and optic lobe having the largest number of expressed genes. Further analyses revealed higher copy numbers of erbb3b (a gene involved in melanocyte development) in A. clarkii compared with other anemonefish, thus suggesting a possible link between erbb3b and the natural melanism polymorphism observed in A. clarkii. The publication of this high-quality genome, along with A. clarkii's many unique traits, position this species as an ideal model organism for addressing scientific questions across a range of disciplines.

Sorghum WRKY transcription factor SbWRKY45 enhanced seed germination under drought stress in transgenic Arabidopsis

The WRKY transcription factors (TFs) family is an important family of plant-specific TFs, playing vital roles in various abiotic and biotic stress responses. WRKY TFs are gaining considerable attention due to their significant roles in stress responses. However, their functions in sorghum (Sorghum bicolor) are lagging behind. In this study, a WRKY gene designated as SbWRKY45 was isolated and characterized from sorghum. SbWRKY45, belonging to Group IIa, consists of one intron and two exons and encodes 430 amino acids. SbWRKY45 is located in chromosome 4. The cis-element prediction analysis showed that the promoter region of SbWRKY45 has several abiotic stress-associated elements. The qRT-PCR results showed that SbWRKY45 was significantly up-regulated in response to drought and cold under salt treatments, was notably changed, and was induced weakly under heat stress. SbWRKY45 exhibited a response to stress in different sorghum tissues, including leaves, stems, and roots. A tissue-specific expression pattern showed that SbWRKY45 was highly expressed in roots compared with leaves and stems, suggesting that SbWRKY45 may play an important role in roots. Overexpression of SbWRKY45 increased germination rates and promoted root growth in transgenic Arabidopsis under drought stress. Taken together, our findings indicate that SbWRKY45 may be involved in mediating the response to drought stress and play a vital role in the abiotic stress response of sorghum

Identification of <i>WUSCHEL</i>-related homeobox (<i>WOX</i>) gene family members and determination of their expression profiles during somatic embryogenesis in <i>Phoebe bournei</i>

WUSCHEL-related homeobox (<italic>WOX</italic>) transcription factor (TF)-encoding genes play crucial roles during embryo development. The function of <italic>WOX</italic> genes in embryonic development has been thoroughly studied in <italic>Arabidopsis thaliana</italic>, but little is known about their function in woody species, especially <italic>Phoebe bournei</italic>, an endemic and endangered species in China. In the present study, a total of 15 <italic>WOX</italic> genes were identified in<italic> P. bournei</italic>, and phylogenetic analysis resulted in their assignment to three typical clades: an ancient clade, an intermediate clade, and a modern/WUS clade. The gene structure and sequence characteristics and the physicochemical properties of WOX proteins were also analyzed. Promoter prediction indicated that <italic>WOX</italic> genes are likely involved in plant growth and development and hormone responses. Subsequently, we evaluated the expression patterns of <italic>WOX</italic> genes in response to auxin (IAA), abscisic acid (ABA), and methyl jasmonate (MeJA) treatments. According to tissue-specific expression patterns, we screened nine <italic>WOX</italic> genes that were present in embryonic calli and that might participate in the somatic embryogenesis (SE) of <italic>P. bournei.</italic> Furthermore, the expression profiles of these nine <italic>WOX</italic> genes during three phases of embryogenic calli development and three phases of somatic embryo development, namely, spheroid embryogenesis, immature cotyledon-producing embryogenesis and mature cotyledon-producing embryogenesis, were monitored. Overall, we systematically analyzed the expression patterns of <italic>WOX</italic> genes in <italic>P. bournei</italic> during SE, the information of which provides a basis for further elucidating the molecular mechanism through which <italic>WOX</italic> TFs function in <italic>P. bournei</italic> embryo development.

Tuning the Transcriptional Activity of the CaMV 35S Promoter in Plants by Single-Nucleotide Changes in the TATA Box.

Synthetic biology uses genetically encoded devices and circuits to implement novel complex functions in living cells and organisms. A hallmark of these genetic circuits is the interaction among their individual parts, according to predefined rules, to process cellular information and produce a circuit output or response. As the number of individual components in a genetic circuit increases, so does the number of interactions needed to achieve the correct behavior, and hence, a greater need to fine-tune the levels of expression of each component. Transcriptional promoters play a key regulatory role in genetic circuits, as they influence the levels of RNA and proteins produced. In multicellular organisms, such as plants, they can also determine developmental, spatial, and tissue-specific patterns of gene expression. The 35S promoter from the Cauliflower Mosaic Virus (CaMV 35S) is widely used in plant biotechnology to direct high levels of gene expression in a variety of plant species. We produced a library of 21 variants of the CaMV 35S promoter by introducing all single nucleotide substitutions to the promoter's TATA box sequence. We then characterized the activity of all variants in homozygous transgenic plants and showed that some of these variants have lower activity than the wild type in plants. These promoter variants could be used to fine-tune the behavior of synthetic genetic circuits in plants.

The Plant Fatty Acyl Reductases.

Fatty acyl reductase (FAR) is a crucial enzyme that catalyzes the NADPH-dependent reduction of fatty acyl-CoA or acyl-ACP substrates to primary fatty alcohols, which in turn acts as intermediate metabolites or metabolic end products to participate in the formation of plant extracellular lipid protective barriers (e.g., cuticular wax, sporopollenin, suberin, and taproot wax). FARs are widely present across plant evolution processes and play conserved roles during lipid synthesis. In this review, we provide a comprehensive view of FAR family enzymes, including phylogenetic analysis, conserved structural domains, substrate specificity, subcellular localization, tissue-specific expression patterns, their varied functions in lipid biosynthesis, and the regulation mechanism of FAR activity. Finally, we pose several questions to be addressed, such as the roles of FARs in tryphine, the interactions between transcription factors (TFs) and FARs in various environments, and the identification of post-transcriptional, translational, and post-translational regulators.

Single-Molecule Real-Time Sequencing of Full-Length Transcriptome and Identification of Genes Related to Male Development in Cannabis sativa.

Female Cannabis sativa plants have important therapeutic properties. The sex ratio of the dioecious cannabis is approximately 1:1. Cultivating homozygous female plants by inducing female plants to produce male flowers is of great practical significance. However, the mechanism underlying cannabis male development remains unclear. In this study, single-molecule real-time (SMRT) sequencing was performed using a mixed sample of female and induced male flowers from the ZYZM1 cannabis variety. A total of 15,241 consensus reads were identified, and 13,657 transcripts were annotated across seven public databases. A total of 48 lncRNAs with an average length of 986.54 bp were identified. In total, 8202 transcripts were annotated as transcription factors, the most common of which were bHLH transcription factors. Moreover, tissue-specific expression pattern analysis showed that 13 MADS transcription factors were highly expressed in male flowers. Furthermore, 232 reads of novel genes were predicted and enriched in lipid metabolism, and qRT-PCR results showed that CER1 may be involved in the development of cannabis male flowers. In addition, 1170 AS events were detected, and two AS events were further validated. Taken together, these results may improve our understanding of the complexity of full-length cannabis transcripts and provide a basis for understanding the molecular mechanism of cannabis male development.

PHD-finger family genes in wheat (Triticum aestivum L.): Evolutionary conservatism, functional diversification, and active expression in abiotic stress.

Plant homeodomain (PHD) transcription factors (TFs) are a class of proteins with conserved Cys4-His-Cys3 domains that play important roles in plant growth and development and in response to abiotic stresses. Although characterization of PHDs has been performed in plants, little is known about their function in wheat (Triticum aestivum L.), especially under stress conditions. In the present study, 244 TaPHDs were identified in wheat using comparative genomics. We renamed them TaPHD1-244 based on their chromosomal distribution, and almost all PHD proteins were predicted to be located in the nucleus. According to the unrooted neighbor-joining phylogenetic tree, gene structure, and motif analyses, PHD genes were divided into four clades. A total of 149 TaPHD genes were assigned to arise from duplication events. Furthermore, 230 gene pairs came from wheat itself, and 119, 186, 168, 7, 2, and 6 gene pairs came from six other species (Hordeum vulgareto, Zea mays, Oryza sativa, Arabidopsis thaliana, Brassica rapa, and Gossypium raimondii, respectively). A total of 548 interacting protein branches were identified to be involved in the protein interaction network. Tissue-specific expression pattern analysis showed that TaPHDs were highly expressed in the stigma and ovary during flowering, suggesting that the TaPHD gene plays an active role in the reproductive growth of wheat. In addition, the qRT-PCR results further confirmed that these TaPHD genes are involved in the abiotic stress response of wheat. In conclusion, our study provides a theoretical basis for deciphering the molecular functions of TaPHDs, particularly in response to abiotic stress.

Genome-wide identification and functional analysis of Dof transcription factor family in Camelina sativa

BackgroundDof transcription factors (TFs) containing C2-C2 zinc finger domains are plant-specific regulatory proteins, playing crucial roles in a variety of biological processes. However, little is known about Dof in Camelina sativa, an important oil crop worldwide, with high stress tolerance. In this study, a genome-wide characterization of Dof proteins is performed to examine their basic structural characteristics, phylogenetics, expression patterns, and functions to identify the regulatory mechanism underlying lipid/oil accumulation and the candidate Dofs mediating stress resistance regulation in C. sativa.ResultsTotal of 103 CsDof genes unevenly distributed on 20 chromosomes were identified from the C. sativa genome, and they were classified into four groups (A, B, C and D) based on the classification of Arabidopsis Dof gene family. All of the CsDof proteins contained the highly-conserved typic CX2C-X21-CX2C structure. Segmental duplication and purifying selection were detected for CsDof genes. 61 CsDof genes were expressed in multiple tissues, and 20 of them showed tissue-specific expression patterns, suggesting that CsDof genes functioned differentially in different tissues of C. sativa. Remarkably, a set of CsDof members were detected to be possible involved in regulation of oil/lipid biosynthesis in C. sativa. Six CsDof genes exhibited significant expression changes in seedlings under salt stress treatment.ConclusionsThe present data reveals that segmental duplication is the key force responsible for the expansion of CsDof gene family, and a strong purifying pressure plays a crucial role in CsDofs’ evolution. Several CsDof TFs may mediate lipid metabolism and stress responses in C. sativa. Several CsDof TFs may mediate lipid metabolism and stress responses in C. sativa. Collectively, our findings provide a foundation for deep understanding the roles of CsDofs and genetic improvements of oil yield and salt stress tolerance in this species and the related crops.

Functional Divergence of the N-Lobe and C-Lobe of Transferrin Gene in Pungitius sinensis (Amur Stickleback).

Transferrin is an important iron-binding glycosylated protein and plays key roles in iron-binding and immune response. Here, a 2037-bp open reading frame was obtained from our previous transcriptome sequencing data of Amur stickleback, which encoded a 679 amino acid putative transferrin protein harbored obvious N-lobe and C-lobe domains. The tissue-specific expression pattern showed that the transcript was detected in a variety of tissues, with the highest signal in liver. Moreover, Streptococcus iniae pathogen stimulation can increase the expression level of this transcript, implying important immune properties for organisms. Next, N-lobes and C-lobes were obtained from 45 fish species. The phylogenetic tree showed that N-lobes and C-lobes were in two different evolutionary branches, and they had different motif composition. Functional divergence indicated a higher evolutionary rate or site-specific alteration among the N-lobe and C-lobe groups. Ka/Ks value of C-lobe group was relatively higher than that of N-lobe group, indicating a faster change rate of C-lobe sequences in evolution. Moreover, some sites experiencing positive selection were also found, which may be involved in the iron- or anion-binding, pathogen resistance and diversification of transferrin protein. Differential iron-binding activity was also detected between N-lobe and C-lobe of Amur stickleback transferrin protein with Chrome Azurol S assay. Compared with the C-lobe, the N-lobe showed stronger growth inhibitory activity of Escherichia coli, implying their potential antibacterial properties. This study will give a reference for subsequent research of transferrin proteins.

Examination of a novel expression-based gene-SNP annotation strategy to identify tissue-specific contributions to heritability in multiple traits.

Complex traits show clear patterns of tissue-specific expression influenced by single nucleotide polymorphisms (SNPs), yet current strategies aggregate SNP effects to genes by employing simple physical proximity-based windows. Here, we examined whether incorporating SNPs with effects on tissue-specific cis-expression would improve our ability to detect trait-relevant tissues across 31 complex traits using stratified linkage disequilibrium score regression (S-LDSC). We found that a physical proximity annotation produced more significant tissue enrichments and larger S-LDSC regression coefficients, as compared to an expression-based annotation. Furthermore, we showed that our expression-based annotation did not outperform an annotation strategy in which an equal number of randomly chosen SNPs were annotated to genes within the same genomic window, suggesting extensive redundancy among SNP effect estimates due to linkage disequilibrium. That said, current sample sizes limit estimation of cis-genetic SNP effects; therefore, we recommend reexamination of the expression-based annotation when larger tissue-specific expression datasets become available. To examine the influence of sample size, we used a large whole blood eQTL reference panel (N = 31,684) applying a similar expression-based annotation strategy. We found that significant cis-expression QTLs in whole blood did not outperform the physical proximity annotation when estimating tissue-specific SNP heritability enrichment for either high- or low-density lipoprotein phenotypes but performed similarly for inflammatory bowel disease. Finally, we report new and updated tissue enrichment estimates across 31 complex traits, such as significant heritability enrichment of the frontal cortex for cognitive performance, educational attainment, and intelligence, providing further evidence of this structure's importance in higher cognitive function.

Genome-wide identification, molecular evolution and expression analysis of the non-specific lipid transfer protein (nsLTP) family in Setaria italica

BackgroundFoxtail millet (Setaria italica L.) is a millet species with high tolerance to stressful environments. Plant non-specific lipid transfer proteins (nsLTPs) are a kind of small, basic proteins involved in many biological processes. So far, the genome of S. italica has been fully sequenced, and a comprehensive understanding of the evolution and expression of the nsLTP family is still lacking in foxtail millet.ResultsForty-five nsLTP genes were identified in S. italica and clustered into 5 subfamilies except three single genes (SinsLTP38, SinsLTP7, and SinsLTP44). The proportion of SinsLTPs was different in each subfamily, and members within the same subgroup shared conserved exon–intron structures. Besides, 5 SinsLTP duplication events were investigated. Both tandem and segmental duplication contributed to nsLTP expansion in S. italica, and the duplicated SinsLTPs had mainly undergone purifying selection pressure, which suggested that the function of the duplicated SinsLTPs might not diverge much. Moreover, we identified the nsLTP members in 5 other monocots, and 41, 13, 10, 4, and 1 orthologous gene pairs were identified between S. italica and S. viridis, S. bicolor, Z. mays, O. sativa, and B. distachyon, respectively. The functional divergence within the nsLTP orthologous genes might be limited. In addition, the tissue-specific expression patterns of the SinsLTPs were investigated, and the expression profiles of the SinsLTPs in response to abiotic stress were analyzed, all the 10 selected SinsLTPs were responsive to drought, salt, and cold stress. Among the selected SinsLTPs, 2 paired duplicated genes shared almost equivalent expression profiles, suggesting that these duplicated genes might retain some essential functions during subsequent evolution.ConclusionsThe present study provided the first systematic analysis for the phylogenetic classification, conserved domain and gene structure, expansion pattern, and expression profile of the nsLTP family in S. italica. These findings could pave a way for further comparative genomic and evolution analysis of nsLTP family in foxtail millet and related monocots, and lay the foundation for the functional analysis of the nsLTPs in S. italica.

Genome-Wide Identification and Expression Analysis of the Aquaporin Gene Family in Lycium barbarum during Fruit Ripening and Seedling Response to Heat Stress.

Plant-water relations mediated by aquaporins (AQPs) play vital roles in both key plant growth processes and responses to environmental challenges. As a well-known medicinal and edible plant, the harsh natural growth habitat endows Lycium plants with ideal materials for stress biology research. However, the details of their molecular switch for water transport remain unclear. In the present work, we first identified and characterized AQP family genes from Lycium (L.) barbarum at the genome scale and conducted systemic bioinformatics and expression analyses. The results showed that there were 38 Lycium barbarum AQPs (LbAQPs) in L.&amp;nbsp;barbarum, which were classified into four subfamilies, including 17 LbPIP, 9 LbTIP, 10 LbNIP, and 2 LbXIP. Their encoded genes were unevenly distributed on all 12 chromosomes, except chromosome 10. Three of these genes encoded truncated proteins and three genes underwent clear gene duplication events. Cis-acting element analysis indicated that the expression of LbAQPs may be mainly regulated by biotic/abiotic stress, phytohormones and light. The qRT-PCR assay indicated that this family of genes presented a clear tissue-specific expression pattern, in which most of the genes had maximal transcript levels in roots, stems, and leaves, while there were relatively lower levels in flowers and fruits. Most of the LbAQP genes were downregulated during L. barbarum fruit ripening and presented a negative correlation with the fruit relative water content (RWC). Most of their transcripts presented a quick and sharp upregulation response to heat stress following exposure of the 2-month-old seedlings to a 42 °C temperature for 0, 1, 3, 12, or 24 h. Our results proposed that LbAQPs were involved in L. barbarum key development events and abiotic stress responses, which may lay a foundation for further studying the molecular mechanism of the water relationship of Lycium plants, especially in harsh environments.