Ensembl Genomes 2013: scaling up access to genome-wide data

Genetic Variants in the ABCB1 and ABCG2 Gene Drug Transporters Involved in Gefitinib-Associated Adverse Reaction: A Systematic Review and Meta-Analysis

This systematic review and meta-analysis aimed to verify the association between the genetic variants of adenosine triphosphate (ATP)-binding cassette subfamily B member 1 (ABCB1) and ATP-binding cassette subfamily G member 2 (ABCG2) genes and the presence and severity of gefitinib-associated adverse reactions. We systematically searched PubMed, Virtual Health Library/Bireme, Scopus, Embase, and Web of Science databases for relevant studies published up to February 2024. In total, five studies were included in the review. Additionally, eight genetic variants related to ABCB1 (rs1045642, rs1128503, rs2032582, and rs1025836) and ABCG2 (rs2231142, rs2231137, rs2622604, and 15622C>T) genes were analyzed. Meta-analysis showed a significant association between the ABCB1 gene rs1045642 TT genotype and presence of diarrhea (OR = 5.41, 95% CI: 1.38-21.14, I2 = 0%), the ABCB1 gene rs1128503 TT genotype and CT + TT group and the presence of skin rash (OR = 4.37, 95% CI: 1.51-12.61, I2 = 0% and OR = 6.99, 95%CI: 1.61-30.30, I2= 0%, respectively), and the ABCG2 gene rs2231142 CC genotype and presence of diarrhea (OR = 3.87, 95% CI: 1.53-9.84, I2 = 39%). No ABCB1 or ABCG2 genes were positively associated with the severity of adverse reactions associated with gefitinib. In conclusion, this study showed that ABCB1 and ABCG2 variants are likely to exhibit clinical implications in predicting the presence of adverse reactions to gefitinib.

Genetic Association of Toll-Like Receptor 4 (TLR4) Gene Polymorphism (rs4986790) With Oral Squamous Cell Carcinoma (OSCC): A Pilot Case-Control Study

Introduction Oral squamous cell carcinoma (OSCC) is a highly prevalent and most common form of oral malignancy in the Indian population. Toll-like receptors belong to an important family of receptors that are involved in the process of pathogen recognition and mounting immune response. The expression of this receptor is dysregulated on the tumor cells as reported across several cancer types. The genetic variants in this gene could have a profound impact on the expression of the Toll-like receptor 4 (TLR4) gene.  Objective This study aimed to understand the association of TLR4 gene polymorphism (rs4986790) with OSCC. The objective of this study was to compare the allele and genotype frequencies between the two groups, viz., OSCC and normal healthy subjects, recruited in the study. Materials and methods The blood samples were collected from normal healthy subjects (N = 25) and OSCC patients (N = 25). Genomic DNA was isolated from all samples, and genotyping was performed for the TLR4 gene polymorphism (rs4986790) employing the polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) approach. The frequency distribution of genotypes and alleles across the study groups was determined by the Chi-square test.  Results The allele frequency for TLR4 gene polymorphism (rs4986790) in the case group was found to be 60% (A allele) and 40% (G allele), respectively. The study population in both groups were found to agree with the Hardy-Weinberg equilibrium (HWE). The genotype frequency did not differ significantly among the two study groups which was evident from the p-value = 0.8285.  Conclusion The present study did not report any significant association of the TLR4 polymorphic marker rs4986790 with OSCC. Further investigations into the association of other polymorphic markers in the TLR4 gene, among the larger population of OSCC patients, could provide evidence of their association with OSCC.

MYB24 orchestrates terpene and flavonol metabolism as light responses to anthocyanin depletion in variegated grape berries

Variegation is a rare type of mosaicism not fully studied in plants, especially fruits. We examined red and white sections of grape (Vitis vinifera cv. 'Béquignol') variegated berries and found that accumulation of products from branches of the phenylpropanoid and isoprenoid pathways showed an opposite tendency. Light-responsive flavonol and monoterpene levels increased in anthocyanin-depleted areas in correlation with increasing MYB24 expression. Cistrome analysis suggested that MYB24 binds to the promoters of 22 terpene synthase (TPS) genes, as well as 32 photosynthesis/light-related genes, including carotenoid pathway members, the flavonol regulator HY5 HOMOLOGUE (HYH), and other radiation response genes. Indeed, TPS35, TPS09, the carotenoid isomerase gene CRTISO2, and HYH were activated in the presence of MYB24 and MYC2. We suggest that MYB24 modulates ultraviolet and high-intensity visible light stress responses that include terpene and flavonol synthesis and potentially affects carotenoids. The MYB24 regulatory network is developmentally triggered after the onset of berry ripening, while the absence of anthocyanin sunscreens accelerates its activation, likely in a dose-dependent manner due to increased radiation exposure. Anthocyanins and flavonols in variegated berry skins act as effective sunscreens but for different wavelength ranges. The expression patterns of stress marker genes in red and white sections of 'Béquignol' berries strongly suggest that MYB24 promotes light stress amelioration but only partly succeeds during late ripening.

Genome-wide analysis of glutamate receptor gene family in allopolyploid Brassica napus and its diploid progenitors

Ionotropic glutamate receptors are ligand-gated nonselective cation channels that mediate neurotransmission in the central nervous system of animals. Plants possess homologous proteins called glutamate receptor-like channels (GLRs) which are involved in vital physiological processes including seed germination, long-distance signaling, chemotaxis, Ca2+ signaling etc. Till now, a comprehensive genome-wide analysis of the GLR gene family members in different economically important species of Brassica is missing. Considering the origin of allotetraploid Brassica napus from the hybridization between the diploid Brassica oleracea and Brassica rapa, we have identified 11, 27 and 65 GLR genes in B. oleracea, B. rapa and B. napus, respectively showing an expansion of this gene family in B. napus. Chromosomal locations revealed several tandemly duplicated GLR genes in all the three species. Moreover, the gene family expanded in B. napus after allopolyploidization. The phylogenetic analysis showed that the 103 GLRs are classified into three main groups. The exon-intron structures of these genes are not very conserved and showed wide variation in intron numbers. However, protein sequences are much conserved as shown by the presence of ten short amino acid sequence motifs. Predicted cis-acting elements in 1 kb promoters of GLR genes are mainly involved in light, stress and hormone responses. RNA-seq analysis showed that in B. oleracea and B. rapa, some GLRs are more tissue specific than others. In B. napus, some GLRs are downregulated under cold stress, while others are upregulated. In summary, this bioinformatic study of the GLR gene family of the three Brassica species provides evidence for the expansion of this gene family in B. napus and also provided useful information for in-depth studies of their biological functions in Brassica.

Genome-Wide Analysis of the Membrane Attack Complex and Perforin Genes and Their Expression Pattern under Stress in the Solanaceae

The Membrane Attack Complex and Perforin (MACPF) proteins play a crucial role in plant development and adaptation to environmental stresses. Heretofore, few MACPF genes have been functionally identified, leaving gaps in our understanding of MACPF genes in other plants, particularly in the Solanaceae family, which includes economically and culturally significant species, such as tomato, potato, and pepper. In this study, we have identified 26 MACPF genes in three Solanaceae species and in the water lily, which serves as the base group for angiosperms. Phylogenetic analysis indicates that angiosperm MACPF genes could be categorized into three distinct groups, with another moss and spikemoss lineage-specific group, which is further supported by the examination of gene structures and domain or motif organizations. Through inter-genome collinearity analysis, it is determined that there are 12 orthologous SolMACPF gene pairs. The expansion of SolMACPF genes is primarily attributed to dispersed duplications, with purifying selection identified as the principal driving force in their evolutionary process, as indicated by the ω values. Furthermore, the analysis of expression patterns revealed that Solanaceae genes are preferentially expressed in reproductive tissues and regulated by various environmental stimuli, particularly induced by submergence. Taken together, these findings offer valuable insights into and a fresh perspective on the evolution and function of SolMACPF genes, thereby establishing a foundation for further investigations into their phenotypic and functional characteristics.

