Tissue-specific Expression Data Research Articles

Genome-wide identification and expression profiling of growth‑regulating factor (GRF) and GRF‑interacting factor (GIF) gene families in chickpea and pigeonpea

The growth-regulating factor (GRF) and GRF-interacting factor (GIF) families encode plant-specific transcription factors and play vital roles in plant development and stress response processes. Although GRF and GIF genes have been identified in various plant species, there have been no reports of the analysis and identification of the GRF and GIF transcription factor families in chickpea (Cicer arietinum) and pigeonpea (Cajanus cajan). The present study identified seven CaGRFs, eleven CcGRFs, four CaGIFs, and four CcGIFs. The identified proteins were grouped into eight and three clades for GRFs and GIFs, respectively based on their phylogenetic relationships. A comprehensive in-silico analysis was performed to determine chromosomal location, sub-cellular localization, and types of regulatory elements present in the putative promoter region. Synteny analysis revealed that GRF and GIF genes showed diploid-polyploid topology in pigeonpea, but not in chickpea. Tissue-specific expression data at the vegetative and reproductive stages of the plant showed that GRFs and GIFs were strongly expressed in tissues like embryos, pods, and seeds, indicating that GRFs and GIFs play vital roles in plant growth and development. This research characterized GRF and GIF families and hints at their primary roles in the chickpea and pigeonpea growth and developmental process. Our findings provide potential gene resources and vital information on GRF and GIF gene families in chickpea and pigeonpea, which will help further understand the regulatory role of these gene families in plant growth and development.

Expanding the phenotypic spectrum of CC2D2A-related ciliopathies: a rare homozygous nonsense variant in a patient with suspected nephronophthisis.

Biallelic pathogenic variants in the gene CC2D2A cause a spectrum of ciliopathies, including Joubert and Meckel syndrome, which frequently involve the kidney; however, no cases of isolated renal disease (i.e., nephronophthisis) have yet been reported. In an adult with a clinical presentation consistent with nephronophthisis, next-generation sequencing identified a rare homozygous nonsense variant in CC2D2A (c.100 C > T; p.(Arg34*)). Tissue-specific expression data and promoter activity analysis demonstrates that this variant primarily affects a transcript isoform predominant in the kidneys but does not affect the transcript isoforms predominant in other tissues typically involved in CC2D2A-related ciliopathies (e.g., cerebellum, liver). Expression analysis of patient-specific cDNA in MDCK cells demonstrates partial translation re-initiation downstream of p.(Arg34*) as a possible escape mechanism from nonsense mediated decay. These data provide mechanistic insights in support of this novel genotype-phenotype association.

GTExome: Modeling commonly expressed missense mutations in the human genome.

A web application, GTExome, is described that quickly identifies, classifies, and models missense mutations in commonly expressed human proteins. GTExome can be used to categorize genomic mutation data with tissue specific expression data from the Genotype-Tissue Expression (GTEx) project. Commonly expressed missense mutations in proteins from a wide range of tissue types can be selected and assessed for modeling suitability. Information about the consequences of each mutation is provided to the user including if disulfide bonds, hydrogen bonds, or salt bridges are broken, buried prolines introduced, buried charges are created or lost, charge is swapped, a buried glycine is replaced, or if the residue that would be removed is a proline in the cis configuration. Also, if the mutation site is in a binding pocket the number of pockets and their volumes are reported. The user can assess this information and then select from available experimental or computationally predicted structures of native proteins to create, visualize, and download a model of the mutated protein using Fast and Accurate Side-chain Protein Repacking (FASPR). For AlphaFold modeled proteins, confidence scores for native proteins are provided. Using this tool, we explored a set of 9,666 common missense mutations from a variety of tissues from GTEx and show that most mutations can be modeled using this tool to facilitate studies of protein-protein and protein-drug interactions. The open-source tool is freely available at https://pharmacogenomics.clas.ucdenver.edu/gtexome/.

