A sum of 48 accessions of physic nut, Jatropha curcas L. were analyzed to determine the genetic diversity and association between geographical origin using RAPD-PCR markers. Eight primers generated a total of 92 fragments with an average of 11.5 amplicons per primer. Polymorphism percentages of J. curcas accessions for Selangor, Kelantan, and Terengganu states were 80.4, 50.0, and 58.7%, respectively, with an average of 63.04%. Jaccard's genetic similarity co-efficient indicated the high level of genetic variation among the accessions which ranged between 0.06 and 0.81. According to UPGMA dendrogram, 48 J. curcas accessions were grouped into four major clusters at coefficient level 0.3 and accessions from same and near states or regions were found to be grouped together according to their geographical origin. Coefficient of genetic differentiation (G(st)) value of J. curcas revealed that it is an outcrossing species.
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.
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.
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.
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.
Type 2 diabetes mellitus (T2DM) is associated with a high incidence of nephropathy. The aim of this study was to investigate the association of a genetic polymorphism of carnosinase (CNDP1-D18S880 and -rs2346061), endothelial nitric oxide synthase (NOS3-rs1799983), and manganese superoxide dismutase (MnSOD-rs4880) genes with the development of diabetic nephropathy among Malaysian type 2 diabetic patients. A case-control association study was performed using 652 T2DM patients comprising 227 Malays (without nephropathy = 96 and nephropathy = 131), 203 Chinese (without nephropathy = 95 and nephropathy = 108), and 222 Indians (without nephropathy = 136 and nephropathy = 86). DNA sequencing was performed for the D18S880 of CNDP1, while the rest were tested using DNA Sequenom MassARRAY to identify the polymorphisms. DNA was extracted from the secondary blood samples taken from the T2DM patients. The alleles and genotypes were tested using four genetic models, and the best mode of inheritance was chosen based on the least p value. The rs2346061 of CNDP1 was significantly associated with diabetic nephropathy among the Indians only with OR = 1.94 and 95% CI = (1.76-3.20) and fitted best the multiplicative model, while D18S880 was associated among all the three major races with the Malays having the strongest association with OR = 2.46 and 95% CI = (1.48-4.10), Chinese with OR = 2.26 and 95% CI = (1.34-3.83), and Indians with OR = 1.77 and 95% CI = (1.18-2.65) in the genotypic multiplicative model. The best mode of inheritance for both MnSOD and NOS3 was the additive model. For MnSOD-rs4880, the Chinese had OR = 2.8 and 95% CI = (0.53-14.94), Indians had OR = 2.4 and 95% CI = (0.69-2.84), and Malays had OR = 2.16 and 95% CI = (0.54-8.65), while for NOS3-rs1799983, the Indians had the highest risk with OR = 3.16 and 95% CI = (0.52-17.56), followed by the Chinese with OR = 3.55 and 95% CI = (0.36-35.03) and the Malays with OR = 2.89 and 95% CI = (0.29-28.32). The four oxidative stress-related polymorphisms have significant effects on the development of nephropathy in type 2 diabetes patients. The genes may, therefore, be considered as risk factors for Malaysian subjects who are predisposed to T2DM nephropathy.
Disease resistance (R) genes from wild relatives could be used to engineer broad-spectrum resistance in domesticated crops. We combined association genetics with R gene enrichment sequencing (AgRenSeq) to exploit pan-genome variation in wild diploid wheat and rapidly clone four stem rust resistance genes. AgRenSeq enables R gene cloning in any crop that has a diverse germplasm panel.
Most B cell precursor acute lymphoblastic leukemia (BCP ALL) can be classified into known major genetic subtypes, while a substantial proportion of BCP ALL remains poorly characterized in relation to its underlying genomic abnormalities. We therefore initiated a large-scale international study to reanalyze and delineate the transcriptome landscape of 1,223 BCP ALL cases using RNA sequencing. Fourteen BCP ALL gene expression subgroups (G1 to G14) were identified. Apart from extending eight previously described subgroups (G1 to G8 associated with MEF2D fusions, TCF3-PBX1 fusions, ETV6-RUNX1-positive/ETV6-RUNX1-like, DUX4 fusions, ZNF384 fusions, BCR-ABL1/Ph-like, high hyperdiploidy, and KMT2A fusions), we defined six additional gene expression subgroups: G9 was associated with both PAX5 and CRLF2 fusions; G10 and G11 with mutations in PAX5 (p.P80R) and IKZF1 (p.N159Y), respectively; G12 with IGH-CEBPE fusion and mutations in ZEB2 (p.H1038R); and G13 and G14 with TCF3/4-HLF and NUTM1 fusions, respectively. In pediatric BCP ALL, subgroups G2 to G5 and G7 (51 to 65/67 chromosomes) were associated with low-risk, G7 (with ≤50 chromosomes) and G9 were intermediate-risk, whereas G1, G6, and G8 were defined as high-risk subgroups. In adult BCP ALL, G1, G2, G6, and G8 were associated with high risk, while G4, G5, and G7 had relatively favorable outcomes. This large-scale transcriptome sequence analysis of BCP ALL revealed distinct molecular subgroups that reflect discrete pathways of BCP ALL, informing disease classification and prognostic stratification. The combined results strongly advocate that RNA sequencing be introduced into the clinical diagnostic workup of BCP ALL.
DNA methylation is instrumental for gene regulation. Global changes in the epigenetic landscape have been recognized as a hallmark of cancer. However, the role of DNA methylation in epithelial ovarian cancer (EOC) remains unclear. In this study, high-density genetic and DNA methylation data in white blood cells from the Framingham Heart Study (N = 1,595) were used to build genetic models to predict DNA methylation levels. These prediction models were then applied to the summary statistics of a genome-wide association study (GWAS) of ovarian cancer including 22,406 EOC cases and 40,941 controls to investigate genetically predicted DNA methylation levels in association with EOC risk. Among 62,938 CpG sites investigated, genetically predicted methylation levels at 89 CpG were significantly associated with EOC risk at a Bonferroni-corrected threshold of P < 7.94 × 10-7. Of them, 87 were located at GWAS-identified EOC susceptibility regions and two resided in a genomic region not previously reported to be associated with EOC risk. Integrative analyses of genetic, methylation, and gene expression data identified consistent directions of associations across 12 CpG, five genes, and EOC risk, suggesting that methylation at these 12 CpG may influence EOC risk by regulating expression of these five genes, namely MAPT, HOXB3, ABHD8, ARHGAP27, and SKAP1. We identified novel DNA methylation markers associated with EOC risk and propose that methylation at multiple CpG may affect EOC risk via regulation of gene expression. SIGNIFICANCE: Identification of novel DNA methylation markers associated with EOC risk suggests that methylation at multiple CpG may affect EOC risk through regulation of gene expression.
Fin whales (Balaenoptera physalus) and blue whales (B. musculus) are the two largest species on Earth and are widely distributed across the world's oceans. Hybrids between these species appear to be relatively widespread and have been reported in both the North Atlantic and North Pacific; they are also relatively common, and have been proposed to occur once in every thousand fin whales. However, despite known hybridization, fin and blue whales are not sibling species. Rather, the closest living relative of fin whales are humpback whales (Megaptera novaeangliae). To improve the quality of fin whale data available for analysis, we assembled and annotated a fin whale nuclear genome using in-silico mate pair libraries and previously published short-read data. Using this assembly and genomic data from a humpback, blue, and bowhead whale, we investigated whether signatures of introgression between the fin and blue whale could be found. We find no signatures of contemporary admixture in the fin and blue whale genomes, although our analyses support ancestral gene flow between the species until 2.4-1.3 Ma. We propose the following explanations for our findings; i) fin/blue whale hybridization does not occur in the populations our samples originate from, ii) contemporary hybrids are a recent phenomenon and the genetic consequences have yet to become widespread across populations, or iii) fin/blue whale hybrids are under large negative selection, preventing them from backcrossing and contributing to the parental gene pools.
The gekkonid genus Cyrtodactylus is a highly diverse group of lizards (280 + species), which covers an expansive geographic range. Although this genus has been the focus of many taxonomic and molecular systematic studies, species on the Southeast Asian island of Borneo have remained understudied, leading to an unclear evolutionary history with cascading effects on taxonomy and biogeographic inferences. We assembled the most comprehensive multilocus Bornean dataset (one mitochondrial and three nuclear loci) that included 129 novel sequences and representatives from each known Cyrtodactylus species on the island to validate taxonomic status, assess species diversity, and elucidate biogeographic patterns. Our results uncovered a high proportion of cryptic diversity and revealed numerous taxonomic complications, especially within the C. consobrinus, C. malayanus, and C. pubisulcus groups. Comparisons of pairwise genetic distances and a preliminary species delimitation analysis using the Automatic Barcode Gap Discovery (ABGD) method demonstrated that some wide-ranging species on Borneo likely comprise multiple distinct and deeply divergent lineages, each with more restricted distributional ranges. We also tested the prevailing biogeographic hypothesis of a single invasion from Borneo into the Philippines. Our analyses revealed that Philippine taxa were not monophyletic, but were likely derived from multiple separate invasions into the geopolitical areas comprising the Philippines. Although our investigation of Bornean Cyrtodactylus is the most comprehensive to-date, it highlights the need for expanded taxonomic sampling and suggests that our knowledge of the evolutionary history, systematics, and biogeography of Bornean Cyrtodactylus is far from complete.
Schizophrenics (n = 250) and normal controls (n = 90) were studied to investigate and compare their dermatoglyphic patterns. Their fingerprint patterns were studied. The frequency of arches in the patient and control groups was similar. The frequency of loops in the control group was higher than in the patient group, and the trend was consistent in all the digits. The whorls in the patient group showed an increase over the control group in all the digits, although this finding was not statistically significant.
Migration among local populations classically has been seen as the principal process retarding genetic microdifferentiation. However, as Sewall Wright pointed out long ago, migration may also act as a random differentiating force. In fact, when migrants comprise a biological kin group, migration may be considered a component of genetic drift. The causes of kin-structured migration (KSM) lie in the common, if not universal, tendency for kin to associate and cooperate. However, similar to genetic drift, KSM has its greatest effect in smaller populations and is most apparent in low-density fission-fusion societies such as the Yanomamo of South America and the Semai of Malaysia, and less salient in higher density, low-mobility populations such as those of the New Guinea Highlands. The evolutionary consequences of KSM begin with increased genetic variation among populations. Such intergroup variation provides a basis for group selection. The origin of larger-scale geographic differentiation can arise from kin-structured migrant groups colonizing new regions. Waves of colonizing kin-structured founder groups may produce gene frequency clines, mimicking demic diffusion and natural selection. Finally, because kin structuring reduces the effective size of a population, it may be speculated that the extremely small effective size inferred for ancestral populations of Homo sapiens may be an artifact of kin-structured demographically larger populations.
Two areas in the chicken anemia virus (CAV) genome have high G:C content with secondary structures. These two G:C rich areas could not be sequenced with Perkin Elmer's Big Dye Terminator Cycle Sequencing Kit. Several modifications were carried out to solve the problem. Finally, a package of modified method was developed to sequence the high G:C areas. The result showed that the Perkin Elmer's Big Dye Terminator Cycle Sequencing Kit with the normal procedures are not suitable for sequencing the high G:C regions of the CAV genome. The present developed method made the Perkin Elmer's Kit useful for the first time to sequence the G:C rich hairpin structures of the CAV genome. The system may be useful to sequence all other G:C rich DNA templates.
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.
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.
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.
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.
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.
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.