Key genetic variants associated with variation of milk oligosaccharides from diverse human populations

Human milk oligosaccharides (HMO), the third most abundant component of human milk, are thought to be important contributors to infant health. Studies have provided evidence that geography, stage of lactation, and Lewis and secretor blood groups are associated with HMO profile. However, little is known about how variation across the genome may influence HMO composition among women in various populations. In this study, we performed genome-wide association analyses of 395 women from 8 countries to identify genetic regions associated with 19 different HMO. Our data support FUT2 as the most significantly associated (P < 4.23⁻⁹ to P < 4.5⁻⁷⁰) gene with seven HMO and provide evidence of balancing selection for FUT2. Although polymorphisms in FUT3 were also associated with variation in lacto-N-fucopentaose II and difucosyllacto-N-tetrose, we found little evidence of selection on FUT3. To our knowledge, this is the first report of the use of genome-wide association analyses on HMO.     

Virus-derived variation in diverse human genomes

Acquisition of genetic material from viruses by their hosts can generate inter-host structural genome variation. We developed computational tools enabling us to study virus-derived structural variants (SVs) in population-scale whole genome sequencing (WGS) datasets and applied them to 3,332 humans. Although SVs had already been cataloged in these subjects, we found previously-overlooked virus-derived SVs. We detected non-germline SVs derived from squirrel monkey retrovirus (SMRV), human immunodeficiency virus 1 (HIV-1), and human T lymphotropic virus (HTLV-1); these variants are attributable to infection of the sequenced lymphoblastoid cell lines (LCLs) or their progenitor cells and may impact gene expression results and the biosafety of experiments using these cells. In addition, we detected new heritable SVs derived from human herpesvirus 6 (HHV-6) and human endogenous retrovirus-K (HERV-K). We report the first solo-direct repeat (DR) HHV-6 likely to reflect DR rearrangement of a known full-length endogenous HHV-6. We used linkage disequilibrium between single nucleotide variants (SNVs) and variants in reads that align to HERV-K, which often cannot be mapped uniquely using conventional short-read sequencing analysis methods, to locate previously-unknown polymorphic HERV-K loci. Some of these loci are tightly linked to trait-associated SNVs, some are in complex genome regions inaccessible by prior methods, and some contain novel HERV-K haplotypes likely derived from gene conversion from an unknown source or introgression. These tools and results broaden our perspective on the coevolution between viruses and humans, including ongoing virus-to-human gene transfer contributing to genetic variation between humans.     

Patterns of genetic variation in anthropogenically impacted populations

Genetic variation is considered critical for allowing natural populations to adapt to their changing environment, and yet the effects of human disturbance on genetic variation in the wild are poorly understood. Different types of human disturbances may genetically impact natural populations in a predictable manner and so the aim of this study was to provide an overview of these changes using a quantitative literature review approach. I examined both allozyme and microsatellite estimates of genetic variation from peer-reviewed journals, using the mean number of alleles per locus and expected heterozygosity as standardized metrics. Populations within each study were categorized according to the type of human disturbance experienced (“hunting/harvest”, “habitat fragmentation”, or “pollution”), and taxon-specific, as well as time- and context-dependent disturbance effects were considered. I found that human disturbances are associated with weak, but consistent changes in neutral genetic variation within natural populations. The direction of change was dependent on the type of human disturbance experienced, with some forms of anthropogenic challenges consistently decreasing genetic variation from background patterns (e.g., habitat fragmentation), whereas others had no effect (e.g., hunting/harvest) or even slightly increased genetic variation (e.g., pollution). These same measures appeared sensitive to both the time of origin and duration of the disturbance as well. This suggests that the presence or absence, strength, type, as well as the spatial and temporal scale of human disturbance experienced may warrant careful consideration when conservation management plans are formulated for natural populations, with particular attention paid to the effects of habitat fragmentation.     

Analysis of the variation in chromosome size among diverse human populations by bivariate flow karyotyping

Chromosomes sampled from seven human populations were analyzed by flow cytometry to survey normal variation in chromosome size. The populations include two African Pygmy groups, two Amerindian tribes, Druze, Khmer Cambodians, and Melanesians. Mitotic chromosomes were isolated from cultured cells and stained with Hoechst 33 258 and chromomycin A3. The relative DNA content and base-pair composition of each homolog was quantified by bivariate flow karyotyping. Significant variation in DNA content, ranging from 10–40%, was observed for chromosomes 1, 13–16, 19, 21, 22, and Y. The measurements for each population appeared to be a random sampling of the total set of 33 individuals for the majority of chromosomes. A few significant differences in the distributions of chromosomal DNA content were observed among the populations, however. The data, when combined with an earlier study of 33 unrelated individuals of unknown ethnic origin, provide a good representation of the variation in chromosome size among humans. Received: 3 September 1996 / Accepted: 10 January 1997     