Evolution and Expression of the Meprin and TRAF Homology Domain-Containing Gene Family in Solanaceae

Meprin and TRAF homology (MATH)-domain-containing proteins are pivotal in modulating plant development and environmental stress responses. To date, members of the MATH gene family have been identified only in a few plant species, including Arabidopsis thaliana, Brassica rapa, maize, and rice, and the functions of this gene family in other economically important crops, especially the Solanaceae family, remain unclear. The present study identified and analyzed 58 MATH genes from three Solanaceae species, including tomato (Solanum lycopersicum), potato (Solanum tuberosum), and pepper (Capsicum annuum). Phylogenetic analysis and domain organization classified these MATH genes into four groups, consistent with those based on motif organization and gene structure. Synteny analysis found that segmental and tandem duplication might have contributed to MATH gene expansion in the tomato and the potato, respectively. Collinearity analysis revealed high conservation among Solanaceae MATH genes. Further cis-regulatory element prediction and gene expression analysis showed that Solanaceae MATH genes play essential roles during development and stress response. These findings provide a theoretical basis for other functional studies on Solanaceae MATH genes.

CaREM1.4 interacts with CaRIN4 to regulate Ralstonia solanacearum tolerance by triggering cell death in pepper

Remorins, plant-specific proteins, have a significant role in conferring on plants the ability to adapt to adverse environments. However, the precise function of remorins in resistance to biological stress remains largely unknown. Eighteen CaREM genes were identified in pepper genome sequences based on the C-terminal conserved domain that is specific to remorin proteins in this research. Phylogenetic relations, chromosomal localization, motif, gene structures, and promoter regions of these remorins were analyzed and a remorin gene, CaREM1.4, was cloned for further study. The transcription of CaREM1.4 in pepper was induced by infection with Ralstonia solanacearum. Knocking down CaREM1.4 in pepper using virus-induced gene silencing (VIGS) technologies reduced the resistance of pepper plants to R. solanacearum and downregulated the expression of immunity-associated genes. Conversely, transient overexpression of CaREM1.4 in pepper and Nicotiana benthamiana plants triggered hypersensitive response-mediated cell death and upregulated expression of defense-related genes. In addition, CaRIN4-12, which interacted with CaREM1.4 at the plasma membrane and cell nucleus, was knocked down with VIGS, decreasing the susceptibility of Capsicum annuum to R. solanacearum. Furthermore, CaREM1.4 reduced ROS production by interacting with CaRIN4-12 upon co-injection in pepper. Taken together, our findings suggest that CaREM1.4 may function as a positive regulator of the hypersensitive response, and it interacts with CaRIN4-12, which negatively regulates plant immune responses of pepper to R. solanacearum. Our study provides new evidence for comprehending the molecular regulatory network of plant cell death.

PlantCADB: A Comprehensive Plant Chromatin Accessibility Database

Chromatin accessibility landscapes are essential for detecting regulatory elements, illustrating the corresponding regulatory networks, and, ultimately, understanding the molecular basis underlying key biological processes. With the advancement of sequencing technologies, a large volume of chromatin accessibility data has been accumulated and integrated for humans and other mammals. These data have greatly advanced the study of disease pathogenesis, cancer survival prognosis, and tissue development. To advance the understanding of molecular mechanisms regulating plant key traits and biological processes, we developed a comprehensive plant chromatin accessibility database (PlantCADB) from 649 samples of 37 species. These samples are abiotic stress-related (such as heat, cold, drought, and salt; 159 samples), development-related (232 samples), and/or tissue-specific (376 samples). Overall, 18,339,426 accessible chromatin regions (ACRs) were compiled. These ACRs were annotated with genomic information, associated genes, transcription factor footprint, motif, and single-nucleotide polymorphisms (SNPs). Additionally, PlantCADB provides various tools to visualize ACRs and corresponding annotations. It thus forms an integrated, annotated, and analyzed plant-related chromatin accessibility resource, which can aid in better understanding genetic regulatory networks underlying development, important traits, stress adaptations, and evolution.PlantCADB is freely available at https://bioinfor.nefu.edu.cn/PlantCADB/.

MicroRNA172b-5p/trehalose-6-phosphate synthase module stimulates trehalose synthesis and microRNA172b-3p/AP2-like module accelerates flowering in barley upon drought stress

MicroRNAs (miRNAs) are major regulators of gene expression during plant development under normal and stress conditions. In this study, we analyzed the expression of 150 conserved miRNAs during drought stress applied to barley ready to flower. The dynamics of miRNAs expression was also observed after rewatering. Target messenger RNA (mRNAs) were experimentally identified for all but two analyzed miRNAs, and 41 of the targets were not reported before. Drought stress applied to barley induced accelerated flowering coordinated by a pair of two differently expressed miRNAs originating from a single precursor: hvu-miR172b-3p and hvu-miR172b-5p. Increased expression of miRNA172b-3p during drought leads to the downregulation of four APETALA2(AP2)-like genes by their mRNA cleavage. In parallel, the downregulation of the miRNA172b-5p level results in an increased level of a newly identified target, trehalose-6-phosphate synthase, a key enzyme in the trehalose biosynthesis pathway. Therefore, drought-treated plants have higher trehalose content, a known osmoprotectant, whose level is rapidly dropping after watering. In addition, trehalose-6-phosphate, an intermediate of the trehalose synthesis pathway, is known to induce flowering. The hvu-miRNA172b-5p/trehalose-6-phosphate synthase and hvu-miRNA172b-3p/AP2-like create a module leading to osmoprotection and accelerated flowering induction during drought.

Was the Last Bacterial Common Ancestor a Monoderm after All?

The very nature of the last bacterial common ancestor (LBCA), in particular the characteristics of its cell wall, is a critical issue to understand the evolution of life on earth. Although knowledge of the relationships between bacterial phyla has made progress with the advent of phylogenomics, many questions remain, including on the appearance or disappearance of the outer membrane of diderm bacteria (also called Gram-negative bacteria). The phylogenetic transition between monoderm (Gram-positive bacteria) and diderm bacteria, and the associated peptidoglycan expansion or reduction, requires clarification. Herein, using a phylogenomic tree of cultivated and characterized bacteria as an evolutionary framework and a literature review of their cell-wall characteristics, we used Bayesian ancestral state reconstruction to infer the cell-wall architecture of the LBCA. With the same phylogenomic tree, we further revisited the evolution of the division and cell-wall synthesis (dcw) gene cluster using homology- and model-based methods. Finally, extensive similarity searches were carried out to determine the phylogenetic distribution of the genes involved with the biosynthesis of the outer membrane in diderm bacteria. Quite unexpectedly, our analyses suggest that all cultivated and characterized bacteria might have evolved from a common ancestor with a monoderm cell-wall architecture. If true, this would indicate that the appearance of the outer membrane was not a unique event and that selective forces have led to the repeated adoption of such an architecture. Due to the lack of phenotypic information, our methodology cannot be applied to all extant bacteria. Consequently, our conclusion might change once enough information is made available to allow the use of an even more diverse organism selection.

Characterization of endogenous promoters of GapC1 and GS for recombinant protein expression in Phaeodactylum tricornutum

