Identification of Candidate Genes Associated With Hypoxia Tolerance in Trachinotus blochii Using Bulked Segregant Analysis and RNA-Seq

Golden Pompano (Trachinotus blochii) has rapidly developed into the one of the main valuable fish species in Chinese marine aquaculture. Due to its rapid growth, active metabolism, and high oxygen consumption, hypoxia will increase its mortality and cause serious economic losses. We constructed two experimental groups of fish with different degrees of tolerance to hypoxia, used BSR-Seq analysis based on genome and genetic linkage groups to locate SNPs and genes that were related to the differences in hypoxia tolerance. The results showed that hypoxia tolerance SNPs of golden pompano may be jointly determined by multiple linkage groups, especially linkage groups 18 and 22. There were 768 and 348 candidate genes located in the candidate regions of the brain and liver, respectively. These genes were mainly involved in anaerobic energy metabolism, stress response, immune response, waste discharge, and cell death. The prostaglandin-endoperoxide synthase 2 (PTGS2) on LG8, which is involved in the metabolism of arachidonic acid, has a G/A nonsynonymous mutation at position 20641628, and the encoded amino acid was changed from hydrophobic aspartic acid to asparaginate. The specific pathway of the RIG-I-like receptor signaling pathway in the liver may mediate the metabolic system and the immune system, linking glucose metabolism with immune regulation. The death of the hypoxia-intolerant group may be due to the accumulation of lactic acid caused by the activation of anaerobic glycolysis during the early stage of hypoxia stress, and the activation of type I interferon was inhibited, which resulted in decreased immunity. Among the genes involved in the RIG-I-like receptor signaling pathway, the CYLD Lysine 63 Deubiquitinase (CYLD) located on LG16 had a G/T nonsynonymous mutation at position 13629651, and the encoded amino acid was changed from alanine acid to valine. The interferon induced with helicase C domain 1 (Ifih1) located on LG18 has a G/C nonsynonymous mutation at position 16153700, and the encoded hydrophilic glycine was changed to hydrophobic alanine. Our findings suggest these SNPs may assist in the molecular breeding of hypoxia-tolerant golden pompano, and speculate that the balance of glucose and lipid metabolism plays a key role in Trachinotus blochii under acute hypoxia.

Selection on BrFLC1 Is Related to Intraspecific Diversity of Brassica rapa Vegetables

Flowering time is important for Brassica rapa vegetables because premature bolting before harvest can lower yield and quality. FLOWERING LOCUS C (FLC) acts as a key repressor of flowering. In this study, we identified a nonsynonymous mutation at the 58th nucleotide of exon1 in BrFLC1 (named as Pe1+58 (A/C)) by screening resequencing data of 199 B. rapa accessions and verified this mutation as being related to flowering time variation. Strong linkage inheritance was detected between this locus and a previously reported splicing site mutation at intron 6 of BrFLC1 (Pi6+1 (G/A)), showing as co-occurrence of BrFLC1Pe1+58(A) and BrFLC1Pi6+1(G), named as haplotype H1: AG, or co-occurrence of BrFLC1Pe1+58(C) and BrFLC1Pi6+1(A), named as haplotype H2: CA. The frequency distribution of BrFLC1 haplotypes skewed to the haplotype H1 in turnip, broccoletto, mizuna, komatsuna, and taicai, while it was skewed to the haplotype H2 in caixin, pak choi, zicaitai, and wutacai. The frequencies of the two haplotypes were comparable in Chinese cabbage. This indicated that BrFLC1 haplotypes were related to B. rapa intraspecific diversification. Further analysis of a Chinese cabbage collection revealed that accessions from the spring ecotype preferred to keep H1: AG and almost all accessions from the summer ecotype were H2: CA. The early flowering haplotype of BrFLC1 was purified in summer Chinese cabbage, indicating that BrFLC1 had been strongly selected during genetic improvement of summer Chinese cabbages. A significant difference in flowering time of F2 individuals with the homologous BrFLC1Pi6+1(G) allele but different BrFLC1Pe1+58 (A/C) alleles, indicated that this locus had independent genetic effects on flowering time. The newly identified allelic diversity of BrFLC1 can be used for breeding of resistance to premature bolting in B. rapa vegetables.

A single amino acid residue substitution in BraA04g017190.3C, a histone methyltransferase, results in premature bolting in Chinese cabbage (Brassica rapa L. ssp. Pekinensis)