Patterns of genetic variation in populations of infectious agents

Background  

The analysis of genetic variation in populations of infectious agents may help us understand their epidemiology and evolution. Here we study a model for assessing the levels and patterns of genetic diversity in populations of infectious agents. The population is structured into many small subpopulations, which correspond to their hosts, that are connected according to a specific type of contact network. We considered different types of networks, including fully connected networks and scale free networks, which have been considered as a model that captures some properties of real contact networks. Infectious agents transmit between hosts, through migration, where they grow and mutate until elimination by the host immune system.     

Bayesian association of haplotypes and non-genetic factors to regulatory and phenotypic variation in human populations

MOTIVATION: With the recent availability of large-scale data sets profiling single nucleotide polymorphisms (SNPs) and quantitative traits data across different human subpopulations, there has been much attention directed towards discovering patterns of genetic variation and their connection to gene regulation and the onset/progression of disease. While previous work has focused primarily on correlating individual SNP markers with gene expression and disease, it has been suggested that using haplotype blocks instead of individual markers can significantly increase statistical power. RESULTS: We present BlockMapper, a probabilistic generative model for genotype data and quantitative traits data, such as gene expression or phenotype measurements. BlockMapper discovers the block structure of genotype data and associates these inferred blocks to patterns of variation in quantitative traits data, whilst accounting for non-genetic factors. Our model achieves high accuracy for predicting Crohn's disease phenotype in Chromosome 5q31 and reveals novel cis-associations between two haplotype blocks in the ENm006 genomic region and GDI1, a gene implicated in X-linked mental retardation. Our results underscore the importance of accounting for the influence of large sets of SNPs on patterns of regulatory/phenotypic variation and represent a step towards an understanding of human genetic variation.     

Patterns of molecular genetic variation among African elephant populations

The highly threatened African elephants have recently been subdivided into two species, Loxodonta africana (savannah or bush elephant) and L. cyclotis (forest elephant) based on morphological and molecular studies. A molecular genetic assessment of 16 microsatellite loci across 20 populations (189 individuals) affirms species level genetic differentiation and provides robust genotypic assessment of species affiliation. Savannah elephant populations show modest levels of phylogeographic subdivision based on composite microsatellite genotype, an indication of recent population isolation and restricted gene flow between locales. The savannah elephants show significantly lower genetic diversity than forest elephants, probably reflecting a founder effect in the recent history of the savannah species.     

Parallel adaptive divergence among geographically diverse human populations

Few genetic differences between human populations conform to the classic model of positive selection, in which a newly arisen mutation rapidly approaches fixation in one lineage, suggesting that adaptation more commonly occurs via moderate changes in standing variation at many loci. Detecting and characterizing this type of complex selection requires integrating individually ambiguous signatures across genomically and geographically extensive data. Here, we develop a novel approach to test the hypothesis that selection has favored modest divergence at particular loci multiple times in independent human populations. We find an excess of SNPs showing non-neutral parallel divergence, enriched for genic and nonsynonymous polymorphisms in genes encompassing diverse and often disease related functions. Repeated parallel evolution in the same direction suggests common selective pressures in disparate habitats. We test our method with extensive coalescent simulations and show that it is robust to a wide range of demographic events. Our results demonstrate phylogenetically orthogonal patterns of local adaptation caused by subtle shifts at many widespread polymorphisms that likely underlie substantial phenotypic diversity.     

Craniometric variation among modern human populations

Previous studies of genetic markers and mitochondrial DNA have found that the amount of variation among major geographic groupings of Homo sapiens is relatively low, accounting for roughly 10% of total variation. This conclusion has had implications for the study of human variation and consideration of alternative models for the origin of modern humans. By contrast, it has often been assumed that the level of among-group variation for morphological traits is much higher. This study examines the level of among-group variation based on Craniometric data from a large sample of modern humans originally collected by W. W. Howells. A multivariate method based on quantitative genetics theory was used to provide an estimate of F_ST — a measure of among-group variation that can be compared with results from studies of genetic markers. Data for 57 Craniometric variables on 1,734 crania were analyzed. These data represent six core areas: Europe, Sub-Saharan Africa, Australasia, Polynesia, the Americas, and the Far East. An additional set of analyses was performed using a three-region subset (Europe, Sub-Saharan Africa, and the Far East) to provide comparability with several genetic studies. The minimum F_ST (assuming complete heritability) for the three-region analysis is 0.065, and the minimum F_ST for the six-region analysis is 0.085. Both of these are less than the average F_ST from genetic studies (average estimates of 0.10–0.11). The smaller value of the minimum F_ST estimates is expected since it provides an estimate of F_ST expected under complete heritability. Using an estimate of average craniometric heritability from the literature provides an estimate of F_ST of 0.112 for the three-region analysis and 0.144 for the six-region analysis. These results show that genetic and craniometric data are in agreement, qualitatively and quantitatively, and that there is limited variation in modern humans among major geographic regions. © 1994 Wiley-Liss, Inc.     