Although diatoms have been utilized as a cellular factory to produce biopharmaceuticals, recombinant proteins, and biofuels, only a few numbers of gene promoters are available. Therefore, the development of novel endogenous promoters is essential for the production of a range of bioactive substances. Here, we characterized the activities of endogenous promoters glyceraldehyde-3-phosphate dehydrogenase (GapC1) and glutamine synthetase (GS) of Phaeodactylum tricornutum using green fluorescent protein (GFP) under different culture conditions. Compared with the widely used fucoxanthin chlorophyll-binding protein A (fcpA) promoter, the GS promoter constitutively drove the expression of GFP throughout all growth phases of P. tricornutum, regardless of culture conditions. Additionally, the GFP level driven by the GapC1 promoter was the highest at the log phase, similar to the fcpA promoter, and increased light and nitrogen-starvation conditions reduced GFP levels by inhibiting promoter activity. These results suggested that the GS promoter could be utilized as a strong endogenous promoter for the genetic engineering of P. tricornutum.

Performance Evaluation of Warfarin Dose Prediction Algorithms and Effects of Clinical Factors on Warfarin Dose in Chinese Patients

BACKGROUND

The clinical utility of warfarin dose prediction algorithms remains controversial, our purpose is to evaluate the performance of warfarin dose prediction algorithms and the effects of clinical factors on warfarin dose in Chinese patients.

METHODS

Clinical data of 217 patients who received warfarin treatment were used to assess 6 warfarin dose prediction algorithms (OHNO, IWPC [International Warfarin Pharmacogenetics Consortium], HUANG, KIM, BRESS, and MIAO). The predicted dose (PD) was compared with the warfarin optimal dose (WOD, defined as the dose that maintains the international normalized ratio within the target range of 2.0-3.0). A multiple regression analysis with WOD as the dependent variable was performed to evaluate the effects of clinical factors on warfarin dose.

RESULTS

The mean absolute error analysis ranked the predictive accuracies of the algorithms as OHNO > IWPC > HUANG > KIM > BRESS > MIAO. Stratified analysis indicated that HUANG most accurately predicted that patients required lower WODs (≤3 mg/d), whereas OHNO was the most effective in predicting medium WODs (3-5 mg/d). KIM was effective in predicting high WODs (>5 mg/d). Multiple linear regression analysis showed that VKORC1 (rs9923231) and body mass index were significantly positively correlated with WOD, whereas concurrent atrial fibrillation status, CYP2C9*3 (rs1057910), and sex were significantly negatively correlated with WOD.

CONCLUSIONS

In Chinese patients, OHNO should be given priority during the prediction and selection of warfarin dose. When using OHNO to predict warfarin dose (≤3 mg/d or >5 mg/d), HUANG or KIM algorithms can provide precise predictions. At the same time, physicians should pay close attention to clinical factors, such as VKORC1 (rs9923231), concurrent atrial fibrillation status, CYP2C9*3 (rs1057910), body mass index, and sex, to improve warfarin dose adjustment strategies in Chinese patients.

Comprehensive molecular dissection of TIFY Transcription factors reveal their dynamic responses to biotic and abiotic stress in wheat (Triticum aestivum L.)

The plant specific TIFY (previously known as ZIM) transcription factor (TF) family plays crucial roles in cross talk between Jasmonic Acid and other phytohormones like gibberellins, salicylic acid, abscisic acid, auxin, and ethylene signaling pathways. Wheat yield is severely affected by rust diseases and many abiotic stresses, where different phytohormone signaling pathways are involved. TIFYs have been studied in many plants yet reports describing their molecular structure and function in wheat are lacking. In the present study, we have identified 23 novel TIFY genes in wheat genome using in silico approaches. The identified proteins were characterized based on their conserved domains and phylogenetically classified into nine subfamilies. Chromosomal localization of the identified TIFY genes showed arbitrary distribution. Forty cis-acting elements including phytohormone, stress and light receptive elements were detected in the upstream regions of TIFY genes. Seventeen wheat microRNAs targeted the identified wheat TIFY genes. Gene ontological studies revealed their major contribution in defense response and phytohormone signaling. Secondary structure of TIFY proteins displayed the characteristic alpha-alpha-beta fold. Synteny analyses indicated all wheat TIFY genes had orthologous sequences in sorghum, rice, maize, barley and Brachypodium indicating presence of similar TIFY domains in monocot plants. Six TIFY genes had been cloned from wheat genomic and cDNA. Sequence characterization revealed similar characteristics as the in silico identified novel TIFY genes. Tertiary structures predicted the active sites in these proteins to play critical roles in DNA binding. Expression profiling of TIFY genes showed their contribution during incompatible and compatible leaf rust infestation. TIFY genes were also highly expressed during the initial hours of phytohormone induced stress. This study furnishes fundamental information on characterization and putative functions of TIFY genes in wheat.

Necrotrophic Fungus Pyrenophora tritici-repentis Triggers Expression of Multiple Resistance Components in Resistant and Susceptible Wheat Cultivars

Tan spot of wheat, caused by Pyrenophora tritici-repentis (Ptr), results in a yield loss through chlorosis and necrosis of healthy leaf tissue. The major objective of this study was to compare gene expression in resistant and susceptible wheat cultivars after infection with Ptr ToxA-producing race 2 and direct infiltration with Ptr ToxA proteins. Greenhouse experiments included exposure of the wheat cultivars to pathogen inoculum or direct infiltration of leaf tissue with Ptr-ToxA protein isolate. Samples from the experiments were subjected to RNA sequencing. Results showed that ToxA RNA sequences were first detected in samples collected eight hours after treatments indicating that upon Ptr contact with wheat tissue, Ptr started expressing ToxA. The resistant wheat cultivar, in response to Ptr inoculum, expressed genes associated with plant resistance responses that were not expressed in the susceptible cultivar; genes of interest included five chitinases, eight transporters, five pathogen-detecting receptors, and multiple classes of signaling factors. Resistant and susceptible wheat cultivars therefore differed in their response in the expression of genes that encode chitinases, transporters, wall-associated kinases, permeases, and wound-induced proteins, among others. Plants exposed to Ptr inoculum expressed transcription factors, kinases, receptors, and peroxidases, which are not expressed as highly in the control samples or samples infiltrated with ToxA. Several of the differentially expressed genes between cultivars were found in the Ptr resistance QTLs on chromosomes 1A, 2D, 3B, and 5A. Future studies should elucidate the specific roles these genes play in the wheat response to Ptr.

Genome-wide transcriptome and proteome profiles indicate an active role of alternative splicing during de-etiolation of maize seedlings

AS events affect genes encoding protein domain composition and make the single gene produce more proteins with a certain number of genes to satisfy the establishment of photosynthesis during de-etiolation. The drastic switch from skotomorphogenic to photomorphogenic development is an excellent system to elucidate rapid developmental responses to environmental stimuli in plants. To decipher the effects of different light wavelengths on de-etiolation, we illuminated etiolated maize seedlings with blue, red, blue-red mixed and white light, respectively. We found that blue light alone has the strongest effect on photomorphogenesis and that this effect can be attributed to the higher number and expression levels of photosynthesis and chlorosynthesis proteins. Deep sequencing-based transcriptome analysis revealed gene expression changes under different light treatments and a genome-wide alteration in alternative splicing (AS) profiles. We discovered 41,188 novel transcript isoforms for annotated genes, which increases the percentage of multi-exon genes with AS to 63% in maize. We provide peptide support for all defined types of AS, especially retained introns. Further in silico prediction revealed that 58.2% of retained introns have changes in domains compared with their most similar annotated protein isoform. This suggests that AS acts as a protein function switch allowing rapid light response through the addition or removal of functional domains. The richness of novel transcripts and protein isoforms also demonstrates the potential and importance of integrating proteomics into genome annotation in maize.

Evolution and Expression of the Membrane Attack Complex and Perforin Gene Family in the Poaceae

