All-in-one sequencing: an improved library preparation method for cost-effective and high-throughput next-generation sequencing

Conservation genomics analysis reveals recent population decline and possible causes in bumblebee Bombus opulentus

Bumblebees are a genus of pollinators (Bombus) that play important roles in natural ecosystem and agricultural production. Several bumblebee species have been recorded as under population decline, and the proportion of species experiencing population decline within subgenus Thoracobombus is higher than average. Bombus opulentus is 1 species in Thoracobombus, but little is known about its recent population dynamics. Here, we employed conservation genomics methods to investigate the population dynamics of B. opulentus during the recent past and identify the likely environmental factors that may cause population decline. Firstly, we placed the scaffold-level of B. opulentus reference genome sequence onto chromosome-level using Hi-C technique. Then, based on this reference genome and whole-genome resequencing data for 51 B. opulentus samples, we reconstructed the population structure and effective population size (Ne) trajectories of B. opulentus and identified genes that were under positive selection. Our results revealed that the collected B. opulentus samples could be divided into 2 populations, and 1 of them experienced a recent population decline; the declining population also exhibited lower genetic diversity and higher inbreeding levels. Genes related to high-temperature tolerance, immune response, and detoxication showed signals of positive selection in the declining population, suggesting that climate warming and pathogen/pesticide exposures may contribute to the decline of this B. opulentus population. Taken together, our study provided insights into the demography of B. opulentus populations and highlighted that populations of the same bumblebee species could have contrasting Ne trajectories and population decline could be caused by a combination of various stressors.

Streamlined whole-genome genotyping through NGS-enhanced thermal asymmetric interlaced (TAIL)-PCR

Whole-genome genotyping (WGG) stands as a pivotal element in genomic-assisted plant breeding. Nevertheless, sequencing-based approaches for WGG continue to be costly, primarily owing to the high expenses associated with library preparation and the laborious protocol. During prior development of foreground and background integrated genotyping by sequencing (FBI-seq), we discovered that any sequence-specific primer (SP) inherently possesses the capability to amplify a massive array of stable and reproducible non-specific PCR products across the genome. Here, we further improved FBI-seq by replacing the adapter ligated by Tn5 transposase with an arbitrary degenerate (AD) primer. The protocol for the enhanced FBI-seq unexpectedly mirrors a simplified thermal asymmetric interlaced (TAIL)-PCR, a technique that is widely used for isolation of flanking sequences. However, the improved TAIL-PCR maximizes the primer-template mismatched annealing capabilities of both SP and AD primers. In addition, leveraging of next-generation sequencing enhances the ability of this technique to assay tens of thousands of genome-wide loci for any species. This cost-effective, user-friendly, and powerful WGG tool, which we have named TAIL-PCR by sequencing (TAIL-peq), holds great potential for widespread application in breeding programs, thereby facilitating genome-assisted crop improvement.

Population Structure, Demographic History, and Adaptation of Giant Honeybees in China Revealed by Population Genomic Data

There have been many population-based genomic studies on human-managed honeybees (Apis mellifera and Apis cerana), but there has been a notable lack of analysis with regard to wild honeybees, particularly in relation to their evolutionary history. Nevertheless, giant honeybees have been found to occupy distinct habitats and display remarkable characteristics, which are attracting an increased amount of attention. In this study, we de novo sequenced and then assembled the draft genome sequence of the Himalayan giant honeybee, Apis laboriosa. Phylogenetic analysis based on genomic information indicated that A. laboriosa and its tropical sister species Apis dorsata diverged ∼2.61 Ma, which supports the speciation hypothesis that links A. laboriosa to geological changes throughout history. Furthermore, we re-sequenced A. laboriosa and A. dorsata samples from five and six regions, respectively, across their population ranges in China. These analyses highlighted major genetic differences for Tibetan A. laboriosa as well as the Hainan Island A. dorsata. The demographic history of most giant honeybee populations has mirrored glacial cycles. More importantly, contrary to what has occurred among human-managed honeybees, the demographic history of these two wild honeybee species indicates a rapid decline in effective population size in the recent past, reflecting their differences in evolutionary histories. Several genes were found to be subject to selection, which may help giant honeybees to adapt to specific local conditions. In summary, our study sheds light on the evolutionary and adaptational characteristics of two wild giant honeybee species, which was useful for giant honeybee conservation.

