To date, there are 16 types of CRF01_AE/B circulating recombinant forms identified, and most of them are distributed in Asian countries such as China, Malaysia, and Singapore. Previous HIV molecular epidemiological surveys showed that CRF01_AE (27.6%) and B (9.6%) subtypes are predominant strains in mainland of China. At the same time, the HIV-1 virus spreads faster in the men who have sex with men (MSM) population than in other risk groups. In Shanghai district, ∼66.0% of newly reported cases were infected through homosexual transmission. In this study, we report a novel recombinant strain of CRF01_AE/B. The near full-length genome phylogenetic tree showed that the strain clustered with the CRF01_AE reference sequence and placed in the peripheral position within the branch of the CRF01_AE strain. Subregional evolutionary results indicated that the CRF01_AE subtype was derived from cluster 4 of CRF01_AE, which is mainly distributed in northern China. The subtype B was correlated with the U.S./Europe B, which are widely prevalent in the Chinese MSM population. In recent years, a large number of recombinant forms between CRF01_AE and B strains are continuously emerging in China. Therefore, understanding the current epidemic recombinant forms will have significant implications for prevention and treatment of HIV/AIDS.
In the present study, we evaluated the association between the TP53BP1 Glu353Asp and T-885G polymorphisms and breast cancer risk as well as with the clinicopathological characteristics of the patients. Genotyping of these polymorphisms was performed on 387 breast cancer patients and 252 normal and healthy women who had no history of any malignancy using PCR-RFLP method in a hospital-based Malaysian population. Breast cancer risk was not observed among women who were heterozygous (OR(adj) = 0.887; 95% CI, 0.632-1.245) or homozygous (OR(adj) = 1.083; 95% CI, 0.595-1.969) for Asp allele, and those carriers of Asp allele (OR(adj) = 0.979; 95% CI, 0.771-1.243). Similarly, women who were TG heterozygotes (OR(adj) = 1.181; 95% CI, 0.842-1.658) or GG homozygotes (OR(adj) = 1.362; 95% CI, 0.746-2.486) and carriers of G allele (OR(adj) = 1.147; 95% CI, 0.903-1.458) were not associated with increased risk of breast cancer. Asp allele genotype was significantly associated with ER negativity (p = 0.0015) and poorly differentiated tumours (p = 0.008), but G allele genotype was not associated with the clinicopathological characteristics. In conclusion, Glu353Asp and T-885G polymorphic variants might not have an influence on breast cancer risk, thus might not be potential candidates for cancer susceptibility. Glu353Asp variant might be associated with tumour aggressiveness as defined by its association with ER negativity and poorly differentiated tumours.
Angiostrongylus cantonensis is a zoonotic parasite that causes eosinophilic meningitis in humans. Earlier work on its mitochondrial genome was based on long polymerase chain reaction method. To date, only the mitogenome of the isolates from China has been studied. We report here the complete mitogenome of the Thailand isolate based on next generation sequencing and compare the genetic diversity with other isolates. The mitogenome of the Thailand isolate (13,519bp) is longer than those of the China isolates (13,497-13,502bp). Five protein-coding genes (atp6, cox1, cox2, cob, nad2) show variations in length among the isolates. The stop codon of the Thailand isolate differs from the China and Taiwan isolates in 4 genes (atp6, cob, nad2, nad6). Additionally, the Thailand isolate has 4 incomplete T stop codon compared to 3 in the China and Taiwan isolates. The control region is longer in the Thailand isolate (258bp) than the China (230-236bp) and Taiwan (237bp) isolates. The intergenic sequence between nad4 and cox1 genes in the Thailand isolate lacks 2bp (indels) at the 5'-end of the sequence as well as differs at 7 other sites compared to the China and Taiwan isolates. In the Thailand isolate, 18 tRNAs lack the entire TΨC-arm, compared to 17 in the China isolate and 16 in the Taiwan isolate. Phylogenetic analyses based on 36 mt-genes, 12 PCGs, 2 rRNA genes, 22 tRNA genes and control region all indicate closer genetic affinity between the China and Taiwan isolates compared to the Thailand isolate. Based on 36 mt-genes, the inter-isolate genetic distance varies from p=3.2% between China and Taiwan isolates to p=11.6% between Thailand and China isolates. The mitogenome will be useful for population, phylogenetics and phylogeography studies.
Angiostrongylus malaysiensis is a nematode parasite of various rat species. When first documented in Malaysia, it was referred to as A. cantonensis. Unlike A. cantonensis, the complete mitochondrial genome of A. malaysiensis has not been documented. We report here its complete mitogenome, its differentiation from A. cantonensis, and the phylogenetic relationships with its congeners and other Metastrongyloid taxa. The whole mitogenome of A. malaysiensis had a total length of 13,516bp, comprising 36 genes (12 PCGs, 2 rRNA and 22 tRNA genes) and a control region. It is longer than that of A. cantonensis (13,509bp). Its control region had a long poly T-stretch of 12bp which was not present in A. cantonensis. A. malaysiensis and A. cantonensis had identical start codon for the 12 PCGs, but four PCGs (atp6, cob, nad2, nad6) had different stop codon. The cloverleaf structure for the 22 tRNAs was similar in A. malaysiensis and A. cantonensis except the TΨC-arm was absent in trnV for A. malaysiensis but present in A. cantonensis. The Angiostrongylus genus was monophyletic, with A. malaysiensis and A. cantonensis forming a distinct lineage from that of A. costaricensis and A. vasorum. The genetic distance between A. malaysiensis and A. cantonensis was p=11.9% based on 12 PCGs, p=9.5% based on 2 rRNA genes, and p=11.6% based on 14 mt-genes. The mitogenome will prove useful for studies on phylogenetics and systematics of Angiostrongylus lungworms and other Metastrongyloid nematodes.
Angiostrongylus cantonensis is a bursate nematode parasite that causes eosinophilic meningitis (or meningoencephalitis) in humans in many parts of the world. The genomic data from A. cantonensis will form a useful resource for comparative genomic and chemogenomic studies to aid the development of diagnostics and therapeutics. We have sequenced, assembled and annotated the genome of A. cantonensis. The genome size is estimated to be ∼260 Mb, with 17,280 genomic scaffolds, 91X coverage, 81.45% for complete and 93.95% for partial score based on CEGMA analysis of genome completeness. The number of predicted genes of ≥300 bp was 17,482. A total of 7737 predicted protein-coding genes of ≥50 amino acids were identified in the assembled genome. Among the proteins of known function, kinases are the most abundant followed by transferases. The draft genome contains 34 excretory-secretory proteins (ES), a minimum of 44 Nematode Astacin (NAS) metalloproteases, 12 Homeobox (HOX) genes, and 30 neurotransmitters. The assembled genome size (260 Mb) is larger than those of Pristionchus pacificus, Caenorhabditis elegans, Necator americanus, Caenorhabditis briggsae, Trichinella spiralis, Brugia malayi and Loa loa, but smaller than Haemonchus contortus and Ascaris suum. The repeat content (25%) is similar to H. contortus. The GC content (41.17%) is lower compared to P. pacificus (42.7%) and H. contortus (43.1%) but higher compared to C. briggsae (37.69%), A. suum (37.9%) and N. americanus (40.2%) while the scaffold N50 is 42,191. This draft genome will facilitate the understanding of many unresolved issues on the parasite and the disorder it causes.
This chapter introduces different aspects of bioinformatics with a brief discussion in the systems biology context. Example applications in network pharmacology of traditional Chinese medicine, systems metabolic engineering, and plant genome-scale modelling are described. Lastly, this chapter concludes on how bioinformatics helps to integrate omics data derived from various studies described in previous chapters for a holistic understanding of secondary metabolite production in P. minus.
Hypertension is a silent killer worldwide, caused by both genetic and environmental factors. Until now, genetic and genomic association studies of hypertension are reporting different degree of association on hypertension. Hence, it is essential to gather all the available information on the reported genetic loci and to determine if any biomarker(s) is/are significantly associated with hypertension. Current review concluded the potential biomarkers for hypertension, with regards to electrolyte and fluid transports, as well as sodium/potassium ions homeostasis, which are supported by the results of case-controls and meta-analyses.
Genome-wide association studies (GWAS) have identified over 100 single nucleotide polymorphisms (SNPs) associated with prostate cancer. However, information on the mechanistic basis for some associations is limited. Recent research has been directed towards the potential association of vitamin D concentrations and prostate cancer, but little is known about whether the aforementioned genetic associations are modified by vitamin D. We investigated the associations of 46 GWAS-identified SNPs, circulating concentrations of 25-hydroxyvitamin D (25(OH)D), and prostate cancer (3,811 cases, 511 of whom died from the disease, compared with 2,980 controls-from 5 cohort studies that recruited participants over several periods beginning in the 1980s). We used logistic regression models with data from the National Cancer Institute Breast and Prostate Cancer Cohort Consortium (BPC3) to evaluate interactions on the multiplicative and additive scales. After allowing for multiple testing, none of the SNPs examined was significantly associated with 25(OH)D concentration, and the SNP-prostate cancer associations did not differ by these concentrations. A statistically significant interaction was observed for each of 2 SNPs in the 8q24 region (rs620861 and rs16902094), 25(OH)D concentration, and fatal prostate cancer on both multiplicative and additive scales (P ≤ 0.001). We did not find strong evidence that associations between GWAS-identified SNPs and prostate cancer are modified by circulating concentrations of 25(OH)D. The intriguing interactions between rs620861 and rs16902094, 25(OH)D concentration, and fatal prostate cancer warrant replication.
Matched MeSH terms: Genome-Wide Association Study*
Whole-genome sequencing across multiple samples in a population provides an unprecedented opportunity for comprehensively characterizing the polymorphic variants in the population. Although the 1000 Genomes Project (1KGP) has offered brief insights into the value of population-level sequencing, the low coverage has compromised the ability to confidently detect rare and low-frequency variants. In addition, the composition of populations in the 1KGP is not complete, despite the fact that the study design has been extended to more than 2,500 samples from more than 20 population groups. The Malays are one of the Austronesian groups predominantly present in Southeast Asia and Oceania, and the Singapore Sequencing Malay Project (SSMP) aims to perform deep whole-genome sequencing of 100 healthy Malays. By sequencing at a minimum of 30× coverage, we have illustrated the higher sensitivity at detecting low-frequency and rare variants and the ability to investigate the presence of hotspots of functional mutations. Compared to the low-pass sequencing in the 1KGP, the deeper coverage allows more functional variants to be identified for each person. A comparison of the fidelity of genotype imputation of Malays indicated that a population-specific reference panel, such as the SSMP, outperforms a cosmopolitan panel with larger number of individuals for common SNPs. For lower-frequency (<5%) markers, a larger number of individuals might have to be whole-genome sequenced so that the accuracy currently afforded by the 1KGP can be achieved. The SSMP data are expected to be the benchmark for evaluating the value of deep population-level sequencing versus low-pass sequencing, especially in populations that are poorly represented in population-genetics studies.
Previous research has shown that polygenic risk scores (PRSs) can be used to stratify women according to their risk of developing primary invasive breast cancer. This study aimed to evaluate the association between a recently validated PRS of 313 germline variants (PRS313) and contralateral breast cancer (CBC) risk. We included 56,068 women of European ancestry diagnosed with first invasive breast cancer from 1990 onward with follow-up from the Breast Cancer Association Consortium. Metachronous CBC risk (N = 1,027) according to the distribution of PRS313 was quantified using Cox regression analyses. We assessed PRS313 interaction with age at first diagnosis, family history, morphology, ER status, PR status, and HER2 status, and (neo)adjuvant therapy. In studies of Asian women, with limited follow-up, CBC risk associated with PRS313 was assessed using logistic regression for 340 women with CBC compared with 12,133 women with unilateral breast cancer. Higher PRS313 was associated with increased CBC risk: hazard ratio per standard deviation (SD) = 1.25 (95%CI = 1.18-1.33) for Europeans, and an OR per SD = 1.15 (95%CI = 1.02-1.29) for Asians. The absolute lifetime risks of CBC, accounting for death as competing risk, were 12.4% for European women at the 10th percentile and 20.5% at the 90th percentile of PRS313. We found no evidence of confounding by or interaction with individual characteristics, characteristics of the primary tumor, or treatment. The C-index for the PRS313 alone was 0.563 (95%CI = 0.547-0.586). In conclusion, PRS313 is an independent factor associated with CBC risk and can be incorporated into CBC risk prediction models to help improve stratification and optimize surveillance and treatment strategies.
Matched MeSH terms: Genome, Human*; Genome-Wide Association Study
By combining data from 160,500 individuals with breast cancer and 226,196 controls of Asian and European ancestry, we conducted genome- and transcriptome-wide association studies of breast cancer. We identified 222 genetic risk loci and 137 genes that were associated with breast cancer risk at a p
Matched MeSH terms: Genome-Wide Association Study*
In polygenic score (PGS) analysis, the coefficient of determination (R2) is a key statistic to evaluate efficacy. R2 is the proportion of phenotypic variance explained by the PGS, calculated in a cohort that is independent of the genome-wide association study (GWAS) that provided estimates of allelic effect sizes. The SNP-based heritability (hSNP2, the proportion of total phenotypic variances attributable to all common SNPs) is the theoretical upper limit of the out-of-sample prediction R2. However, in real data analyses R2 has been reported to exceed hSNP2, which occurs in parallel with the observation that hSNP2 estimates tend to decline as the number of cohorts being meta-analyzed increases. Here, we quantify why and when these observations are expected. Using theory and simulation, we show that if heterogeneities in cohort-specific hSNP2 exist, or if genetic correlations between cohorts are less than one, hSNP2 estimates can decrease as the number of cohorts being meta-analyzed increases. We derive conditions when the out-of-sample prediction R2 will be greater than hSNP2 and show the validity of our derivations with real data from a binary trait (major depression) and a continuous trait (educational attainment). Our research calls for a better approach to integrating information from multiple cohorts to address issues of between-cohort heterogeneity.
Matched MeSH terms: Genome-Wide Association Study*
Genome-wide association studies (GWASs) have revealed SNP rs889312 on 5q11.2 to be associated with breast cancer risk in women of European ancestry. In an attempt to identify the biologically relevant variants, we analyzed 909 genetic variants across 5q11.2 in 103,991 breast cancer individuals and control individuals from 52 studies in the Breast Cancer Association Consortium. Multiple logistic regression analyses identified three independent risk signals: the strongest associations were with 15 correlated variants (iCHAV1), where the minor allele of the best candidate, rs62355902, associated with significantly increased risks of both estrogen-receptor-positive (ER(+): odds ratio [OR] = 1.24, 95% confidence interval [CI] = 1.21-1.27, ptrend = 5.7 × 10(-44)) and estrogen-receptor-negative (ER(-): OR = 1.10, 95% CI = 1.05-1.15, ptrend = 3.0 × 10(-4)) tumors. After adjustment for rs62355902, we found evidence of association of a further 173 variants (iCHAV2) containing three subsets with a range of effects (the strongest was rs113317823 [pcond = 1.61 × 10(-5)]) and five variants composing iCHAV3 (lead rs11949391; ER(+): OR = 0.90, 95% CI = 0.87-0.93, pcond = 1.4 × 10(-4)). Twenty-six percent of the prioritized candidate variants coincided with four putative regulatory elements that interact with the MAP3K1 promoter through chromatin looping and affect MAP3K1 promoter activity. Functional analysis indicated that the cancer risk alleles of four candidates (rs74345699 and rs62355900 [iCHAV1], rs16886397 [iCHAV2a], and rs17432750 [iCHAV3]) increased MAP3K1 transcriptional activity. Chromatin immunoprecipitation analysis revealed diminished GATA3 binding to the minor (cancer-protective) allele of rs17432750, indicating a mechanism for its action. We propose that the cancer risk alleles act to increase MAP3K1 expression in vivo and might promote breast cancer cell survival.
Genome-wide association studies have identified SNPs near ZNF365 at 10q21.2 that are associated with both breast cancer risk and mammographic density. To identify the most likely causal SNPs, we fine mapped the association signal by genotyping 428 SNPs across the region in 89,050 European and 12,893 Asian case and control subjects from the Breast Cancer Association Consortium. We identified four independent sets of correlated, highly trait-associated variants (iCHAVs), three of which were located within ZNF365. The most strongly risk-associated SNP, rs10995201 in iCHAV1, showed clear evidence of association with both estrogen receptor (ER)-positive (OR = 0.85 [0.82-0.88]) and ER-negative (OR = 0.87 [0.82-0.91]) disease, and was also the SNP most strongly associated with percent mammographic density. iCHAV2 (lead SNP, chr10: 64,258,684:D) and iCHAV3 (lead SNP, rs7922449) were also associated with ER-positive (OR = 0.93 [0.91-0.95] and OR = 1.06 [1.03-1.09]) and ER-negative (OR = 0.95 [0.91-0.98] and OR = 1.08 [1.04-1.13]) disease. There was weaker evidence for iCHAV4, located 5' of ADO, associated only with ER-positive breast cancer (OR = 0.93 [0.90-0.96]). We found 12, 17, 18, and 2 candidate causal SNPs for breast cancer in iCHAVs 1-4, respectively. Chromosome conformation capture analysis showed that iCHAV2 interacts with the ZNF365 and NRBF2 (more than 600 kb away) promoters in normal and cancerous breast epithelial cells. Luciferase assays did not identify SNPs that affect transactivation of ZNF365, but identified a protective haplotype in iCHAV2, associated with silencing of the NRBF2 promoter, implicating this gene in the etiology of breast cancer.
Genome-wide polygenic risk scores (GW-PRSs) have been reported to have better predictive ability than PRSs based on genome-wide significance thresholds across numerous traits. We compared the predictive ability of several GW-PRS approaches to a recently developed PRS of 269 established prostate cancer-risk variants from multi-ancestry GWASs and fine-mapping studies (PRS269). GW-PRS models were trained with a large and diverse prostate cancer GWAS of 107,247 cases and 127,006 controls that we previously used to develop the multi-ancestry PRS269. Resulting models were independently tested in 1,586 cases and 1,047 controls of African ancestry from the California Uganda Study and 8,046 cases and 191,825 controls of European ancestry from the UK Biobank and further validated in 13,643 cases and 210,214 controls of European ancestry and 6,353 cases and 53,362 controls of African ancestry from the Million Veteran Program. In the testing data, the best performing GW-PRS approach had AUCs of 0.656 (95% CI = 0.635-0.677) in African and 0.844 (95% CI = 0.840-0.848) in European ancestry men and corresponding prostate cancer ORs of 1.83 (95% CI = 1.67-2.00) and 2.19 (95% CI = 2.14-2.25), respectively, for each SD unit increase in the GW-PRS. Compared to the GW-PRS, in African and European ancestry men, the PRS269 had larger or similar AUCs (AUC = 0.679, 95% CI = 0.659-0.700 and AUC = 0.845, 95% CI = 0.841-0.849, respectively) and comparable prostate cancer ORs (OR = 2.05, 95% CI = 1.87-2.26 and OR = 2.21, 95% CI = 2.16-2.26, respectively). Findings were similar in the validation studies. This investigation suggests that current GW-PRS approaches may not improve the ability to predict prostate cancer risk compared to the PRS269 developed from multi-ancestry GWASs and fine-mapping.
Dromedary camels are a domestic species characterized by various adaptive traits. Limited efforts have been employed toward identifying genetic regions and haplotypes under selection that might be related to such adaptations. These genetic elements are considered valuable sources that should be conserved to maintain the dromedaries' adaptability. Here, we have analyzed whole genome sequences of 40 dromedary camels from different Arabian Peninsula populations to assess their genetic relationship and define regions with signatures of selection. Genetic distinction based on geography was observed, classifying the populations into four groups: (1) North and Central, (2) West, (3) Southwest, and (4) Southeast, with substantial levels of genetic admixture. Using the de-correlated composite of multiple signal approach, which combines four intra-population analyses (Tajima's D index, nucleotide diversity, integrated haplotype score, and number of segregating sites by length), a total of 36 candidate regions harboring 87 genes were identified to be under positive selection. These regions overlapped with 185 haplotype blocks encompassing 1 340 haplotypes, of which 30 (∼2%) were found to be approaching fixation. The defined candidate genes are associated with different biological processes related to the dromedaries' adaptive physiologies, including neurological pathways, musculoskeletal development, fertility, fat distribution, immunity, visual development, and kidney physiology. The results of this study highlight opportunities for further investigations at the whole-genome level to enhance our understanding of the evolutionary pressures shaping the dromedary genome.
Understanding dietary diversity is a fundamental task in the study of stump-tailed macaque, Macaca arctoides in its natural habitat. However, direct feeding observation and morphological identification using fecal samples are not effective and nearly impossible to obtain in natural habitats because this species is sensitive to human presence. As ecological methods are challenging and time-consuming, DNA metabarcoding offers a more powerful assessment of the diet. We used a chloroplast tRNL DNA metabarcoding approach to identify the diversity of plants consumed by free-ranging M. arctoides in the Malaysia-Thailand border region located in Perlis State Park, Peninsular Malaysia. DNA was extracted from three fecal samples, and chloroplast tRNL DNA was amplified and sequenced using the Illumina MiniSeq platform. Sequences were analyzed using the CLC Genomic Workbench software. A total of 145 plant species from 46 families were successfully identified as being consumed by M. arctoides. The most abundant species were yellow saraca, Saraca thaipingensis (11.70%), common fig, Ficus carica (9.33%), aramata, Clathrotropis brachypetala (5.90%), sea fig, Ficus superba (5.44%), and envireira, Malmea dielsiana (1.70%). However, Clathrotropis and Malmea are not considered Malaysian trees because of limited data available from Malaysian plant DNA. Our study is the first to identify plant taxa up to the species level consumed by stump-tailed macaques based on a DNA metabarcoding approach. This result provides an important understanding on diet of wild M. arctoides that only reside in Perlis State Park, Malaysia.
The firefly genus Luciola sensu McDermott contains 282 species that are distributed across major parts of Asia, Europe, Africa, Australia, and the Pacific islands. Due to phenotypic similarities, species identification using external morphological characters can be unreliable for this group. Consequently, decades of piecemeal taxonomic treatments have resulted in numerous erroneous and contentious classifications. Furthermore, our understanding of the group's evolutionary history is limited due to the lack of a robust phylogenetic framework that has also impeded efforts to stabilize its taxonomy. Here, we constructed molecular phylogenies of Luciola and its allies based on combined mitogenomes and Cytochrome c oxidase subunit 1 (COX1) sequences including a newly sequenced mitogenome of an unidentified taxon from Singapore. Our results showed that this taxon represents a distinct and hitherto undescribed evolutionary lineage that forms a clade with L. filiformis from Japan and L. curtithorax from China. Additionally, the Singaporean lineage can be differentiated from other congeners through several external and internal diagnostic morphological characters, and is thus described herein as a new species. Our phylogeny also strongly supported the paraphyly of Luciola with regard to L. cruciata and L. owadai, which were inferred to be more closely related to the genus Aquatica as opposed to other members of Luciola sensu stricto. The genus Hotaria was inferred as a derived clade within Luciola (sister to L. italica), supporting its status as a subgenus of Luciola instead of a distinct genus. This is the first time since 1909 that a new species of luminous firefly has been discovered in Singapore, highlighting the need for continued biodiversity research, even in small, well-studied and highly developed countries, such as Singapore.