In order to clarify the phylogenetic relationships among the main marine myxosporean clades including newly established Ceratonova clade and scrutinizing their evolutionary origins, we performed large-scale phylogenetic analysis of all myxosporean species from the marine myxosporean lineage based on three gene analyses and statistical topology tests. Furthermore, we obtained new molecular data for Ceratonova shasta, C. gasterostea, eight Ceratomyxa species and one Myxodavisia species. We described five new species: Ceratomyxa ayami n. sp., C. leatherjacketi n. sp., C. synaphobranchi n. sp., C. verudaensis n. sp. and Myxodavisia bulani n. sp.; two of these formed a new, basal Ceratomyxa subclade. We identified that the Ceratomyxa clade is basal to all other marine myxosporean lineages, and Kudoa with Enteromyxum are the most recently branching clades. Topologies were least stable at the nodes connecting the marine urinary clade, the marine gall bladder clade and the Ceratonova clade. Bayesian inference analysis of SSU rDNA and the statistical tree topology tests suggested that Ceratonova is closely related to the Enteromyxum and Kudoa clades, which represent a large group of histozoic species. A close relationship between Ceratomyxa and Ceratonova was not supported, despite their similar myxospore morphologies. Overall, the site of sporulation in the vertebrate host is a more accurate predictor of phylogenetic relationships than the morphology of the myxospore.
Subtilisin-like serine proteases (EC 3.4.21) consist of a widespread family of enzymes that is involved in various processes including in plants. The full-length cDNA (CpSUB1) and the corresponding genomic DNA for papaya subtilase have been obtained using rapid amplification of cDNA ends (RACEs) and PCR primer walking techniques, respectively. The cDNA clone contains an open reading frame of 2316bp encoding 772 amino acids with a calculated molecular mass of 82.6kDa and an isoelectric point (pI) of 8.97. The CpSUB1 gene is composed of nine exons and eight introns. The amino acid sequence encoded by CpSUB1 shared high identity (>60%) with the amino acid sequence of other plant subtilisin-like proteases. Sequence analysis of CpSUB1 revealed the presence of a possible signal peptide (25 amino acid residues) and an NH(2)-terminal prosequence (88 amino acid residues). In addition, papaya subtilase possesses the characteristic subtilisin catalytic triad amino acids namely Asp, His and Ser, together with the substrate-binding site, Asn. DNA hybridization analysis showed that subtilase gene exists as a single copy in the papaya genome. RNA hybridization analyses showed that expression of the subtilase transcripts was only detected in mesocarp but not in non-fruit tissues. Gene expression in fruit tissues reached the highest level during the ripening stage at which the fruits undergo dramatic softening process. Subsequently, pro-subtilase ( approximately 80kDa) was expressed as recombinant pro-enzyme ( approximately 97kDa), which was used to generate antiserum against papaya subtilase, anti-sub. Protein gel blot analysis using anti-sub towards total protein extracted from all ripening stages revealed that a protein with a molecular mass of approximately 70kDa reacted with the antiserum. Hence both RNA hybridization and protein gel blot analyses confirmed the presence of subtilase during papaya fruit ripening, pointing to its possible involvement in this important process.
Using Genetic Algorithm, this paper presents a modelling method to generate novel logical-based features from DNA sequences enriched with H3K4mel histone signatures. Current histone signature is mostly represented using k-mers content features incapable of representing all the possible complex interactions of various DNA segments. The main contributions are, among others: (a) demonstrating that there are complex interactions among sequence segments in the histone regions; (b) developing a parse tree representation of the logical complex features. The proposed novel feature is compared to the k-mers content features using datasets from the mouse (mm9) genome. Evaluation results show that the new feature improves the prediction performance as shown by f-measure for all datasets tested. Also, it is discovered that tree-based features generated from a single chromosome can be generalized to predict histone marks in other chromosomes not used in the training. These findings have a great impact on feature design considerations for histone signatures as well as other classifier design features.
Although the spliceogenic nature of the BRCA2 c.68-7T > A variant has been demonstrated, its association with cancer risk remains controversial. In this study, we accurately quantified by real-time PCR and digital PCR (dPCR), the BRCA2 isoforms retaining or missing exon 3. In addition, the combined odds ratio for causality of the variant was estimated using genetic and clinical data, and its associated cancer risk was estimated by case-control analysis in 83,636 individuals. Co-occurrence in trans with pathogenic BRCA2 variants was assessed in 5,382 families. Exon 3 exclusion rate was 4.5-fold higher in variant carriers (13%) than controls (3%), indicating an exclusion rate for the c.68-7T > A allele of approximately 20%. The posterior probability of pathogenicity was 7.44 × 10-115 . There was neither evidence for increased risk of breast cancer (OR 1.03; 95% CI 0.86-1.24) nor for a deleterious effect of the variant when co-occurring with pathogenic variants. Our data provide for the first time robust evidence of the nonpathogenicity of the BRCA2 c.68-7T > A. Genetic and quantitative transcript analyses together inform the threshold for the ratio between functional and altered BRCA2 isoforms compatible with normal cell function. These findings might be exploited to assess the relevance for cancer risk of other BRCA2 spliceogenic variants.
Bactrocera dorsalis s.s. (Hendel) and B. papayae Drew & Hancock, are invasive pests belonging to the B. dorsalis complex. Their species status, based on morphology, is sometimes arguable. Consequently, the existence of cryptic species and/or population isolation may decrease the effectiveness of the sterile insect technique (SIT) due to an unknown degree of sexual isolation between released sterile flies and wild counterparts. To evaluate the genetic relationship and current demography in wild populations for guiding the application of area-wide integrated pest management using SIT, seven microsatellite-derived markers from B. dorsalis s.s. and another five from B. papayae were used for surveying intra- and inter-specific variation, population structure, and recent migration among sympatric and allopatric populations of the two morphological forms across Southern Thailand and West Malaysia.
Oil palm is the most productive oil-bearing crop. Although it is planted on only 5% of the total world vegetable oil acreage, palm oil accounts for 33% of vegetable oil and 45% of edible oil worldwide, but increased cultivation competes with dwindling rainforest reserves. We report the 1.8-gigabase (Gb) genome sequence of the African oil palm Elaeis guineensis, the predominant source of worldwide oil production. A total of 1.535 Gb of assembled sequence and transcriptome data from 30 tissue types were used to predict at least 34,802 genes, including oil biosynthesis genes and homologues of WRINKLED1 (WRI1), and other transcriptional regulators, which are highly expressed in the kernel. We also report the draft sequence of the South American oil palm Elaeis oleifera, which has the same number of chromosomes (2n = 32) and produces fertile interspecific hybrids with E. guineensis but seems to have diverged in the New World. Segmental duplications of chromosome arms define the palaeotetraploid origin of palm trees. The oil palm sequence enables the discovery of genes for important traits as well as somaclonal epigenetic alterations that restrict the use of clones in commercial plantings, and should therefore help to achieve sustainability for biofuels and edible oils, reducing the rainforest footprint of this tropical plantation crop.
HIV-1 epidemics among men who have sex with men (MSM) continue to expand in developed and developing countries. Although HIV infection in MSM is amongst the highest of the key affected populations in many countries in Southeast Asia, comprehensive molecular epidemiological study of HIV-1 among MSM remains inadequate in the region including in Malaysia. Here, we reported the phylodynamic profiles of HIV-1 genotypes circulating among MSM population in Kuala Lumpur, Malaysia. A total of n = 459 newly-diagnosed treatment-naïve consenting subjects were recruited between March 2006 and August 2012, of whom 87 (18.9%) were self-reported MSM. Transmitted drug resistance mutations were absent in these isolates. Cumulatively, phylogenetic reconstructions of the pro-rt gene (HXB2∶2253-3275) showed that HIV-1 subtype B and CRF01_AE were predominant and contributed to approximately 80% of the total HIV-1 infection among MSM. In addition to numerous unique transmission lineages within these genotypes, twelve monophyletic transmission clusters of different sizes (2-7 MSM sequences, supported by posterior probability value of 1) were identified in Malaysia. Bayesian coalescent analysis estimated that the divergence times for these clusters were mainly dated between 1995 and 2005 with four major transmission clusters radiating at least 12 years ago suggesting that active spread of multiple sub-epidemic clusters occurred during this period. The changes in effective population size of subtype B showed an exponential growth within 5 years between 1988 and 1993, while CRF01_AE lineage exhibited similar expansion between 1993 and 2003. Our study provides the first insight of the phylodynamic profile of HIV-1 subtype B and CRF01_AE circulating among MSM population in Kuala Lumpur, Malaysia, unravelling the importance of understanding transmission behaviours as well as evolutionary history of HIV-1 in assessing the risk of outbreak or epidemic expansion.
Multivariate analyses were performed using 13 morphological traits and 13 molecular markers (10 SSRs and three ISSRs) to assess the phylogenetic relationship among tungro resistant genotypes. For morphological traits, the genotypes were grouped into six clusters, according to D(2) statistic and Canonical vector analysis. Plant height, days to flowering, days to maturity, panicle length, number of spikelet per panicle, number of unfilled grain per panicle and yield were important contributors to genetic divergence in 14 rice genotypes. Based on Nei's genetic distance for molecular studies, seven clusters were formed among the tungro resistant and susceptible genotypes. Mantel's test revealed a significant correlation (r = 0.834*) between the morphological and molecular data. To develop high yielding tungro resistant varieties based on both morphological and molecular analyses, crosses could be made with susceptible (BR10 and BR11) genotypes with low yielding but highly resistant genotypes, Sonahidemota, Kumragoir, Nakuchimota, Khaiyamota, Khairymota and Kachamota. The chi-square analysis for seven alleles (RM11, RM17, RM20, RM23, RM80, RM108 and RM531) of SSR and five loci (RY1, MR1, MR2, MR4 and GF5) of three ISSR markers in F2 population of cross, BR11×Sonahidemota, showed a good fit to the expected segregation ratio (1:2:1) for a single gene model.
Andrographis paniculata (AP) is a medicinal plant species introduced into Malaysia. To address the genetic structure and evolutionary connectedness of the Malaysian AP with the Indian AP, a DNA sequence analysis was conducted based on 24 microsatellite markers. Out of the 24 primer sets, seven novel microsatellite primers were designed and amplified intra-specifically according to the available Indian AP sequences at the National Centre for Biotechnology Information (NCBI), where 17 of them were amplified using the cross-species strategy by employing the primers belonging to Acanthus ilicifolius Linn (Acanthaceae) and Lumnitzera racemosa Wild (Combretaceae). The primers were then applied on the Malaysian AP accessions. Sixteen of the new microsatellite loci were amplified successfully. Analysis of these microsatellite sequences, revealed some significant differences between the Indian and Malaysian AP accessions in terms of the size and type of the repeat motifs. These findings depicted the cryptic feature of this species. Despite identifying several heterozygous alleles no polymorphism was observed in the detected loci of the selected accessions. This situation was in concordance with the presence of "fixed heterozygosity" phenomenon in the mentioned loci. Accordingly, this was fully consistent with the occurrence of the genetic bottleneck and founder effect within Malaysian AP population. Apart from the amplification of new microsatellites in this species, our observations could be in agreement with the risk of genetic depletion and consequently extinction of this precious herb in Malaysia. This issue should be taken into consideration in the future studies.
Understanding the mechanisms of gene regulation during breast cancer is one of the most difficult problems among oncologists because this regulation is likely comprised of complex genetic interactions. Given this complexity, a computational study using the Bayesian network technique has been employed to construct a gene regulatory network from microarray data. Although the Bayesian network has been notified as a prominent method to infer gene regulatory processes, learning the Bayesian network structure is NP hard and computationally intricate. Therefore, we propose a novel inference method based on low-order conditional independence that extends to the case of the Bayesian network to deal with a large number of genes and an insufficient sample size. This method has been evaluated and compared with full-order conditional independence and different prognostic indices on a publicly available breast cancer data set. Our results suggest that the low-order conditional independence method will be able to handle a large number of genes in a small sample size with the least mean square error. In addition, this proposed method performs significantly better than other methods, including the full-order conditional independence and the St. Gallen consensus criteria. The proposed method achieved an area under the ROC curve of 0.79203, whereas the full-order conditional independence and the St. Gallen consensus criteria obtained 0.76438 and 0.73810, respectively. Furthermore, our empirical evaluation using the low-order conditional independence method has demonstrated a promising relationship between six gene regulators and two regulated genes and will be further investigated as potential breast cancer metastasis prognostic markers.
It has recently been shown that ancestors of New Guineans and Bougainville Islanders have inherited a proportion of their ancestry from Denisovans, an archaic hominin group from Siberia. However, only a sparse sampling of populations from Southeast Asia and Oceania were analyzed. Here, we quantify Denisova admixture in 33 additional populations from Asia and Oceania. Aboriginal Australians, Near Oceanians, Polynesians, Fijians, east Indonesians, and Mamanwa (a "Negrito" group from the Philippines) have all inherited genetic material from Denisovans, but mainland East Asians, western Indonesians, Jehai (a Negrito group from Malaysia), and Onge (a Negrito group from the Andaman Islands) have not. These results indicate that Denisova gene flow occurred into the common ancestors of New Guineans, Australians, and Mamanwa but not into the ancestors of the Jehai and Onge and suggest that relatives of present-day East Asians were not in Southeast Asia when the Denisova gene flow occurred. Our finding that descendants of the earliest inhabitants of Southeast Asia do not all harbor Denisova admixture is inconsistent with a history in which the Denisova interbreeding occurred in mainland Asia and then spread over Southeast Asia, leading to all its earliest modern human inhabitants. Instead, the data can be most parsimoniously explained if the Denisova gene flow occurred in Southeast Asia itself. Thus, archaic Denisovans must have lived over an extraordinarily broad geographic and ecological range, from Siberia to tropical Asia.
Long columns of migrating larval sciarid armyworms were discovered in central and northern Japan, specifically Kanagawa, Gunma, Miyagi and Akita prefectures, as well as Hokkaido. This is the first examination of armyworms in East Asia. In Europe, armyworms have been identified as Sciara militaris, belonging to the family Sciaridae (sciarid flies or black fungus gnats), by rearing them to adulthood. In Japan, we were unable to obtain live samples for rearing; therefore, DNA barcodes were obtained from the samples of armyworms collected in the Gunma and Miyagi prefectures. The DNA barcodes were compared with those obtained from the following samples: pupae of S. militaris from UK, adults of Sciara kitakamiensis, Sciara humeralis, Sciara hemerobioides, Sciara thoracica, Sciara helvola and Sciara melanostyla from Japan, and adults of one undescribed Sciara species from Malaysia. Neighbour-joining, maximum parsimony, and maximum likelihood analyses revealed that the armyworms discovered in Japan are S. kitakamiensis. Although adults of this species have been recorded in several locations in Japan, this is the first report of migrating larval armyworms. DNA barcodes were effectively used to link different life stages of this species. The average intraspecific and interspecific pairwise genetic distances of the genus Sciara were 0.3% and 12.6%, respectively. The present study illustrates that DNA barcodes are an effective means of identifying sciarid flies in Japan.
Signal transduction through protein-protein interactions and protein modifications are the main mechanisms controlling many biological processes. Here we described the implementation of MedScan information extraction technology and Pathway Studio software (Ariadne Genomics Inc.) to create a Salmonella specific molecular interaction database. Using the database, we have constructed several signal transduction pathways in Salmonella enterica serovar Typhi which causes Typhoid Fever, a major health threat especially in developing countries. S. Typhi has several pathogenicity islands that control rapid switching between different phenotypes including adhesion and colonization, invasion, intracellular survival, proliferation, and biofilm formation in response to environmental changes. Understanding of the detailed mechanism for S. Typhi survival in host cells is necessary for development of efficient detection and treatment of this pathogen. The constructed pathways were validated using publically available gene expression microarray data for Salmonella.
Nasopharyngeal carcinoma (NPC) is one of the most common cancers in Malaysia, mainly occurring among the Chinese population. To detect common genetic alterations in NPC, we screened seven cases of NPC using the comparative genomic hybridization (CGH) technique. Before proceeding to the CGH technique, the tumors were first confirmed to consist of 75% tumor cells or more. In brief, the technique consists of binding tumor DNA with normal DNA and human Cot-1 DNA, which is then hybridized to normal metaphase spreads. The slides were then counterstained with 4,6 diamino-2-phenylindole (DAPI II) for detection. Analyses were performed using CGH software (Cytovision). We found genetic alterations in all seven NPC samples. The common chromosomal gains (57%, four cases) were found on chromosome arms 1q, 4p, 5, 7q, 11, 14p, 15q, 18p, and 21p, and common chromosomal losses (43%, three cases) were found on chromosome arm 16p. Our results showed chromosomal alterations in all seven NPC cases in the Malaysian population. This result provides the platform for further investigations to locate tumor suppressor genes and oncogenes at specific chromosomal regions in Malaysian NPC patients.
By investigating genealogical relationships, we estimated the phylogenetic history and biogeography in the megophryid genus Leptobrachium (sensu lato, including Vibrissaphora) from southern China, Indochina, Thailand and the Sundaland. The genealogical relationships among the 30 named and unnamed taxa were estimated using 2009 bp of sequences from the mitochondrial DNA genes 12S rRNA, tRNA(val), and 16S rRNA using maximum parsimony, maximum likelihood, and Bayesian inference methods. The genus Leptobrachium was a well-supported monophyletic group that contained two major clades. One clade had three subclades primarily from disjunct regions including Borneo, Peninsular Malaysia and Java, and Thailand. The Bornean subclade included one species each from the Philippines and Sumatra. The other major clade consisted of two subclades, one from Indochina and the other from southern China (Vibrissaphora). Divergence times estimated an old evolutionary history of each subclade, one that could not be explained by the geohistory of Southeast Asian major landmasses.
Several anuran groups of Laurasian origin are each co-distributed in four isolated regions of the Northern Hemisphere: central/southern Europe and adjacent areas, Korean Peninsula and adjacent areas, Indo-Malaya, and southern North America. Similar distribution patterns have been observed in diverse animal and plant groups. Savage [Savage, J.M., 1973. The geographic distribution of frogs: patterns and predictions. In: Vial, J.L. (Ed.), Evolutionary Biology of the Anurans. University of Missouri Press, Columbia, pp. 351-445] hypothesized that the Miocene global cooling and increasing aridities in interiors of Eurasia and North America caused a southward displacement and range contraction of Laurasian frogs (and other groups). We use the frog genus Bombina to test Savage's biogeographical hypothesis. A phylogeny of Bombina is reconstructed based on three mitochondrial and two nuclear gene fragments. The genus is divided into three major clades: an Indo-Malaya clade includes B. fortinuptialis, B. lichuanensis, B. maxima, and B. microdeladigitora; a European clade includes B. bombina, B. pachypus, and B. variegata; and a Korean clade contains B. orientalis. The European and Korean clades form sister-group relationship. Molecular dating of the phylogenetic tree using the penalized likelihood and Bayesian analyses suggests that the divergence between the Indo-Malaya clade and other Bombina species occurred 5.9-28.6 million years ago. The split time between the European clade and the Korean clade is estimated at 5.1-20.9 million years ago. The divergence times of these clades are not significantly later than the timing of Miocene cooling and drying, and therefore can not reject Savage's hypothesis. Some other aspects of biogeography of Bombina also are discussed. The Korean Peninsula and the Shandong Peninsula might have supplied distinct southern refugia for B. orientalis during the Pleistocene glacial maxima. In the Indo-Malaya clade, the uplift of the Tibetan Plateau might have promoted the split between B. maxima and the other species.
The extent to which response to environmental change is mediated by species-specific ecology is an important aspect of the population histories of tropical taxa. During the Pleistocene glacial cycles and associated sea level fluctuations, the Sunda region in Southeast Asia experienced concurrent changes in landmass area and the ratio of forest to open habitat, providing an ideal setting to test the expectation that habitat associations played an important role in determining species' response to the opportunity for geographic expansion. We used mitochondrial control region sequences and six microsatellite loci to compare the phylogeographic structure and demographic histories of four broadly sympatric species of Old World fruit bats in the genus, Cynopterus. Two forest-associated species and two open-habitat generalists were sampled along a latitudinal transect in Singapore, peninsular Malaysia, and southern Thailand. Contrary to expectations based on habitat associations, the geographic scale of population structure was not concordant across ecologically similar species. We found evidence for long and relatively stable demographic history in one forest and one open-habitat species, and inferred non-coincident demographic expansions in the second forest and open-habitat species. Thus, while these results indicate that Pleistocene climate change did not have a single effect on population structure across species, a correlation between habitat association and response to environmental change was supported in only two of four species. We conclude that interactions between multiple factors, including historical and contemporary environmental change, species-specific ecology and interspecific interactions, have shaped the recent evolutionary histories of Cynopterus fruit bats in Southeast Asia.
Directed evolution of nucleotide libraries using recombination or mutagenesis is an important technique for customizing catalytic or biophysical traits of proteins. Conventional directed evolution methods, however, suffer from cumbersome digestion and ligation steps. Here, we describe a simple method to increase nucleotide diversity using single-stranded DNA (ssDNA) as a starting template. An initial PCR amplification using phosphorylated primers with overlapping regions followed by treatment with lambda exonuclease generates ssDNA templates that can then be annealed via the overlap regions. Double-stranded DNA (dsDNA) is then generated through extension with Klenow fragment. To demonstrate the applicability of this methodology for directed evolution of nucleotide libraries, we generated both gene shuffled and regional mutagenesis synthetic antibody libraries with titers of 2×108 and 6×107, respectively. We conclude that our method is an efficient and convenient approach to generate diversity in nucleic acid based libraries, especially recombinant antibody libraries.