Membrane Attack Complex and Perforin (MACPF) proteins play crucial roles in plant development and plant responses to environmental stresses. To date, only four MACPF genes have been identified in Arabidopsis thaliana, and the functions of the MACPF gene family members in other plants, especially in important crop plants, such as the Poaceae family, remain largely unknown. In this study, we identified and analyzed 42 MACPF genes from six completely sequenced and well annotated species representing the major Poaceae clades. A phylogenetic analysis of MACPF genes resolved four groups, characterized by shared motif organizations and gene structures within each group. MACPF genes were unevenly distributed along the Poaceae chromosomes. Moreover, segmental duplications and dispersed duplication events may have played significant roles during MACPF gene family expansion and functional diversification in the Poaceae. In addition, phylogenomic synteny analysis revealed a high degree of conservation among the Poaceae MACPF genes. In particular, Group I, II, and III MACPF genes were exposed to strong purifying selection with different evolutionary rates. Temporal and spatial expression analyses suggested that Group III MACPF genes were highly expressed relative to the other groups. In addition, most MACPF genes were highly expressed in vegetative tissues and up-regulated by several biotic and abiotic stresses. Taken together, these findings provide valuable information for further functional characterization and phenotypic validation of the Poaceae MACPF gene family.

Transcriptome Analysis Reveals Long Intergenic Non-Coding RNAs Contributed to Intramuscular Fat Content Differences between Yorkshire and Wei Pigs

Intramuscular fat (IMF) content is closely related to various meat traits, such as tenderness, juiciness, and flavor. The IMF content varies considerably among pig breeds with different genetic backgrounds. Long intergenic non-coding RNAs (lincRNAs) have been widely identified in many species and found to be an important class of regulators that can participate in multiple biological processes. However, the mechanism behind lincRNAs regulation of pig IMF content remains unknown and requires further study. In our study, we identified a total of 156 lincRNAs in the longissimus dorsi muscle of Wei (fat-type) and Yorkshire (lean-type) pigs using previously published data. These identified lincRNAs have shorter transcript length, longer exon length, lower exon number, and lower expression level as compared with protein-coding transcripts. We predicted potential target genes (PTGs) that are potentially regulated by lincRNAs in cis or trans regulation. Gene ontology and pathway analyses indicated that many potential lincRNAs target genes are involved in IMF-related processes or pathways, such as fatty acid catabolic process and adipocytokine signaling pathway. In addition, we analyzed quantitative trait locus (QTL) sites that differentially expressed lincRNAs (DE lincRNAs) between Wei and Yorkshire pigs co-localized. The QTL sites where DE lincRNAs co-localize are mostly related to IMF content. Furthermore, we constructed a co-expressed network between DE lincRNAs and their differentially expressed PTGs (DEPTGs). On the basis of their expression levels, we suggest that many DE lincRNAs can affect IMF development by positively or negatively regulating their PTGs. This study identified and analyzed some lincRNAs- and PTGs-related IMF development of the two pig breeds and provided new insight into research on the roles of lincRNAs in the two types of breeds.

Alternative splicing of a barley gene results in an excess-tillering and semi-dwarf mutant

An excess-tillering semi-dwarf gene Hvhtd was identified from an EMS-induced mutant in barley and alternative splicing results in excess-tillering semi-dwarf traits. Tillering and plant height are important traits determining plant architecture and grain production in cereal crops. This study identified an excess-tillering semi-dwarf mutant (htd) from an EMS-treated barley population. Genetic analysis of the F₁, F₂, and F_2:3 populations showed that a single recessive gene controlled the excess-tillering semi-dwarf in htd. Using BSR-Seq and gene mapping, the Hvhtd gene was delimited within a 1.8 Mb interval on chromosome 2HL. Alignment of the RNA-Seq data with the functional genes in the interval identified a gene HORVU2Hr1G098820 with alternative splicing between exon2 and exon3 in the mutant, due to a G to A single-nucleotide substitution at the exon and intron junction. An independent mutant with a similar phenotype confirmed the result, with alternative splicing between exon3 and exon4. In both cases, the alternative splicing resulted in a non-functional protein. And the gene HORVU2Hr1G098820 encodes a trypsin family protein and may be involved in the IAA signaling pathway and differs from the mechanism of Green Revolution genes in the gibberellic acid metabolic pathway.

Impact of demography on linked selection in two outcrossing Brassicaceae species

Genetic diversity is shaped by mutation, genetic drift, gene flow, recombination, and selection. The dynamics and interactions of these forces shape genetic diversity across different parts of the genome, between populations and species. Here, we have studied the effects of linked selection on nucleotide diversity in outcrossing populations of two Brassicaceae species, Arabidopsis lyrata and Capsella grandiflora, with contrasting demographic history. In agreement with previous estimates, we found evidence for a modest population size expansion thousands of generations ago, as well as efficient purifying selection in C. grandiflora. In contrast, the A. lyrata population exhibited evidence for very recent strong population size decline and weaker efficacy of purifying selection. Using multiple regression analyses with recombination rate and other genomic covariates as explanatory variables, we can explain 47% of the variance in neutral diversity in the C. grandiflora population, while in the A. lyrata population, only 11% of the variance was explained by the model. Recombination rate had a significant positive effect on neutral diversity in both species, suggesting that selection at linked sites has an effect on patterns of neutral variation. In line with this finding, we also found reduced neutral diversity in the vicinity of genes in the C. grandiflora population. However, in A. lyrata no such reduction in diversity was evident, a finding that is consistent with expectations of the impact of a recent bottleneck on patterns of neutral diversity near genes. This study thus empirically demonstrates how differences in demographic history modulate the impact of selection at linked sites in natural populations.

Arabidopsis EMSY-like (EML) histone readers are necessary for post-fertilization seed development, but prevent fertilization-independent seed formation

Seed development is a complex regulatory process that includes a transition from gametophytic to sporophytic program. The synchronized development of different seed compartments (seed coat, endosperm and embryo) depends on a balance in parental genome contributions. In the most severe cases, the disruption of the balance leads to seed abortion. This represents one of the main obstacles for the engineering of asexual reproduction through seeds (apomixis), and for generating new interspecies hybrids. The repression of auxin synthesis by the Polycomb Repressive Complex 2 (PRC2) is a major mechanism contributing to sensing genome balance. However, current efforts focusing on decreasing PRC2 or elevating auxin levels have not yet resulted in the production of apomictic seed. Here, we show that EMSY-Like Tudor/Agenet H3K36me3 histone readers EML1 and EML3 are necessary for early stages of seed development to proceed at a normal rate in Arabidopsis. We further demonstrate that both EML1 and EML3 are required to prevent the initiation of seed development in the absence of fertilization. Based on the whole genome expression analysis, we identify auxin transport and signaling genes as the most enriched downstream targets of EML1 and EML3. We hypothesize that EML1 and EML3 function to repress and soften paternal gene expression by fine-tuning auxin responses. Discovery of this pathway may contribute to the engineering of apomixis and interspecies hybrids.

Hybridization preceded radiation in diploid wheats

Evolutionary relationships among the Aegilops-Triticum relatives of cultivated wheats have been difficult to resolve owing to incomplete lineage sorting and reticulate evolution. Recent studies have suggested that the wheat D-genome lineage (progenitor of Ae. tauschii) originated through homoploid hybridization between the A-genome lineage (progenitor of Triticum s.str.) and the B-genome lineage (progenitor of Ae. speltoides). This scenario of reticulation has been debated, calling for adequate phylogenetic analyses based on comprehensive sampling. To reconstruct the evolution of Aegilops-Triticum diploids, we here combined high-throughput sequencing of 38 nuclear low-copy loci of multiple accessions of all 13 species with inferences of the species phylogeny using the full-parameterized MCMC_SEQ method. Phylogenies recovered a monophyletic Aegilops-Triticum lineage that began diversifying ~6.6 Ma ago and gave rise to four sublineages, i.e. the A- (2 species), B- (1 species), D- (9 species) and T- (Ae. mutica) genome lineage. Full-parameterized phylogenies as well as patterns of tree dilation and tree compression supported a hybrid origin of the D-genome lineage from A and B ~3.0-4.0 Ma ago, and did not indicate additional hybridization events. Conflicting ABBA-BABA tests suggestive of further reticulation were shown here to result from ancestral population structure rather than hybridization. This comprehensive and dated phylogeny of wheat relatives indicates that the origin of the hybrid D-genome was followed by intense diversification into the majority of extant diploid as well as allopolyploid wild wheats.

