The evolutionary gain of spliceosomal introns: sequence and phase preferences

Theories regarding the evolution of spliceosomal introns differ in the extent to which the distribution of introns reflects either a formative role in the evolution of protein-coding genes or the adventitious gain of genetic elements. Here, systematic methods are used to assess the causes of the present-day distribution of introns in 10 families of eukaryotic protein-coding genes comprising 1,868 introns in 488 distinct alignment positions. The history of intron evolution inferred using a probabilistic model that allows ancestral inheritance of introns, gain of introns, and loss of introns reveals that the vast majority of introns in these eukaryotic gene families were not inherited from the most recent common ancestral genes, but were gained subsequently. Furthermore, among inferred events of intron gain that meet strict criteria of reliability, the distribution of sites of gain with respect to reading-frame phase shows a 5:3:2 ratio of phases 0, 1 and 2, respectively, and exhibits a nucleotide preference for MAG GT (positions -3 to +2 relative to the site of gain). The nucleotide preferences of intron gain may prove to be the ultimate cause for the phase bias. The phase bias of intron gain is sufficient to account quantitatively for the well-known 5:3:2 bias in phase frequencies among extant introns, a conclusion that holds even when taxonomic heterogeneity in phase patterns is considered. Thus, intron gain accounts for the vast majority of extant introns and for the bias toward phase 0 introns that previously was interpreted as evidence for ancient formative introns.     

Whole genome sequencing of a snailfish from the Yap Trench (~7,000 m) clarifies the molecular mechanisms underlying adaptation to the deep sea

Hadal environments (depths below 6,000 m) are characterized by extremely high hydrostatic pressures, low temperatures, a scarce food supply, and little light. The evolutionary adaptations that allow vertebrates to survive in this extreme environment are poorly understood. Here, we constructed a high-quality reference genome for Yap hadal snailfish (YHS), which was captured at a depth of ~7,000 m in the Yap Trench. The final YHS genome assembly was 731.75 Mb, with a contig N50 of 0.75 Mb and a scaffold N50 of 1.26 Mb. We predicted 24,329 protein-coding genes in the YHS genome, and 24,265 of these genes were successfully functionally annotated. Phylogenetic analyses suggested that YHS diverged from a Mariana Trench snailfish approximately 0.92 million years ago. Many genes associated with DNA repair show evidence of positive selection and have expanded copy numbers in the YHS genome, possibly helping to maintain the integrity of DNA under increased hydrostatic pressure. The levels of trimethylamine N-oxide (TMAO), a potent protein stabilizer, are much higher in the muscles of YHS than in those of shallow-water fish. This difference is perhaps due to the five copies of the TMAO-generating enzyme flavin-containing monooxygenase-3 gene (fmo3) in the YHS genome and the abundance of trimethylamine (TMA)-generating bacteria in the YHS gut. Thus, the high TMAO content might help YHS adapt to high hydrostatic pressure by improving protein stability. Additionally, the evolutionary features of the YHS genes encoding sensory-related proteins are consistent with the scarce food supply and darkness in the hadal environments. These results clarify the molecular mechanisms underlying the adaptation of hadal organisms to the deep-sea environment and provide valuable genomic resources for in-depth investigations of hadal biology.     

Ultrasound Based Assessment of Coronary Artery Flow and Coronary Flow Reserve Using the Pressure Overload Model in Mice