Illuminate the Functions of Dark Proteins Using the Reactome-IDG Web Portal.

Understudied or dark proteins have the potential to shed light on as-yet undiscovered molecular mechanisms that underlie phenotypes and suggest innovative therapeutic approaches for many diseases. The Reactome-IDG (Illuminating the Druggable Genome) project aims to place dark proteins in the context of manually curated, highly reliable pathways in Reactome, the most comprehensive, open-source biological pathway knowledgebase, facilitating the understanding functions and predicting therapeutic potentials of dark proteins. The Reactome-IDG web portal, deployed at https://idg.reactome.org, provides a simple, interactive web page for users to search pathways that may functionally interact with dark proteins, enabling the prediction of functions of dark proteins in the context of Reactome pathways. Enhanced visualization features implemented at the portal allow users to investigate the functional contexts for dark proteins based on tissue-specific gene or protein expression, drug-target interactions, or protein or gene pairwise relationships in the original Reactome's systems biology graph notation (SBGN) diagrams or the new simplified functional interaction (FI) network view of pathways. The protocols in this chapter describe step-by-step procedures to use the web portal to learn biological functions of dark proteins in the context of Reactome pathways. © 2023 Wiley Periodicals LLC. Basic Protocol 1: Search for interacting pathways of a protein Support Protocol: Interacting pathway results for an annotated protein Alternate Protocol: Use individual pairwise relationships to predict interacting pathways of a protein Basic Protocol 2: Using the IDG pathway browser to study interacting pathways Basic Protocol 3: Overlaying tissue-specific expression data Basic Protocol 4: Overlaying protein/gene pairwise relationships in the pathway context Basic Protocol 5: Visualizing drug/target interactions.

Reference genome of lumpfish Cyclopterus lumpus Linnaeus provides evidence of male heterogametic sex determination through the AMH pathway.

Teleosts exhibit extensive diversity of sex determination (SD) systems and mechanisms, providing the opportunity to study the evolution of SD and sex chromosomes. Here we sequenced the genome of the common lumpfish (Cyclopterus lumpus Linnaeus), a species of increasing importance to aquaculture, and identified the SD region and master SD locus using a 70K single nucleotide polymorphism array and tissue-specific expression data. The chromosome-level assembly identified 25 diploid chromosomes with a total size of 572.89Mb, a scaffold N50 of 23.86Mb and genome annotation-predicted 21,480 protein-coding genes. Genome-wide association analysis located a highly sex-associated region on chromosome 13, suggesting that anti-Müllerian hormone (AMH) is the putative SD factor. Linkage disequilibrium and heterozygosity across chromosome 13support a proto-XX/XY system, with an absence of widespread chromosome divergence between sexes. We identified three copies of AMH in the lumpfish primary and alternate haplotype assemblies localized in the SD region. Comparison to sequences from other teleosts suggested a monophyletic relationship and conservation within the Cottioidei. One AMH copy showed similarity to AMH/AMHY in a related species and was also the only copy with expression in testis tissue, suggesting this copy may be the functional copy of AMH in lumpfish. The two other copies arranged in tandem inverted duplication were highly similar, suggesting a recent duplication event. This study provides a resource for the study of early sex chromosome evolution and novel genomic resources that benefits lumpfish conservation management and aquaculture.

Genome-Wide Identification and Analysis of Cell Cycle Genes in Birch

