New developments in the InterPro database

Computational tools for plant genomics and breeding

Plant genomics and crop breeding are at the intersection of biotechnology and information technology. Driven by a combination of high-throughput sequencing, molecular biology and data science, great advances have been made in omics technologies at every step along the central dogma, especially in genome assembling, genome annotation, epigenomic profiling, and transcriptome profiling. These advances further revolutionized three directions of development. One is genetic dissection of complex traits in crops, along with genomic prediction and selection. The second is comparative genomics and evolution, which open up new opportunities to depict the evolutionary constraints of biological sequences for deleterious variant discovery. The third direction is the development of deep learning approaches for the rational design of biological sequences, especially proteins, for synthetic biology. All three directions of development serve as the foundation for a new era of crop breeding where agronomic traits are enhanced by genome design.

Genome-Wide Identification and Characterization of CDPK Gene Family in Cultivated Peanut (Arachis hypogaea L.) Reveal Their Potential Roles in Response to Ca Deficiency

This study identified 45 calcium-dependent protein kinase (CDPK) genes in cultivated peanut (Arachis hypogaea L.), which are integral in plant growth, development, and stress responses. These genes, classified into four subgroups based on phylogenetic relationships, are unevenly distributed across all twenty peanut chromosomes. The analysis of the genetic structure of AhCDPKs revealed significant similarity within subgroups, with their expansion primarily driven by whole-genome duplications. The upstream promoter sequences of AhCDPK genes contained 46 cis-acting regulatory elements, associated with various plant responses. Additionally, 13 microRNAs were identified that target 21 AhCDPK genes, suggesting potential post-transcriptional regulation. AhCDPK proteins interacted with respiratory burst oxidase homologs, suggesting their involvement in redox signaling. Gene ontology and KEGG enrichment analyses affirmed AhCDPK genes' roles in calcium ion binding, protein kinase activity, and environmental adaptation. RNA-seq data revealed diverse expression patterns under different stress conditions. Importantly, 26 AhCDPK genes were significantly induced when exposed to Ca deficiency during the pod stage. During the seedling stage, four AhCDPKs (AhCDPK2/-25/-28/-45) in roots peaked after three hours, suggesting early signaling roles in pod Ca nutrition. These findings provide insights into the roles of CDPK genes in plant development and stress responses, offering potential candidates for predicting calcium levels in peanut seeds.

Functional characterization of all-trans retinoic acid-induced differentiation factor (ATRAID)

All-trans retinoic acid-induced differentiation (ATRAID) factor was first identified in HL60 cells. Several mRNA isoforms exist, but the respective proteins have not been fully characterized. In transfected cells expressing Myc-Flag-tagged ATRAID Isoform (Iso) A, B, and C, Iso C was found to be expressed at high levels, Iso A was found to be expressed at low levels due to rapid degradation, and the predicted protein expressed from Iso B was not detected. Iso C was present mainly in an N-glycosylated form. In subcellular fractionation experiments, Iso C localized to the membranous and nuclear fractions, while immunofluorescence analysis revealed that Iso C is located close to the plasma membrane, mainly in cytoplasmic vesicles and in the Golgi area. We confirm that Iso C colocalizes to some extent with endosomal/lysosomal markers LAMP1 and LAMP2. Furthermore, we show that ATRAID co-localizes with RAB11, a GTPase associated with recycling endosomes and implicated in regulating vesicular trafficking.

Teleost Piscidins-In Silico Perspective of Natural Peptide Antibiotics from Marine Sources

Fish, like all other animals, are exposed to constant contact with microbes, both on their skin and on the surfaces of their respiratory and digestive systems. Fish have a system of non-specific immune responses that provides them with initial protection against infection and allows them to survive under normal conditions despite the presence of these potential invaders. However, fish are less protected against invading diseases than other marine vertebrates because their epidermal surface, composed primarily of living cells, lacks the keratinized skin that serves as an efficient natural barrier in other marine vertebrates. Antimicrobial peptides (AMPs) are one type of innate immune protection present in all life forms. AMPs have been shown to have a broader range of biological effects than conventional antibiotics, including antibacterial, antiviral, antiprotozoal, and antifungal effects. Although other AMPs, such as defensins and hepcidins, are found in all vertebrates and are relatively well conserved, piscidins are found exclusively in Teleost fish and are not found in any other animal. Therefore, there is less information on the expression and bioactivity of piscidins than on other AMPs. Piscidins are highly effective against Gram-positive and Gram-negative bacteria that cause disease in fish and humans and have the potential to be used as pharmacological anti-infectives in biomedicine and aquaculture. To better understand the potential benefits and limitations of using these peptides as therapeutic agents, we are conducting a comprehensive study of the Teleost piscidins included in the "reviewed" category of the UniProt database using bioinformatics tools. They all have amphipathic alpha-helical structures. The amphipathic architecture of piscidin peptides and positively charged residues influence their antibacterial activity. These alpha-helices are intriguing antimicrobial drugs due to their stability in high-salt and metal environments. New treatments for multidrug-resistant bacteria, cancer, and inflammation may be inspired by piscidin peptides.

Recent trends in next generation immunoinformatics harnessed for universal coronavirus vaccine design

The ongoing pandemic of coronavirus disease 2019 (COVID-19), caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has globally devastated public health, the economies of many countries and quality of life universally. The recent emergence of immune-escaped variants and scenario of vaccinated individuals being infected has raised the global concerns about the effectiveness of the current available vaccines in transmission control and disease prevention. Given the high rate mutation of SARS-CoV-2, an efficacious vaccine targeting against multiple variants that contains virus-specific epitopes is desperately needed. An immunoinformatics approach is gaining traction in vaccine design and development due to the significant reduction in time and cost of immunogenicity studies and increasing reliability of the generated results. It can underpin the development of novel therapeutic methods and accelerate the design and production of peptide vaccines for infectious diseases. Structural proteins, particularly spike protein (S), along with other proteins have been studied intensively as promising coronavirus vaccine targets. Numbers of promising online immunological databases, tools and web servers have widely been employed for the design and development of next generation COVID-19 vaccines. This review highlights the role of immunoinformatics in identifying immunogenic peptides as potential vaccine targets, involving databases, and prediction and characterization of epitopes which can be harnessed for designing future coronavirus vaccines.

Wild emmer wheat, the progenitor of modern bread wheat, exhibits great diversity in the VERNALIZATION1 gene

Wild emmer wheat is an excellent reservoir of genetic variability that can be utilized to improve cultivated wheat to address the challenges of the expanding world population and climate change. Bearing this in mind, we have collected a panel of 263 wild emmer wheat (WEW) genotypes across the Fertile Crescent. The genotypes were grown in different locations and phenotyped for heading date. Genome-wide association mapping (GWAS) was carried out, and 16 SNPs were associated with the heading date. As the flowering time is controlled by photoperiod and vernalization, we sequenced the VRN1 gene, the most important of the vernalization response genes, to discover new alleles. Unlike most earlier attempts, which characterized known VRN1 alleles according to a partial promoter or intron sequences, we obtained full-length sequences of VRN-A1 and VRN-B1 genes in a panel of 95 wild emmer wheat from the Fertile Crescent and uncovered a significant sequence variation. Phylogenetic analysis of VRN-A1 and VRN-B1 haplotypes revealed their evolutionary relationships and geographic distribution in the Fertile Crescent region. The newly described alleles represent an attractive resource for durum and bread wheat improvement programs.

