PubMed Health⌕ Search

Biomedical subjects

Yun-Xin Fu

Publications and source records attributed to Yun-Xin Fu.

At least 19 recordsLinked to original sources

Molecular motions in HIV-1 gp120 mutants reveal their preferences for different conformations.

Both the crystal structures of the HIV-1 gp120 core bound by the CD4 and antigen 17b, and the SIV gp120 core pre-bound by CD4 are known. We built the homology models of gp120 with loops V3 and V4 in the CD4-complex, CD4-free and CD4-unliganded states, and models of the 375 S/W and 423 I/P mutants in the CD4-free and unliganded states, respectively. CONCOORD was utilized for generating ensembles of the seven gp120 models that were analyzed by essential dynamics analyses to identify their preferred concerted motions. The revealed large-scale concerted motions are related to either the receptor association/release or the conformational transition between different conformational states. Essential subspace overlap analyses were performed to quantitatively distinguish the preference for conformational transitions between states of the gp120 mutants and further to ascertain what kind of conformational state that the mutants prefer to adopt. Results indicate that the 375 S/W mutant, in which the tryptophan indole group is predicted to occupy the phe43 pocket in the gp120 interior, favors a conformation close to the CD4-bound state. However, the other mutant 423 I/P inclines to prevent the formation of bridging sheet and stabilize the conformation in the unliganded state. Our theoretical analyses are in agreement with experimentally determined mutation effects, and can be extended to a new approach to design or screen mutants that have effects on conformation/function of a protein.

Amino Acid Sequence↗

Dynamic domains and geometrical properties of HIV-1 gp120 during conformational changes induced by CD4 binding.

The HIV-1 gp120 exterior envelope glycoprotein undergoes a series of conformational rearrangements while sequentially interacting with the receptor CD4 and coreceptor CCR5 or CXCR4 on the surface of host cells to initiate virus entry. Both the crystal structures of the HIV-1 gp120 core bound by the CD4 and antigen 17b, and the SIV gp120 core pre-bound by the CD4 are known. We have performed dynamic domain studies on the homology models of the CD4-bound and unliganded HIV-1 gp120 with modeled V3 and V4 loops to explore details of conformational changes, hinge axes, and hinge bending regions in the gp120 structures upon CD4 binding. Four dynamic domains were clustered and intricately motional modes for domain pairs were discovered. Together with the detailed comparative analyses of geometrical properties between the unliganded and liganded gp120 models, an induced fit model was proposed to explain events accompanying the CD4 engagement to the gp120, which provided new insight into the dynamics of the molecular induced binding mechanism that complements the molecular dynamics and crystallographic studies.

Binding Sites↗

A new strategy for estimating recombination fractions between dominant markers from an F2 population.

Although most high-density linkage maps have been constructed from codominant markers such as single-nucleotide polymorphisms (SNPs) and microsatellites due to their high linkage information, dominant markers can be expected to be even more significant as proteomic technique becomes widely applicable to generate protein polymorphism data from large samples. However, for dominant markers, two possible linkage phases between a pair of markers complicate the estimation of recombination fractions between markers and consequently the construction of linkage maps. The low linkage information of the repulsion phase and high linkage information of coupling phase have led geneticists to construct two separate but related linkage maps. To circumvent this problem, we proposed a new method for estimating the recombination fraction between markers, which greatly improves the accuracy of estimation through distinction between the coupling phase and the repulsion phase of the linked loci. The results obtained from both real and simulated F2 dominant marker data indicate that the recombination fractions estimated by the new method contain a large amount of linkage information for constructing a complete linkage map. In addition, the new method is also applicable to data with mixed types of markers (dominant and codominant) with unknown linkage phase.

Algorithms↗

Evidence for recombination in Mycobacterium tuberculosis.

Due to its mostly isolated living environment, Mycobacterium tuberculosis is generally believed to be highly clonal, and thus recombination between different strains must be rare and is not critical for the survival of the species. To investigate the roles recombination could have possibly played in the evolution of M. tuberculosis, an analysis was conducted on previously determined genotypes of 36 synonymous single nucleotide polymorphisms (SNPs) in 3,320 M. tuberculosis isolates. The results confirmed the predominant clonal structure of the M. tuberculosis population. However, recombination between different strains was also suggested. To further resolve the issue, 175 intergenic SNPs and 234 synonymous SNPs were genotyped in 37 selected representative strains. A clear mosaic polymorphic pattern ahead of the MT0105 locus encoding a PPE (Pro-Pro-Glu) protein was obtained, which is most likely a result of recombination hot spot. Given that PPE proteins are thought to be critical in host-pathogen interactions, we hypothesize that recombination has been influential in the history of M. tuberculosis and possibly a major contributor to the diversity observed ahead of the MT0105 locus.

