We have sequenced miRNA libraries from human embryonic, neural and foetal mesenchymal stem cells. We report that the majority of miRNA genes encode mature isomers that vary in size by one or more bases at the 3' and/or 5' end of the miRNA. Northern blotting for individual miRNAs showed that the proportions of isomiRs expressed by a single miRNA gene often differ between cell and tissue types. IsomiRs were readily co-immunoprecipitated with Argonaute proteins in vivo and were active in luciferase assays, indicating that they are functional. Bioinformatics analysis predicts substantial differences in targeting between miRNAs with minor 5' differences and in support of this we report that a 5' isomiR-9-1 gained the ability to inhibit the expression of DNMT3B and NCAM2 but lost the ability to inhibit CDH1 in vitro. This result was confirmed by the use of isomiR-specific sponges. Our analysis of the miRGator database indicates that a small percentage of human miRNA genes express isomiRs as the dominant transcript in certain cell types and analysis of miRBase shows that 5' isomiRs have replaced canonical miRNAs many times during evolution. This strongly indicates that isomiRs are of functional importance and have contributed to the evolution of miRNA genes.
Tibetan high-altitude adaptation (HAA) has been studied extensively, and many candidate genes have been reported. Subsequent efforts targeting HAA functional variants, however, have not been that successful (e.g., no functional variant has been suggested for the top candidate HAA gene, EPAS1). With WinXPCNVer, a method developed in this study, we detected in microarray data a Tibetan-enriched deletion (TED) carried by 90% of Tibetans; 50% were homozygous for the deletion, whereas only 3% carried the TED and 0% carried the homozygous deletion in 2,792 worldwide samples (p < 10(-15)). We employed long PCR and Sanger sequencing technologies to determine the exact copy number and breakpoints of the TED in 70 additional Tibetan and 182 diverse samples. The TED had identical boundaries (chr2: 46,694,276-46,697,683; hg19) and was 80 kb downstream of EPAS1. Notably, the TED was in strong linkage disequilibrium (LD; r(2) = 0.8) with EPAS1 variants associated with reduced blood concentrations of hemoglobin. It was also in complete LD with the 5-SNP motif, which was suspected to be introgressed from Denisovans, but the deletion itself was absent from the Denisovan sequence. Correspondingly, we detected that footprints of positive selection for the TED occurred 12,803 (95% confidence interval = 12,075-14,725) years ago. We further whole-genome deep sequenced (>60×) seven Tibetans and verified the TED but failed to identify any other copy-number variations with comparable patterns, giving this TED top priority for further study. We speculate that the specific patterns of the TED resulted from its own functionality in HAA of Tibetans or LD with a functional variant of EPAS1.
Western honey bees (Apis mellifera) far exceed the commonly observed 1–2 meiotic recombination events per chromosome and exhibit the highest Metazoan recombination rate (20 cM/Mb) described thus far. However, the reasons for this exceptional rate of recombination are not sufficiently understood. In a comparative study, we report on the newly constructed genomic linkage maps of Apis florea and Apis dorsata that represent the two honey bee lineages without recombination rate estimates so far. Each linkage map was generated de novo, based on SNP genotypes of haploid male offspring of a single female. The A. florea map spans 4,782 cM with 1,279 markers in 16 linkage groups. The A. dorsata map is 5,762 cM long and contains 1,189 markers in 16 linkage groups. Respectively, these map sizes result in average recombination rate estimates of 20.8 and 25.1 cM/Mb. Synteny analyses indicate that frequent intra-chromosomal rearrangements but no translocations among chromosomes accompany the high rates of recombination during the independent evolution of the three major honey bee lineages. Our results imply a common cause for the evolution of very high recombination rates in Apis. Our findings also suggest that frequent homologous recombination during meiosis might increase ectopic recombination and rearrangements within but not between chromosomes. It remains to be investigated whether the resulting inversions may have been important in the evolutionary differentiation between honey bee species.
Simulium (Simulium) phraense sp. nov. (Diptera: Simuliidae) is described from females, males, pupae, and larvae from Thailand. This new species is placed in the Simulium striatum species group and is most similar to Simulium (Simulium) nakhonense Takaoka & Suzuki (Diptera: Simuliidae) from Thailand among species of the same species group but is barely distinguished from the latter species by lacking annular ridges on the surface of the pupal gill filaments. The fast-evolving nuclear big zinc finger (BZF) gene has successfully differentiated this new species from its allies, S. (S.) nakhonense and Simulium (Simulium) chiangmaiense Takaoka & Suzuki (Diptera: Simuliidae) of the S. striatum species group. The BZF gene sequences show that this new species is more closely related to S. (S.) nakhonense than to S. (S.) chiangmaiense, further supporting its morphological classification.
Laccases are industrially attractive enzymes and their applications have expanded to the field of bioremediation. The challenge of today's biotechnology in enzymatic studies is to design enzymes that not only have a higher activity but are also more stable and could fit well with the condition requirements. Laccases are known to oxidize non-natural substrates like polycyclic aromatic hydrocarbons (PAHs). We suppose by increasing the hydrophobicity of laccase, it would increase the chance of the enzyme to meet the hydrophobic substrates in a contamination site, therefore increasing the bioremediation efficacy of PAHs from environment. In this attempt, the applications of evolutionary trace (ET), molecular surface accessibility and hydrophobicity analysis on laccase sequences and laccase's crystal structure (1KYA) are described for optimal design of an enzyme with higher hydrophobicity. Our analysis revealed that Q23A, Q45I, N141A, Q237V, N262L, N301V, N331A, Q360L and Q482A could be promising exchanges to be tested in mutagenesis experiments.
A group of 31 rabies viruses (RABVs), recovered primarily from dogs, one deer and one human case, were collected from various areas in China between 1989 and 2006. Complete G gene sequences determined for these isolates indicated identities of nucleotide and amino acid sequences of >or=87% and 93.8%, respectively. Phylogenetic analysis of these and some additional Chinese isolates clearly supported the placement of all Chinese viruses in Lyssavirus genotype 1 and divided all Chinese isolates between four distinct groups (I-IV). Several variants identified within the most commonly encountered group I were distributed according to their geographical origins. A comparison of representative Chinese viruses with other isolates retrieved world-wide indicated a close evolutionary relationship between China group I and II viruses and those of Indonesia while China group III viruses formed an outlying branch to variants from Malaysia and Thailand. China group IV viruses were closely related to several vaccine strains. The predicted glycoprotein sequences of these RABVs variants are presented and discussed with respect to the utility of the anti-rabies biologicals currently employed in China.
A new species of Ansonia is described based on genetic and morphological differentiation. Ansonia lumut sp. nov. is most closely related to three other Peninsular Malaysian species, A. penangensis, A. malayana, and A. jeetsukumarani but differs from these and other congeners by at least 6.9% sequence divergence at the 12S, 16S rRNA and t-RNA-val genes and the following combination of morphological characters: (1) SVL 21.0-23.6 mm in males, 27.7-31.6 mm in females; (2) first finger shorter than second; (3) interorbital and tarsal ridges absent; (4) light interscapular spot absent; (5) presence of large, yellow rictal tubercle; (6) dorsum black with greenish-yellow reticulations; (7) flanks with small yellow spots; (8) fore and hind limbs with yellow cross-bars; and (9) venter light gray with fine, white spotting.
The HIV epidemic is primarily characterised by the circulation of HIV-1 group M (main) comprising of 11 subtypes and sub-subtypes (A1, A2, B-D, F1, F2, G, H, J, and K) and to date 55 circulating recombinant forms (CRFs). In Southeast Asia, active inter-subtype recombination involving three main circulating genotypes--subtype B (including subtype B', the Thai variant of subtype B), CRF01_AE, and CRF33_01B--have contributed to the emergence of novel unique recombinant forms. In the present study, we conducted the molecular epidemiological surveillance of HIV-1 gag-RT genes among 258 people who inject drugs (PWIDs) in Kuala Lumpur, Malaysia, between 2009 and 2011 whereby a novel CRF candidate was recently identified. The near full-length genome sequences obtained from six epidemiologically unlinked individuals showed identical mosaic structures consisting of subtype B' and CRF01_AE, with six unique recombination breakpoints in the gag-RT, pol, and env regions. Among the high-risk population of PWIDs in Malaysia, which was predominantly infected by CRF33_01B (>70%), CRF58_01B circulated at a low but significant prevalence (2.3%, 6/258). Interestingly, the CRF58_01B shared two unique recombination breakpoints with other established CRFs in the region: CRF33_01B, CRF48_01B, and CRF53_01B in the gag gene, and CRF15_01B (from Thailand) in the env gene. Extended Bayesian Markov chain Monte Carlo sampling analysis showed that CRF58_01B and other recently discovered CRFs were most likely to have originated in Malaysia, and that the recent spread of recombinant lineages in the country had little influence from neighbouring countries. The isolation, genetic characterization, and evolutionary features of CRF58_01B among PWIDs in Malaysia signify the increasingly complex HIV-1 diversity in Southeast Asia that may hold an implication on disease treatment, control, and prevention.
As a simian species, the langurs are not known to harbor simian retroviruses, except for one report on a simian Type D endogenous retrovirus from the spectacled langur (Trachypithecus obscurus) from Malaysia. The present report describes for the first time natural infection of the common Hanuman langur (Semnopithecus entellus) from India by a novel simian retrovirus (SRV). The new SRV is phylogenetically related to but distinct from the three molecularly characterized serotypes, SRV 1-3, of the five known serotypes of SRVs, based on sequence analyses from the 3'orf and env regions of the viral genome. The novel SRV isolated from the Indian Hanuman langur is provisionally named SRV-6.
African-born Australians are a recognised "priority population" in Australia's Sixth National HIV/AIDS Strategy. We compared exposure location and route for African-born people living with HIV (PLHIV) in Victoria, Australia, with HIV-1 pol subtype from drug resistance assays and geographical origin suggested by phylogenetic analysis of env gene. Twenty adult HIV positive African-born Victorian residents were recruited via treating doctors. HIV exposure details were obtained from interviews and case notes. Viral RNA was extracted from participant stored plasma or whole blood. The env V3 region was sequenced and compared to globally representative reference HIV-1 sequences in the Los Alamos National Library HIV Database. Twelve participants reported exposure via heterosexual sex and two via iatrogenic blood exposures; four were men having sex with men (MSM); two were exposed via unknown routes. Eight participants reported exposure in their countries of birth, seven in Australia, three in other countries and two in unknown locations. Genotype results (pol) were available for ten participants. HIV env amplification was successful in eighteen cases. HIV-1 subtype was identified in all participants: eight both pol and env; ten env alone and two pol alone. Twelve were subtype C, four subtype B, three subtype A and one subtype CRF02_AG. Reported exposure location was consistent with the phylogenetic clustering of env sequences. African Australians are members of multiple transnational social and sexual networks influencing their exposure to HIV. Phylogenetic analysis may complement traditional surveillance to discern patterns of HIV exposure, providing focus for HIV prevention programs in mobile populations.
Southeast Asia houses various culturally and linguistically diverse ethnic groups. In Malaysia, where the Malay, Chinese, and Indian ethnic groups form the majority, there exist minority groups such as the "negritos" who are believed to be descendants of the earliest settlers of Southeast Asia. Here we report patterns of genetic substructure and admixture in two Malaysian negrito populations (Jehai and Kensiu), using ~50,000 genome-wide single-nucleotide polymorphism (SNP) data. We found traces of recent admixture in both the negrito populations, particularly in the Jehai, with the Malay through principal component analysis and STRUCTURE analysis software, which suggested that the admixture was as recent as one generation ago. We also identified significantly differentiated nonsynonymous SNPs and haplotype blocks related to intracellular transport, metabolic processes, and detection of stimulus. These results highlight the different levels of admixture experienced by the two Malaysian negritos. Delineating admixture and differentiated genomic regions should be of importance in designing and interpretation of molecular anthropology and disease association studies.
Beta-amyloid precursor protein cleavage enzyme 1 (BACE1) and beta-amyloid precursor protein cleavage enzyme 2 (BACE2), members of aspartyl protease family, are close homologues and have high similarity in their protein crystal structures. However, their enzymatic properties differ leading to disparate clinical consequences. In order to identify the residues that are responsible for such differences, we used evolutionary trace (ET) method to compare the amino acid conservation patterns of BACE1 and BACE2 in several mammalian species. We found that, in BACE1 and BACE2 structures, most of the ligand binding sites are conserved which indicate their enzymatic property of aspartyl protease family members. The other conserved residues are more or less randomly localized in other parts of the structures. Four group-specific residues were identified at the ligand binding site of BACE1 and BACE2. We postulated that these residues would be essential for selectivity of BACE1 and BACE2 biological functions and could be sites of interest for the design of selective inhibitors targeting either BACE1 or BACE2.
Lentiviruses (genus Lentivirus) are complex retroviruses that infect a broad range of mammals, including humans. Unlike many other retrovirus genera, lentiviruses have only rarely been incorporated into the mammalian germline. However, a small number of endogenous retrovirus (ERV) lineages have been identified, and these rare genomic "fossils" can provide crucial insights into the long-term history of lentivirus evolution. Here, we describe a previously unreported endogenous lentivirus lineage in the genome of the South African springhare (Pedetes capensis), demonstrating that the host range of lentiviruses has historically extended to rodents (order Rodentia). Furthermore, through comparative and phylogenetic analysis of lentivirus and ERV genomes, considering the biogeographic and ecological characteristics of host species, we reveal broader insights into the long-term evolutionary history of the genus.
Since its discovery in 1969, enterovirus 71 (EV71) has been recognised as a frequent cause of epidemics of hand-foot-and-mouth disease (HFMD) associated with severe neurological sequelae in a small proportion of cases. There has been a significant increase in EV71 epidemic activity throughout the Asia-Pacific region since 1997. Recent HFMD epidemics in this region have been associated with a severe form of brainstem encephalitis associated with pulmonary oedema and high case-fatality rates. The emergence of large-scale epidemic activity in the Asia-Pacific region has been associated with the circulation of three genetic lineages that appear to be undergoing rapid evolutionary change. Two of these lineages (B3 and B4) have not been described previously and appear to have arisen from an endemic focus in equatorial Asia, which has served as a source of virus for HFMD epidemics in Malaysia, Singapore and Australia. The third lineage (C2) has previously been identified [Brown, B.A. et al. (1999) J. Virol. 73, 9969-9975] and was primarily responsible for the large HFMD epidemic in Taiwan during 1998. As EV71 appears not to be susceptible to newly developed antiviral agents and a vaccine is not currently available, control of EV71 epidemics through high-level surveillance and public health intervention needs to be maintained and extended throughout the Asia-Pacific region. Future research should focus on (1) understanding the molecular genetics of EV71 virulence, (2) identification of the receptor(s) for EV71, (3) development of antiviral agents to ameliorate the severity of neurological disease and (4) vaccine development to control epidemics. Following the successful experience of the poliomyelitis control programme, it may be possible to control EV71 epidemics if an effective live-attenuated vaccine is developed.
The concept of long-lived (ancient) lakes has had a great influence on the development of evolutionary biogeography. According to this insight, a number of lakes on Earth have existed for several million years (e.g., Baikal and Tanganyika) and represent unique evolutionary hotspots with multiple intra-basin radiations. In contrast, rivers are usually considered to be variable systems, and the possibility of their long-term existence during geological epochs has never been tested. In this study, we reconstruct the history of freshwater basin interactions across continents based on the multi-locus fossil-calibrated phylogeny of freshwater mussels (Unionidae). These mussels most likely originated in Southeast and East Asia in the Jurassic, with the earliest expansions into North America and Africa (since the mid-Cretaceous) following the colonization of Europe and India (since the Paleocene). We discovered two ancient monophyletic mussel radiations (mean age ~51-55 Ma) within the paleo-Mekong catchment (i.e., the Mekong, Siam, and Malacca Straits paleo-river drainage basins). Our findings reveal that the Mekong may be considered a long-lived river that has existed throughout the entire Cenozoic epoch.
Bamboo and rattan are widely grown for manufacturing, horticulture, and agroforestry. Bamboo and rattan production might help reduce poverty, boost economic growth, mitigate climate change, and protect the natural environment. Despite progress in research, sufficient molecular and genomic resources to study these species are lacking. We launched the Genome Atlas of Bamboo and Rattan (GABR) project, a comprehensive, coordinated international effort to accelerate understanding of bamboo and rattan genetics through genome analysis. GABR includes 2 core subprojects: Bamboo-T1K (Transcriptomes of 1000 Bamboos) and Rattan-G5 (Genomes of 5 Rattans), and several other subprojects. Here we describe the organization, directions, and status of GABR.
Aquatic ecosystems that form major biodiversity hotspots are critically threatened due to environmental and anthropogenic stressors. We believe that, in this genomic era, computational methods can be applied to promote aquatic biodiversity conservation by addressing questions related to the evolutionary history of aquatic organisms at the molecular level. However, huge amounts of genomics data generated can only be discerned through the use of bioinformatics. Here, we examine the applications of next-generation sequencing technologies and bioinformatics tools to study the molecular evolution of aquatic animals and discuss the current challenges and future perspectives of using bioinformatics toward aquatic animal conservation efforts.
Alternative splicing can expand the diversity of proteomes. Homologous mutually exclusive exons (MXEs) originate from the same ancestral exon and result in polypeptides with similar structural properties but altered sequence. Why would some genes switch homologous exons and what are their biological impact? Here, we analyse the extent of sequence, structural and functional variability in MXEs and report the first large scale, structure-based analysis of the biological impact of MXE events from different genomes. MXE-specific residues tend to map to single domains, are highly enriched in surface exposed residues and cluster at or near protein functional sites. Thus, MXE events are likely to maintain the protein fold, but alter specificity and selectivity of protein function. This comprehensive resource of MXE events and their annotations is available at: http://gene3d.biochem.ucl.ac.uk/mxemod/. These findings highlight how small, but significant changes at critical positions on a protein surface are exploited in evolution to alter function.
A growing number of metagenomics-based approaches have been used for the discovery of viruses in insects, cultivated plants, and water in agricultural production systems. In this study, sixteen blueberry root transcriptomes from eight clonally propagated blueberry plants of cultivar 'Emerald' (interspecific hybrid of Vaccinium corymbosum and V. darrowi) generated as part of a separate study on varietal tolerance to soil salinity were analyzed for plant viral sequences. The objective was to determine if the asymptomatic plants harbored the latent blueberry red ringspot virus (BRRV) in their roots. The only currently known mechanism of transmission of BRRV is through vegetative propagation; however, the virus can remain latent for years with some plants of 'Emerald' never developing red ringspot symptoms. Bioinformatic analyses of 'Emerald' transcriptomes using de novo assembly and reference-based mapping approaches yielded eight complete viral genomes of BRRV (genus Soymovirus, family Caulimoviridae). Validation in vitro by PCR confirmed the presence of BRRV in 100% of the 'Emerald' root samples. Sequence and phylogenetic analyses showed 94% to 97% nucleotide identity between BRRV genomes from Florida and sequences from Czech Republic, Japan, Poland, Slovenia, and the United States. Taken together, this study documented the first detection of a complete BRRV genome from roots of asymptomatic blueberry plants and in Florida through in silico analysis of plant transcriptomes.
Ectodysplasin A1 receptor (EDAR) is a TNF receptor family member with roles in the development and growth of hair, teeth and glands. A derived allele of EDAR, single-nucleotide variant rs3827760, encodes EDAR:p.(Val370Ala), a receptor with more potent signalling effects than the ancestral EDAR370Val. This allele of rs3827760 is at very high frequency in modern East Asian and Native American populations as a result of ancient positive selection and has been associated with straighter, thicker hair fibres, alteration of tooth and ear shape, reduced chin protrusion and increased fingertip sweat gland density. Here we report the characterisation of another SNV in EDAR, rs146567337, encoding EDAR:p.(Ser380Arg). The derived allele of this SNV is at its highest global frequency, of up to 5%, in populations of southern China, Vietnam, the Philippines, Malaysia and Indonesia. Using haplotype analyses, we find that the rs3827760 and rs146567337 SNVs arose on distinct haplotypes and that rs146567337 does not show the same signs of positive selection as rs3827760. From functional studies in cultured cells, we find that EDAR:p.(Ser380Arg) displays increased EDAR signalling output, at a similar level to that of EDAR:p.(Val370Ala). The existence of a second SNV with partly overlapping geographic distribution, the same in vitro functional effect and similar evolutionary age as the derived allele of rs3827760, but of independent origin and not exhibiting the same signs of strong selection, suggests a northern focus of positive selection on EDAR function in East Asia.