Comparative genomics of Sarcoptes scabiei provide new insights into adaptation to permanent parasitism and within-host species divergence

Sarcoptes scabiei is the causative agent of a highly contagious skin disease in humans and more than 100 mammals. Here, we reported the first chromosome-level reference genome of S. scabiei isolated from rabbits, with a contig N50 size of 5.92 Mb, a total assembled length of 57.30 Mb, ∼12.65% repetitive sequences and 9333 predicted protein-coding genes. The phylogenetic tree based on 1338 shared high-confidence single-copy orthologous genes estimated that the mammalian ectoparasite S. scabiei and the plant-feeding mite Tetranychus urticae separated ∼340 million years ago. Both neighbour-joining tree and principal component analysis of 20 mite populations isolated from four hosts (humans, pigs, dogs and rabbits) distributed in three countries (China, Australia and the USA) consistently supported genetic subdivisions according to host species rather than geographical location. The demographic history of S. scabiei reconstructed by multiple sequentially Markovian coalescent analysis suggested that S. scabiei isolated from rabbits, humans, dogs and pigs diverged ∼5000 years ago. Investigation of the homeobox (Hox) genes revealed that S. scabiei contains 8 of 10 canonical Hox genes that are present in the arthropod ancestor, and the absence of the Abd-A gene may correlate with the long gap between their front and back legs. Comparative genomics demonstrated that genes specific to scabies mites were mainly enriched in nutrition digestive systems, whereas genes in the families that involved detoxification (cytochrome P450, carboxyl/cholinesterases and the ATP-binding cassette transporter C group) were extremely contracted compared with that of other mites analysed in this study. Selective sweep analysis of mite populations from various pairs of two out of the four host species revealed that the strongest selective sweep signals were mainly enriched in cysteine-type peptidase activity and apoptosis. The results provided clues for the mechanisms of S. scabiei adaptation to a permanent parasitic lifestyle and knowledge that would enable further control of this highly contagious skin disease.

1H, 13C and 15N resonance assignments and solution structures of the two RRM domains of Matrin-3

Matrin-3 is a multifunctional protein that binds to both DNA and RNA. Its DNA-binding activity is linked to the formation of the nuclear matrix and transcriptional regulation, while its RNA-binding activity is linked to mRNA metabolism including splicing, transport, stabilization, and degradation. Correspondingly, Matrin-3 has two zinc finger domains for DNA binding and two consecutive RNA recognition motif (RRM) domains for RNA binding. Matrin-3 has been reported to cause amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD) when its disordered region contains pathogenic mutations. Simultaneously, it has been shown that the RNA-binding activity of Matrin-3 mediated by its RRM domains, affects the formation of insoluble cytoplasmic granules, which are related to the pathogenic mechanism of ALS/FTD. Thus, the effect of the RRM domains on the phase separation of condensed protein/RNA mixtures has to be clarified for a comprehensive understanding of ALS/FTD. Here, we report the ¹H, ¹⁵N, and ¹³C resonance assignments of the two RNA binding domains and their solution structures. The resonance assignments and the solution structures obtained in this work will contribute to the elucidation of the molecular basis of Matrin-3 in the pathogenic mechanism of ALS and/or FTD.

Expansion of bilin-based red light sensors in the subaerial desert cyanobacterium Nostoc flagelliforme

Light is the crucial environmental signal for desiccation-tolerant cyanobacteria to activate photosynthesis and prepare for desiccation at dawn. However, the photobiological characteristics of desert cyanobacteria adaptation to one of the harshest habitats on Earth remain unresolved. In this study, we surveyed the genome of a subaerial desert cyanobacterium Nostoc flagelliforme and identified two phytochromes and seven cyanobacteriochromes (CBCRs) with one or more bilin-binding GAF (cGMP phosphodiesterase/adenylyl cyclase/FhlA) domains. Biochemical and spectroscopic analyses of 69 purified GAF-containing proteins from recombinant phycocyanobilin (PCB), biliverdin or phycoerythrobilin-producing Escherichia coli indicated that nine of these proteins bind chromophores. Further investigation revealed that 11 GAFs form covalent adducts responsive to near-UV and visible light: eight GAFs contained PCB chromophores, three GAFs contained biliverdin chromophores and one contained the PCB isomer, phycoviolobilin. Interestingly, COO91_03972 is the first-ever reported GAF-only CBCR capable of sensing five wavelengths of light. Bioinformatics and biochemical analyses revealed that residue P132 of COO91_03972 is essential for chromophore binding to dual-cysteine CBCRs. Furthermore, the complement of N. flagelliforme CBCRs is enriched in red light sensors. We hypothesize that these sensors are critical for the acclimatization of N. flagelliforme to weak light environments at dawn.

Network metrics, structural dynamics and density functional theory calculations identified a novel Ursodeoxycholic Acid derivative against therapeutic target Parkin for Parkinson's disease

•

GIN analysis revealed PARK2, LRRK2, PARK7, PINK1 and SNCA as hub-genes.

•

Topologically favoured Parkin was considered as a therapeutic target.

•

ADMET screening identified a novel UDCA derivative as potential lead candidate.

•

Chemical reactivity and ligand stability were analysed through DFT simulation.

•

Docking and MDS established novel lead as potential Parkin inhibitor.

Parkinson's disease (PD) has been designated as one of the priority neurodegenerative disorders worldwide. Although diagnostic biomarkers have been identified, early onset detection and targeted therapy are still limited. An integrated systems and structural biology approach were adopted to identify therapeutic targets for PD. From a set of 49 PD associated genes, a densely connected interactome was constructed. Based on centrality indices, degree of interaction and functional enrichments, LRRK2, PARK2, PARK7, PINK1 and SNCA were identified as the hub-genes. PARK2 (Parkin) was finalized as a potent theranostic candidate marker due to its strong association (score > 0.99) with α-synuclein (SNCA), which directly regulates PD progression. Besides, modeling and validation of Parkin structure, an extensive virtual-screening revealed small (commercially available) inhibitors against Parkin. Molecule-258 (ZINC5022267) was selected as a potent candidate based on pharmacokinetic profiles, Density Functional Theory (DFT) energy calculations (ΔE = 6.93 eV) and high binding affinity (Binding energy = -6.57 ± 0.1 kcal/mol; Inhibition constant = 15.35 µM) against Parkin. Molecular dynamics simulation of protein-inhibitor complexes further strengthened the therapeutic propositions with stable trajectories (low structural fluctuations), hydrogen bonding patterns and interactive energies (>0kJ/mol). Our study encourages experimental validations of the novel drug candidate to prevent the auto-inhibition of Parkin mediated ubiquitination in PD.

Chromosome-scale assembly and high-density genetic map of the yellow drum, Nibea albiflora

The yellow drum (Nibea albiflora) is an economically important sciaenid fish in East Asian countries. In this study, we sequenced and assembled a near-complete gynogenetic yellow drum genome. We generated 45.63 Gb of Illumina short-reads and 80.27 Gb of PacBio long-reads and assembled them into a 628.01-Mb genome with a contig N50 of 4.42 Mb. Twenty-four chromosomes with a scaffold N50 of 26.73 Mb were obtained using the Hi-C analysis. We predicted a set of 27,069 protein-coding genes, of which 1,581 and 2,583 were expanded and contracted gene families, respectively. The most expanded genes were categorised into the protein binding, zinc-ion binding and ATP binding functional pathways. We built a high-density genetic linkage map that spanned 4,300.2 cM with 24 linkage groups and a resolution of 0.69 cM. The high-quality reference genome and annotated profiles that we produced will not only increase our understanding of the genetic architecture of economic traits in the yellow drum, but also help us explore the evolution and unique biological characteristics of sciaenid fishes.