A prolific and robust whole-genome genotyping method using PCR amplification via primer-template mismatched annealing

Whole-genome genotyping methods are important for breeding. However, it has been challenging to develop a robust method for simultaneous foreground and background genotyping that can easily be adapted to different genes and species. In our study, we accidently discovered that in adapter ligation-mediated PCR, the amplification by primer-template mismatched annealing (PTMA) along the genome could generate thousands of stable PCR products. Based on this observation, we consequently developed a novel method for simultaneous foreground and background integrated genotyping by sequencing (FBI-seq) using one specific primer, in which foreground genotyping is performed by primer-template perfect annealing (PTPA), while background genotyping employs PTMA. Unlike DNA arrays, multiple PCR, or genome target enrichments, FBI-seq requires little preliminary work for primer design and synthesis, and it is easily adaptable to different foreground genes and species. FBI-seq therefore provides a prolific, robust, and accurate method for simultaneous foreground and background genotyping to facilitate breeding in the post-genomics era.

Transcriptome Analysis to Identify Genes Related to Flowering Reversion in Tomato

Flowering reversion is a common phenomenon in plant development in which differentiated floral organs switch from reproductive growth to vegetative growth and ultimately form abnormal floral organs or vegetative organs. This greatly reduces tomato yield and quality. Research on this phenomenon has recently increased, but there is a lack of research at the molecular and gene expression levels. Here, transcriptomic analyses of the inflorescence meristem were performed in two kinds of materials at different developmental stages, and a total of 3223 differentially expressed genes (DEGs) were screened according to the different developmental stages and trajectories of the two materials. The analysis of database annotations showed that these DEGs were closely related to starch and sucrose metabolism, DNA replication and modification, plant hormone synthesis and signal transduction. It was further speculated that tomato flowering reversion may be related to various biological processes, such as cell signal transduction, energy metabolism and protein post-transcriptional regulation. Combined with the results of previous studies, our work showed that the gene expression levels of CLE9, FA, PUCHI, UF, CLV3, LOB30, SFT, S-WOX9 and SVP were significantly different in the two materials. Endogenous hormone analysis and exogenous hormone treatment revealed a variety of plant hormones involved in flowering reversion in tomato. Thus, tomato flowering reversion was studied comprehensively by transcriptome analysis for the first time, providing new insights for the study of flower development regulation in tomato and other plants.

ERH Interacts With EIF2α and Regulates the EIF2α/ATF4/CHOP Pathway in Bladder Cancer Cells

Background

There is a lack of research on the molecular interaction of the enhancers of rudimentary homolog (ERH) in bladder cancer (BC) cells. This study aimed to determine the interacting proteins of ERH in human T24 cells.

Methods

First, the ERH gene was overexpressed in human T24 cells. Coimmunoprecipitation (co-IP) and shotgun mass spectrometry (MS) analyses were performed to obtain a list of proteins that interact with ERH. Subsequently, bioinformatic analyses with Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) and protein–protein interaction (PPI) studies were performed to analyze the ERH-interactive protein list (ERH-IPL). Then, we selected one of the interacting proteins, EIF2α for verification. An immunofluorescence colocalization assay was performed to validate the co-expression of the selected protein, and the binding sites of the two proteins were predicted by ZDOCK technology. Finally, PCR analysis on the downstream molecules of the interacting protein was performed for verification.

Results

ERH protein was successfully overexpressed in human T24 cells. We obtained a list of 205 proteins that might directly or indirectly interact with the ERH protein by mass spectrometric analysis. The bioinformatic analysis showed that ERH-interacting proteins were related to “ribonucleoprotein complex”, “ATPase activity”, “nuclear speck”, and “translation factor activity, RNA binding”. We further identified one of the key genes, EIF2S1, and confirmed that the corresponding protein EIF2α is co-expressed and may bind with ERH in human T24 cells. The mRNA levels of molecules ATF4 and CHOP were found to be upregulated by ERH.

Conclusion

