The Genotype-Tissue Expression (GTEx) project was established to characterize genetic effects on the transcriptome across human tissues and to link these regulatory mechanisms to trait and disease associations. Here, we present analyses of the version 8 data, examining 15,201 RNA-sequencing samples from 49 tissues of 838 postmortem donors. We comprehensively characterize genetic associations for gene expression and splicing in cis and trans, showing that regulatory associations are found for almost all genes, and describe the underlying molecular mechanisms and their contribution to allelic heterogeneity and pleiotropy of complex traits. Leveraging the large diversity of tissues, we provide insights into the tissue specificity of genetic effects and show that cell type composition is a key factor in understanding gene regulatory mechanisms in human tissues.

The capacity of conventional breeding to simultaneously improve the yield and quality of cotton fiber is limited. The accumulation of the plant hormone indole-3-acetic acid (IAA) in cotton fiber initials prompted us to investigate the effects of genetically engineering increased IAA levels in the ovule epidermis. Targeted expression of the IAA biosynthetic gene iaaM, driven by the promoter of the petunia MADS box gene Floral Binding protein 7 (FBP7), increased IAA levels in the epidermis of cotton ovules at the fiber initiation stage. This substantially increased the number of lint fibers, an effect that was confirmed in a 4-year field trial. The lint percentage of the transgenic cotton, an important component of fiber yield, was consistently higher in our transgenic plants than in nontransgenic controls, resulting in a >15% increase in lint yield. Fiber fineness was also notably improved.

Highlights d Human germline cell specification begins from a transitional pluripotent state d Human primordial germ cells are specified from lineageprimed progenitors d Lineage-primed TFAP2A progenitors have gastrulating and amnion cell identity d TFAP2C regulates SOX17 at the point of human primordial germ cell specification

Gene co-expression networks capture biologically important patterns in gene expression data, enabling functional analyses of genes, discovery of biomarkers, and interpretation of genetic variants. Most network analyses to date have been limited to assessing correlation between total gene expression levels in a single tissue or small sets of tissues. Here, we built networks that additionally capture the regulation of relative isoform abundance and splicing, along with tissue-specific connections unique to each of a diverse set of tissues. We used the Genotype-Tissue Expression (GTEx) project v6 RNA sequencing data across 50 tissues and 449 individuals. First, we developed a framework called Transcriptome-Wide Networks (TWNs) for combining total expression and relative isoform levels into a single sparse network, capturing the interplay between the regulation of splicing and transcription. We built TWNs for 16 tissues and found that hubs in these networks were strongly enriched for splicing and RNA binding genes, demonstrating their utility in unraveling regulation of splicing in the human transcriptome. Next, we used a Bayesian biclustering model that identifies network edges unique to a single tissue to reconstruct Tissue-Specific Networks (TSNs) for 26 distinct tissues and 10 groups of related tissues. Finally, we found genetic variants associated with pairs of adjacent nodes in our networks, supporting the estimated network structures and identifying 20 genetic variants with distant regulatory impact on transcription and splicing. Our networks provide an improved understanding of the complex relationships of the human transcriptome across tissues.

DNA methylation plays important roles in genome protection, regulation of gene expression and is associated with plants development. Plant DNA methylation pattern was mediated by cytosine-5 DNA methyltransferase and demethylase. Although the genomes of AA and BB wild peanuts have been fully sequenced, these two gene families have not been studied. In this study we report the identification and analysis of putative cytosine-5 DNA methyltransferases (C5-MTases) and demethylases in AA and BB wild peanuts. Cytosine-5 DNA methyltransferases in AA and BB wild peanuts could be classified in MET, CMT, and DRM2 groups based on their domain organization. This result was supported by the gene and protein structural characteristics and phylogenetic analysis. We found that some wild peanut DRM2 members didn't contain UBA domain which was different from other plants such as Arabidopsis, maize and soybean. Five DNA demethylase encoding genes were found in AA genome and five in BB genome. The selective pressure analysis showed that wild peanut C5-MTase genes mainly underwent purifying selection but many positive selection sites can be detected. Conversely, DNA demethylase genes mainly underwent positive selection during evolution. Additionally, the expression dynamic of cytosine-5 DNA methyltransferase and demethylase genes in different cultivated peanut tissues were analyzed. Expression result showed that cold, heat or PEG stress could influence the expression level of C5-MTase and DNA demethylase genes in cultivated peanut. These results are useful for better understanding the complexity of these two gene families, and will facilitate epigenetic studies in peanut in the future.