A chromosome-scale draft genome sequence of horsegram (Macrotyloma uniflorum)

Horsegram (Macrotyloma uniflorum [Lam.] Verdc.) is an underutilized warm-season diploid legume (2n = 20, 22). Because of its ability to grow under water-deficient and marginal soil conditions, horsegram is a preferred choice in the era of global climate change. In recognition of its potential as a crop species, we generated and analyzed a draft genome sequence for a horsegram variety, HPK-4. Ten chromosome-scale pseudomolecules were created by aligning Illumina scaffold sequences onto a linkage map. The total length of the ten pseudomolecules was 259.2 Mbp, covering 89% of the total length of the assembled sequences. A total of 36,105 genes were predicted on the assembled sequences. Diversity analysis of 89 horsegram accessions by dd-RAD-Seq identified 277 single nucleotide polymorphisms (SNPs), suggesting narrow genetic diversity among the horsegram accessions. This is the first attempt to generate a draft genome sequence of horsegram and will provide a reference for sequence-based analysis of horsegram germplasm.

Understanding the Role of Populus ECERIFERUM2-Likes in the Biosynthesis of Very-Long-Chain Fatty Acids for Cuticular Waxes

Cuticular waxes are derived from very-long-chain fatty acid (VLCFA) precursors made by the concerted action of four enzymes that form the fatty acid (FA) elongation complex. The condensing enzyme of the complex confers specificity to substrates of different chain lengths, yet on its own cannot account for the biosynthesis of VLCFAs longer than 28 carbons (C28). Recent evidence from Arabidopsis thaliana points to a synergistic role of clade II BAHD acyltransferases and condensing enzymes in the elongation of VLCFAs beyond C28. In Populus trichocarpa, clade II is composed of seven uncharacterized paralogous genes (PtCER2-like1-7). In the present study, five of these genes were heterologously expressed in yeast and their respective FA profiles were determined. PtCER2-likes differentially altered the accumulation of C28 and C30 FAs when expressed in the presence of the condensing enzyme AtCER6. Among these, PtCER2-like5 produced the highest levels of C28 FAs in yeast and its expression was localized to the epidermis in β-glucuronidase-reporter poplar lines, consistent with a role in cuticular wax biosynthesis. Complementation of the A. thaliana cer2-5 mutant with PtCER2-like5 increased the levels of C28-derived cuticular waxes at the expense of C30-derived components. Together, these results demonstrate that the role of CER2-likes in cuticular wax biosynthesis is conserved in Populus clade II BAHD acyltransferases.

Genomic and transcriptomic insights into the habitat adaptation of the diazotrophic paddy-field cyanobacterium Nostoc sphaeroides

Nitrogen-fixing cyanobacteria are common in paddy fields, one of the most productive wetland ecosystems. Here, we present the complete genome of Nostoc sphaeroides, a paddy-field diazotroph used for food and medicine for more than 1700 years and deciphered the transcriptional regulation during the developmental transition from hormogonia to vegetative filaments with heterocysts. The genome of N. sphaeroides consists of one circular chromosome (6.48 Mb), one of the largest ever reported megaplasmids (2.34 Mb), and seven plasmids. Multiple gene families involved in the adaption to high solar radiation and water fluctuation conditions were found expanded, while genes involved in anoxic adaptation and phosphonate utilization are located on the megaplasmid, suggesting its indispensable role in environmental adaptation. Distinct gene expression patterns were observed during the light-intensity-regulated transition from hormogonia to vegetative filaments, specifically, genes encoding proteins involved in photosynthetic light reaction, carbon fixation, nitrogen metabolism and heterocyst differentiation were significantly upregulated, whereas genes related to cell motility were down-regulated. Our results provide genomic and transcriptomic insights into the adaptation of a filamentous nitrogen-fixing cyanobacterium to the highly dynamic paddy-field habitat, suggesting N. sphaeroides as an excellent system to understand the transition from aquatic to terrestrial habitats and to support sustainable rice production.

1H, 13C and 15N resonance assignment of the YTH domain of YTHDC2

In humans, YTH (YT521-B homology) domain containing protein 2 (YTHDC2) plays a crucial role in the phase-shift from mitosis to meiosis. YTH domains bind to methylated adenosine nucleotides such as m⁶A. In a phylogenic tree, the YTH domain of YTHDC2 (YTH2) and that of the YTH containing protein YTHDC1 (YTH1) belong to the same sub-group. However, the binding affinity of m⁶A differs between these proteins. Here, we report ¹H, ¹³C and ¹⁵N resonance assignment of YTH2 and its solution structure to examine the difference of the structural architecture and the dynamic properties of YTH1 and YTH2. YTH2 adopts a β1-α1-β2-α2-β3-β4-β5-α3-β6-α4 topology, which was also observed in YTH1. However, the β4-β5 loops of YTH1 and YTH2 are distinct in length and amino acid composition. Our data revealed that, unlike in YTH1, the structure of m⁶A-binding pocket of YTH2 formed by the β4-β5 loop is stabilized by electrostatic interaction. This assignment and the structural information for YTH2 will provide the insight on the further functional research of YTHDC2.

Genetic variants of small airways and interstitial pulmonary disease in children

Genetic variants of small airways and interstitial pulmonary disease have not been comprehensively studied. This cluster of respiratory disorders usually manifests from early infancy (‘lung disease in utero’). In this study, 24 variants linked to these entities are described. The variants involved two genes associated with surfactant metabolism dysfunction (ABCA3 and CSF2RB), two with pulmonary fibrosis (MUC5B and SFTP), one with bronchiectasis (SCNN1B), and one with alpha-1-antitrypsin deficiency (SERPINA1). A nonsense variant, MUC5B:c.16861G > T, p.Glu5621*, was found in homozygous state in two siblings with severe respiratory disease from birth. One of the siblings also had heterozygous SFTPA1:c.675C > G, p.Asn225Lys, which resulted in a more severe respiratory disease. The sibling with only the homozygous MUC5B variant had lung biopsy, which showed alveolar simplification, interstitial fibrosis, intra-alveolar lipid-laden macrophages, and foci of foreign body giant cell reaction in distal airspaces. Two missense variants, MUC5B:c.14936 T > C, p.Ile4979Thr (rs201287218) and MUC5B:c.16738G > A, p.Gly5580Arg (rs776709402), were also found in compound heterozygous state in two siblings with severe respiratory disease from birth. Overall, the results emphasize the need for genetic studies for patients with complex respiratory problems. Identifying pathogenic variants, such as those presented here, assists in effective family counseling aimed at genetic prevention. In addition, results of genetic studies improve the clinical care and provide opportunities for participating in clinical trials, such as those involving molecularly-targeted therapies.

Donkey genomes provide new insights into domestication and selection for coat color