Research Highlights: This study identified the cell cycle genes in birch that likely play important roles during the plant’s growth and development. This analysis provides a basis for understanding the regulatory mechanism of various cell cycles in Betula pendula Roth. Background and Objectives: The cell cycle factors not only influence cell cycles progression together, but also regulate accretion, division, and differentiation of cells, and then regulate growth and development of the plant. In this study, we identified the putative cell cycle genes in the B. pendula genome, based on the annotated cell cycle genes in Arabidopsis thaliana (L.) Heynh. It can be used as a basis for further functional research. Materials and Methods: RNA-seq technology was used to determine the transcription abundance of all cell cycle genes in xylem, roots, leaves, and floral tissues. Results: We identified 59 cell cycle gene models in the genome of B. pendula, with 17 highly expression genes among them. These genes were BpCDKA.1, BpCDKB1.1, BpCDKB2.1, BpCKS1.2, BpCYCB1.1, BpCYCB1.2, BpCYCB2.1, BpCYCD3.1, BpCYCD3.5, BpDEL1, BpDpa2, BpE2Fa, BpE2Fb, BpKRP1, BpKRP2, BpRb1, and BpWEE1. Conclusions: By combining phylogenetic analysis and tissue-specific expression data, we identified 17 core cell cycle genes in the Betulapendula genome.

LDexpress: an online tool for integrating population-specific linkage disequilibrium patterns with tissue-specific expression data

Genome-wide association studies have identified thousands of genetic susceptibility loci associated with cancer as well as other traits and diseases. Mapping germline variation in identified genetic susceptibility regions to alterations in nearby gene expression nominates candidate genes potentially related to disease risk for further functional investigation. We developed LDexpress as an online resource that integrates population-specific linkage disequilibrium data from the 1000 Genomes (1000G) project and tissue-specific expression data from the Genotype-Tissue Expression project to better study regional germline variation impacting gene expression. LDexpress is a publicly available web tool designed to be easy to use, flexible to conduct a wide range of variant queries, and quick to efficiently investigate dozens of query variants across multiple tissue types. We demonstrate the utility of LDexpress using example genomic queries and anticipate this tool will accelerate understanding of disease etiology by uncovering associations of regional germline variation to nearby gene expression.

Genome-wide identification of FRK genes in Populus trichocarpa and their expression under different nitrogen treatments.

The online version contains supplementary material available at 10.1007/s12298-021-01055-6.

Variable expression quantitative trait loci analysis of breast cancer risk variants

Genome wide association studies (GWAS) have identified more than 180 variants associated with breast cancer risk, however the underlying functional mechanisms and biological pathways which confer disease susceptibility remain largely unknown. As gene expression traits are under genetic regulation we hypothesise that differences in gene expression variability may identify causal breast cancer susceptibility genes. We performed variable expression quantitative trait loci (veQTL) analysis using tissue-specific expression data from the Genotype-Tissue Expression (GTEx) Common Fund Project. veQTL analysis identified 70 associations (p < 5 × 10–8) consisting of 60 genes and 27 breast cancer risk variants, including 55 veQTL that were observed in breast tissue only. Pathway analysis of genes associated with breast-specific veQTL revealed an enrichment of four genes (CYP11B1, CYP17A1 HSD3B2 and STAR) involved in the C21-steroidal biosynthesis pathway that converts cholesterol to breast-related hormones (e.g. oestrogen). Each of these four genes were significantly more variable in individuals homozygous for rs11075995 (A/A) breast cancer risk allele located in the FTO gene, which encodes an RNA demethylase. The A/A allele was also found associated with reduced expression of FTO, suggesting an epi-transcriptomic mechanism may underlie the dysregulation of genes involved in hormonal biosynthesis leading to an increased risk of breast cancer. These findings provide evidence that genetic variants govern high levels of expression variance in breast tissue, thus building a more comprehensive insight into the underlying biology of breast cancer risk loci.

Yellow nutsedge WRI4-like gene improves drought tolerance in Arabidopsis thaliana by promoting cuticular wax biosynthesis