SummaryCotton (Gossypium hirsutum L.) fibers, one of the most important natural raw materials for textile industry, are highly elongated trichomes from epidermal cells of cotton ovules. DET2, an Arabidopsis steroid 5d-reductase, is considered to catalyze a major rate-limiting in brassinosteroid (BR) biosynthesis. To understand the role of BRs in cotton fiber development, GhDET2, which putatively encodes a steroid 5a-reductase by sequence comparison, was cloned from developing fiber cells. In vitro assessment of GhDET2 protein activity confirmed that GhDET2 encodes a functional steroid 5a-redutase. High levels of GhDET2 transcript were detected during the fiber initiation stage and the fiber rapid elongation stage. Antisense-mediated suppression of GhDET2 inhibited both fiber initiation and fiber elongation. Similarly, treating cultured ovules with finasteride, a steroid 5a-reductase inhibitor, reduced fiber elongation. Inhibition of fiber cell elongation by expression of antisense GhDET2 or the finasteride treatment could be reversed by epibrassinolide, a biologically active BR. Furthermore, seed coat-specific expression of GhDET2 increased fiber number and length. Therefore, GhDET2 and BRs play a crucial role in the initiation and elongation of cotton fiber cells, suggesting that modulation of BR biosynthesis factors may improve fiber quality or yield.

Ammonia oxidation to nitrite and its subsequent oxidation to nitrate provides energy to the two populations of nitrifying chemoautotrophs in the energy-starved dark ocean, driving a coupling between reduced inorganic nitrogen (N) pools and production of new organic carbon (C) in the dark ocean. However, the relationship between the flux of new C production and the fluxes of N of the two steps of oxidation remains unclear. Here, we show that, despite orders-of-magnitude difference in cell abundances between ammonia oxidizers and nitrite oxidizers, the two populations sustain similar bulk N-oxidation rates throughout the deep waters with similarly high affinities for ammonia and nitrite under increasing substrate limitation, thus maintaining overall homeostasis in the oceanic nitrification pathway. Our observations confirm the theoretical predictions of a redox-informed ecosystem model. Using balances from this model, we suggest that consistently low ammonia and nitrite concentrations are maintained when the two populations have similarly high substrate affinities and their loss rates are proportional to their maximum growth rates. The stoichiometric relations between the fluxes of C and N indicate a threefold to fourfold higher C-fixation efficiency per mole of N oxidized by ammonia oxidizers compared to nitrite oxidizers due to nearly identical apparent energetic requirements for C fixation of the two populations. We estimate that the rate of chemoautotrophic C fixation amounts to ∼1 × 1013to ∼2 × 1013mol of C per year globally through the flux of ∼1 × 1014to ∼2 × 1014mol of N per year of the two steps of oxidation throughout the dark ocean.

WRKY proteins are plant specific transcription factors involved in various developmental and physiological processes, especially in biotic and abiotic stress resistance. Although previous studies suggested that WRKY proteins in soybean (Glycine max var. Williams 82) involved in both abiotic and biotic stress responses, the global information of WRKY proteins in the latest version of soybean genome (Wm82.a2v1) and their response to dehydration and salt stress have not been reported. In this study, we identified 176 GmWRKY proteins from soybean Wm82.a2v1 genome. These proteins could be classified into three groups, namely group I (32 proteins), group II (120 proteins), and group III (24 proteins). Our results showed that most GmWRKY genes were located on Chromosome 6, while chromosome 11, 12, and 20 contained the least number of this gene family. More GmWRKY genes were distributed on the ends of chromosomes to compare with other regions. The cis-acting elements analysis suggested that GmWRKY genes were transcriptionally regulated upon dehydration and salt stress. RNA-seq data analysis indicated that three GmWRKY genes responded negatively to dehydration, and 12 genes positively responded to salt stress at 1, 6, and 12 h, respectively. We confirmed by qRT-PCR that the expression of GmWRKY47 and GmWRKY 58 genes was decreased upon dehydration, and the expression of GmWRKY92, 144 and 165 genes was increased under salt treatment.