Current knowledge about the evolutionary history of donkeys is still incomplete due to the lack of archeological and whole-genome diversity data. To fill this gap, we have de novo assembled a chromosome-level reference genome of one male Dezhou donkey and analyzed the genomes of 126 domestic donkeys and seven wild asses. Population genomics analyses indicate that donkeys were domesticated in Africa and conclusively show reduced levels of Y chromosome variability and discordant paternal and maternal histories, possibly reflecting the consequences of reproductive management. We also investigate the genetic basis of coat color. While wild asses show diluted gray pigmentation (Dun phenotype), domestic donkeys display non-diluted black or chestnut coat colors (non-Dun) that were probably established during domestication. Here, we show that the non-Dun phenotype is caused by a 1 bp deletion downstream of the TBX3 gene, which decreases the expression of this gene and its inhibitory effect on pigment deposition.

A new donkey reference genome and comparisons with wild asses yields insights into the evolutionary history of donkey domestication and identifies a genetic variant that results in the non-Dun coat colours of domestic donkeys.

Engineering and application of a biosensor with focused ligand specificity

Cell factories converting bio-based precursors to chemicals present an attractive avenue to a sustainable economy, yet screening of genetically diverse strain libraries to identify the best-performing whole-cell biocatalysts is a low-throughput endeavor. For this reason, transcriptional biosensors attract attention as they allow the screening of vast libraries when used in combination with fluorescence-activated cell sorting (FACS). However, broad ligand specificity of transcriptional regulators (TRs) often prohibits the development of such ultra-high-throughput screens. Here, we solve the structure of the TR LysG of Corynebacterium glutamicum, which detects all three basic amino acids. Based on this information, we follow a semi-rational engineering approach using a FACS-based screening/counterscreening strategy to generate an l-lysine insensitive LysG-based biosensor. This biosensor can be used to isolate l-histidine-producing strains by FACS, showing that TR engineering towards a more focused ligand spectrum can expand the scope of application of such metabolite sensors.

Transcriptional biosensors represent powerful tools for the screening of vast strain libraries, but the broad ligand specificity of some transcriptional regulators (TRs) can prohibit such applications. Here authors present the engineering of a LysG-based biosensor with a focused ligand specificity to isolate L-histidine-producing strains.

Genome sequence and comparative analysis of reindeer (Rangifer tarandus) in northern Eurasia

Reindeer are semi-domesticated ruminants that have adapted to the challenging northern Eurasian environment characterized by long winters and marked annual fluctuations in daylight. We explored the genetic makeup behind their unique characteristics by de novo sequencing the genome of a male reindeer and conducted gene family analyses with nine other mammalian species. We performed a population genomics study of 23 additional reindeer representing both domestic and wild populations and several ecotypes from various geographic locations. We assembled 2.66 Gb (N50 scaffold of 5 Mb) of the estimated 2.92 Gb reindeer genome, comprising 27,332 genes. The results from the demographic history analysis suggested marked changes in the effective population size of reindeer during the Pleistocene period. We detected 160 reindeer-specific and expanded genes, of which zinc finger proteins (n = 42) and olfactory receptors (n = 13) were the most abundant. Comparative genome analyses revealed several genes that may have promoted the adaptation of reindeer, such as those involved in recombination and speciation (PRDM9), vitamin D metabolism (TRPV5, TRPV6), retinal development (PRDM1, OPN4B), circadian rhythm (GRIA1), immunity (CXCR1, CXCR2, CXCR4, IFNW1), tolerance to cold-triggered pain (SCN11A) and antler development (SILT2). The majority of these characteristic reindeer genes have been reported for the first time here. Moreover, our population genomics analysis suggested at least two independent reindeer domestication events with genetic lineages originating from different refugial regions after the Last Glacial Maximum. Taken together, our study has provided new insights into the domestication, evolution and adaptation of reindeer and has promoted novel genomic research of reindeer.

The Grass Carp Genomic Visualization Database (GCGVD): an informational platform for genome biology of grass carp

With the release of the draft genome of the grass carp, researches on the grass carp from the genetic level and the further molecular mechanisms of economically valuable physiological behaviors have gained great attention. In this paper, we integrated a large number of genomic, genetic and some other data resources and established a web-based grass carp genomic visualization database (GCGVD). To view these data more effectively, we visualized grass carp and zebrafish gene collinearity and genetic linkage map using Scalable Vector Graphics (SVG) format in the browser, and genomic annotations by JBrowse. Furthermore, we carried out some preliminary study on a whole-genome alternative splicing (AS)of the grass carp. The RNA-seq reads of 15 samples were aligned to the reference genome of the grass carp by Bowtie2 software. RNA-seq reads of each sample and density map of reads were also exhibited in JBrowse. Additionally, we designed a universal grass carp genome annotation data model to improve the retrieval speed and scalability. Compared with the published database GCGD previously, we newly added the visualization of some more genomic annotations, conserved domain and RNA-seq reads aligned to the reference genome. GCGVD can be accessed at http://122.112.216.104.

Expansion and Functional Diversification of SKP1-Like Genes in Wheat (Triticum aestivum L.)

The ubiquitin proteasome 26S system (UPS), involving monomeric and multimeric E3 ligases is one of the most important signaling pathways in many organisms, including plants. The SCF (SKP1/Cullin/F-box) multimeric complex is particularly involved in response to development and stress signaling. The SKP1 protein (S-phase kinase-associated protein 1) is the core subunit of this complex. In this work, we firstly identified 92 and 87 non-redundant Triticum aestivumSKP1-like (TaSKP) genes that were retrieved from the latest release of the wheat genome database (International Wheat Genome Sequencing Consortium (IWGSC) RefSeq v1.0) and the genome annotation of the TGAC v1 respectively. We then investigated the structure, phylogeny, duplication events and expression patterns of the SKP1-like gene family in various tissues and environmental conditions using a wheat expression platform containing public data. TaSKP1-like genes were expressed differentially in response to stress conditions, displaying large genomic variations or short insertions/deletions which suggests functional specialization within TaSKP1-like genes. Finally, interactions between selected wheat FBX (F-box) proteins and putative ancestral TaSKP1-like proteins were tested using the yeast two-hybrid (Y2H) system to examine the molecular interactions. These observations suggested that six Ta-SKP1 genes are likely to be ancestral genes, having similar functions as ASK1 and ASK2 in Arabidopsis, OSK1 and OSK20 in rice and PpSKP1 and PpSKP2 in Physcomitrella patens.

Identifying differentially expressed proteins in sorghum cell cultures exposed to osmotic stress

Drought stress triggers remarkable physiological changes and growth impediments, which significantly diminish plant biomass and crop yield. However, certain plant species show notable resilience, maintaining nearly normal yields under severe water deficits. For example, sorghum is a naturally drought-tolerant crop, which is ideal for studying plant adaptive responses to drought. Here we used sorbitol treatments to simulate drought-induced osmotic stress in sorghum cell suspension cultures and analysed fractions enriched for extracellular matrix proteins using isobaric tags for relative and absolute quantification technology. Sorbitol induced an overall increase in protein secretion, with putative redox proteins, proteases, and glycosyl hydrolases featuring prominently among the responsive proteins. Gene expression analysis of selected candidates revealed regulation at the transcriptional level. There was a notable differential gene expression between drought-tolerant and drought-sensitive sorghum varieties for some of the candidates. This study shows that protein secretion is a major component of the sorghum response to osmotic stress. Additionally, our data provide candidate genes, which may have putative functions in sorghum drought tolerance, and offer a pool of genes that could be developed as potential biomarkers for rapid identification of drought tolerant lines in plant breeding programs.

Strategies for high-altitude adaptation revealed from high-quality draft genome of non-violacein producing Janthinobacterium lividum ERGS5:01