BackgroundCuticular wax plays important role in protecting plants from drought stress. In Arabidopsis WRI4 improves drought tolerance by regulating the biosynthesis of fatty acids and cuticular wax. Cyperus esculentus (yellow nutsedge) is a tough weed found in tropical and temperate zones as well as in cooler regions. In the current study, we report the molecular cloning of a WRI4-like gene from Cyperus esculentus and its functional characterization in Arabidopsis.ResultsUsing RACE PCR, full-length WRI-like gene was amplified from yellow nutsedge. Phylogenetic analyses and amino acid comparison suggested it to be a WRI4-like gene. According to the tissue-specific expression data, the highest expression of WRI4-like gene was found in leaves, followed by roots and tuber. Transgenic Arabidopsis plants expressing nutsedge WRI4-like gene manifested improved drought stress tolerance. Transgenic lines showed significantly reduced stomatal conductance, transpiration rate, chlorophyll leaching, water loss and improved water use efficiency (WUE). In the absence of drought stress, expression of key genes for fatty acid biosynthesis was not significantly different between transgenic lines and WT while that of cuticular wax biosynthesis genes was significantly higher in transgenic lines than WT. The PEG-simulated drought stress significantly increased expression of key genes for fatty acid as well as wax biosynthesis in transgenic Arabidopsis lines but not in WT plants. Consistent with the gene expression data, cuticular wax load and deposition was significantly higher in stem and leaves of transgenic lines compared with WT under control as well as drought stress conditions.ConclusionsWRI4-like gene from Cyperus esculentus improves drought tolerance in Arabidopsis probably by promoting cuticular wax biosynthesis and deposition. This in turn lowers chlorophyll leaching, stomatal conductance, transpiration rate, water loss and improves water use efficiency under drought stress conditions. Therefore, CeWRI4-like gene could be a good candidate for improving drought tolerance in crops.

Jasmonate and aluminum crosstalk in tomato: Identification and expression analysis of WRKYs and ALMTs during JA/Al-regulated root growth

The WRKY transcription factors (TFs) are involved in aluminum (Al) stress and jasmonic acid (JA)-regulated resistance responses. WRKYs act as regulators of Al-activated malate transporter (ALMT) proteins (anion channels) by directly binding to their promoters and altering malate efflux, thereby regulating Al ion toxicity in plant roots. JA enhances Al-induced root growth inhibition in Arabidopsis. However, the relationship between WRKY and ALMT genes and their involvement in JA-mediated root growth inhibition during Al stress in tomato remain unknown. Here, we demonstrate a similar phenomenon that JA enhances Al-induced root growth inhibition in tomato (Solanum lycopersicum). By analyzing RNA-seq data and tissue-specific expression data from public databases, we selected 17 WRKY and 6 ALMT family genes to identify the genes participated in this process. The promoters of many of the selected genes contained MeJA responsive element, G-box (target site of MYC2, a core TF of JA signaling), and W-box (target site for WRKY). Quantitative real-time PCR was performed to evaluate the expression levels of selected WRKY and ALMT genes under AlCl3 and Methyl jasmonate (MeJA) treatment. SlMYC2-VIGS seedlings and jasmonic acid-insensitive1 (jai1) mutant were also employed to analyze the expression patterns of selected genes. We find that SlALMT3 is responsible for the crosstalk regulatory mechanism between Al and JA in root growth inhibition, and 6 SlWRKYs may act as the upstream regulators of SlALMT3 in this crosstalk response. This study is initial and informative in exploring the crosstalk regulatory mechanism between JA and Al in tomato.

Identification of putative lignin biosynthesis genes in Betula pendula

We identified 15 genes encoding enzymes for phenylpropanoid biosynthesis in the genome of birch by combining phylogenetic analysis and tissue-specific expression data. Lignin is one of the most abundant terrestrial biopolymers and is essential for plant structure and defense. An essential step in lignin formation is phenylpropanoid synthesis, we identified 120 gene models in 10 protein families encoding enzymes for phenylpropanoid biosynthesis in the genome of Betula pendula. The transcript abundance was determined for all 120 genes in xylem, root, leaf, and flower tissues using RNA-seq technology. We identified 15 genes that likely encode phenylpropanoid biosynthesis enzymes during wood formation. Ten of these genes are evolutionarily conserved compared to the lignin genes in Populus trichocarpa.