Background Flowering is an important inflection point in the transformation from vegetative to reproductive growth, and premature bolting severely decreases crop yield and quality. Results In this study, a stable early-bolting mutant, ebm3, was identified in an ethyl methanesulfonate (EMS)-mutagenized population of a Chinese cabbage doubled haploid (DH) line ‘FT’. Compared with ‘FT’, ebm3 showed early bolting under natural cultivation in autumn, and curled leaves. Genetic analysis showed that the early-bolting phenotype was controlled by a single recessive nuclear gene. Modified MutMap sequencing, genotyping analyses and allelism test provide strong evidence that BrEBM3 (BraA04g017190.3 C), encoding the histone methyltransferase CURLY LEAF (CLF), was the strongly candidate gene of the emb3. A C to T base substitution in the 14th exon of BrEBM3 resulted in an amino acid change (S to F) and the early-bolting phenotype of emb3. The mutation occurred in the SET domain (Suppressor of protein-effect variegation 3–9, Enhancer-of-zeste, Trithorax), which catalyzes site- and state-specific lysine methylation in histones. Tissue-specific expression analysis showed that BrEBM3 was highly expressed in the flower and bud. Promoter activity assay confirmed that BrEBM3 promoter was active in inflorescences. Subcellular localization analysis revealed that BrEBM3 localized in the nucleus. Transcriptomic studies supported that BrEBM3 mutation might repress H3K27me3 deposition and activate expression of the AGAMOUS (AG) and AGAMOUS-like (AGL) loci, resulting in early flowering. Conclusions Our study revealed that an EMS-induced early-bolting mutant ebm3 in Chinese cabbage was caused by a nonsynonymous mutation in BraA04g017190.3 C, encoding the histone methyltransferase CLF. These results improve our knowledge of the genetic and genomic resources of bolting and flowering, and may be beneficial to the genetic improvement of Chinese cabbage.

Dissection of PIK3CA Aberration for Cervical Adenocarcinoma Outcomes

Personalized treatment of genetically stratified subgroups has the potential to improve outcomes in many malignant tumors. This study distills clinically meaningful prognostic/predictive genomic marker for cervical adenocarcinoma using signature genomic aberrations and single-point nonsynonymous mutation-specific droplet digital PCR (ddPCR). Mutations in PIK3CA E542K, E545K, or H1047R were detected in 41.7% of tumors. PIK3CA mutation detected in the patient’s circulating DNA collected before treatment or during follow-up was significantly associated with decreased progression-free survival or overall survival. PIK3CA mutation in the circulating DNA during follow-up after treatment predicted recurrence with 100% sensitivity and 64.29% specificity. It is the first indication of the predictive power of PIK3CA mutations in cervical adenocarcinoma. The work contributes to the development of liquid biopsies for follow up surveillance and a possibility of tailoring management of this particular women’s cancer.

Comparative transcriptomic analysis reveals the regulatory mechanism of the gibberellic acid pathway of Tartary buckwheat (Fagopyrum tataricum (L.) Gaertn.) dwarf mutants

BackgroundTartary buckwheat is an important minor crop species with high nutritional and medicinal value and is widely planted worldwide. Cultivated Tartary buckwheat plants are tall and have hollow stems that lodge easily, which severely affects their yield and hinders the development of the Tartary buckwheat industry.MethodsHeifeng No. 1 seeds were treated with ethylmethanesulfonate (EMS) to generate a mutant library. The dwarf mutantftdmwas selected from the mutagenized population, and the agronomic characteristics giving rise to the dwarf phenotype were evaluated. Ultra-fast liquid chromatography-electrospray ionization tandem mass spectrometry (UFLC-ESI–MS/MS) was performed to determine the factors underlying the different phenotypes between the wild-type (WT) andftdmplants. In addition, RNA sequencing (RNA-seq) was performed via the HiSeq 2000 platform, and the resulting transcriptomic data were analysed to identify differentially expressed genes (DEGs). Single-nucleotide polymorphism (SNP) variant analysis revealed possible sites associated with dwarfism. The expression levels of the potential DEGs between the WT andftdmmutant were then measured via qRT-PCR and fragments per kilobase of transcript per million mapped reads (FPKM).ResultThe plant height (PH) of theftdmmutant decreased to 42% of that of the WT, and compared with the WT, the mutant and had a higher breaking force (BF) and lower lodging index (LI). Lower GA4 and GA7 contents and higher contents of jasmonic acid (JA), salicylic acid (SA) and brassinolactone (BR) were detected in the stems of theftdmmutant compared with the WT. Exogenous application of GAs could not revert the dwarfism of theftdmmutant. On the basis of the transcriptomic analysis, 146 homozygous SNP loci were identified. In total, 12 DEGs with nonsynonymous mutations were ultimately identified, which were considered potential candidate genes related to the dwarf trait. When the sequences of eight genes whose expression was downregulated and four genes whose expression was upregulated were compared, SKIP14, an F-box protein whose sequence is 85% homologous to that of SLY1 in Arabidopsis, presented an amino acid change (from Ser to Asn) and was expressed at a lower level in the stems of theftdmmutant compared with the WT. Hence, we speculated that this amino acid change in SKIP14 resulted in a disruption in GA signal transduction, indirectly decreasing the GA content and downregulating the expression of genes involved in GA biosynthesis or the GA response. Further studies are needed to determine the molecular basis underlying the dwarf phenotype of theftdmmutant.ConclusionWe report a Tartary buckwheat EMS dwarf mutant,ftdm, suitable for high-density planting and commercial farming. A significant decrease in GA4 and GA7 levels was detected in theftdmmutant, and 12 DEGs expressed in the stems of theftdmmutant were selected as candidates of the dwarfing gene. One nonsynonymous mutation was detected in theSKIP14gene in theftdmmutant, and this gene had a lower transcript level compared with that in the WT.