A light pink coloured bacterial strain ERGS5:01 isolated from glacial stream water of Sikkim Himalaya was affiliated to Janthinobacterium lividum based on 16S rRNA gene sequence identity and phylogenetic clustering. Whole genome sequencing was performed for the strain to confirm its taxonomy as it lacked the typical violet pigmentation of the genus and also to decipher its survival strategy at the aquatic ecosystem of high elevation. The PacBio RSII sequencing generated genome of 5,168,928 bp with 4575 protein-coding genes and 118 RNA genes. Whole genome-based multilocus sequence analysis clustering, in silico DDH similarity value of 95.1% and, the ANI value of 99.25% established the identity of the strain ERGS5:01 (MCC 2953) as a non-violacein producing J. lividum. The genome comparisons across genus Janthinobacterium revealed an open pan-genome with the scope of the addition of new orthologous cluster to complete the genomic inventory. The genomic insight provided the genetic basis of freezing and frequent freeze-thaw cycle tolerance and, for industrially important enzymes. Extended insight into the genome provided clues of crucial genes associated with adaptation in the harsh aquatic ecosystem of high altitude.

Predicting Alpha Helical Transmembrane Proteins Using HMMs

Alpha helical transmembrane (TM) proteins constitute an important structural class of membrane proteins involved in a wide variety of cellular functions. The prediction of their transmembrane topology, as well as their discrimination in newly sequenced genomes, is of great importance for the elucidation of their structure and function. Several methods have been applied for the prediction of the transmembrane segments and the topology of alpha helical transmembrane proteins utilizing different algorithmic techniques. Hidden Markov Models (HMMs) have been efficiently used in the development of several computational methods used for this task. In this chapter we give a brief review of different available prediction methods for alpha helical transmembrane proteins pointing out sequence and structural features that should be incorporated in a prediction method. We then describe the procedure of the design and development of a Hidden Markov Model capable of predicting the transmembrane alpha helices in proteins and discriminating them from globular proteins.

The genome sequence of sweet cherry (Prunus avium) for use in genomics-assisted breeding