ERH protein affects “ribonucleoprotein complex”, “ATPase activity”, “nuclear speck”, and “translation factor activity, RNA binding”. The ERH protein can interact with EIF2α and regulate the EIF2α-ATF4/CHOP signaling pathway in human T24 cells.

STAT3 Is the Master Regulator for the Forming of 3D Spheroids of 3T3-L1 Preadipocytes

To elucidate the currently unknown mechanisms responsible for the diverse biological aspects between two-dimensional (2D) and three-dimensional (3D) cultured 3T3-L1 preadipocytes, RNA-sequencing analyses were performed. During a 7-day culture period, 2D- and 3D-cultured 3T3-L1 cells were subjected to lipid staining by BODIPY, qPCR for adipogenesis related genes, including peroxisome proliferator-activated receptor γ (Pparγ), CCAAT/enhancer-binding protein alpha (Cebpa), Ap2 (fatty acid-binding protein 4; Fabp4), leptin, and AdipoQ (adiponectin), and RNA-sequencing analysis. Differentially expressed genes (DEGs) were detected by next-generation RNA sequencing (RNA-seq) and validated by a quantitative reverse transcription–polymerase chain reaction (qRT–PCR). Bioinformatic analyses were performed on DEGs using a Gene Ontology (GO) enrichment analysis and an Ingenuity Pathway Analysis (IPA). Significant spontaneous adipogenesis was observed in 3D 3T3-L1 spheroids, but not in 2D-cultured cells. The mRNA expression of Pparγ, Cebpa, and Ap2 among the five genes tested were significantly higher in 3D spheroids than in 2D-cultured cells, thus providing support for this conclusion. RNA analysis demonstrated that a total of 826 upregulated and 725 downregulated genes were identified as DEGs. GO enrichment analysis and IPA found 50 possible upstream regulators, and among these, 6 regulators—transforming growth factor β1 (TGFβ1), signal transducer and activator of transcription 3 (STAT3), interleukin 6 (IL6), angiotensinogen (AGT), FOS, and MYC—were, in fact, significantly upregulated. Further analyses of these regulators by causal networks of the top 14 predicted diseases and functions networks (IPA network score indicated more than 30), suggesting that STAT3 was the most critical upstream regulator. The findings presented herein suggest that STAT3 has a critical role in regulating the unique biological properties of 3D spheroids that are produced from 3T3-L1 preadipocytes.

Identification of Genes Related to Cold Tolerance and Novel Genetic Markers for Molecular Breeding in Taiwan Tilapia (Oreochromis spp.) via Transcriptome Analysis

Taiwan tilapia is one of the primary species used in aquaculture practices in Taiwan. However, as a tropical fish, it is sensitive to cold temperatures that can lead to high mortality rates during winter months. Genetic and broodstock management strategies using marker-assisted selection and breeding are the best tools currently available to improve seed varieties for tilapia species. The purpose of this study was to develop molecular markers for cold stress-related genes using digital gene expression analysis of next-generation transcriptome sequencing in Taiwan tilapia (Oreochromis spp.). We constructed and sequenced cDNA libraries from the brain, gill, liver, and muscle tissues of cold-tolerance (CT) and cold-sensitivity (CS) strains. Approximately 35,214,833,100 nucleotides of raw sequencing reads were generated, and these were assembled into 128,147 unigenes possessing a total length of 185,382,926 bp and an average length of 1446 bp. A total of 25,844 unigenes were annotated using five protein databases and Venny analysis, and 38,377 simple sequence repeats (SSRs) and 65,527 single nucleotide polymorphisms (SNPs) were identified. Furthermore, from the 38-cold tolerance-related genes that were identified using differential gene expression analysis in the four tissues, 13 microsatellites and 37 single nucleotide polymorphism markers were identified. The results of the genotype analysis revealed that the selected markers could be used for population genetics. In addition to the diversity assessment, one of the SNP markers was determined to be significantly related to cold-tolerance traits and could be used as a molecular marker to assist in the selection and verification of cold-tolerant populations. The specific genetic markers explored in this study can be used for the identification of genetic polymorphisms and cold tolerance traits in Taiwan tilapia, and they can also be used to further explore the physiological and biochemical molecular regulation pathways of fish that are involved in their tolerance to environmental temperature stress.