Universal constraints on protein evolution in the long-term evolution experiment with Escherichia coli

Although it is well known that abundant proteins evolve slowly across the tree of life, there is little consensus for why this is true. Here, I report that abundant proteins evolve slowly in the hypermutator populations of Lenski’s long-term evolution experiment with Escherichia coli (LTEE). Specifically, the density of all observed mutations per gene, as measured in metagenomic time series covering 60,000 generations of the LTEE, significantly anti-correlates with mRNA abundance, protein abundance, and degree of protein-protein interaction. The same pattern holds for nonsynonymous mutation density. However, synonymous mutation density, measured across the LTEE hypermutator populations, positively correlates with protein abundance. These results show that universal constraints on protein evolution are visible in data spanning three decades of experimental evolution. Therefore, it should be possible to design experiments to answer why abundant proteins evolve slowly.

The role of synonymous mutations in the evolution of TEM β-lactamase genes

Nonsynonymous mutations are well documented in TEM β-lactamases. The resulting amino acid changes often alter the conferred phenotype from broad spectrum (2b) conferred by TEM-1 to extended spectrum (2be), inhibitor resistant (2br), or both extended spectrum and inhibitor resistant (2ber). The encoding blaTEM genes also deviate in numerous synonymous mutations, which are not well understood. blaTEM-3 (2be), blaTEM-33 (2br), and blaTEM-109 (2ber) were studied in comparison to blaTEM-1. blaTEM-33 was chosen for more detailed studies, because it deviates from blaTEM-1 by a single nonsynonymous mutation and three additional, synonymous mutations. Genes encoding the enzymes with only nonsynonymous or all, including synonymous, mutations plus all permutations between blaTEM-1 and blaTEM-33 were expressed in Escherichia coli cells. In disc diffusion assays, genes encoding TEM-3, TEM-33, and TEM-109 with all synonymous mutations resulted in higher resistance levels than genes without synonymous mutations. Disc diffusion assays with the 16 genes carrying all possible nucleotide change combinations between blaTEM-1 and blaTEM-33 indicated different susceptibilities for different variants. Nucleotide BLAST searches did not identify genes without synonymous mutations but some without nonsynonymous mutations. Energies of possible secondary mRNA structures calculated with mfold are generally higher with synonymous mutations, suggesting that their role could be to destabilize the mRNA and facilitate its unfolding for efficient translation. In summary, our data indicate that transitions from blaTEM-1 to other variant genes by simply acquiring the nonsynonymous mutations is not favored. Instead, synonymous mutations seem to support the transition to other variant genes with nonsynonymous mutations leading to different phenotypes.

Tumor hepatitis B virus RNA identifies a clinically and molecularly distinct subset of hepatocellular carcinoma