We determined the genome sequence of sweet cherry (Prunus avium) using next-generation sequencing technology. The total length of the assembled sequences was 272.4 Mb, consisting of 10,148 scaffold sequences with an N50 length of 219.6 kb. The sequences covered 77.8% of the 352.9 Mb sweet cherry genome, as estimated by k-mer analysis, and included >96.0% of the core eukaryotic genes. We predicted 43,349 complete and partial protein-encoding genes. A high-density consensus map with 2,382 loci was constructed using double-digest restriction site–associated DNA sequencing. Comparing the genetic maps of sweet cherry and peach revealed high synteny between the two genomes; thus the scaffolds were integrated into pseudomolecules using map- and synteny-based strategies. Whole-genome resequencing of six modern cultivars found 1,016,866 SNPs and 162,402 insertions/deletions, out of which 0.7% were deleterious. The sequence variants, as well as simple sequence repeats, can be used as DNA markers. The genomic information helps us to identify agronomically important genes and will accelerate genetic studies and breeding programs for sweet cherries. Further information on the genomic sequences and DNA markers is available in DBcherry (http://cherry.kazusa.or.jp (8 May 2017, date last accessed)).

Exoproteome Analysis of the Seaweed Pathogen Nautella italica R11 Reveals Temperature-Dependent Regulation of RTX-Like Proteins

Climate fluctuations have been linked to an increased prevalence of disease in seaweeds, including the red alga Delisea pulchra, which is susceptible to a bleaching disease caused by the bacterium Nautella italica R11 under elevated seawater temperatures. To further investigate the role of temperature in the induction of disease by N. italica R11, we assessed the effect of temperature on the expression of the extracellular proteome (exoproteome) in this bacterium. Label-free quantitative mass spectrometry was used to identify 207 proteins secreted into supernatant fraction, which is equivalent to 5% of the protein coding genes in the N. italica R11 genome. Comparative analysis demonstrated that expression of over 30% of the N. italica R11 exoproteome is affected by temperature. The temperature-dependent proteins include traits that could facilitate the ATP-dependent transport of amino acid and carbohydrate, as well as several uncharacterized proteins. Further, potential virulence determinants, including two RTX-like proteins, exhibited significantly higher expression in the exoproteome at the disease inducing temperature of 24°C relative to non-inducing temperature (16°C). This is the first study to demonstrate that temperature has an influence exoproteome expression in a macroalgal pathogen. The results have revealed several temperature regulated candidate virulence factors that may have a role in macroalgal colonization and invasion at elevated sea-surface temperatures, including novel RTX-like proteins.

Melanisation of Aspergillus terreus-Is Butyrolactone I Involved in the Regulation of Both DOPA and DHN Types of Pigments in Submerged Culture?

Pigments and melanins of fungal spores have been investigated for decades, revealing important roles in the survival of the fungus in hostile environments. The key genes and the encoded enzymes for pigment and melanin biosynthesis have recently been found in Ascomycota, including Aspergillus spp. In Aspergillus terreus, the pigmentation has remained mysterious with only one class of melanin biogenesis being found. In this study, we examined an intriguing, partially annotated gene cluster of A. terreus strain NIH2624, utilizing previously sequenced transcriptome and improved gene expression data of strain MUCL 38669, under the influence of a suggested quorum sensing inducing metabolite, butyrolactone I. The core polyketide synthase (PKS) gene of the cluster was predicted to be significantly longer on the basis of the obtained transcriptional data, and the surrounding cluster was positively regulated by butyrolactone I at the late growth phase of submerged culture, presumably during sporulation. Phylogenetic analysis of the extended PKS revealed remarkable similarity with a group of known pigments of Fusarium spp., indicating a similar function for this PKS. We present a hypothesis of this PKS cluster to biosynthesise a 1,8-dihydroxynaphthalene (DHN)-type of pigment during sporulation with the influence of butyrolactone I under submerged culture.

Transcriptomic Complexity of Aspergillus terreus Velvet Gene Family under the Influence of Butyrolactone I

Filamentous fungi of the Ascomycota phylum are known to contain a family of conserved conidiation regulating proteins with distinctive velvet domains. In Aspergilli, this velvet family includes four proteins, VeA, VelB, VelC and VosA, and is involved in conidiation and secondary metabolism along with a global regulator LaeA. In A. terreus, the overexpression of LaeA has been observed to increase the biogenesis of the pharmaceutically-important secondary metabolite, lovastatin, while the role of the velvet family has not been studied. The secondary metabolism and conidiation of A. terreus have also been observed to be increased by butyrolactone I in a quorum-sensing manner. An enlightenment of the interplay of these regulators will give potential advancement to the industrial use of this fungus, as well as in resolving the pathogenic features. In this study, the Aspergillus terreus MUCL 38669 transcriptome was strand-specifically sequenced to enable an in-depth gene expression analysis to further investigate the transcriptional role of butyrolactone I in these processes. The sequenced transcriptome revealed intriguing properties of the velvet family transcripts, including the regulator laeA, and uncovered the velC gene in A. terreus. The reliability refining microarray gene expression analysis disclosed a positive regulatory role for butyrolactone I in laeA expression, as well as an influence on the expression of the canonical conidiation-regulating genes under submerged culture. All of this supports the suggested regulative role of butyrolactone I in A. terreus secondary metabolism, as well as conidiation.

A three-way approach for protein function classification

The knowledge of protein functions plays an essential role in understanding biological cells and has a significant impact on human life in areas such as personalized medicine, better crops and improved therapeutic interventions. Due to expense and inherent difficulty of biological experiments, intelligent methods are generally relied upon for automatic assignment of functions to proteins. The technological advancements in the field of biology are improving our understanding of biological processes and are regularly resulting in new features and characteristics that better describe the role of proteins. It is inevitable to neglect and overlook these anticipated features in designing more effective classification techniques. A key issue in this context, that is not being sufficiently addressed, is how to build effective classification models and approaches for protein function prediction by incorporating and taking advantage from the ever evolving biological information. In this article, we propose a three-way decision making approach which provides provisions for seeking and incorporating future information. We considered probabilistic rough sets based models such as Game-Theoretic Rough Sets (GTRS) and Information-Theoretic Rough Sets (ITRS) for inducing three-way decisions. An architecture of protein functions classification with probabilistic rough sets based three-way decisions is proposed and explained. Experiments are carried out on Saccharomyces cerevisiae species dataset obtained from Uniprot database with the corresponding functional classes extracted from the Gene Ontology (GO) database. The results indicate that as the level of biological information increases, the number of deferred cases are reduced while maintaining similar level of accuracy.

Epigenetic map and genetic map basis of complex traits in cassava population

Cassava (Manihot esculenta Crantz) is an important tropical starchy root crop that is adapted to drought but extremely cold sensitive. A cold-tolerant, high-quality, and robust supply of cassava is urgently needed. Here, we clarify genome-wide distribution and classification of CCGG hemi-methylation and full-methylation, and detected 77 much candidate QTLs^epi for cold stress and 103 much candidate QTLs^epi for storage root quality and yield in 186 cassava population, generated by crossing two non-inbred lines with female parent KU50 and male parent SC124 (KS population). We developed amplified-fragment single nucleotide polymorphism and methylation (AFSM) genetic map in this population. We also constructed the AFSM QTL map, identified 260 much candidate QTL genes for cold stress and 301 much candidate QTL genes for storage root quality and yield, based on the years greenhouse and field trials. This may accounted for a significant amount of the variation in the key traits controlling cold tolerance and the high quality and yield of cassava.

Role of LRRTMs in synapse development and plasticity

Leucine-rich-repeat transmembrane neuronal proteins (LRRTMs) are a family of four synapse organizing proteins critical for the development and function of excitatory synapses. The genes encoding LRRTMs and their binding partners, neurexins and HSPGs, are strongly associated with multiple psychiatric disorders. Here, we review the literature covering their structural features, expression patterns in the developing and adult brains, evolutionary origins, and discovery as synaptogenic proteins. We also discuss their role in the development and plasticity of excitatory synapses as well as their disease associations.

Retention and Molecular Evolution of Lipoxygenase Genes in Modern Rosid Plants

Whole-genome duplication events have occurred more than once in the genomes of some rosids and played a significant role over evolutionary time. Lipoxygenases (LOXs) are involved in many developmental and resistance processes in plants. Our study concerns the subject of the LOX gene family; we tracked the evolutionary process of ancestral LOX genes in four modern rosids. Here we show that some members of the LOX gene family in the Arabidopsis genome are likely to be lost during evolution, leading to a smaller size than that in Populus, Vitis, and Carica. Strong purifying selection acted as a critical role in almost all of the paralogous and orthologous genes. The structure of LOX genes in Carica and Populus are relatively stable, whereas Vitis and Arabidopsis have a difference. By searching conserved motifs of LOX genes, we found that each sub-family shared similar components. Research on intraspecies gene collinearity show that recent duplication holds an important position in Populus and Arabidopsis. Gene collinearity analysis within and between these four rosid plants revealed that all LOX genes in each modern rosid were the offspring from different ancestral genes. This study traces the evolution of LOX genes which have been differentially retained and expanded in rosid plants. Our results presented here may aid in the selection of special genes retained in the rosid plants for further analysis of biological function.

Draft genome sequence of subterranean clover, a reference for genus Trifolium

Clovers (genus Trifolium) are widely cultivated across the world as forage legumes and make a large contribution to livestock feed production and soil improvement. Subterranean clover (T. subterraneum L.) is well suited for genomic and genetic studies as a reference species in the Trifolium genus, because it is an annual with a simple genome structure (autogamous and diploid), unlike the other economically important perennial forage clovers, red clover (T. pratense) and white clover (T. repens). This report represents the first draft genome sequence of subterranean clover. The 471.8 Mb assembled sequence covers 85.4% of the subterranean clover genome and contains 42,706 genes. Eight pseudomolecules of 401.1 Mb in length were constructed, based on a linkage map consisting of 35,341 SNPs. The comparative genomic analysis revealed that different clover chromosomes showed different degrees of conservation with other Papilionoideae species. These results provide a reference for genetic and genomic analyses in the genus Trifolium and new insights into evolutionary divergence in Papilionoideae species.

A Factor Graph Approach to Automated GO Annotation

As volume of genomic data grows, computational methods become essential for providing a first glimpse onto gene annotations. Automated Gene Ontology (GO) annotation methods based on hierarchical ensemble classification techniques are particularly interesting when interpretability of annotation results is a main concern. In these methods, raw GO-term predictions computed by base binary classifiers are leveraged by checking the consistency of predefined GO relationships. Both formal leveraging strategies, with main focus on annotation precision, and heuristic alternatives, with main focus on scalability issues, have been described in literature. In this contribution, a factor graph approach to the hierarchical ensemble formulation of the automated GO annotation problem is presented. In this formal framework, a core factor graph is first built based on the GO structure and then enriched to take into account the noisy nature of GO-term predictions. Hence, starting from raw GO-term predictions, an iterative message passing algorithm between nodes of the factor graph is used to compute marginal probabilities of target GO-terms. Evaluations on Saccharomyces cerevisiae, Arabidopsis thaliana and Drosophila melanogaster protein sequences from the GO Molecular Function domain showed significant improvements over competing approaches, even when protein sequences were naively characterized by their physicochemical and secondary structure properties or when loose noisy annotation datasets were considered. Based on these promising results and using Arabidopsis thaliana annotation data, we extend our approach to the identification of most promising molecular function annotations for a set of proteins of unknown function in Solanum lycopersicum.

Complete genome sequencing of Dehalococcoides sp. strain UCH007 using a differential reads picking method

A novel Dehalococcoides sp. strain UCH007 was isolated from the groundwater polluted with chlorinated ethenes in Japan. This strain is capable of dechlorinating trichloroethene, cis-1,2-dichloroethene and vinyl chloride to ethene. Dehalococcoides bacteria are hardly cultivable, so genome sequencing has presented a challenge. In this study, we developed a differential reads picking method for mixed genomic DNA obtained from a co-culture, and applied it to the sequencing of strain UCH007. The genome of strain UCH007 consists of a 1,473,548-bp chromosome that encodes 1509 coding sequences including 29 putative reductive dehalogenase genes. Strain UCH007 is the first strain in the Victoria subgroup found to possess the pceA, tceA and vcrA genes.

Solute Carrier Family 26 Member a2 (slc26a2) Regulates Otic Development and Hair Cell Survival in Zebrafish

Hearing loss is one of the most prevalent human birth defects. Genetic factors contribute to the pathogenesis of deafness. It is estimated that one-third of deafness genes have already been identified. The current work is an attempt to find novel genes relevant to hearing loss using guilt-by-profiling and guilt-by-association bioinformatics analyses of approximately 80 known non-syndromic hereditary hearing loss (NSHL) genes. Among the 300 newly identified candidate deafness genes, slc26a2 were selected for functional studies in zebrafish. The slc26a2 gene was knocked down using an antisense morpholino (MO), and significant defects were observed in otolith patterns, semicircular canal morphology, and lateral neuromast distributions in morphants. Loss-of-function defects are caused primarily by apoptosis, and morphants are insensitive to sound stimulation and imbalanced swimming behaviours. Morphant defects were found to be partially rescued by co-injection of human SLC26A2 mRNA. All the results suggest that bioinformatics is capable of predicting new deafness genes and this showed slc26a2 is to be a critical otic gene whose dysfunction may induce hearing impairment.

toxoMine: an integrated omics data warehouse for Toxoplasma gondii systems biology research

Toxoplasma gondii (T. gondii) is an obligate intracellular parasite that must monitor for changes in the host environment and respond accordingly; however, it is still not fully known which genetic or epigenetic factors are involved in regulating virulence traits of T. gondii. There are on-going efforts to elucidate the mechanisms regulating the stage transition process via the application of high-throughput epigenomics, genomics and proteomics techniques. Given the range of experimental conditions and the typical yield from such high-throughput techniques, a new challenge arises: how to effectively collect, organize and disseminate the generated data for subsequent data analysis. Here, we describe toxoMine, which provides a powerful interface to support sophisticated integrative exploration of high-throughput experimental data and metadata, providing researchers with a more tractable means toward understanding how genetic and/or epigenetic factors play a coordinated role in determining pathogenicity of T. gondii. As a data warehouse, toxoMine allows integration of high-throughput data sets with public T. gondii data. toxoMine is also able to execute complex queries involving multiple data sets with straightforward user interaction. Furthermore, toxoMine allows users to define their own parameters during the search process that gives users near-limitless search and query capabilities. The interoperability feature also allows users to query and examine data available in other InterMine systems, which would effectively augment the search scope beyond what is available to toxoMine. toxoMine complements the major community database ToxoDB by providing a data warehouse that enables more extensive integrative studies for T. gondii. Given all these factors, we believe it will become an indispensable resource to the greater infectious disease research community.

Database URL:http://toxomine.org

The resurrection genome of Boea hygrometrica: A blueprint for survival of dehydration

"Drying without dying" is an essential trait in land plant evolution. Unraveling how a unique group of angiosperms, the Resurrection Plants, survive desiccation of their leaves and roots has been hampered by the lack of a foundational genome perspective. Here we report the ∼1,691-Mb sequenced genome of Boea hygrometrica, an important resurrection plant model. The sequence revealed evidence for two historical genome-wide duplication events, a compliment of 49,374 protein-coding genes, 29.15% of which are unique (orphan) to Boea and 20% of which (9,888) significantly respond to desiccation at the transcript level. Expansion of early light-inducible protein (ELIP) and 5S rRNA genes highlights the importance of the protection of the photosynthetic apparatus during drying and the rapid resumption of protein synthesis in the resurrection capability of Boea. Transcriptome analysis reveals extensive alternative splicing of transcripts and a focus on cellular protection strategies. The lack of desiccation tolerance-specific genome organizational features suggests the resurrection phenotype evolved mainly by an alteration in the control of dehydration response genes.

Complete genome sequence of Selenomonas ruminantium subsp. lactilytica will accelerate further understanding of the nature of the class Negativicutes

Selenomonas ruminantium subsp. lactilytica, a strictly anaerobic ruminal bacterium, possesses typical Gram-negative cell surface structure comprising cytoplasmic membrane, peptidoglycan layer and outer membrane, whereas its 16S rRNA-based taxonomy shows that the bacteria belongs to Gram-positive Firmicutes. Complete genome analysis showed that genes or gene clusters involved in Gram-negative cell structure were scattered in the S. ruminantium genome, and might provide the new insight of phylogenetic relationship between the bacterium and other bacterial species.

Survey of genome sequences in a wild sweet potato, Ipomoea trifida (H. B. K.) G. Don

Ipomoea trifida (H. B. K.) G. Don. is the most likely diploid ancestor of the hexaploid sweet potato, I. batatas (L.) Lam. To assist in analysis of the sweet potato genome, de novo whole-genome sequencing was performed with two lines of I. trifida, namely the selfed line Mx23Hm and the highly heterozygous line 0431-1, using the Illumina HiSeq platform. We classified the sequences thus obtained as either ‘core candidates’ (common to the two lines) or ‘line specific’. The total lengths of the assembled sequences of Mx23Hm (ITR_r1.0) was 513 Mb, while that of 0431-1 (ITRk_r1.0) was 712 Mb. Of the assembled sequences, 240 Mb (Mx23Hm) and 353 Mb (0431-1) were classified into core candidate sequences. A total of 62,407 (62.4 Mb) and 109,449 (87.2 Mb) putative genes were identified, respectively, in the genomes of Mx23Hm and 0431-1, of which 11,823 were derived from core sequences of Mx23Hm, while 28,831 were from the core candidate sequence of 0431-1. There were a total of 1,464,173 single-nucleotide polymorphisms and 16,682 copy number variations (CNVs) in the two assembled genomic sequences (under the condition of log₂ ratio of >1 and CNV size >1,000 bases). The results presented here are expected to contribute to the progress of genomic and genetic studies of I. trifida, as well as studies of the sweet potato and the genus Ipomoea in general.

Comparative analysis of the human and zebrafish kinomes: focus on the development of kinase inhibitors

Targeting kinases with semi-selective kinase inhibitors is one of the most successful drug development strategies of the 21^st century. Zebrafish have become an increasingly useful model for pharmaceutical development. Water-soluble compounds can be screened for zebrafish phenotypes in a high throughput format against a living vertebrate, and cell-signaling events can be imaged in transparent living fish. Despite zebrafish being a more relevant model than more distantly related systems such as the well-annotated kinome of yeast and drosophila, there is no comparative analysis of the human and zebrafish kinome. Furthermore most approved kinase inhibitors, often called 'DFG in' ATP competitive inhibitors, act on conserved active site residues in the kinase. Since the active site residues can be identified by examining the primary sequence, primary sequence identity can be a rough guide as to whether a particular inhibitor will have activity against another kinase. There is a need to evaluate the utility of zebrafish as a drug development model for active site inhibitors of kinases. Here we offer a systematic comparison of the catalytic domains of classical human kinases with the catalytic domains of all annotated zebrafish kinases. We found a high degree of identity between the catalytic domains of most human kinases and their zebrafish homologs, and we ranked 504 human kinase catalytic domains by order of similarity. We found only 23 human kinases with no easily recognizable homologous zebrafish catalytic domain. On the other hand we found 78 zebrafish kinase catalytic domains with no close human counterpart. These 'additional kinase active sites' could represent potential mediators of zebrafish toxicity that may not be relevant to human kinase inhibitors. We used two clinically approved human kinase inhibitors, one targeting a highly homologous target and one targeting a lesser homologous target, and we compared the known human kinase target structures with modeled zebrafish target structures. As expected, the homologous target had high structural identity, but even the less homologous target had high structural identity in residues contacted by the inhibitor. Overall this analysis should help guide researchers interested in studying human kinases and their inhibitors in more tractable systems.

Functional features and protein network of human sperm-egg interaction

Elucidation of the sperm-egg interaction at the molecular level is one of the unresolved problems in sexual reproduction, and understanding the molecular mechanism is crucial in solving problems in infertility and failed in vitro fertilization (IVF). Many molecular interactions in the form of protein-protein interactions (PPIs) mediate the sperm-egg membrane interaction. Due to the complexity of the problem such as difficulties in analyzing in vivo membrane PPIs, many efforts have failed to comprehensively elucidate the fusion mechanism and the molecular interactions that mediate sperm-egg membrane fusion. The main purpose of this study was to reveal possible protein interactions and associated molecular function during sperm-egg interaction using a protein interaction network approach. Different databases have been used to construct the human sperm-egg interaction network. The constructed network revealed new interactions. These included CD151 and CD9 in human oocyte that interact with CD49 in sperm, and CD49 and ITGA4 in sperm that interact with CD63 and CD81, respectively, in the oocyte. These results showed that the different integrins in sperm may be involved in human sperm-egg interaction. It was also suggested that sperm ADAM2 plays a role as a protein candidate involved in sperm-egg membrane interaction by interacting with CD9 in the oocyte. Interleukin-4 receptor activity, receptor signaling protein tyrosine kinase activity, and manganese ion transmembrane transport activity are the major molecular functions in sperm-egg interaction protein network. The disease association analysis indicated that sperm-egg interaction defects are also reflected in other disease networks such as cardiovascular, hematological, and breast cancer diseases. By analyzing the network, we identified the major molecular functions and disease association genes in sperm-egg interaction protein. Further experimental studies will be required to confirm the significance of these new computationally resolved interactions and the genetic links between sperm-egg interaction abnormalities and the associated disease.

Differential retention and expansion of the ancestral genes associated with the paleopolyploidies in modern rosid plants, as revealed by analysis of the extensins super-gene family

BACKGROUND

All modern rosids originated from a common hexapolyploid ancestor, and the genomes of some rosids have undergone one or more cycles of paleopolyploidy. After the duplication of the ancient genome, wholesale gene loss and gene subfunctionalization has occurred. Using the extensin super-gene family as an example, we tracked the differential retention and expansion of ancestral extensin genes in four modern rosids, Arabidopsis, Populus, Vitis and Carica, using several analytical methods.

RESULTS

The majority of extensin genes in each of the modern rosids were found to originate from different ancestral genes. In Arabidopsis and Populus, almost half of the extensins were paralogous duplicates within the genome of each species. By contrast, no paralogous extensins were detected in Vitis and Carica, which have only undergone the common γ-triplication event. It was noteworthy that a group of extensins containing the IPR006706 domain had actively duplicated in Arabidopsis, giving rise to a neo-extensin around every 3 million years. However, such extensins were absent from, or rare in, the other three rosids. A detailed examination revealed that this group of extensins had proliferated significantly in the genomes of a number of species in the Brassicaceae. We propose that this group of extensins might play important roles in the biology and in the evolution of the Brassicaceae. Our analyses also revealed that nearly all of the paralogous and orthologous extensin-pairs have been under strong purifying selection, leading to the strong conservation of the function of extensins duplicated from the same ancestral gene.

CONCLUSIONS

Our analyses show that extensins originating from a common ancestor have been differentially retained and expanded among four modern rosids. Our findings suggest that, if Arabidopsis is used as the model plant, we can only learn a limited amount about the functions of a particular gene family. These results also provide an example of how it is essential to learn the origination of a gene when analyzing its function across different plant species.

Differential retention and expansion of the ancestral genes associated with the paleopolyploidies in modern rosid plants, as revealed by analysis of the extensins super-gene family

Background

All modern rosids originated from a common hexapolyploid ancestor, and the genomes of some rosids have undergone one or more cycles of paleopolyploidy. After the duplication of the ancient genome, wholesale gene loss and gene subfunctionalization has occurred. Using the extensin super-gene family as an example, we tracked the differential retention and expansion of ancestral extensin genes in four modern rosids, Arabidopsis, Populus, Vitis and Carica, using several analytical methods.

Results

The majority of extensin genes in each of the modern rosids were found to originate from different ancestral genes. In Arabidopsis and Populus, almost half of the extensins were paralogous duplicates within the genome of each species. By contrast, no paralogous extensins were detected in Vitis and Carica, which have only undergone the common γ-triplication event. It was noteworthy that a group of extensins containing the IPR006706 domain had actively duplicated in Arabidopsis, giving rise to a neo-extensin around every 3 million years. However, such extensins were absent from, or rare in, the other three rosids. A detailed examination revealed that this group of extensins had proliferated significantly in the genomes of a number of species in the Brassicaceae. We propose that this group of extensins might play important roles in the biology and in the evolution of the Brassicaceae. Our analyses also revealed that nearly all of the paralogous and orthologous extensin-pairs have been under strong purifying selection, leading to the strong conservation of the function of extensins duplicated from the same ancestral gene.

Conclusions

Our analyses show that extensins originating from a common ancestor have been differentially retained and expanded among four modern rosids. Our findings suggest that, if Arabidopsis is used as the model plant, we can only learn a limited amount about the functions of a particular gene family. These results also provide an example of how it is essential to learn the origination of a gene when analyzing its function across different plant species.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-612) contains supplementary material, which is available to authorized users.