Transthoracic Doppler echocardiography (TTDE) is a clinically useful, noninvasive tool for studying coronary artery flow velocity and coronary flow reserve (CFR) in humans. Reduced CFR is accompanied by marked intramyocardial and pericoronary fibrosis and is used as an indication of the severity of dysfunction. This study explores, step-by-step, the real-time changes measured in the coronary flow velocity, CFR and systolic to diastolic peak velocity (S/D) ratio in the setting of an aortic banding model in mice. By using a Doppler transthoracic imaging technique that yields reproducible and reliable data, the method assesses changes in flow in the septal coronary artery (SCA), for a period of over two weeks in mice, that previously either underwent aortic banding or thoracotomy. During imaging, hyperemia in all mice was induced by isoflurane, an anesthetic that increased coronary flow velocity when compared with resting flow. All images were acquired by a single imager. Two ratios, (1) CFR, the ratio between hyperemic and baseline flow velocities, and (2) systolic (S) to diastolic (D) flow were determined, using a proprietary software and by two independent observers. Importantly, the observed changes in coronary flow preceded LV dysfunction as evidenced by normal LV mass and fractional shortening (FS). The method was benchmarked against the current gold standard of coronary assessment, histopathology. The latter technique showed clear pathologic changes in the coronary artery in the form of peri-coronary fibrosis that correlated to the flow changes as assessed by echocardiography. The study underscores the value of using a non-invasive technique to monitor coronary circulation in mouse hearts. The method minimizes redundant use of research animals and demonstrates that advanced ultrasound-based indices, such as CFR and S/D ratios, can serve as viable diagnostic tools in a variety of investigational protocols including drug studies and the study of genetically modified strains.     

Anaplastic Lymphoma Kinase Rearrangement in Digestive Tract Cancer: Implication for Targeted Therapy in Chinese Population

Background

Anaplastic lymphoma kinase (ALK) rearrangements define a subgroup of lung cancer which is eligible to targeted kinase inhibition. The aim of this study is to observe the incidence rate of ALK fusion in a large cohort of Chinese digestive tract cancer patients.

Patients and Methods

Tissue microarray (TMA) was constructed from 808 digestive tract cancer cases, including 169 esophageal squamous cell carcinoma, 182 gastric cancer and 457 colorectal cancer (CRC) cases. We tested all cases for ALK expression via a fully automated immunohistochemistry (IHC) assay. The IHC-positive cases were subjected to fluorescence in situ hybridization (FISH), real-time polymerase chain reaction (qRT-PCR), target gene enrichment and sequencing for confirmation of ALK gene rearrangement and discovery of novel fusion partner.

Results

Among the tested cases, 2 (0.44%) CRC cases showed positive both by IHC and FISH. By qRT-PCR, EML4–ALK fusion was found in one IHC-positive CRC case. In another IHC-positive CRC case, target gene enrichment and sequencing revealed ALK was fused to a novel partner, spectrin beta non-erythrocytic 1 (SPTBN1). One gastric cancer case showed partially positive IHC result, but no fusion was found by FISH and gene sequencing.

Conclusions

The incidence rate of ALK gene fusion in Chinese CRC patients was 0.44%,but not detectable in gastric and esophageal cancers. The novel SPTBN1 -ALK fusion, together with other ALK fusion genes, may become a potential target for anti-ALK therapy.     

Triterpenoids and Flavonoids from the Leaves of Astragalus membranaceus and Their Inhibitory Effects on Nitric Oxide Production

Four new cycloartane triterpenes, named huangqiyegenins V and VI and huangqiyenins K and L ( 1 – 4 , resp.), together with nine known triterpenoids, 5 – 13 , and eight flavonoids, 14 – 21 , were isolated from a 70%‐EtOH extract of Astragalus membranaceus leaves. The structures of the new compounds were elucidated by detailed spectroscopic analyses, and the compounds were identified as (9β,11α,16β,20R,24S)‐11,16,25‐trihydroxy‐20,24‐epoxy‐9,19‐cyclolanostane‐3,6‐dione ( 1 ), (9β,16β,24S)‐16,24,25‐trihydroxy‐9,19‐cyclolanostane‐3,6‐dione ( 2 ), (3β,6α,9β,16β,20R,24R)‐16,25‐dihydroxy‐3‐(β‐D ‐xylopyranosyloxy)‐20,24‐epoxy‐9,19‐cyclolanostan‐6‐yl acetate ( 3 ), and (3β,6α,9β,16β,24E)‐26‐(β‐D ‐glucopyranosyloxy)‐16‐hydroxy‐3‐(β‐D ‐xylopyranosyloxy)‐9,19‐cyclolanost‐24‐en‐6‐yl acetate ( 4 ). All isolated compounds were evaluated for their inhibitory activities against LPS‐induced NO production in RAW264.7 macrophage cells. Compounds 1 – 3, 14, 15 , and 18 exhibited strong inhibition on LPS‐induced NO release by macrophages with IC₅₀ values of 14.4–27.1 μM .     