Hepatitis B virus (HBV) infection contributes to hepatocellular carcinoma (HCC) initiation and is associated with worse outcomes. Many prior studies of HBV-related HCC have not accounted for potential heterogeneity among HBV-related tumors by assessing whether HBV activity is present in tumor tissue. Here, we measured tumor HBV RNA, a proxy for viral activity, and investigated the association between HBV RNA status and several clinicogenomic characteristics. We obtained clinical, mutation, RNA-Seq and survival data for 439 HCC tumors from The Cancer Genome Atlas and International Cancer Genome Consortium. Tumors were classified as HBV RNA positive if they harbored >1 HBV RNA read per million human reads. We investigated the association between HBV RNA status and nonsynonymous somatic mutations, gene set expression, homologous recombination deficiency (HRD) score and mutation-specific survival. HBV RNA positive status was associated with higher nonsynonymous mutation rates of multiple genes, including TP53 and CDKN2A, while HBV RNA negative status was associated with higher nonsynonymous BAP1 mutation rate. HBV RNA positive status was also associated with increased transcription of genes involved in multiple DNA damage repair pathways, genes upregulated by MYC and mTORC1, and genes overexpressed in several HCC subclasses associated with a proliferative phenotype. Further, HBV RNA positive status was associated with increased three-biomarker HRD score (22.2 for HBV RNA+ vs. 16.0 for HBV RNA-). Finally, HBV RNA status was associated with multiple mutation-specific survival differences, including decreased survival for HBV RNA positive patients with nonsynonymous KEAP1 mutations compared to those without (hazard ratio 4.26). HCC tumors harboring genomic evidence of HBV activity therefore constitute a distinct HCC subset characterized by specific differences in nonsynonymous mutations, gene set expression, three-biomarker HRD score and mutation-specific survival.

A Single Amino Acid Change In Histone Methyltransferase CURLY LEAF Results In Premature Bolting In Chinese Cabbage (Brassica Rapa L. Ssp. Pekinensis)

Background: Flowering is an important inflection point in the transformation from vegetative to reproductive growth, and premature bolting severely decreases crop yield and quality. Results: In this study, a stable early-bolting mutant, ebm3, was identified in an ethyl methanesulfonate (EMS)-mutagenized population of a Chinese cabbage doubled haploid (DH) line ‘FT’. Compared with ‘FT’, ebm3 showed early bolting under natural cultivation in autumn, and curled leaves. Genetic analysis showed that the early-bolting phenotype was controlled by a single recessive nuclear gene. Modified MutMap and genotyping analyses revealed that Brebm3 (BraA04g017190.3C), encoding the histone methyltransferase CURLY LEAF (CLF), was the causal gene of the emb3. A C to T base substitution in the 14th exon of Brebm3 resulted in an amino acid change (S to F) and the early-bolting phenotype of emb3. The mutation occurred in the SET domain (Suppressor of protein-effect variegation 3-9, Enhancer-of-zeste, Trithorax), which catalyzes site- and state-specific lysine methylation in histones. Tissue-specific expression analysis showed that Brebm3 was highly expressed in the flower and bud. Promoter activity assay confirmed that Brebm3 promoter was active in inflorescences. Subcellular localization analysis revealed that Brebm3 localized in the nucleus. Transcriptomic studies supported that Brebm3 mutation might repress H3K27me3 deposition and activate expression of the AGAMOUS (AG) and AGAMOUS-like (AGL) loci, resulting in early flowering.Conclusions: Our study revealed that an EMS-induced early-bolting mutant ebm3 in Chinese cabbage was caused by a nonsynonymous mutation in BraA04g017190.3C, encoding the histone methyltransferase CLF.These results improve our knowledge of the genetic and genomic resources of bolting and flowering, and may be beneficial to the genetic improvement of Chinese cabbage.

Whole Genome Sequencing Reveals the Effects of Recent Artificial Selection on Litter Size of Bamei Mutton Sheep

Bamei mutton sheep is a Chinese domestic sheep breed developed by crossing German Mutton Merino sheep and indigenous Mongolian sheep for meat production. Here, we focused on detecting candidate genes associated with the increasing of the litter size in this breeds under recent artificial selection to improve the efficiency of mutton production. We selected five high- and five low-fecundity Bamei mutton sheep for whole-genome resequencing to identify candidate genes for sheep prolificacy. We used the FST and XP-EHH statistical approach to detect the selective sweeps between these two groups. Combining the two selective sweep methods, the reproduction-related genes JUN, ITPR3, PLCB2, HERC5, and KDM4B were detected. JUN, ITPR3, and PLCB2 play vital roles in GnRH (gonadotropin-releasing hormone), oxytocin, and estrogen signaling pathway. Moreover, KDM4B, which had the highest FST value, exhibits demethylase activity. It can affect reproduction by binding the promoters of estrogen-regulated genes, such as FOXA1 (forkhead box A1) and ESR1 (estrogen receptor 1). Notably, one nonsynonymous mutation (p.S936A) specific to the high-prolificacy group was identified at the TUDOR domain of KDM4B. These observations provide a new opportunity to research the genetic variation influencing fecundity traits within a population evolving under artificial selection. The identified genomic regions that are responsible for litter size can in turn be used for further selection.