Genome-Wide Identification and Characterization of the Wheat Remorin (TaREM) Family during Cold Acclimation

Remorins (REMs) are plant-specific proteins that play an essential role in plant-microbe interactions. However, their roles in vernalization and abiotic stress responses remain speculative. Most remorins have a variable proline-rich -half and a more conserved -half that is predicted to form coils. A search of the wheat ( L.) database revealed the existence of 20 different genes, which we classified into six groups on the basis of whether they shared a common phylogenetic and structural origin. Analysis of the physical genomic distributions demonstrated that genes are dispersed in the wheat genome and have one to seven introns. Promoter analysis of genes revealed the presence of putative -elements related to diverse functions like development, hormonal regulation, and biotic and abiotic stress responsiveness. Expression levels of genes were measured in plants grown under field and controlled conditions and in response to hormone treatment. Our analyses revealed that 12 members of the REM family are regulated during cold acclimation in wheat in four different tissues (roots, crowns, stems, and leaves), with the highest expression in roots. Differential gene expression was found between wheat cultivars with contrasting degrees of cold tolerance, suggesting the implication of genes in cold response and tolerance. Additionally, eight genes were induced in response to abscisic acid and methyl jasmonate treatment. This genome-wide analysis of genes provides valuable resources for functional analysis aimed at understanding their role in stress adaptation.

Identification and Functional Prediction of Long Intergenic Non-coding RNAs Related to Subcutaneous Adipose Development in Pigs

An increasing number of studies have shown that long intergenic non-coding RNAs (lincRNAs) are a very important class of non-coding RNAs that plays a vital role in many biological processes. Adipose tissue is an important place for storing energy, but few studies on lincRNAs were related to pig subcutaneous fat development. Here, we used published RNA-seq data from subcutaneous adipose tissue of Italian Large White pigs and identified 252 putative lincRNAs, wherein 34 were unannotated. These lincRNAs had relatively shorter length, lower number of exons, and lower expression level compared with protein-coding transcripts. Gene ontology and pathway analysis indicated that the adjacent genes of lincRNAs were involved in lipid metabolism. In addition, differentially expressed lincRNAs (DELs) between low and high backfat thickness pigs were identified. Through the detection of quantitative trait locus (QTL), DELs were mainly located in QTLs related to adipose development. Based on the expression correlation of DEL genes and their differentially expressed potential target genes, we constructed a co-expression network and a potential pathway of DEL's effect on lipid metabolism. Our study identified and analyzed lincRNAs in subcutaneous adipose tissue, and results suggested that lincRNAs may be involved in the regulation of subcutaneous fat development. Our findings provided new insights into the biological function of porcine lincRNAs.

Synthetic lethality guiding selection of drug combinations in ovarian cancer

BACKGROUND

Synthetic lethality describes a relationship between two genes where single loss of either gene does not trigger significant impact on cell viability, but simultaneous loss of both gene functions results in lethality. Targeting synthetic lethal interactions with drug combinations promises increased efficacy in tumor therapy.

MATERIALS AND METHODS

We established a set of synthetic lethal interactions using publicly available data from yeast screens which were mapped to their respective human orthologs using information from orthology databases. This set of experimental synthetic lethal interactions was complemented by a set of predicted synthetic lethal interactions based on a set of protein meta-data like e.g. molecular pathway assignment. Based on the combined set, we evaluated drug combinations used in late stage clinical development (clinical phase III and IV trials) or already in clinical use for ovarian cancer with respect to their effect on synthetic lethal interactions. We furthermore identified a set of drug combinations currently not being tested in late stage ovarian cancer clinical trials that however have impact on synthetic lethal interactions thus being worth of further investigations regarding their therapeutic potential in ovarian cancer.

RESULTS

Twelve of the tested drug combinations addressed a synthetic lethal interaction with the anti-VEGF inhibitor bevacizumab in combination with paclitaxel being the most studied drug combination addressing the synthetic lethal pair between VEGFA and BCL2. The set of 84 predicted drug combinations for example holds the combination of the PARP inhibitor olaparib and paclitaxel, which showed efficacy in phase II clinical studies.

CONCLUSION

A set of drug combinations currently not tested in late stage ovarian cancer clinical trials was identified having impact on synthetic lethal interactions thus being worth of further investigations regarding their therapeutic potential in ovarian cancer.

Regulation of Root Angle and Gravitropism

Regulation of plant root angle is critical for obtaining nutrients and water and is an important trait for plant breeding. A plant’s final, long-term root angle is the net result of a complex series of decisions made by a root tip in response to changes in nutrient availability, impediments, the gravity vector and other stimuli. When a root tip is displaced from the gravity vector, the short-term process of gravitropism results in rapid reorientation of the root toward the vertical. Here, we explore both short- and long-term regulation of root growth angle, using natural variation in tomato to identify shared and separate genetic features of the two responses. Mapping of expression quantitative trait loci mapping and leveraging natural variation between and within species including Arabidopsis suggest a role for PURPLE ACID PHOSPHATASE 27 and CELL DIVISION CYCLE 73 in determining root angle.

Synthetic lethal combinations of low-toxicity drugs for breast cancer identified in silico by genetic screens in yeast

In recent years, the concept of synthetic lethality, describing a cellular state where loss of two genes leads to a non-viable phenotype while loss of one gene can be compensated, has emerged as a novel strategy for cancer therapy. Various compounds targeting synthetic lethal pathways are either under clinical investigation or are already routinely used in multiple cancer entities such as breast cancer. Most of them target the well-described synthetic lethal interplay between PARP1 and BRCA1/2. In our study, we investigated, using an in silico methodological approach, clinically utilized drug combinations for breast cancer treatment, by correlating their known molecular targets with known homologous interaction partners that cause synthetic lethality in yeast. Further, by creating a machine-learning algorithm, we were able to suggest novel synthetic lethal drug combinations of low-toxicity drugs in breast cancer and showed their negative effects on cancer cell viability in vitro. Our findings foster the understanding of evolutionarily conserved synthetic lethality in breast cancer cells and might lead to new drug combinations with favorable toxicity profile in this entity.

Comparative genomics reveals the distinct evolutionary trajectories of the robust and complex coral lineages

BACKGROUND

Despite the biological and economic significance of scleractinian reef-building corals, the lack of large molecular datasets for a representative range of species limits understanding of many aspects of their biology. Within the Scleractinia, based on molecular evidence, it is generally recognised that there are two major clades, Complexa and Robusta, but the genomic bases of significant differences between them remain unclear.

RESULTS

Draft genome assemblies and annotations were generated for three coral species: Galaxea fascicularis (Complexa), Fungia sp., and Goniastrea aspera (Robusta). Whilst phylogenetic analyses strongly support a deep split between Complexa and Robusta, synteny analyses reveal a high level of gene order conservation between all corals, but not between corals and sea anemones or between sea anemones. HOX-related gene clusters are, however, well preserved across all of these combinations. Differences between species are apparent in the distribution and numbers of protein domains and an apparent correlation between number of HSP20 proteins and stress tolerance. Uniquely amongst animals, a complete histidine biosynthesis pathway is present in robust corals but not in complex corals or sea anemones. This pathway appears to be ancestral, and its retention in the robust coral lineage has important implications for coral nutrition and symbiosis.