Genetic Determinants of Reduced Arsenic Metabolism Efficiency in the 10q24.32 Region Are Associated With Reduced AS3MT Expression in Multiple Human Tissue Types.

Approximately 140 million people worldwide are exposed to inorganic arsenic through contaminated drinking water. Chronic exposure increases risk for cancers as well as cardiovascular, respiratory, and neurologic diseases. Arsenic metabolism involves the AS3MT (arsenic methyltransferase) gene, and arsenic metabolism efficiency (AME, measured as relative concentrations of arsenic metabolites in urine) varies among individuals. Inherited genetic variation in the 10q24.32 region, containing AS3MT, influences AME, but the mechanisms remain unclear. To better understand these mechanisms, we use tissue-specific expression data from GTEx (Genotype-tissue Expression project) to identify cis-eQTLs (expression quantitative trait loci) for AS3MT and other nearby genes. We combined these data with results from a genome-wide association study of AME using "colocalization analysis," to determine if 10q24.32 SNPs (single nucleotide polymorphisms) that affect AME also affect expression of AS3MT or nearby genes. These analyses identified cis-eQTLs for AS3MT in 38 tissue types. Colocalization results suggest that the casual variant represented by AME lead SNP rs4919690 impacts expression of AS3MT in 13 tissue types (> 80% probability). Our results suggest this causal SNP also regulates/coregulates expression of nearby genes: BORCS7 (43 tissues), NT5C2 (2 tissues), CYP17A1-AS1 (1 tissue), and RP11-724N1.1 (1 tissue). The rs4919690 allele associated with decreased AME is associated with decreased expression of AS3MT (and other coregulated genes). Our study provides a potential biological mechanism for the association between 10q24.32 variation and AME and suggests that the causal variant, represented by rs4919690, may impact AME (as measured in urine) through its effects on arsenic metabolism occurring in multiple tissue types.

Genome-wide identification of BXL genes in Populus trichocarpa and their expression under different nitrogen treatments.

β-d-xylosidase (BXL) hydrolyzes xylobiose and xylo-oligosaccharides into xylose monomers, and is a rate-limiting enzyme in the degradation of hemicellulose in the cell wall. In this study, ten genes encoding putative BXL proteins were identified in the Populus trichocarpa genome by bioinformatics methods. In the phylogenetic analysis, the PtBXLs formed two subfamilies. PtBXL8 and PtBXL9 were closely related to AtBXL1, an important enzyme in the normal development of the Arabidopsis cell wall structure. Chromosomal distribution and genome synteny analyses revealed two tandem-duplicated gene pairs PtBXL3/4 and PtBXL6/7 on chromosomes II and V, respectively, and six segmental-duplicated gene pairs on chromosomes II, V, VIII, X, and XIV among the PtBXL gene family. Tissue-specific expression data from PlantGenIE indicated that PtBXL2, 4, 5, and 10 were highly expressed in stems. Quantitative real-time RT-PCR analyses revealed that PtBXL4, 5, and 9 were up-regulated in the upper stem in response to the low and high ammonium and nitrate treatments. The influence of nitrogen on the expression of PtBXL4, 5, and 9 genes may affect the formation of the plant secondary cell wall. This comprehensive analysis of the BXL family in poplar provides new insights into their regulation by nitrogen and increases our understanding of the roles of BXLs in hemicellulose metabolism in the secondary cell wall and during plant development.

Genome-wide identification of GRF transcription factors in soybean and expression analysis of GmGRF family under shade stress