Sequence variation in ZFX introns in human populations

DNA variation in human populations was studied by examining the last intron of the ZFX gene (about 1, 151 bp) with a worldwide sample of 29 individuals. Only one polymorphic site was found, which is located in an Alu sequence. This polymorphism is present at an intermediate frequency in all populations studied, and could be a shared polymorphism or due to migration among populations in Asia, Europe, and Africa. The nucleotide diversity is 0.04%, supporting the view that the level of nucleotide variation in nuclear DNA is very low in humans. From the sequence data, the age (T) of the most recent common ancestor of the sampled sequences is estimated: the mode of T is about 306,000 years, and the 95% confidence interval of T is 162,000-952,000 years. This mode estimate is considerably older than the estimates from Y- linked sequences.      

Patterns of ribosomal gene variation in elite commercial chicken pure line populations

The nucleolus organizer region (NOR) encodes the tandemly repeated 18S, 5.8S and 28S ribosomal (r) RNA genes. The NORs of broiler and layer commercial chicken pure lines were studied to establish the type and extent of genetic variation at this important locus. The parameters studied were gene copy number, repeat size, and diversity of NOR-types. The populations were organized into three groups for analysis including brown-egg broiler (13 lines), brown-egg layer (six lines), and white-egg layer (eight lines). The ribosomal gene copy number average of the white-egg layer populations was significantly lower (329 genes) than that of the brown-egg layers (372 genes); the brown-egg broiler ribosomal gene average was intermediate (350 genes). The white-egg layer populations exhibited a ribosomal repeat unit average size of 36 kb, significantly different from the brown-egg layer and brown-egg broiler average repeat unit size of 32.5 and 33.9 kb, respectively. NOR array size was similar among the three groups (6 mb). The brown-egg broiler populations exhibited polymorphic NOR patterns, intra- and interline, whereas the white-egg layer populations were essentially monomorphic for NOR-type; brown-egg layers exhibited an intermediate level of NOR diversity. Some NOR array characteristics may be a function of breed origin as brown-egg commercial populations, both broilers and layers, have similar breed origins and exhibited similarities for predominant repeat unit size as compared with white-egg layer populations. However, the finding that brown-egg broiler lines typically exhibit a greater number of segregating NOR-types than brown-egg layer lines suggests that the selection schemes of broiler vs. layer pure line populations may also have influenced the degree of variation at this gene complex.     

Metric dental variation of major human populations

Mesiodistal and buccolingual crown diameters of all teeth recorded in 72 major human population groups and seven geographic groups were analyzed. The results obtained are fivefold. First, the largest teeth are found among Australians, followed by Melanesians, Micronesians, sub-Saharan Africans, and Native Americans. Philippine Negritos, Jomon/Ainu, and Western Eurasians have small teeth, while East/Southeast Asians and Polynesians are intermediate in overall tooth size. Second, in terms of odontometric shape factors, world extremes are Europeans, aboriginal New World populations, and to a lesser extent, Australians. Third, East/Southeast Asians share similar dental features with sub-Saharan Africans, and fall in the center of the phenetic space occupied by a wide array of samples. Fourth, the patterning of dental variation among major geographic populations is more or less consistent with those obtained from genetic and craniometric data. Fifth, once differences in population size between sub-Saharan Africa, Europe, South/West Asia, Australia, and Far East, and genetic drift are taken into consideration, the pattern of sub-Saharan African distinctiveness becomes more or less comparable to that based on genetic and craniometric data. As such, worldwide patterning of odontometric variation provides an additional avenue in the ongoing investigation of the origin(s) of anatomically modern humans.     

Identification of copy number variation hotspots in human populations