Draft sequences of the radish (Raphanus sativus L.) genome

Radish (Raphanus sativus L., n = 9) is one of the major vegetables in Asia. Since the genomes of Brassica and related species including radish underwent genome rearrangement, it is quite difficult to perform functional analysis based on the reported genomic sequence of Brassica rapa. Therefore, we performed genome sequencing of radish. Short reads of genomic sequences of 191.1 Gb were obtained by next-generation sequencing (NGS) for a radish inbred line, and 76,592 scaffolds of ≥300 bp were constructed along with the bacterial artificial chromosome-end sequences. Finally, the whole draft genomic sequence of 402 Mb spanning 75.9% of the estimated genomic size and containing 61,572 predicted genes was obtained. Subsequently, 221 single nucleotide polymorphism markers and 768 PCR-RFLP markers were used together with the 746 markers produced in our previous study for the construction of a linkage map. The map was combined further with another radish linkage map constructed mainly with expressed sequence tag-simple sequence repeat markers into a high-density integrated map of 1,166 cM with 2,553 DNA markers. A total of 1,345 scaffolds were assigned to the linkage map, spanning 116.0 Mb. Bulked PCR products amplified by 2,880 primer pairs were sequenced by NGS, and SNPs in eight inbred lines were identified.

Population genomics reveal recent speciation and rapid evolutionary adaptation in polar bears

Polar bears are uniquely adapted to life in the High Arctic and have undergone drastic physiological changes in response to Arctic climates and a hyper-lipid diet of primarily marine mammal prey. We analyzed 89 complete genomes of polar bear and brown bear using population genomic modeling and show that the species diverged only 479-343 thousand years BP. We find that genes on the polar bear lineage have been under stronger positive selection than in brown bears; nine of the top 16 genes under strong positive selection are associated with cardiomyopathy and vascular disease, implying important reorganization of the cardiovascular system. One of the genes showing the strongest evidence of selection, APOB, encodes the primary lipoprotein component of low-density lipoprotein (LDL); functional mutations in APOB may explain how polar bears are able to cope with life-long elevated LDL levels that are associated with high risk of heart disease in humans.