Antigens, Bacterial↗

Ranking analysis of microarray data: a powerful method for identifying differentially expressed genes.

Microarray technology provides a powerful tool for the expression profile of thousands of genes simultaneously, which makes it possible to explore the molecular and metabolic etiology of the development of a complex disease under study. However, classical statistical methods and technologies fail to be applicable to microarray data. Therefore, it is necessary and motivating to develop powerful methods for large-scale statistical analyses. In this paper, we described a novel method, called Ranking Analysis of Microarray Data (RAM). RAM, which is a large-scale two-sample t-test method, is based on comparisons between a set of ranked T statistics and a set of ranked Z values (a set of ranked estimated null scores) yielded by a "randomly splitting" approach instead of a "permutation" approach and a two-simulation strategy for estimating the proportion of genes identified by chance, i.e., the false discovery rate (FDR). The results obtained from the simulated and observed microarray data show that RAM is more efficient in identification of genes differentially expressed and estimation of FDR under undesirable conditions such as a large fudge factor, small sample size, or mixture distribution of noises than Significance Analysis of Microarrays.

Animals↗

Statistical tests for detecting positive selection by utilizing high-frequency variants.

By comparing the low-, intermediate-, and high-frequency parts of the frequency spectrum, we gain information on the evolutionary forces that influence the pattern of polymorphism in population samples. We emphasize the high-frequency variants on which positive selection and negative (background) selection exhibit different effects. We propose a new estimator of theta (the product of effective population size and neutral mutation rate), thetaL, which is sensitive to the changes in high-frequency variants. The new thetaL allows us to revise Fay and Wu's H-test by normalization. To complement the existing statistics (the H-test and Tajima's D-test), we propose a new test, E, which relies on the difference between thetaL and Watterson's thetaW. We show that this test is most powerful in detecting the recovery phase after the loss of genetic diversity, which includes the postselective sweep phase. The sensitivities of these tests to (or robustness against) background selection and demographic changes are also considered. Overall, D and H in combination can be most effective in detecting positive selection while being insensitive to other perturbations. We thus propose a joint test, referred to as the DH test. Simulations indicate that DH is indeed sensitive primarily to directional selection and no other driving forces.

Animals↗

Growth hormone locus expands and diverges after the separation of New and Old World Monkeys.

While most mammals including the prosimians have a single copy of the growth hormone (GH) gene, anthropoids possess a cluster of GH-related genes. Throughout the evolution of the main anthropoid groups [New World Monkeys (NWM), Old World Monkeys (OWM), and apes], two features stand out of the GH loci. The first is the appearance of chorionic somatommamotropin hormone (CSH) genes within the OWM lineage and the second is the expansion of the loci intergenic regions in the OWM and apes. In relation with this loci expansion, the NWM possess intergenic regions of homogeneous lengths (3.5 kb). In contrast, heterogeneous lengths (6 and 13 kb) have been reported for species of the OWM. At the present, none of the OWM genomic GH loci organizations have been described. Here, we report the genomic organization of the GH locus in the rhesus monkey, this locus has six GH-related genes separated by five intergenic regions. The 5' end gene (GH-1) encodes for the pituitary GH and is followed by CSH-1, GH-2, CSH-2, CSH-3 and CSH-4 genes. The five intergenic regions have heterogeneous lengths and also present more or less the same Alu distribution as the human GH locus. To analyze the events that contributed to the extension of the intergenic regions of the GH locus and the emergence of the regulatory elements, the five GH locus intergenic regions of the spider monkey (NWM) were sequenced. The results of comparing the loci from both species suggest that the long intergenic regions (13 kb) of the rhesus GH locus share a common ancestor with the 3.5 kb intergenic regions of the spider monkey. However, the observed increased length of the former is due to an insertion (approximately 8.7 kb) at their 3' end. Interestingly in this insert, we discovered a DNA element resembling the enhancer of the CSH genes of the human GH locus. On the other hand, we observed that the short intergenic regions (6 kb) increased by a different recombination event.

Alu Elements↗

A novel method for estimating linkage maps.