BackgroundThe Growth-regulating factor (GRF) family encodes plant-specific transcription factors which contain two conserved domains, QLQ and WRC. Members of this family play vital roles in plant development and stress response processes. Although GRFs have been identified in various plant species, we still know little about the GRF family in soybean (Glycine max).ResultsIn the present study, 22 GmGRFs distributed on 14 chromosomes and one scaffold were identified by searching soybean genome database and were clustered into five subgroups according to their phylogenetic relationships. GmGRFs belonging to the same subgroup shared a similar motif composition and gene structure. Synteny analysis revealed that large-scale duplications played key roles in the expansion of the GmGRF family. Tissue-specific expression data showed that GmGRFs were strongly expressed in growing tissues, including the shoot apical meristems, developing seeds and flowers, indicating that GmGRFs play critical roles in plant growth and development. On the basis of expression analysis of GmGRFs under shade conditions, we found that all GmGRFs responded to shade stress. Most GmGRFs were down-regulated in soybean leaves after shade treatment.ConclusionsTaken together, this research systematically analyzed the characterization of the GmGRF family and its primary roles in soybean development and shade stress response. Further studies of the function of the GmGRFs in the growth, development and stress tolerance of soybean, especially under shade stress, will be valuable.

LncRRIsearch: A Web Server for lncRNA-RNA Interaction Prediction Integrated With Tissue-Specific Expression and Subcellular Localization Data.

Long non-coding RNAs (lncRNAs) play critical roles in various biological processes, but the function of the majority of lncRNAs is still unclear. One approach for estimating a function of a lncRNA is the identification of its interaction target because functions of lncRNAs are expressed through interaction with other biomolecules in quite a few cases. In this paper, we developed “LncRRIsearch,” which is a web server for comprehensive prediction of human and mouse lncRNA-lncRNA and lncRNA-mRNA interaction. The prediction was conducted using RIblast, which is a fast and accurate RNA-RNA interaction prediction tool. Users can investigate interaction target RNAs of a particular lncRNA through a web interface. In addition, we integrated tissue-specific expression and subcellular localization data for the lncRNAs with the web server. These data enable users to examine tissue-specific or subcellular localized lncRNA interactions. LncRRIsearch is publicly accessible at http://rtools.cbrc.jp/LncRRIsearch/.

Tissue- and sex-specific small RNAomes reveal sex differences in response to the environment.

RNA interference (RNAi) related pathways are essential for germline development and fertility in metazoa and can contribute to inter- and trans-generational inheritance. In the nematode Caenorhabditis elegans, environmental double-stranded RNA provided by feeding can lead to heritable changes in phenotype and gene expression. Notably, transmission efficiency differs between the male and female germline, yet the underlying mechanisms remain elusive. Here we use high-throughput sequencing of dissected gonads to quantify sex-specific endogenous piRNAs, miRNAs and siRNAs in the C. elegans germline and the somatic gonad. We identify genes with exceptionally high levels of secondary 22G RNAs that are associated with low mRNA expression, a signature compatible with silencing. We further demonstrate that contrary to the hermaphrodite germline, the male germline, but not male soma, is resistant to environmental RNAi triggers provided by feeding, in line with previous work. This sex-difference in silencing efficacy is associated with lower levels of gonadal RNAi amplification products. Moreover, this tissue- and sex-specific RNAi resistance is regulated by the germline, since mutant males with a feminized germline are RNAi sensitive. This study provides important sex- and tissue-specific expression data of miRNA, piRNA and siRNA as well as mechanistic insights into sex-differences of gene regulation in response to environmental cues.

Meta-Analysis of Genetic Influences on Initial Alcohol Sensitivity.