CONCLUSIONS

The availability of three new coral genomes enabled recognition of a de novo histidine biosynthesis pathway in robust corals which is only the second identified biosynthetic difference between corals. These datasets provide a platform for understanding many aspects of coral biology, particularly the interactions of corals with their endosymbionts.

Whole genome screen reveals a novel relationship between Wolbachia levels and Drosophila host translation

Wolbachia is an intracellular bacterium that infects a remarkable range of insect hosts. Insects such as mosquitos act as vectors for many devastating human viruses such as Dengue, West Nile, and Zika. Remarkably, Wolbachia infection provides insect hosts with resistance to many arboviruses thereby rendering the insects ineffective as vectors. To utilize Wolbachia effectively as a tool against vector-borne viruses a better understanding of the host-Wolbachia relationship is needed. To investigate Wolbachia-insect interactions we used the Wolbachia/Drosophila model that provides a genetically tractable system for studying host-pathogen interactions. We coupled genome-wide RNAi screening with a novel high-throughput fluorescence in situ hybridization (FISH) assay to detect changes in Wolbachia levels in a Wolbachia-infected Drosophila cell line JW18. 1117 genes altered Wolbachia levels when knocked down by RNAi of which 329 genes increased and 788 genes decreased the level of Wolbachia. Validation of hits included in depth secondary screening using in vitro RNAi, Drosophila mutants, and Wolbachia-detection by DNA qPCR. A diverse set of host gene networks was identified to regulate Wolbachia levels and unexpectedly revealed that perturbations of host translation components such as the ribosome and translation initiation factors results in increased Wolbachia levels both in vitro using RNAi and in vivo using mutants and a chemical-based translation inhibition assay. This work provides evidence for Wolbachia-host translation interaction and strengthens our general understanding of the Wolbachia-host intracellular relationship.

Systematic miRNome profiling reveals differential microRNAs in transgenic maize metabolism

BACKGROUND

While some genetically modified organisms (GMOs) are created to produce new double-stranded RNA molecules (dsRNA), in others, such molecules may occur as an unintended effect of the genetic engineering process. Furthermore, GMOs might produce naturally occurring dsRNA molecules in higher or lower quantities than its non-transgenic counterpart. This study is the first to use high-throughput technology to characterize the miRNome of commercialized GM maize events and to investigate potential alterations in miRNA regulatory networks.

RESULTS

Thirteen different conserved miRNAs were found to be dys-regulated in GM samples. The insecticide Bt GM variety had the most distinct miRNome. These miRNAs target a range of endogenous transcripts, such as transcription factors and nucleic acid binding domains, which play key molecular functions in basic genetic regulation. In addition, we have identified 20 potential novel miRNAs with target transcripts involved in lipid metabolism in maize. isomiRs were also found in 96 conserved miRNAs sequences, as well as potential transgenic miRNA sequences, which both can be a source of potential off-target effects in the plant genome. We have also provided information on technical limitations and when to carry on additional in vivo experimental testing.

CONCLUSIONS

These findings do not reveal hazards per se but show that robust and reproducible miRNA profiling technique can strengthen the assessment of risk by detecting any new intended and unintended dsRNA molecules, regardless of the outcome, at any stage of GMO development.

Schizophrenia-associated rs4702 G allele-specific downregulation of FURIN expression by miR-338-3p reduces BDNF production

Genome-wide association studies (GWAS) reveal numerous schizophrenia (SCZ)-associated single-nucleotide polymorphisms (SNPs); however, functional characterizations of the risk variants remain to be established. Using data from 108 SCZ GWAS loci, we performed systematic miRNA binding site screening of 128 SCZ-associated SNPs and found that 2 out of 3 SNPs located in the 3'UTR were predicted to alter 3 miRNAs' binding sites in 2 target genes. Of the identified SNPs, the most genome-wide significant SNP rs4702 (A/G) in the FURIN 3'UTR, previously identified as an SCZ-associated cis-expression quantitative trait loci (downregulated by the risk G allele), is located in the binding site of miR-338-3p in the presence of the risk G allele. Allele-specific downregulation of FURIN by miR-338-3p was validated with a luciferase reporter assay. Furthermore, we demonstrated that miR-338-3p-mediated FURIN inhibition reduced brain-derived neurotrophic factor (BDNF) maturation and secretion in human embryonic kidney 293T cells. Our data reveal that schizophrenia-associated rs4702 G allele-specific downregulation of FURIN by miR-338-3p reduces mature BDNF production. These data help elucidate the mechanism of genetic predisposition toward schizophrenia or other neurodevelopmental diseases.

Establishment of a Han Chinese-specific pharmacogenetic-guided warfarin dosing algorithm

Warfarin is the most common oral anticoagulant. Because of a narrow therapeutic range, interindividual differences in drug responses, and the risk of bleeding, there are many challenges in using warfarin. We need to predict the warfarin maintenance dose. However, ethnic-specific algorithms may be required, and some Chinese algorithms do not perform adequately. Therefore, we aimed to establish a Han Chinese appropriate algorithm.We recruited a study group consisting of 361 Han Chinese patients receiving warfarin treatment who had heart valve replacements. Genotyping of 38 single nucleotide polymorphisms (SNPs) in 13 candidate genes was carried out using the MassARRAY. In the derivation cohort, a multiple linear regression model was constructed to predict the warfarin dosage. We evaluated the accuracy of our algorithm in the validation cohort and compared it with the other 5 algorithms based on Han Chinese and other races.We established a Han Chinese-specific pharmacogenetic-guided warfarin dosing algorithm. Warfarin maintenance dosage (mg/day) = 1.787 - 0.023 × (Age) + 1.151 × (BSA [m]) + 0.917 × (VKORC1 AG) + 4.619 × (VKORC1 GG) + 0.595 × (CYP4F2 TT) + 0.707 × (CYP2C19 CC). It explained 58.3% of the variance in warfarin doses in Han Chinese patients and was superior to the other 5 algorithms. The ability of the 6 algorithms which estimate the required dose correctly was tested. Our model had a mean absolute error of 0.74 mg/day, the other 5 models have mean absolute error of 0.81 mg/day,1.05 mg/day, 1.24 mg/day, 1.18 mg/day, and 0.85 mg/day, respectively. Our model had a mean percentage error of 26.9%, the other 5 models have the mean percentage error of 27.7%, 27.2%, 52.3%, 45.7%, and 29.3%, respectively.Physicians can not adopt algorithm from other race directly to predict warfarin dose in patients with heart valve replacements, they should establish a new algorithm or adjust another algorithm to fit their patients. The algorithm established in this study has the potential to assist physicians in determining warfarin doses that are close to the appropriate doses.

Solid phase chemistry to covalently and reversibly capture thiolated RNA

Here, we describe an approach to enrich newly transcribed RNAs from primary mouse neurons using 4-thiouridine (s4U) metabolic labeling and solid phase chemistry. This one-step enrichment procedure captures s4U-RNA by using highly efficient methane thiosulfonate (MTS) chemistry in an immobilized format. Like solution-based methods, this solid-phase enrichment can distinguish mature RNAs (mRNA) with differential stability, and can be used to reveal transient RNAs such as enhancer RNAs (eRNAs) and primary microRNAs (pri-miRNAs) from short metabolic labeling. Most importantly, the efficiency of this solid-phase chemistry made possible the first large scale measurements of RNA polymerase II (RNAPII) elongation rates in mouse cortical neurons. Thus, our approach provides the means to study regulation of RNA metabolism in specific tissue contexts as a means to better understand gene expression in vivo.