The goal of linkage mapping is to find the true order of loci from a chromosome. Since the number of possible orders is large even for a modest number of loci, the problem of finding the optimal solution is known as a NP-hard problem or traveling salesman problem (TSP). Although a number of algorithms are available, many either are low in the accuracy of recovering the true order of loci or require tremendous amounts of computational resources, thus making them difficult to use for reconstructing a large-scale map. We developed in this article a novel method called unidirectional growth (UG) to help solve this problem. The UG algorithm sequentially constructs the linkage map on the basis of novel results about additive distance. It not only is fast but also has a very high accuracy in recovering the true order of loci according to our simulation studies. Since the UG method requires n-1 cycles to estimate the ordering of n loci, it is particularly useful for estimating linkage maps consisting of hundreds or even thousands of linked codominant loci on a chromosome.

Algorithms↗

Nucleotide variation and haplotype diversity in a 10-kb noncoding region in three continental human populations.

Noncoding regions are usually less subject to natural selection than coding regions and so may be more useful for studying human evolution. The recent surveys of worldwide DNA variation in four 10-kb noncoding regions revealed many interesting but also some incongruent patterns. Here we studied another 10-kb noncoding region, which is in 6p22. Sixty-six single-nucleotide polymorphisms were found among the 122 worldwide human sequences, resulting in 46 genotypes, from which 48 haplotypes were inferred. The distribution patterns of DNA variation, genotypes, and haplotypes suggest rapid population expansion in relatively recent times. The levels of polymorphism within human populations and divergence between humans and chimpanzees at this locus were generally similar to those for the other four noncoding regions. Fu and Li's tests rejected the neutrality assumption in the total sample and in the African sample but Tajima's test did not reject neutrality. A detailed examination of the contributions of various types of mutations to the parameters used in the neutrality tests clarified the discrepancy between these test results. The age estimates suggest a relatively young history in this region. Combining three autosomal noncoding regions, we estimated the long-term effective population size of humans to be 11,000 +/- 2800 using Tajima's estimator and 17,600 +/- 4700 using Watterson's estimator and the age of the most recent common ancestor to be 860,000 +/- 258,000 years ago.

Africa↗

Exact coalescent for the Wright-Fisher model.

The Kingman coalescent, which has become the foundation for a wide range of theoretical as well as empirical studies, was derived as an approximation of the Wright-Fisher (WF) model. The approximation heavily relies on the assumption that population size is large and sample size is much smaller than the population size. Whether the sample size is too large compared to the population size is rarely questioned in practice when applying statistical methods based on the Kingman coalescent. Since WF model is the most widely used population genetics model for reproduction, it is desirable to develop a coalescent framework for the WF model, which can be used whenever there are concerns about the accuracy of the Kingman coalescent as an approximation. This paper described the exact coalescent theory for the WF model and develops a simulation algorithm, which is then used, together with an analytical approach, to study the properties of the exact coalescent as well as its differences to the Kingman coalescent. We show that the Kingman coalescent differs from the exact coalescent by: (1) shorter waiting time between successive coalescent events; (2) different probability of observing a topological relationship among sequences in a sample; and (3) slightly smaller tree length in the genealogy of a large sample. On the other hand, there is little difference in the age of the most recent common ancestor (MRCA) of the sample. The exact coalescent makes up the longer waiting time between successive coalescent events by having multiple coalescence at the same time. The most significant difference among various summary statistics of a coalescent examined is the sum of lengths of external branches, which can be more than 10% larger for exact coalescent than that for the Kingman coalescent. As a whole, the Kingman coalescent is a remarkably accurate approximation to the exact coalescent for sample and population sizes falling considerably outside the region that was originally anticipated.

Algorithms↗

American founder mutation for Lynch syndrome. Prevalence estimates and implications.

BACKGROUND: Recently, a new founder mutation, an exon 1-6 deletion in a mismatch repair gene (MMR), MSH2, in nine kindreds with Lynch syndrome was reported. In 3 of the kindreds this mutation was traced by genealogy through 11-12 generations to a common founder, and thus termed the American Founder Mutation (AFM). Since then, 13 additional 'unrelated' kindreds with AFM were detected by a recently designed single polymerase chain reaction. This test might serve as first-line screening for Lynch syndrome mutations, provided AFM was prevalent, which is assessed in the current study. METHODS: The number of current AFM carriers and the incidence of Lynch syndrome caused by AFM were estimated based on population growth of mutation carriers derived from genealogy data. For cross-checking, its annual incidence was also estimated based on published epidemiology data. RESULTS: There are 18,981 (5th and 95th percentiles, 6038 and 34,466, respectively) expected current AFM carriers, or 160 (range 51-290) Lynch syndrome cases diagnosed per year due to AFM estimated based on genealogy data. The incidence estimate closely overlaps with that based on published epidemiology data, which is 114-400 cases per year. CONCLUSIONS: A large number of AFM carriers are likely to exist in the U.S., which harbors significant implications for cancer control. Given the ease of detection, testing for AFM not only among members of the existing AFM families, but also among all patients with Lynch syndrome in the U.S. is proposed.