Previous studies indicate that low initial sensitivity to alcohol may be a risk factor for later alcohol misuse. Evidence suggests that initial sensitivity is influenced by genetic factors, but few molecular genetic studies have been reported. We conducted a meta-analysis of 2 population-based genome-wide association studies of the Self-Rating of the Effects of Alcohol scale. Our final sample consisted of 7,339 individuals (82.3% of European descent; 59.2% female) who reported having used alcohol at least 5 times. In addition, we estimated single nucleotide polymorphism (SNP)-based heritability and conducted a series of secondary aggregate genetic analyses. No individual locus reached genome-wide significance. Gene and set based analyses, both overall and using tissue-specific expression data, yielded largely null results, and genes previously implicated in alcohol problems and consumption were overall not associated with initial sensitivity. Only 1 gene set, related to hormone signaling and including core clock genes, survived correction for multiple testing. A meta-analysis of SNP-based heritability resulted in a modest estimate of =0.19 (SE=0.10), though this was driven by 1 sample (N=3,683, =0.36, SE=0.14, p=0.04). No significant genetic correlations with other relevant outcomes were observed. Findings yielded only modest support for a genetic component underlying initial alcohol sensitivity. Results suggest that its biological underpinnings may diverge somewhat from that of other alcohol outcomes and may be related to core clock genes or other aspects of hormone signaling. Larger samples, ideally of prospectively assessed samples, are likely necessary to improve gene identification efforts and confirm the current findings.

Genome-wide characterization, evolution, and expression analysis of the leucine-rich repeat receptor-like protein kinase (LRR-RLK) gene family in Rosaceae genomes

BackgroundLeucine-rich repeat receptor-like protein kinase (LRR-RLK) is the largest gene family of receptor-like protein kinases (RLKs) and actively participates in regulating the growth, development, signal transduction, immunity, and stress responses of plants. However, the patterns of LRR-RLK gene family evolution in the five main Rosaceae species for which genome sequences are available have not yet been reported. In this study, we performed a comprehensive analysis of LRR-RLK genes for five Rosaceae species: Fragaria vesca (strawberry), Malus domestica (apple), Pyrus bretschneideri (Chinese white pear), Prunus mume (mei), and Prunus persica (peach), which contained 201, 244, 427, 267, and 258 LRR-RLK genes, respectively.ResultsAll LRR-RLK genes were further grouped into 23 subfamilies based on the hidden Markov models approach. RLK-Pelle_LRR-XII-1, RLK-Pelle_LRR-XI-1, and RLK-Pelle_LRR-III were the three largest subfamilies. Synteny analysis indicated that there were 236 tandem duplicated genes in the five Rosaceae species, among which subfamilies XII-1 (82 genes) and XI-1 (80 genes) comprised 68.6%.ConclusionsOur results indicate that tandem duplication made a large contribution to the expansion of the subfamilies. The gene expression, tissue-specific expression, and subcellular localization data revealed that LRR-RLK genes were differentially expressed in various organs and tissues, and the largest subfamily XI-1 was highly expressed in all five Rosaceae species, suggesting that LRR-RLKs play important roles in each stage of plant growth and development. Taken together, our results provide an overview of the LRR-RLK family in Rosaceae genomes and the basis for further functional studies.

Genome-wide association mapping of QTL underlying seed oil and protein contents of a diverse panel of soybean accessions

To investigate the genetic basis of variation in oil and protein contents in soybean seeds, a diverse collection of 421 mainly Chinese soybean cultivars was genotyped using 1536 SNPs, mostly from candidate genes related to acyl-lipid metabolism and from regions harboring known QTL. Six significant associations were identified for each of seed oil and protein contents which individually explained 2.7–5.9% of the phenotypic variance. Six associations occurred in or near known QTL and the remaining are putative novel QTL. Ten significant associations influenced the oil content without decreasing protein content, and vice versa. One SNP was pleiotropic, with opposite effects on oil and protein contents. The genetic region covering Map-6076 and-6077 was shown to be involved in controlling oil content in soybean by integrating the results of association mapping with information on known QTL and tissue-specific expression data. This region was subject to strong selection during the genetic improvement of soybean. Our results not only confirm and refine the map positions of known QTL but also contribute to a further elucidation of the genetic architecture of protein and oil contents in soybean seeds by identifying new associations exhibiting pleiotropic effects on seed protein and oil contents.