Copy number variants (CNVs) in the human genome contribute to both Mendelian and complex traits as well as to genomic plasticity in evolution. The investigation of mutational rates of CNVs is critical to understanding genomic instability and the etiology of the copy number variation (CNV)-related traits. However, the evaluation of the CNV mutation rate at the genome level poses an insurmountable practical challenge that requires large samples and accurate typing. In this study, we show that an approximate estimation of the CNV mutation rate could be achieved by using the phylogeny information of flanking SNPs. This allows a genome-wide comparison of mutation rates between CNVs with the use of vast, readily available data of SNP genotyping. A total of 4187 CNV regions (CNVRs) previously identified in HapMap populations were investigated in this study. We showed that the mutation rates for the majority of these CNVRs are at the order of 10⁻⁵ per generation, consistent with experimental observations at individual loci. Notably, the mutation rates of 104 (2.5%) CNVRs were estimated at the order of 10⁻³ per generation; therefore, they were identified as potential hotspots. Additional analyses revealed that genome architecture at CNV loci has a potential role in inciting mutational hotspots in the human genome. Interestingly, 49 (47%) CNV hotspots include human genes, some of which are known to be functional CNV loci (e.g., CNVs of C4 and β-defensin causing autoimmune diseases and CNVs of HYDIN with implication in control of cerebral cortex size), implicating the important role of CNV in human health and evolution, especially in common and complex diseases.     

Genetic variation of recent Alu insertions in human populations

The Alu family of intersperesed repeats is comprised of ovr 500,000 members which may be divided into discrete subfamilies based upon mutations held in common between members. Distinct subfamilies of Alu sequences have amplified within the human genome in recent evolutionary history. Several individual Alu family members have amplified so recently in human evolution that they are variable as to presence and absence at specific loci within different human populations. Here, we report on the distribution of six polymorphic Alu insetions in a survey of 563 individuals from 14 human population groups across several continents. Our results indicate that these polymorphic Alu insertions probably have an African origin and that there is a much smaller amount of genetic variation between European populations than that found between other populations groups. Present address: Department of Pathology, Stanley S. Scott Cancer Center, Louisiana State University Medical Center, 1901 Perdido St., New Orleans, LA 70112 Correspondence to: M.A. Batzer     

Gene-expression variation within and among human populations

Understanding patterns of gene-expression variation within and among human populations will provide important insights into the molecular basis of phenotypic diversity and the interpretation of patterns of expression variation in disease. However, little is known about how gene-expression variation is apportioned within and among human populations. Here, we characterize patterns of natural gene-expression variation in 16 individuals of European and African ancestry. We find extensive variation in gene-expression levels and estimate that approximately 83% of genes are differentially expressed among individuals and that approximately 17% of genes are differentially expressed among populations. By decomposing total gene-expression variation into within- versus among-population components, we find that most expression variation is due to variation among individuals rather than among populations, which parallels observations of extant patterns of human genetic variation. Finally, we performed allele-specific quantitative polymerase chain reaction to demonstrate that cis-regulatory variation in the lymphocyte adaptor protein (SH2B adapter protein 3) contributes to differential expression between European and African samples. These results provide the first insight into how human population structure manifests itself in gene-expression levels and will help guide the search for regulatory quantitative trait loci.     

Patterns of ecological specialization among microbial populations in the Red Sea and diverse oligotrophic marine environments

Large swaths of the nutrient‐poor surface ocean are dominated numerically by cyanobacteria (Prochlorococcus), cyanobacterial viruses (cyanophage), and alphaproteobacteria (SAR11). How these groups thrive in the diverse physicochemical environments of different oceanic regions remains poorly understood. Comparative metagenomics can reveal adaptive responses linked to ecosystem‐specific selective pressures. The Red Sea is well‐suited for studying adaptation of pelagic‐microbes, with salinities, temperatures, and light levels at the extreme end for the surface ocean, and low nutrient concentrations, yet no metagenomic studies have been done there. The Red Sea (high salinity, high light, low N and P) compares favorably with the Mediterranean Sea (high salinity, low P), Sargasso Sea (low P), and North Pacific Subtropical Gyre (high light, low N). We quantified the relative abundance of genetic functions among Prochlorococcus, cyanophage, and SAR11 from these four regions. Gene frequencies indicate selection for phosphorus acquisition (Mediterranean/Sargasso), DNA repair and high‐light responses (Red Sea/Pacific Prochlorococcus), and osmolyte C1 oxidation (Red Sea/Mediterranean SAR11). The unexpected connection between salinity‐dependent osmolyte production and SAR11 C1 metabolism represents a potentially major coevolutionary adaptation and biogeochemical flux. Among Prochlorococcus and cyanophage, genes enriched in specific environments had ecotype distributions similar to nonenriched genes, suggesting that inter‐ecotype gene transfer is not a major source of environment‐specific adaptation. Clustering of metagenomes using gene frequencies shows similarities in populations (Red Sea with Pacific, Mediterranean with Sargasso) that belie their geographic distances. Taken together, the genetic functions enriched in specific environments indicate competitive strategies for maintaining carrying capacity in the face of physical stressors and low nutrient availability.