Colorectal Neoplasms, Hereditary Nonpolyposis↗

Evolutionary implications of Avian Infectious Bronchitis Virus (AIBV) analysis.

For developing efficient vaccines, it is essential to identify which amino acid changes are most important to the survival of the virus. We investigate the amino acid substitution features in the Avian Infectious Bronchitis Virus (AIBV) antigenic domain of a vaccine serotype (DE072) and a virulent viral strain (GA98) to better understand adaptive evolution of AIBV. In addition, the SARS Coronavirus (SARS-CoV) was also analyzed in the same way. It is interesting to find that extreme comparability exists between AIBV and SARS in amino acid substitution pattern. It suggests that amino acid changes that result in overall shift of residue charge and polarity should be paid special attention to during the development of vaccines.

Amino Acid Substitution↗

[Population genetics of Rhinopithecus bieti: a study of the mitochondrial control region].

Yunnan snub-nose monkey (Rhinopithecus bieti) is a famous endangered primate in China. So far, however, studies on its population genetics based on DNA sequences are not available. In this paper, the whole mitochondria control region of the samples from Weixi, Yunan Province as well as the whole cytochrome b gene in some individuals were sequenced. A deep divergence was observed within the Weixi population, which was confirmed after excluding the possibility of it being a nuclear pseudogene. Nonetheless, if the effects of population structure and migration are considered, the true level of polymorphism of the Weixi population may be not as high as observed.

Animals↗

Bone mineral density of recent African immigrants in the United States.

BACKGROUND: Racial/ethnic difference in bone mineral density (BMD) exists. The underlying mechanism is unclear and needs investigation. PURPOSE: To determine BMD and its relation to environmental exposure in recent African immigrants. METHODS: BMD in recent sub-Saharan Sudanese immigrants (55 men and 88 premenopausal women) in the United States was measured. Analysis of covariance (ANCOVA) model was performed, with total body, spine and hip BMD as dependent variables; and sex, age, body weight, the length of stay in the United States, and milk intake as independent variables. RESULTS: BMD Z score in the spine but not total body or hip in the Sudanese immigrants was significantly lower compared with the normative values of African Americans and Caucasians. Total body and hip BMD was positively correlated (p < 0.015) with their length of stay in the United States. Hip BMD was significantly correlated with milk intake (p < 0.02) and marginally (p = 0.052) with their length of stay in the United States, independent of body weight. CONCLUSIONS: Spinal BMD was significantly lower in recent Sudanese immigrants than in African Americans or Caucasians. Their hip and total body BMD was associated with their length of stay in the United States, suggesting a potential environmental factors in the ethnic diversity of BMD.

Absorptiometry, Photon↗

Bone mineral density-affecting genes in Africans.

BACKGROUND: We have recently reported the role of environmental exposure in the ethnic diversity of bone mineral density (BMD). Potential genetic difference has not been adequately assessed. PURPOSE: To determine allele frequencies of BMD-affecting genes and their association with BMD in Africans. METHODS: Allele frequencies at 18 polymorphic sites in 13 genes that affect BMD in Asians and/or Caucasians were determined in 143 recent immigrants (55 men and 88 women, 18-51 years of age) from sub-Saharan Sudan to the United States. Genetic association studies were performed. RESULTS: Among the 14 single-nucleotide polymorphisms (SNPs), 10 were significantly different in allele frequency between Sudanese and Asians, and 10 between Sudanese and Caucasians. Only the osteocalcin gene was not significantly different in allele frequency among Sudanese, Asians and Caucasians. Allele frequencies in the TGFB, COL1A1 and CSR genes were extremely low (<0.04) in the Sudanese. Frequencies of microsatellite alleles in four genes were significantly different among Sudanese, Asians and Caucasians. SNPs in the VDR and ERalpha genes were associated with BMD and/or BMC (bone mineral content) at several bone sites. CONCLUSIONS: Genetic difference may play a role in the ethnic diversity in BMD and/or BMC.

Adolescent↗

Mutation rate variation at human dinucleotide microsatellites.