Automatic identification of optimal marker genes for phenotypic and taxonomic groups of microorganisms

Finding optimal markers for microorganisms important in the medical, agricultural, environmental or ecological fields is of great importance. Thousands of complete microbial genomes now available allow us, for the first time, to exhaustively identify marker proteins for groups of microbial organisms. In this work, we model the biological task as the well-known mathematical “hitting set” problem, solving it based on both greedy and randomized approximation algorithms. We identify unique markers for 17 phenotypic and taxonomic microbial groups, including proteins related to the nitrite reductase enzyme as markers for the non-anammox nitrifying bacteria group, and two transcription regulation proteins, nusG and yhiF, as markers for the Archaea and Escherichia/Shigella taxonomic groups, respectively. Additionally, we identify marker proteins for three subtypes of pathogenic E. coli, which previously had no known optimal markers. Practically, depending on the completeness of the database this algorithm can be used for identification of marker genes for any microbial group, these marker genes may be prime candidates for the understanding of the genetic basis of the group's phenotype or to help discover novel functions which are uniquely shared among a group of microbes. We show that our method is both theoretically and practically efficient, while establishing an upper bound on its time complexity and approximation ratio; thus, it promises to remain efficient and permit the identification of marker proteins that are specific to phenotypic or taxonomic groups, even as more and more bacterial genomes are being sequenced.

The RNASeq-er API-a gateway to systematically updated analysis of public RNA-seq data

Motivation

The exponential growth of publicly available RNA-sequencing (RNA-Seq) data poses an increasing challenge to researchers wishing to discover, analyse and store such data, particularly those based in institutions with limited computational resources. EMBL-EBI is in an ideal position to address these challenges and to allow the scientific community easy access to not just raw, but also processed RNA-Seq data. We present a Web service to access the results of a systematically and continually updated standardized alignment as well as gene and exon expression quantification of all public bulk (and in the near future also single-cell) RNA-Seq runs in 264 species in European Nucleotide Archive, using Representational State Transfer.

Results

The RNASeq-er API (Application Programming Interface) enables ontology-powered search for and retrieval of CRAM, bigwig and bedGraph files, gene and exon expression quantification matrices (Fragments Per Kilobase Of Exon Per Million Fragments Mapped, Transcripts Per Million, raw counts) as well as sample attributes annotated with ontology terms. To date over 270 00 RNA-Seq runs in nearly 10 000 studies (1PB of raw FASTQ data) in 264 species in ENA have been processed and made available via the API.

Availability and Implementation

The RNASeq-er API can be accessed at http://www.ebi.ac.uk/fg/rnaseq/api . The commands used to analyse the data are available in supplementary materials and at https://github.com/nunofonseca/irap/wiki/iRAP-single-library .

Contact

rnaseq@ebi.ac.uk ; rpetry@ebi.ac.uk.

Supplementary information

Supplementary data are available at Bioinformatics online.

miRDis: a Web tool for endogenous and exogenous microRNA discovery based on deep-sequencing data analysis

Small RNA sequencing is the most widely used tool for microRNA (miRNA) discovery, and shows great potential for the efficient study of miRNA cross-species transport, i.e., by detecting the presence of exogenous miRNA sequences in the host species. Because of the increased appreciation of dietary miRNAs and their far-reaching implication in human health, research interests are currently growing with regard to exogenous miRNAs bioavailability, mechanisms of cross-species transport and miRNA function in cellular biological processes. In this article, we present microRNA Discovery (miRDis), a new small RNA sequencing data analysis pipeline for both endogenous and exogenous miRNA detection. Specifically, we developed and deployed a Web service that supports the annotation and expression profiling data of known host miRNAs and the detection of novel miRNAs, other noncoding RNAs, and the exogenous miRNAs from dietary species. As a proof-of-concept, we analyzed a set of human plasma sequencing data from a milk-feeding study where 225 human miRNAs were detected in the plasma samples and 44 show elevated expression after milk intake. By examining the bovine-specific sequences, data indicate that three bovine miRNAs (bta-miR-378, -181* and -150) are present in human plasma possibly because of the dietary uptake. Further evaluation based on different sets of public data demonstrates that miRDis outperforms other state-of-the-art tools in both detection and quantification of miRNA from either animal or plant sources. The miRDis Web server is available at: http://sbbi.unl.edu/miRDis/index.php.

Fungal Genomes and Insights into the Evolution of the Kingdom

The kingdom Fungi comprises species that inhabit nearly all ecosystems. Fungi exist as both free-living and symbiotic unicellular and multicellular organisms with diverse morphologies. The genomes of fungi encode genes that enable them to thrive in diverse environments, invade plant and animal cells, and participate in nutrient cycling in terrestrial and aquatic ecosystems. The continuously expanding databases of fungal genome sequences have been generated by individual and large-scale efforts such as Génolevures, Broad Institute's Fungal Genome Initiative, and the 1000 Fungal Genomes Project (http://1000.fungalgenomes.org). These efforts have produced a catalog of fungal genes and genomic organization. The genomic datasets can be utilized to better understand how fungi have adapted to their lifestyles and ecological niches. Large datasets of fungal genomic and transcriptomic data have enabled the use of novel methodologies and improved the study of fungal evolution from a molecular sequence perspective. Combined with microscopes, petri dishes, and woodland forays, genome sequencing supports bioinformatics and comparative genomics approaches as important tools in the study of the biology and evolution of fungi.

The Sequenced Angiosperm Genomes and Genome Databases

Angiosperms, the flowering plants, provide the essential resources for human life, such as food, energy, oxygen, and materials. They also promoted the evolution of human, animals, and the planet earth. Despite the numerous advances in genome reports or sequencing technologies, no review covers all the released angiosperm genomes and the genome databases for data sharing. Based on the rapid advances and innovations in the database reconstruction in the last few years, here we provide a comprehensive review for three major types of angiosperm genome databases, including databases for a single species, for a specific angiosperm clade, and for multiple angiosperm species. The scope, tools, and data of each type of databases and their features are concisely discussed. The genome databases for a single species or a clade of species are especially popular for specific group of researchers, while a timely-updated comprehensive database is more powerful for address of major scientific mysteries at the genome scale. Considering the low coverage of flowering plants in any available database, we propose construction of a comprehensive database to facilitate large-scale comparative studies of angiosperm genomes and to promote the collaborative studies of important questions in plant biology.

A phylogenomic analysis of lichen-feeding tiger moths uncovers evolutionary origins of host chemical sequestration

Host species utilize a variety of defenses to deter feeding, including secondary chemicals. Some phytophagous insects have evolved tolerance to these chemical defenses, and can sequester secondary defense compounds for use against their own predators and parasitoids. While numerous studies have examined plant-insect interactions, little is known about lichen-insect interactions. Our study focused on reconstructing the evolution of lichen phenolic sequestration in the tiger moth tribe Lithosiini (Lepidoptera: Erebidae: Arctiinae), the most diverse lineage of lichen-feeding moths, with 3000 described species. We built an RNA-Seq dataset and examined the adult metabolome for the presence of lichen-derived phenolics. Using the transcriptomic dataset, we recover a well-resolved phylogeny of the Lithosiini, and determine that the metabolomes within species are more similar than those among species. Results from an initial ancestral state reconstruction suggest that the ability to sequester phenolics produced by a single chemical pathway preceded generalist sequestration of phenolics produced by multiple chemical pathways. We conclude that phenolics are consistently and selectively sequestered within Lithosiini. Furthermore, sequestration of compounds from a single chemical pathway may represent a synapomorphy of the tribe, and the ability to sequester phenolics produced by multiple pathways arose later. These findings expand on our understanding of the interactions between Lepidoptera and their lichen hosts.