Whole-exome sequencing identification of a recurrent CRYBB2 variant in a four-generation Chinese family with congenital nuclear cataracts

Congenital cataracts is the most common cause of visual impairment and blindness in children. Although there have been extensive studies into the pathogenesis of congenital cataracts, the pathogenic mechanism underlying the recurrent variant CRYBB2:c.62T>A(p.I21N) has not been previously reported. Thus, the present study aimed to use whole-exome sequencing (WES) to identify potential genetic variants and investigate how they may have induced the occurrence of cataracts in a four-generation Chinese family with congenital nuclear cataracts. The medical history of this family was recorded and WES was conducted for one proband. Sanger sequencing was used to verify the presence of the putative variant in all participants. PolyPhen-2, SIFT and ProtScale were used to analyze the effect of the identified variants on protein function and hydrophobicity, and Pymol was used to show the structure of the wild-type (Wt) and mutant β-crystallin B2 (CRYBB2) protein. Full-length Wt-CRYBB2 or mutant-CRYBB2 (I21N-CRYBB2) were fused to green fluorescent protein (GFP), and the recombinant plasmids were transfected into HeLa cells. Reverse transcription-quantitative PCR and western blotting were used to detect the expression levels of CRYBB2 mRNA and protein. Immunofluorescence and flow cytometry analyses were used to detect protein localization and apoptosis, respectively. A recurrent variant CRYBB2:c.62T>A(p.I21N) was identified in a four-generation Chinese family with congenital nuclear cataracts. Multiple-sequence alignment of CRYBB2 demonstrated that codon 21 was highly conserved. Pymol revealed that the structure of the I21N-CRYBB2 protein was distinct from that of Wt-CRYBB2. PolyPhen-2 predicted that it had a variant provean score 1.0, suggesting it was 'probably damaging', and SIFT predicted it had a variant provean score of -5.113, indicating it was 'deleterious'. ProtScale indicated that the hydrophobicity of the mutation site was significantly reduced. The protein expression levels of the I21N-CRYBB2 were decreased compared with the Wt-CRYBB2. Immunofluorescence analysis revealed that the variant I21N-CRYBB2 protein tended to accumulate around the nucleus, and flow cytometry analysis indicated that it increased cell apoptosis. Furthermore, I21N-CRYBB2 induced the activation of the unfolded protein response (UPR). In conclusion, a pathogenic variant of CRYBB2:c.62T>A(p.I21N) was identified via WES in a four-generation Chinese family with congenital nuclear cataracts. Through biological analysis, it was found that the variant induced abnormal protein aggregation, activated the UPR and triggered excessive cell apoptosis, which may lead to the occurrence of congenital nuclear cataracts in this family.

High-quality genome assembly of Huazhan and Tianfeng, the parents of an elite rice hybrid Tian-you-hua-zhan

High-quality rice reference genomes have accelerated the comprehensive identification of genome-wide variations and research on functional genomics and breeding. Tian-you-hua-zhan has been a leading hybrid in China over the past decade. Here, de novo genome assembly strategy optimization for the rice indica lines Huazhan (HZ) and Tianfeng (TF), including sequencing platforms, assembly pipelines and sequence depth, was carried out. The PacBio and Nanopore platforms for long-read sequencing were utilized, with the Canu, wtdbg2, SMARTdenovo, Flye, Canu-wtdbg2, Canu-SMARTdenovo and Canu-Flye assemblers. The combination of PacBio and Canu was optimal, considering the contig N50 length, contig number, assembled genome size and polishing process. The assembled contigs were scaffolded with Hi-C data, resulting in two "golden quality" rice reference genomes, and evaluated using the scaffold N50, BUSCO, and LTR assembly index. Furthermore, 42,625 and 41,815 non-transposable element genes were annotated for HZ and TF, respectively. Based on our assembly of HZ and TF, as well as Zhenshan97, Minghui63, Shuhui498 and 9311, comprehensive variations were identified using Nipponbare as a reference. The de novo assembly strategy for rice we optimized and the "golden quality" rice genomes we produced for HZ and TF will benefit rice genomics and breeding research, especially with respect to uncovering the genomic basis of the elite traits of HZ and TF.