Mutation is the ultimate source of genetic variation, and mutation rate is thus an important parameter governing the extent of genetic variation. Microsatellites are highly informative genetic markers that have been widely used in genetic studies. While previous studies showed that the mutation rate differs in di-, tri-, and tetranucleotide repeats, how mutation rate distributes within each class of repeat is poorly understood. This study first revealed the pattern of the mutation rate variation within the dinucleotide repeats. Two data sets were used. The first is the allele frequency data from 115 microsatellites with dinucleotide repeats distributed along the human genome in 10 worldwide populations. The second data set is much larger, consisting of the allele frequency of 5252 dinucleotide repeats from the Genome Database. Mutation rate for each locus is estimated through a new homozygosity-based estimator, which has been shown to be unbiased and highly efficient and is reasonably robust against deviations from the single-step model. The mutation rates among loci can be approximated well by a gamma distribution and its shape parameter can be accurately estimated with this approach. This result provides the basic guidelines for analyzing the large-scale genomic data from microsatellite loci.

Alleles↗

Moderate mutation rate in the SARS coronavirus genome and its implications.

BACKGROUND: The outbreak of severe acute respiratory syndrome (SARS) caused a severe global epidemic in 2003 which led to hundreds of deaths and many thousands of hospitalizations. The virus causing SARS was identified as a novel coronavirus (SARS-CoV) and multiple genomic sequences have been revealed since mid-April, 2003. After a quiet summer and fall in 2003, the newly emerged SARS cases in Asia, particularly the latest cases in China, are reinforcing a wide-spread belief that the SARS epidemic would strike back. With the understanding that SARS-CoV might be with humans for years to come, knowledge of the evolutionary mechanism of the SARS-CoV, including its mutation rate and emergence time, is fundamental to battle this deadly pathogen. To date, the speed at which the deadly virus evolved in nature and the elapsed time before it was transmitted to humans remains poorly understood. RESULTS: Sixteen complete genomic sequences with available clinical histories during the SARS outbreak were analyzed. After careful examination of multiple-sequence alignment, 114 single nucleotide variations were identified. To minimize the effects of sequencing errors and additional mutations during the cell culture, three strategies were applied to estimate the mutation rate by 1) using the closely related sequences as background controls; 2) adjusting the divergence time for cell culture; or 3) using the common variants only. The mutation rate in the SARS-CoV genome was estimated to be 0.80 - 2.38 x 10-3 nucleotide substitution per site per year which is in the same order of magnitude as other RNA viruses. The non-synonymous and synonymous substitution rates were estimated to be 1.16 - 3.30 x 10-3 and 1.67 - 4.67 x 10-3 per site per year, respectively. The most recent common ancestor of the 16 sequences was inferred to be present as early as the spring of 2002. CONCLUSIONS: The estimated mutation rates in the SARS-CoV using multiple strategies were not unusual among coronaviruses and moderate compared to those in other RNA viruses. All estimates of mutation rates led to the inference that the SARS-CoV could have been with humans in the spring of 2002 without causing a severe epidemic.

Coronavirus↗

DNA sequence duplication in Rhodobacter sphaeroides 2.4.1: evidence of an ancient partnership between chromosomes I and II.

The complex genome of Rhodobacter sphaeroides 2.4.1, composed of chromosomes I (CI) and II (CII), has been sequenced and assembled. We present data demonstrating that the R. sphaeroides genome possesses an extensive amount of exact DNA sequence duplication, 111 kb or approximately 2.7% of the total chromosomal DNA. The chromosomal DNA sequence duplications were aligned to each other by using MUMmer. Frequency and size distribution analyses of the exact DNA duplications revealed that the interchromosomal duplications occurred prior to the intrachromosomal duplications. Most of the DNA sequence duplications in the R. sphaeroides genome occurred early in species history, whereas more recent sequence duplications are rarely found. To uncover the history of gene duplications in the R. sphaeroides genome, 44 gene duplications were sampled and then analyzed for DNA sequence similarity against orthologous DNA sequences. Phylogenetic analysis revealed that approximately 80% of the total gene duplications examined displayed type A phylogenetic relationships; i.e., one copy of each member of a duplicate pair was more similar to its orthologue, found in a species closely related to R. sphaeroides, than to its duplicate, counterpart allele. The data reported here demonstrate that a massive level of gene duplications occurred prior to the origin of the R. sphaeroides 2.4.1 lineage. These findings lead to the conclusion that there is an ancient partnership between CI and CII of R. sphaeroides 2.4.1.

Bacterial Proteins↗