Genome-Wide Analysis of Colonization History and Concomitant Selection in Arabidopsis lyrata

The high climatic variability in the past hundred thousand years has affected the demographic and adaptive processes in many species, especially in boreal and temperate regions undergoing glacial cycles. This has also influenced the patterns of genome-wide nucleotide variation, but the details of these effects are largely unknown. Here we study the patterns of genome-wide variation to infer colonization history and patterns of selection of the perennial herb species Arabidopsis lyrata, in locally adapted populations from different parts of its distribution range (Germany, UK, Norway, Sweden, and USA) representing different environmental conditions. Using site frequency spectra based demographic modeling, we found strong reduction in the effective population size of the species in general within the past 100,000 years, with more pronounced effects in the colonizing populations. We further found that the northwestern European A. lyrata populations (UK and Scandinavian) are more closely related to each other than with the Central European populations, and coalescent based population split modeling suggests that western European and Scandinavian populations became isolated relatively recently after the glacial retreat. We also highlighted loci showing evidence for local selection associated with the Scandinavian colonization. The results presented here give new insights into postglacial Scandinavian colonization history and its genome-wide effects.

BRUTUS and its paralogs, BTS LIKE1 and BTS LIKE2, encode important negative regulators of the iron deficiency response in Arabidopsis thaliana

Iron (Fe) is required for plant health, but it can also be toxic when present in excess. Therefore, Fe levels must be tightly controlled. The Arabidopsis thaliana E3 ligase BRUTUS (BTS) is involved in the negative regulation of the Fe deficiency response and we show here that the two A. thaliana BTS paralogs, BTS LIKE1 (BTSL1) and BTS LIKE2 (BTSL2) encode proteins that act redundantly as negative regulators of the Fe deficiency response. Loss of both of these E3 ligases enhances tolerance to Fe deficiency. We further generated a triple mutant with loss of both BTS paralogs and a partial loss of BTS expression that exhibits even greater tolerance to Fe-deficient conditions and increased Fe accumulation without any resulting Fe toxicity effects. Finally, we identified a mutant carrying a novel missense mutation of BTS that exhibits an Fe deficiency response in the root when grown under both Fe-deficient and Fe-sufficient conditions, leading to Fe toxicity when plants are grown under Fe-sufficient conditions.

A Genetic Variant Ameliorates β-Thalassemia Severity by Epigenetic-Mediated Elevation of Human Fetal Hemoglobin Expression

A delayed fetal-to-adult hemoglobin (Hb) switch ameliorates the severity of β-thalassemia and sickle cell disease. The molecular mechanism underlying the epigenetic dysregulation of the switch is unclear. To explore the potential cis-variants responsible for the Hb switching, we systematically analyzed an 80-kb region spanning the β-globin cluster using capture-based next-generation sequencing of 1142 Chinese β-thalassemia persons and identified 31 fetal hemoglobin (HbF)-associated haplotypes of the selected 28 tag regulatory single-nucleotide polymorphisms (rSNPs) in seven linkage disequilibrium (LD) blocks. A Ly1 antibody reactive (LYAR)-binding motif disruptive rSNP rs368698783 (G/A) from LD block 5 in the proximal promoter of hemoglobin subunit gamma 1 (HBG1) was found to be a significant predictor for β-thalassemia clinical severity by epigenetic-mediated variant-dependent HbF elevation. We found this rSNP accounted for 41.6% of β-hemoglobinopathy individuals as an ameliorating factor in a total of 2,738 individuals from southern China and Thailand. We uncovered that the minor allele of the rSNP triggers the attenuation of LYAR and two repressive epigenetic regulators DNA methyltransferase 3 alpha (DNMT3A) and protein arginine methyltransferase 5 (PRMT5) from the HBG promoters, mediating allele-biased γ-globin elevation by facilitating demethylation of HBG core promoter CpG sites in erythroid progenitor cells from β-thalassemia persons. The present study demonstrates that this common rSNP in the proximal ^Aγ-promoter is a major genetic modifier capable of ameliorating the severity of thalassemia major through the epigenetic-mediated regulation of the delayed fetal-to-adult Hb switch and provides potential targets for the treatment of β-hemoglobinopathy.

Data on the genome-wide identification of CNL R-genes in Setaria italica (L.) P. Beauv

We report data associated with the identification of 242 disease resistance genes (R-genes) in the genome of Setaria italica as presented in “Genetic diversity of disease resistance genes in foxtail millet (Setaria italica L.)” (Andersen and Nepal, 2017) [1]. Our data describe the structure and evolution of the Coiled-coil, Nucleotide-binding site, Leucine-rich repeat (CNL) R-genes in foxtail millet. The CNL genes were identified through rigorous extraction and analysis of recently available plant genome sequences using cutting-edge analytical software. Data visualization includes gene structure diagrams, chromosomal syntenic maps, a chromosomal density plot, and a maximum-likelihood phylogenetic tree comparing Sorghum bicolor, Panicum virgatum, Setaria italica, and Arabidopsis thaliana. Compilation of InterProScan annotations, Gene Ontology (GO) annotations, and Basic Local Alignment Search Tool (BLAST) results for the 242 R-genes identified in the foxtail millet genome are also included in tabular format.

Identification of introns harboring functional sequence elements through positional conservation

Many human introns carry out a function, in the sense that they are critical to maintain normal cellular activity. Their identification is fundamental to understanding cellular processes and disease. However, being noncoding elements, such functional introns are poorly predicted based on traditional approaches of sequence and structure conservation. Here, we generated a dataset of human functional introns that carry out different types of functions. We showed that functional introns share common characteristics, such as higher positional conservation along the coding sequence and reduced loss rates, regardless of their specific function. A unique property of the data is that if an intron is unknown to be functional, it still does not mean that it is indeed non-functional. We developed a probabilistic framework that explicitly accounts for this unique property, and predicts which specific human introns are functional. We show that we successfully predict function even when the algorithm is trained on introns with a different type of function. This ability has many implications in studying regulatory networks, gene regulation, the effect of mutations outside exons on human disease, and on our general understanding of intron evolution and their functional exaptation in mammals.

Meta Analysis of Human AlzGene Database: Benefits and Limitations of Using C. elegans for the Study of Alzheimer's Disease and Co-morbid Conditions

Human genome-wide association studies (GWAS) and linkage studies have identified 695 genes associated with Alzheimer's disease (AD), the vast majority of which are associated with late-onset AD. Although orthologs of these AD genes have been studied in several model species, orthologs in the nematode, Caenorhabditis elegans, remain incompletely identified, with orthologs to only 17 AD-related genes identified in the C. elegans database, WormBase. Therefore, we performed a comprehensive search for additional C. elegans orthologs of AD genes using well-established programs, including OrthoList, which utilizes four ontology prediction programs. We also validated 680 of the AD genes as a unique gene from the AlzGene database, including 431 genes (63%) that are predicted to have orthologs in C. elegans. Another 178 human AD genes (26%) were associated with one or more other neurological diseases, including amyotrophic lateral sclerosis, multiple sclerosis, Parkinson's disease, and schizophrenia. Of these, there were 105 genes (59%) with orthologs in C. elegans. Interestingly, three AD genes (ACE, TNF, and MTHFR) were associated with all four of the other neurological diseases. The human AD genes were enriched in three major ontology pathway groups, including lipoprotein metabolism, hemostasis, and extracellular matrix organizations, as well as in pathways that are amyloid related (NOTCH signaling) and associated with neural (neurotransmitter clearance) and immune (advanced glycation end-product receptor signaling and TRAF6-NF-kappaB) systems. Thus, the results from this study provide a potentially useful system for assessing comorbidities that may be associated with late-onset AD and other neurological conditions. The technical advantages and limitations of the ortholog searches are further discussed.