An innovative protein expression system using RNA polymerase I for large-scale screening of high-nucleic-acid content Saccharomyces cerevisiae strains

Saccharomyces cerevisiae is the preferred source of RNA derivatives, which are widely used as supplements for foods and pharmaceuticals. As the most abundant RNAs, the ribosomal RNAs (rRNAs) transcribed by RNA polymerase I (Pol I) have no 5′ caps, thus cannot be translated to proteins. To screen high-nucleic-acid content yeasts more efficiently, a cap-independent protein expression system mediated by Pol I has been designed and established to monitor the regulatory changes of rRNA synthesis by observing the variation in the reporter genes expression. The elements including Pol I-recognized rDNA promoter, the internal ribosome entry site from cricket paralytic virus which can recruit ribosomes internally, reporter genes (URA3 and yEGFP3), oligo-dT and an rDNA terminator were ligated to a yeast episomal plasmid. This system based on the URA3 gene worked well by observing the growth phenotype and did not require the disruption of cap-dependent initiation factors. The fluorescence intensity of strains expressing the yEGFP3 gene increased and drifted after mutagenesis. Combined with flow cytometry, cells with higher GFP level were sorted out. A strain showed 58% improvement in RNA content and exhibited no sequence alteration in the whole expression cassette introduced. This study provides a novel strategy for breeding high-nucleic-acid content yeasts.     

Identification of curcumin derivatives as human glyoxalase I inhibitors: A combination of biological evaluation, molecular docking, 3D-QSAR and molecular dynamics simulation studies

Several recent developments suggest that the human glyoxalase I (GLO I) is a potential target for anti-tumor drug development. In present study, a series of curcumin derivatives with high inhibitory activity against human GLO I were discovered. Inhibition constant (K(i)) values of compounds 8, 9, 10, 11 and 13 to GLO I are 4.600μM, 2.600μM, 3.200μM, 3.600μM and 3.600μM, respectively. To elucidate the structural features of potent inhibitors, docking-based three-dimensional structure-activity relationship (3D-QSAR) analyses were performed. Satisfactory agreement between experiment and theory suggests that comparative molecular similarity index analysis (CoMSIA) modeling exhibit much better correlation and predictive power. The cross-validated q(2) value is 0.638 while no-validation r(2) value is 0.930. Integrated with docking-based 3D-QSAR CoMSIA modeling, molecular surface property (electrostatic and steric) mapping and molecular dynamics simulation, a set of receptor-ligand binding models and bio-affinity predictive models for rational design of more potent inhibitors of GLO I are established.     

Study of strategies for selecting quantitative trait locus mapping procedures by computer simulation

Along with the development and integration of molecular genetics and quantitative genetics, many quantitative trait locus (QTL) mapping studies have been conducted using different mapping populations in various crop species. Existing QTLs can be used for marker-assisted breeding and map-based cloning, whereas the false-positive QTLs are no use. The purpose of this study is to evaluate the suitability of different mapping procedures for data from different genetic models. In this study, four types of recombinant inbred lines (RILs) with different genetic models, viz. additive QTLs (Model I), additive and epistatic QTLs (Model II), additive QTLs and QTL × environment interaction (Model III), additive, epistatic QTLs and QTL × environment interaction (Model IV), were simulated by computer. Six types of QTL mapping procedures, viz. CIM, MIMF, MIMR, ICIM, MQM and NWIM, on four kinds of QTL mapping software, viz. WinQTL Cartographer Version 2.5, IciMapping Version 2.0, MapQTL Version 5.0 and QTLnetwork Version 2.0, were used for screening QTLs of the simulated RILs. The results showed that different mapping procedures have different suitability for different genetic models. CIM and MQM can only screen Model I data. MIMR, MIMF and ICIM can only screen Model I and Model II data. NWIM can screen all four models’ data. It can be concluded that different genetic models’ data have different most suitable mapping procedures. In practical experiments where the genetic model of the data is unknown, a multiple model mapping strategy should be used, that is a full model scanning with complex model procedure followed by verification with other procedures corresponding to the scanning results.