Ethnicity and human genetic linkage maps

2023 lifetime achievement award: "If you want to go fast, go alone; if you want to go far, go together"

This article is based on the address given by the author at the 2023 meeting of The American Society of Human Genetics (ASHG). A video of the original address can be found at the ASHG website.

Genetic Variation, GWAS and Accuracy of Prediction for Host Resistance to Sparicotyle chrysophrii in Farmed Gilthead Sea Bream (Sparus aurata)

Gilthead sea bream (Sparus aurata) belongs to a group of teleost which has high importance in Mediterranean aquaculture industry. However, industrial production is increasingly compromised by an elevated outbreak of diseases in sea cages, especially a disease caused by monogeneans parasite Sparicotyle chrysophrii. This parasite mainly colonizes gill tissues of host and causes considerable economical losses with mortality and reduction in growth. The aim of current study was to explore the genetics of host resistance against S. chrysophrii and investigate the potential for genomic selection to possibly accelerate genetic progress. To achieve the desired goals, a test population derived from the breeding nucleus of Andromeda Group was produced. This experimental population was established by crossing of parents mated in partial factorial crosses of ∼8 × 8 using 58 sires and 62 dams. The progeny obtained from this mating design was challenged with S. chrysophrii using a controllable cohabitation infection model. At the end of the challenge, fish were recorded for parasite count, and all the recorded fish were tissue sampled for genotyping by sequencing using 2b-RAD methodology. The initial (before challenge test) and the final body weight (after challenge test) of the fish were also recorded. The results obtained through the analysis of phenotypic records (n = 615) and the genotypic data (n = 841, 724 offspring and 117 parents) revealed that the resistance against this parasite is lowly heritable (h² = 0.147 with pedigree and 0.137 with genomic information). We observed moderately favorable genetic correlation (R_g = −0.549 to −0.807) between production traits (i.e., body weight and specific growth rate) and parasite count, which signals a possibility of indirect selection. A locus at linkage group 17 was identified that surpassed chromosome-wide Bonferroni threshold which explained 22.68% of the total genetic variance, and might be playing role in producing genetic variation. The accuracy of prediction was improved by 8% with genomic information compared to pedigree.

Gene-dense autosomal chromosomes show evidence for increased selection

Purifying selection tends to reduce nucleotide and haplotype diversity leading to increased linkage disequilibrium. However, detection of evidence for selection is difficult as the signature is confounded by wide variation in the recombination rate which has a complex relationship with selection. The effective bottleneck time (the ratio of the linkage disequilibrium map to the genetic map in Morgans) controls for variability in the recombination rate. Reduced effective bottleneck times indicate stronger residual linkage disequilibrium, consistent with increased selection. Using whole genome sequence data from one European and three Sub-Saharan African human populations we find, in the African samples, strong correlations between high gene densities and reduced effective bottleneck time for autosomal chromosomes. This suggests that gene-dense autosomes have been subject to increased purifying selection reducing effective bottleneck times compared to gene-poor autosomes. Although previous studies have shown unusually strong linkage disequilibrium for the sex chromosomes variation within the autosomes has not been recognised. The strongest relationship is between effective bottleneck time and the density of essential genes, which are likely targets of greater selective pressure (p = 0.006, for the 22 autosomes). The magnitude of the reduction in chromosome-specific effective bottleneck times from the least to the most gene-dense autosomes is ~17-21% for Sub-Saharan African populations. The effect size is greater in Sub-Saharan African populations, compared to a European sample, consistent with increased efficiency of selection in populations with larger effective population sizes which have not been subject to intense population bottlenecks as experienced by populations of European ancestry. The findings highlight the value of deeper analyses of selection within Sub-Saharan African populations.

Genetic Determinants of Telomere Length in African American Youth

Telomere length (TL) is associated with numerous disease states and is affected by genetic and environmental factors. However, TL has been mostly studied in adult populations of European or Asian ancestry. These studies have identified 34 TL-associated genetic variants recently used as genetic proxies for TL. The generalizability of these associations to pediatric populations and racially diverse populations, specifically of African ancestry, remains unclear. Furthermore, six novel variants associated with TL in a population of European children have been identified but not validated. We measured TL from whole blood samples of 492 healthy African American youth (children and adolescents between 8 and 20 years old) and performed the first genome-wide association study of TL in this population. We were unable to replicate neither the 34 reported genetic associations found in adults nor the six genetic associations found in European children. However, we discovered a novel genome-wide significant association between TL and rs1483898 on chromosome 14. Our results underscore the importance of examining genetic associations with TL in diverse pediatric populations such as African Americans.

Racial and Socioeconomic Variation in Genetic Markers of Telomere Length: A Cross-Sectional Study of U.S. Older Adults

Background

Shorter telomere length (TL) has been associated with stress and adverse socioeconomic conditions, yet U.S. blacks have longer TL than whites. The role of genetic versus environmental factors in explaining TL by race and socioeconomic position (SEP) remains unclear.

Methods

We used data from the U.S. Health and Retirement Study (N = 11,934) to test the hypothesis that there are differences in TL-associated SNPs by race and SEP. We constructed a TL polygenic risk score (PRS) and examined its association with race/ethnicity, educational attainment, assets, gender, and age.

Results

U.S. blacks were more likely to have a lower PRS for TL, as were older individuals and men. Racial differences in TL were statistically accounted for when controlling for population structure using genetic principal components. The GWAS-derived SNPs for TL, however, may not have consistent associations with TL across different racial/ethnic groups.

Conclusions

This study showed that associations of race/ethnicity with TL differed when accounting for population stratification. The role of race/ethnicity for TL remains uncertain, however, as the genetic determinants of TL may differ by race/ethnicity. Future GWAS samples should include racially diverse participants to allow for better characterization of the determinants of TL in human populations.

•

Blacks, older individuals, and men are more likely to have a polygenic risk score predisposing them to longer telomeres.

•

There is no association between telomere length and race/ethnicity after controlling for population structure.

•

GWAS studies have not included diverse samples, and genetic associations with telomere length may differ by race/ethnicity.

Telomeres are structures that protect the ends of chromosomes from damage. Shorter telomeres may be a marker of human aging. Shorter telomeres have been associated with higher stress and lower levels of education, but U.S. blacks have longer telomeres than whites. We show that blacks, older individuals, and men have genetic markers that may predispose them to longer telomeres. After accounting for genetic population structure, there is no longer an association between telomere length and self-reported race/ethnicity. Because genetic determinants of telomere length may differ by race/ethnicity, it is critical that future genetic studies include racially/ethnically diverse populations.

Correlations between Synaptic Initiation and Meiotic Recombination: A Study of Humans and Mice

Meiotic recombination is initiated by programmed double strand breaks (DSBs), only a small subset of which are resolved into crossovers (COs). The mechanism determining the location of these COs is not well understood. Studies in plants, fungi, and insects indicate that the same genomic regions are involved in synaptic initiation and COs, suggesting that early homolog alignment is correlated with the eventual resolution of DSBs as COs. It is generally assumed that this relationship extends to mammals, but little effort has been made to test this idea. Accordingly, we conducted an analysis of synaptic initiation sites (SISs) and COs in human and mouse spermatocytes and oocytes. In contrast to our expectation, we observed remarkable sex- and species-specific differences, including pronounced differences between human males and females in both the number and chromosomal location of SISs. Further, the combined data from our studies in mice and humans suggest that the relationship between SISs and COs in mammals is a complex one that is not dictated by the sites of synaptic initiation as reported in other organisms, although it is clearly influenced by them.

Comparative mapping between Coho Salmon (Oncorhynchus kisutch) and three other salmonids suggests a role for chromosomal rearrangements in the retention of duplicated regions following a whole genome duplication event

Whole genome duplication has been implicated in evolutionary innovation and rapid diversification. In salmonid fishes, however, whole genome duplication significantly pre-dates major transitions across the family, and re-diploidization has been a gradual process between genomes that have remained essentially collinear. Nevertheless, pairs of duplicated chromosome arms have diverged at different rates from each other, suggesting that the retention of duplicated regions through occasional pairing between homeologous chromosomes may have played an evolutionary role across species pairs. Extensive chromosomal arm rearrangements have been a key mechanism involved in re-dipliodization of the salmonid genome; therefore, we investigated their influence on degree of differentiation between homeologs across salmon species. We derived a linkage map for coho salmon and performed comparative mapping across syntenic arms within the genus Oncorhynchus, and with the genus Salmo, to determine the phylogenetic relationship between chromosome arrangements and the retention of undifferentiated duplicated regions. A 6596.7 cM female coho salmon map, comprising 30 linkage groups with 7415 and 1266 nonduplicated and duplicated loci, respectively, revealed uneven distribution of duplicated loci along and between chromosome arms. These duplicated regions were conserved across syntenic arms across Oncorhynchus species and were identified in metacentric chromosomes likely formed ancestrally to the divergence of Oncorhynchus from Salmo. These findings support previous studies in which observed pairings involved at least one metacentric chromosome. Re-diploidization in salmon may have been prevented or retarded by the formation of metacentric chromosomes after the whole genome duplication event and may explain lineage-specific innovations in salmon species if functional genes are found in these regions.

Examining variation in recombination levels in the human female: a test of the production-line hypothesis

The most important risk factor for human aneuploidy is increasing maternal age, but the basis of this association remains unknown. Indeed, one of the earliest models of the maternal-age effect--the "production-line model" proposed by Henderson and Edwards in 1968--remains one of the most-cited explanations. The model has two key components: (1) that the first oocytes to enter meiosis are the first ovulated and (2) that the first to enter meiosis have more recombination events (crossovers) than those that enter meiosis later in fetal life. Studies in rodents have demonstrated that the first oocytes to enter meiosis are indeed the first to be ovulated, but the association between the timing of meiotic entry and recombination levels has not been tested. We recently initiated molecular cytogenetic studies of second-trimester human fetal ovaries, allowing us to directly examine the number and distribution of crossover-associated proteins in prophase-stage oocytes. Our observations on over 8,000 oocytes from 191 ovarian samples demonstrate extraordinary variation in recombination within and among individuals but provide no evidence of a difference in recombination levels between oocytes entering meiosis early in fetal life and those entering late in fetal life. Thus, our data provide a direct test of the second tenet of the production-line model and suggest that it does not provide a plausible explanation for the human maternal-age effect, meaning that-45 years after its introduction-we can finally conclude that the production-line model is not the basis for the maternal-age effect on trisomy.

A dense linkage map for Chinook salmon (Oncorhynchus tshawytscha) reveals variable chromosomal divergence after an ancestral whole genome duplication event

Comparisons between the genomes of salmon species reveal that they underwent extensive chromosomal rearrangements following whole genome duplication that occurred in their lineage 58−63 million years ago. Extant salmonids are diploid, but occasional pairing between homeologous chromosomes exists in males. The consequences of re-diploidization can be characterized by mapping the position of duplicated loci in such species. Linkage maps are also a valuable tool for genome-wide applications such as genome-wide association studies, quantitative trait loci mapping or genome scans. Here, we investigated chromosomal evolution in Chinook salmon (Oncorhynchus tshawytscha) after genome duplication by mapping 7146 restriction-site associated DNA loci in gynogenetic haploid, gynogenetic diploid, and diploid crosses. In the process, we developed a reference database of restriction-site associated DNA loci for Chinook salmon comprising 48528 non-duplicated loci and 6409 known duplicated loci, which will facilitate locus identification and data sharing. We created a very dense linkage map anchored to all 34 chromosomes for the species, and all arms were identified through centromere mapping. The map positions of 799 duplicated loci revealed that homeologous pairs have diverged at different rates following whole genome duplication, and that degree of differentiation along arms was variable. Many of the homeologous pairs with high numbers of duplicated markers appear conserved with other salmon species, suggesting that retention of conserved homeologous pairing in some arms preceded species divergence. As chromosome arms are highly conserved across species, the major resources developed for Chinook salmon in this study are also relevant for other related species.

Cytological studies of human meiosis: sex-specific differences in recombination originate at, or prior to, establishment of double-strand breaks

Meiotic recombination is sexually dimorphic in most mammalian species, including humans, but the basis for the male:female differences remains unclear. In the present study, we used cytological methodology to directly compare recombination levels between human males and females, and to examine possible sex-specific differences in upstream events of double-strand break (DSB) formation and synaptic initiation. Specifically, we utilized the DNA mismatch repair protein MLH1 as a marker of recombination events, the RecA homologue RAD51 as a surrogate for DSBs, and the synaptonemal complex proteins SYCP3 and/or SYCP1 to examine synapsis between homologs. Consistent with linkage studies, genome-wide recombination levels were higher in females than in males, and the placement of exchanges varied between the sexes. Subsequent analyses of DSBs and synaptic initiation sites indicated similar male:female differences, providing strong evidence that sex-specific differences in recombination rates are established at or before the formation of meiotic DSBs. We then asked whether these differences might be linked to variation in the organization of the meiotic axis and/or axis-associated DNA and, indeed, we observed striking male:female differences in synaptonemal complex (SC) length and DNA loop size. Taken together, our observations suggest that sex specific differences in recombination in humans may derive from chromatin differences established prior to the onset of the recombination pathway.

The origin, global distribution, and functional impact of the human 8p23 inversion polymorphism

Genomic inversions are an increasingly recognized source of genetic variation. However, a lack of reliable high-throughput genotyping assays for these structures has precluded a full understanding of an inversion's phylogenetic, phenotypic, and population genetic properties. We characterize these properties for one of the largest polymorphic inversions in man (the ∼4.5-Mb 8p23.1 inversion), a structure that encompasses numerous signals of natural selection and disease association. We developed and validated a flexible bioinformatics tool that utilizes SNP data to enable accurate, high-throughput genotyping of the 8p23.1 inversion. This tool was applied retrospectively to diverse genome-wide data sets, revealing significant population stratification that largely follows a clinal “serial founder effect” distribution model. Phylogenetic analyses establish the inversion's ancestral origin within the Homo lineage, indicating that 8p23.1 inversion has occurred independently in the Pan lineage. The human inversion breakpoint was localized to an inverted pair of human endogenous retrovirus elements within the large, flanking low-copy repeats; experimental validation of this breakpoint confirmed these elements as the likely intermediary substrates that sponsored inversion formation. In five data sets, mRNA levels of disease-associated genes were robustly associated with inversion genotype. Moreover, a haplotype associated with systemic lupus erythematosus was restricted to the derived inversion state. We conclude that the 8p23.1 inversion is an evolutionarily dynamic structure that can now be accommodated into the understanding of human genetic and phenotypic diversity.

The landscape of recombination in African Americans

Recombination, together with mutation, gives rise to genetic variation in populations. Here we leverage the recent mixture of people of African and European ancestry in the Americas to build a genetic map measuring the probability of crossing over at each position in the genome, based on about 2.1 million crossovers in 30,000 unrelated African Americans. At intervals of more than three megabases it is nearly identical to a map built in Europeans. At finer scales it differs significantly, and we identify about 2,500 recombination hotspots that are active in people of West African ancestry but nearly inactive in Europeans. The probability of a crossover at these hotspots is almost fully controlled by the alleles an individual carries at PRDM9 (P value < 10(-245)). We identify a 17-base-pair DNA sequence motif that is enriched in these hotspots, and is an excellent match to the predicted binding target of PRDM9 alleles common in West Africans and rare in Europeans. Sites of this motif are predicted to be risk loci for disease-causing genomic rearrangements in individuals carrying these alleles. More generally, this map provides a resource for research in human genetic variation and evolution.

Enhanced genetic maps from family-based disease studies: population-specific comparisons

Background

Accurate genetic maps are required for successful and efficient linkage mapping of disease genes. However, most available genome-wide genetic maps were built using only small collections of pedigrees, and therefore have large sampling errors. A large set of genetic studies genotyped by the NHLBI Mammalian Genotyping Service (MGS) provide appropriate data for generating more accurate maps.

Results

We collected a large sample of uncleaned genotype data for 461 markers generated by the MGS using the Weber screening sets 9 and 10. This collection includes genotypes for over 4,400 pedigrees containing over 17,000 genotyped individuals from different populations. We identified and cleaned numerous relationship and genotyping errors, as well as verified the marker orders. We used this dataset to test for population-specific genetic maps, and to re-estimate the genetic map distances with greater precision; standard errors for all intervals are provided. The map-interval sizes from the European (or European descent), Chinese, and Hispanic samples are in quite good agreement with each other. We found one map interval on chromosome 8p with a statistically significant size difference between the European and Chinese samples, and several map intervals with significant size differences between the African American and Chinese samples. When comparing Palauan with European samples, a statistically significant difference was detected at the telomeric region of chromosome 11p. Several significant differences were also identified between populations in chromosomal and genome lengths.

Conclusions

Our new population-specific screening set maps can be used to improve the accuracy of disease-mapping studies. As a result of the large sample size, the average length of the 95% confidence interval (CI) for a 10 cM map interval is only 2.4 cM, which is considerably smaller than on previously published maps.

Extensive recombination rate variation in the house mouse species complex inferred from genetic linkage maps

The rate of recombination is a key genomic parameter that displays considerable variation among taxa. Species comparisons have demonstrated that the rate of evolution in recombination rate is strongly dependent on the physical scale of measurement. Individual recombination hotspots are poorly conserved among closely related taxa, whereas genomic-scale recombination rate variation bears a strong signature of phylogenetic history. In contrast, the mode and tempo of evolution in recombination rates measured on intermediate physical scales is poorly understood. Here, we conduct a detailed statistical comparison between two whole-genome F₂ genetic linkage maps constructed from experimental intercrosses between closely related house mouse subspecies (Mus musculus). Our two maps profile a common wild-derived inbred strain of M. m. domesticus crossed to distinct wild-derived inbred strains representative of two other house mouse subspecies, M. m. castaneus and M. m. musculus. We identify numerous orthologous genomic regions with significant map length differences between these two crosses. Because the genomes of these recently diverged house mice are highly collinear, observed differences in map length (centimorgans) are suggestive of variation in broadscale recombination rate (centimorgans per megabase) within M. musculus. Collectively, these divergent intervals span 19% of the house mouse genome, disproportionately aggregating on the X chromosome. In addition, we uncover strong statistical evidence for a large effect, sex-linked, site-specific modifier of recombination rate segregating within M. musculus. Our findings reveal considerable variation in the megabase-scale recombination landscape among recently diverged taxa and underscore the continued importance of genetic linkage maps in the post-genome era.

Sex-specific genetic architecture of human fatness in Chinese: the SAPPHIRe Study

To dissect the genetic architecture of sexual dimorphism in obesity-related traits, we evaluated the sex-genotype interaction, sex-specific heritability and genome-wide linkages for seven measurements related to obesity. A total of 1,365 non-diabetic Chinese subjects from the family study of the Stanford Asia-Pacific Program of Hypertension and Insulin Resistance were used to search for quantitative trait loci (QTLs) responsible for the obesity-related traits. Pleiotropy and co-incidence effects from the QTLs were also examined using the bivariate linkage approach. We found that sex-specific differences in heritability and the genotype-sex interaction effects were substantially significant for most of these traits. Several QTLs with strong linkage evidence were identified after incorporating genotype by sex (G × S) interactions into the linkage mapping, including one QTL for hip circumference [maximum LOD score (MLS) = 4.22, empirical p = 0.000033] and two QTLs: for BMI on chromosome 12q with MLS 3.37 (empirical p = 0.0043) and 3.10 (empirical p = 0.0054). Sex-specific analyses demonstrated that these linkage signals all resulted from females rather than males. Most of these QTLs for obesity-related traits replicated the findings in other ethnic groups. Bivariate linkage analyses showed several obesity traits were influenced by a common set of QTLs. All regions with linkage signals were observed in one gender, but not in the whole sample, suggesting the genetic architecture of obesity-related traits does differ by gender. These findings are useful for further identification of the liability genes for these phenotypes through candidate genes or genome-wide association analysis.

Genetic crossovers are predicted accurately by the computed human recombination map

Hotspots of meiotic recombination can change rapidly over time. This instability and the reported high level of inter-individual variation in meiotic recombination puts in question the accuracy of the calculated hotspot map, which is based on the summation of past genetic crossovers. To estimate the accuracy of the computed recombination rate map, we have mapped genetic crossovers to a median resolution of 70 Kb in 10 CEPH pedigrees. We then compared the positions of crossovers with the hotspots computed from HapMap data and performed extensive computer simulations to compare the observed distributions of crossovers with the distributions expected from the calculated recombination rate maps. Here we show that a population-averaged hotspot map computed from linkage disequilibrium data predicts well present-day genetic crossovers. We find that computed hotspot maps accurately estimate both the strength and the position of meiotic hotspots. An in-depth examination of not-predicted crossovers shows that they are preferentially located in regions where hotspots are found in other populations. In summary, we find that by combining several computed population-specific maps we can capture the variation in individual hotspots to generate a hotspot map that can predict almost all present-day genetic crossovers.

In eukaryotes genetic crossovers are responsible for generating genetic diversity and ensuring the proper segregation of chromosomes. Genetic crossovers are tightly clustered in hotspots. Although the existence of hotspots in humans is clearly proven, mechanisms of their formation and the regulation of meiotic recombination in general remain poorly understood. An additional complication in studies of meiotic recombination is the fact that the direct experimental mapping of human hotspots on a genome-wide scale is not feasible with current methods. The best available indirect methods compute the position of hotspots from patterns of historic associations between genetic markers in population samples. In this study we determined the positions of genetic crossovers in ten pedigrees of European origin and then compared the positions of crossovers with the hotspots computed from HapMap data. Importantly, we find that the population-averaged computed map is in close agreement with the observed distribution of genetic crossovers. We also find that cryptic hotspots that are not easily detected in the computed European map can be more effectively identified if other populations are included in the analysis. Our analysis shows that high-resolution recombination profiles are highly similar between distantly related populations and that by including computed hotspots from several populations we can predict nearly all crossovers.

Meiotic recombination in human oocytes

Studies of human trisomies indicate a remarkable relationship between abnormal meiotic recombination and subsequent nondisjunction at maternal meiosis I or II. Specifically, failure to recombine or recombination events located either too near to or too far from the centromere have been linked to the origin of human trisomies. It should be possible to identify these abnormal crossover configurations by using immunofluorescence methodology to directly examine the meiotic recombination process in the human female. Accordingly, we initiated studies of crossover-associated proteins (e.g., MLH1) in human fetal oocytes to analyze their number and distribution on nondisjunction-prone human chromosomes and, more generally, to characterize genome-wide levels of recombination in the human female. Our analyses indicate that the number of MLH1 foci is lower than predicted from genetic linkage analysis, but its localization pattern conforms to that expected for a crossover-associated protein. In studies of individual chromosomes, our observations provide evidence for the presence of "vulnerable" crossover configurations in the fetal oocyte, consistent with the idea that these are subsequently translated into nondisjunctional events in the adult oocyte.

Expression differences by continent of origin point to the immortalization process

Analysis of recently available microarray expression data sets obtained from immortalized cell lines of the individuals represented in the HapMap project have led to inconclusive comparisons across cohorts with different ancestral continent of origin (ACOO). To address this apparent inconsistency, we applied a novel approach to accentuate population-specific gene expression signatures for the CEU [homogeneous US residents with northern and western European ancestry (HapMap samples)] and YRI [homogenous Yoruba people of Ibadan, Nigeria (HapMap samples)] trios. In this report, we describe how four independent data sets point to the differential expression across ACOO of gene networks implicated in transforming the normal lymphoblast into immortalized lymphoblastoid cells. In particular, Werner syndrome helicase and related genes are differentially expressed between the YRI and CEU cohorts. We further demonstrate that these differences correlate with viral titer and that both the titer and expression differences are associated with ACOO. We use the 14 genes most differentially expressed to construct an ACOO-specific ‘immortalization network’ comprised of 40 genes, one of which show significant correlation with genomic variation (eQTL). The extent to which these measured group differences are due to differences in the immortalization procedures used for each group or reflect ACOO-specific biological differences remains to be determined. That the ACOO group differences in gene expression patterns may depend strongly on the process of transforming cells to establish immortalized lines should be considered in such comparisons.

Genome scan of M. tuberculosis infection and disease in Ugandans

Tuberculosis (TB), caused by Mycobacterium tuberculosis (Mtb), is an enduring public health problem globally, particularly in sub-Saharan Africa. Several studies have suggested a role for host genetic susceptibility in increased risk for TB but results across studies have been equivocal. As part of a household contact study of Mtb infection and disease in Kampala, Uganda, we have taken a unique approach to the study of genetic susceptibility to TB, by studying three phenotypes. First, we analyzed culture confirmed TB disease compared to latent Mtb infection (LTBI) or lack of Mtb infection. Second, we analyzed resistance to Mtb infection in the face of continuous exposure, defined by a persistently negative tuberculin skin test (PTST-); this outcome was contrasted to LTBI. Third, we analyzed an intermediate phenotype, tumor necrosis factor-alpha (TNFα) expression in response to soluble Mtb ligands enriched with molecules secreted from Mtb (culture filtrate). We conducted a full microsatellite genome scan, using genotypes generated by the Center for Medical Genetics at Marshfield. Multipoint model-free linkage analysis was conducted using an extension of the Haseman-Elston regression model that includes half sibling pairs, and HIV status was included as a covariate in the model. The analysis included 803 individuals from 193 pedigrees, comprising 258 full sibling pairs and 175 half sibling pairs. Suggestive linkage (p<10⁻³) was observed on chromosomes 2q21-2q24 and 5p13-5q22 for PTST-, and on chromosome 7p22-7p21 for TB; these findings for PTST- are novel and the chromosome 7 region contains the IL6 gene. In addition, we replicated recent linkage findings on chromosome 20q13 for TB (p = 0.002). We also observed linkage at the nominal α = 0.05 threshold to a number of promising candidate genes, SLC11A1 (PTST- p = 0.02), IL-1 complex (TB p = 0.01), IL12BR2 (TNFα p = 0.006), IL12A (TB p = 0.02) and IFNGR2 (TNFα p = 0.002). These results confirm not only that genetic factors influence the interaction between humans and Mtb but more importantly that they differ according to the outcome of that interaction: exposure but no infection, infection without progression to disease, or progression of infection to disease. Many of the genetic factors for each of these stages are part of the innate immune system.

A genome-wide Asian genetic map and ethnic comparison: the GENDISCAN study

BACKGROUND

Genetic maps provide specific positions of genetic markers, which are required for performing genetic studies. Linkage analyses of Asian families have been performed with Caucasian genetic maps, since appropriate genetic maps of Asians were not available. Different ethnic groups may have different recombination rates as a result of genomic variations, which would generate misspecification of the genetic map and reduce the power of linkage analyses.

RESULTS

We constructed the genetic map of a Mongolian population in Asia with CRIMAP software. This new map, called the GENDISCAN map, is based on genotype data collected from 1026 individuals of 73 large Mongolian families, and includes 1790 total and 1500 observable meioses. The GENDISCAN map provides sex-averaged and sex-specific genetic positions of 1039 microsatellite markers in Kosambi centimorgans (cM) with physical positions. We also determined 95% confidence intervals of genetic distances of the adjacent marker intervals. Genetic lengths of the whole genome, chromosomes and adjacent marker intervals are compared with those of Rutgers Map v.2, which was constructed based on Caucasian populations (Centre d'Etudes du Polymorphisme Humain (CEPH) and Icelandic families) by mapping methods identical to those of the GENDISCAN map, CRIMAP software and the Kosambi map function. Mongolians showed approximately 1.9 fewer recombinations per meiosis than Caucasians. As a result, genetic lengths of the whole genome and chromosomes of the GENDISCAN map are shorter than those of Rutgers Map v.2. Thirty-eight marker intervals differed significantly between the Mongolian and Caucasian genetic maps.

CONCLUSION

The new GENDISCAN map is applicable to the genetic study of Asian populations. Differences in the genetic distances between the GENDISCAN and Caucasian maps could facilitate elucidation of genomic variations between different ethnic groups.

Family physicians' beliefs about genetic contributions to racial/ethnic and gender differences in health and clinical decision-making

Greater attention towards genetics as a contributor to group health differences may lead to inappropriate use of race/ethnicity and gender as genetic heuristics and exacerbate health disparities. As part of a web-based survey, 1,035 family physicians (FPs) rated the contribution of genetics and environment to racial/ethnic and gender differences in health outcomes, and the importance of race/ethnicity and gender in their clinical decision-making. FPs attributed racial/ethnic and gender differences in health outcomes equally to environment and genetics. These beliefs were not associated with rated importance of race/ethnicity or gender in clinical decision-making. FPs appreciate the complexity of genetic and environmental influences on health differences by race/ethnicity and gender.

Preeclampsia risk and angiotensinogen polymorphisms M235T and AGT -217 in African American and Caucasian women

INTRODUCTION

Genetic variants of the angiotensinogen gene have been linked to both hypertension and preeclampsia. The M235T polymorphism is more common in hypertension and preeclampsia in some populations. A polymorphism in the angiotensinogen basal promoter region of AGT -217 is more common in African Americans with hypertension. The authors investigated the frequency of M235T and AGT -217 in Caucasian and African American women with and without preeclampsia.

METHODS

The study was a nested case-control study of primiparous women with singleton pregnancies. Genomic DNA from preeclamptic and control subjects underwent polymerase chain reaction amplification and restriction digestion.

RESULTS

The M235T and AGT -217 polymorphisms were both more common in African American women; however, the variants were not more common in preeclampsia.

CONCLUSION

The frequency of angiotensinogen polymorphisms M235T and AGT -217 is different by race; however, these polymorphisms are not associated with an increased risk of preeclampsia.

Genome-wide association studies of cancer

Genome-wide association studies provide a new and powerful approach to investigate the effect of inherited genetic variation on the risk of human disease. These studies rely on high throughput DNA microarray technology to genotype hundreds of thousands of genetic variants across the human genome. The first genome-wide association studies have identified previously unknown genetic risk factors that influence a range of diseases, including prostate cancer, breast cancer, myocardial infarction, age-related macular degeneration, diabetes, Crohn's disease and obesity. Many more studies are currently underway, including a number that will focus on other cancers (e.g., colorectal). Here we discuss the major issues involved in conducting genome-wide association studies and how these studies can be used to examine cancer phenotypes.

Variation in patterns of human meiotic recombination

In the last 30 years it has become evident that patterns of meiotic recombination can be highly variable among individuals. The evidence comes from both low and high resolution analyses of hotspots of recombination in human and other species. In addition, a comparison of the recombination profiles in closely related species such as human and chimpanzee reveals essentially no correlation in the position of hotspots. Although the variation in hotspots of meiotic recombination is clearly documented, the mechanisms responsible for such variation are far from being understood. Here we will review the available evidence of natural variation in meiotic recombination and will discuss potential implications of this variation on the functional mechanisms of crossover formation and control.

The ambiguous meanings of the racial/ethnic categories routinely used in human genetics research

Many researchers are currently studying the distribution of genetic variations among diverse groups, with particular interest in explaining racial/ethnic health disparities. However, the use of racial/ethnic categories as variables in biological research is controversial. Just how racial/ethnic categories are conceptualized, operationalized, and interpreted is a key consideration in determining the legitimacy of their use, but has received little attention. We conducted semi-structured, open-ended interviews with 30 human genetics scientists from the US and Canada who use racial/ethnic variables in their research. They discussed the types of classifications they use, the criteria upon which they are based, and their methods for classifying individual samples and subjects. We found definitions of racial/ethnic variables were often lacking or unclear, the specific categories they used were inconsistent and context specific, and classification practices were often implicit and unexamined. We conclude that such conceptual and practical problems are inherent to routinely used racial/ethnic categories themselves, and that they lack sufficient rigor to be used as key variables in biological research. It is our position that it is unacceptable to persist in the constructing of scientific arguments based on these highly ambiguous variables.

Crossover analysis using immunofluorescent detection of MLH1 foci in frozen–thawed testicular tissue

To date, the effects of freezing on spermatogenesis have not yet been fully investigated at a molecular level. Antibody localization studies have identified the MutL homolog 1 (MLH1) protein, a mis-match repair protein, at the prophase I stage of meiosis, which allows the detection of recombination foci during pachytene. This study investigated the effect of long-term testicular tissue cryopreservation on meiotic prophase I, identified by recombination foci frequency and synaptonemal complex (SC) integrity. Frozen-thawed testicular tissues from 12 males who had each fathered a child were used. Because vasectomy or reverse vasectomy procedures are rare in the locale of the investigation, it was not possible to obtain fresh testicular tissue and use the males as their own controls. Immunocytogenetic analysis of 612 spermatocytes at the pachytene stage was performed. The results indicated a mean number of MLH1 foci of 49.2 (SD +/- 5.9), and no correlation was found between the freezing period, the MLH1 frequency and the SC integrity. The results suggest that freezing of testicular tissue taken post-puberty does not appear to be detrimental to the crossover process as identified by occurrence of MLH1 loci.

Biological research on drug abuse and addiction in Hispanics: current status and future directions

Impressive progress has been made in the understanding of biological contributions to drug abuse and addiction. An area that has only recently begun to receive attention is potential ethnic and racial differences in biological systems that contribute to, or protect from, problem drug use. This article reviews recent research on drug abuse and addiction in Hispanics in which biological questions have been addressed, including work on genes, gene products (proteins), physiology and pharmacotherapy. Taken together, work to date suggests that there are both similarities and differences between Hispanics and other ethnic groups in biological factors related to drug abuse and addiction. Although the results are intriguing, relatively few studies have been done, and those that have been done have often been inconclusive due to low numbers of Hispanic subjects. Moreover, studies have often failed to recognize the complexity and heterogeneity of Hispanic populations in the United States and around the world. After reviewing the current status of the field, recommendations are given for future research in both humans and relevant animal models that will lead to a better understanding of drug abuse and addiction in Hispanics.

Resources for Genetic Variation Studies

The rapid growth of genome-wide diversity databases, as well as ongoing large-scale resequencing projects targeting genes and other functional components of our genome, provide valuable resources of natural variation at the DNA sequence level. In this review, we briefly summarize the wealth of data on DNA polymorphisms in humans, the distribution of this diversity in the genome as well as among individuals, and the consequence of recombination on its organization. These data provide a set of powerful tools that can be used to better understand inherited phenotypic variation in humans. We discuss the implications for the design of studies investigating correlations between genotypes and phenotypes, both at the fundamental level of genome function and regulation, and for the mapping of disease genes.

Accommodating chromosome inversions in linkage analysis

This work develops a population-genetics model for polymorphic chromosome inversions. The model precisely describes how an inversion changes the nature of and approach to linkage equilibrium. The work also describes algorithms and software for allele-frequency estimation and linkage analysis in the presence of an inversion. The linkage algorithms implemented in the software package Mendel estimate recombination parameters and calculate the posterior probability that each pedigree member carries the inversion. Application of Mendel to eight Centre d'Etude du Polymorphisme Humain pedigrees in a region containing a common inversion on 8p23 illustrates its potential for providing more-precise estimates of the location of an unmapped marker or trait gene. Our expanded cytogenetic analysis of these families further identifies inversion carriers and increases the evidence of linkage.

The contribution of individual and pairwise combinations of SNPs in the APOA1 and APOC3 genes to interindividual HDL-C variability

Apolipoproteins (apo) A-I and C-III are components of high-density lipoprotein-cholesterol (HDL-C), a quantitative trait negatively correlated with risk of cardiovascular disease (CVD). We analyzed the contribution of individual and pairwise combinations of single nucleotide polymorphisms (SNPs) in the APOA1/APOC3 genes to HDL-C variability to evaluate (1) consistency of published single-SNP studies with our single-SNP analyses; (2) consistency of single-SNP and two-SNP phenotype-genotype relationships across race-, gender-, and geographical location-dependent contexts; and (3) the contribution of single SNPs and pairs of SNPs to variability beyond that explained by plasma apo A-I concentration. We analyzed 45 SNPs in 3,831 young African-American (N=1,858) and European-American (N=1,973) females and males ascertained by the Coronary Artery Risk Development in Young Adults (CARDIA) study. We found three SNPs that significantly impact HDL-C variability in both the literature and the CARDIA sample. Single-SNP analyses identified only one of five significant HDL-C SNP genotype relationships in the CARDIA study that was consistent across all race-, gender-, and geographical location-dependent contexts. The other four were consistent across geographical locations for a particular race-gender context. The portion of total phenotypic variance explained by single-SNP genotypes and genotypes defined by pairs of SNPs was less than 3%, an amount that is miniscule compared to the contribution explained by variability in plasma apo A-I concentration. Our findings illustrate the impact of context-dependence on SNP selection for prediction of CVD risk factor variability.

Ethnicity may be a reason for lipid changes and high Lp(a) levels in rheumatoid arthritis

There are so many studies that suggest the changes in lipid profiles and lipoprotein (a) [Lp(a)] are associated with early atherosclerosis in rheumatoid arthritis (RA). But there are some opposite studies also. Because of marked ethnicity differences in the distribution of Lp(a), we aimed to investigate the associations of Lp(a) levels and lipid changes in Turkish RA patients. There were 30 women and 20 men, a total of 50 patients with RA (mean age 47.6 +/- 13.2 years), included and 21 healthy women and 14 healthy men (mean age 45.7 +/- 14.5 years) were recruited as a control (C) group. Serum Lp(a), total cholesterol (TC), triglyceride (TG), HDL cholesterol (HDL-C) and LDL cholesterol (LDL-C) levels were analysed for each group. Analysis of six different studies was performed. In the RA and C groups, mean serum Lp(a) levels were 39.7 +/- 64.4 and 10.5 +/- 13.4 mg/dl, respectively (P=0.001). Mean TC levels were 189.2 +/- 142.5 and 174.0 +/- 29.3 mg/dl (P=0.294), mean TG levels were 121.4 +/- 65.4 and 106.5 +/- 80.0 mg/dl (P=0.030), mean HDL-C levels were 44.5 +/- 10.0 and 47.7 +/- 4.8 mg/dl (P=0.014) and mean LDL-C levels were 94.3 +/- 35.3 and 102.0 +/- 24.6 mg/dl (P=0.98), respectively. Analysis of the six studies showed Lp(a) level was higher and HDL level was lower in RA patients than in healthy controls. Patients with RA may have altered lipid profiles from one country to another one. Especially in Turkey, higher serum Lp(a), lower HDL-C and higher TG levels may be found in RA patients instead of some findings of other countries showing different results. Ethnicity may be a reason for these findings.

Large-scale recombination rate patterns are conserved among human populations

In humans, most recombination events occur in a small fraction of the genome. These hotspots of recombination show considerable variation in intensity and/or location across species and, potentially, across human populations. On a larger scale, the patterns of recombination rates have been mostly investigated in individuals of European ancestry, and it remains unknown whether the results obtained can be directly applied to other human populations. Here, we investigate this question using genome-wide polymorphism data. We show that population recombination rates recapitulate a large part of the genetic map information, regardless of the population considered. We also show that the ratio of the population recombination rate estimate of two populations is overall constant along the chromosomes. These two observations support the hypothesis that large-scale recombination patterns are conserved across human populations. Local deviations from the overall pattern of conservation of the recombination rates can be used to select candidate regions with large polymorphic inversions or under local selection.

Coverage and power in genomewide association studies

The ability of genomewide association studies to decipher genetic traits is driven in part by how well the measured single-nucleotide polymorphisms "cover" the unmeasured causal variants. Estimates of coverage based on standard linkage-disequilibrium measures, such as the average maximum squared correlation coefficient (r2), can lead to inaccurate and inflated estimates of the power of genomewide association studies. In contrast, use of the "cumulative r2 adjusted power" measure presented here gives more-accurate estimates of power for genomewide association studies.

Direct detection of null alleles in SNP genotyping data

Pinpointing genetic associations in the human genome relies heavily on the accuracy of the underlying genotype data. Null alleles can generate significant inaccuracies in genotype data and can negatively affect the statistical power of a study. Existing quality control (QC) tests, including tests of Hardy-Weinberg equilibrium, are not sensitive enough to detect the presence of even moderately frequent null alleles in the data. We show that direct analysis of raw data from a quantitative genotyping platform can detect up to 75% of null alleles, even at frequencies below the sensitivity of more traditional methods. Detecting unexpected null alleles not only has benefits in QC of genotype data but may also be valuable in detecting rare, functional null alleles that would otherwise be missed.

Examination of a region showing linkage map discrepancies across sheep breeds

The availability of accurate linkage maps is an important step for the localization of genetic variants of interest. However, most studies in livestock assume the published map is applicable in their population despite the large differences between the breeds of a species. A region of sheep Chromosome 1 was previously identified as providing evidence for a marker order inconsistent with the published linkage map. In this study the identified region was investigated in more detail. Four microsatellite markers covering the central 5 cM of the inconsistent region and two flanking markers were genotyped in three sheep breeds, a commercial population (Charollais), an experimental population (Scottish Blackface), and a feral population (Soay). With the inclusion of the published linkage map, this provided evidence for three different marker orders across four sheep populations. Evidence for selection in this region was investigated using both a single-point allelic competition model and a multipoint allele-sharing statistic. Only the Charollais population provided evidence for selection, with significant transmission bias observed at marker BM7145. The implications of variation in linkage maps on the design and analysis of fine-mapping studies are discussed.

The elusive multiple self-healing squamous epithelioma (MSSE) gene: further mapping, analysis of candidates, and loss of heterozygosity

The MSSE gene predisposes to multiple invasive but self-healing skin tumours (multiple self-healing epitheliomata). MSSE was previously mapped to chromosome 9q22-q31 and a shared haplotype in affected families suggested a founder mutation. We have refined the MSSE critical region (<1 cM, <1 Mb) between the zinc-finger gene ZNF169 and the Fanconi anaemia gene FANCC. By genetic mapping we have excluded ZNF169 and FANCC as well as PTCH (PATCHED) and TGFBR1 (transforming growth factor beta receptor type-1) genes. The CDC14B cell cycle phosphatase gene also lies in the region but screening of the complete coding region revealed no mutation in MSSE patients. Somatic cell hybrids created by haploid conversion of an MSSE patient's cells enabled screening of the MSSE chromosome 9 and showed no CDC14B deletion or mutation that abrogates CDC14B mRNA expression. Thus, CDC14B is unlikely to be the MSSE gene. We also report the first molecular analysis of MSSE tumours showing loss of heterozygosity of the MSSE region, with loss of the normal allele, providing the first evidence that MSSE is a tumour suppressor gene.

The synaptonemal complex and meiotic recombination in humans: new approaches to old questions

Meiotic prophase serves as an arena for the interplay of two important cellular activities, meiotic recombination and synapsis of homologous chromosomes. Synapsis is mediated by the synaptonemal complex (SC), originally characterized as a structure linked to pairing of meiotic chromosomes (Moses (1958) J Biophys Biochem Cytol 4:633-638). In 1975, the first electron micrographs of human pachytene stage SCs were presented (Moses et al. (1975) Science 187:363-365) and over the next 15 years the importance of the SC to normal meiotic progression in human males and females was established (Jhanwar and Chaganti (1980) Hum Genet 54:405-408; Pathak and Elder (1980) Hum Genet 54:171-175; Solari (1980) Chromosoma 81:315-337; Speed (1984) Hum Genet 66:176-180; Wallace and Hulten (1985) Ann Hum Genet 49(Pt 3):215-226). Further, these studies made it clear that abnormalities in the assembly or maintenance of the SC were an important contributor to human infertility (Chaganti et al. (1980) Am J Hum Genet 32:833-848; Vidal et al. (1982) Hum Genet 60:301-304; Bojko (1983) Carlsberg Res Commun 48:285-305; Bojko (1985) Carlsberg Res Commun 50:43-72; Templado et al. (1984) Hum Genet 67:162-165; Navarro et al. (1986) Hum Reprod 1:523-527; Garcia et al. (1989) Hum Genet 2:147-53). However, the utility of these early studies was limited by lack of information on the structural composition of the SC and the identity of other SC-associated proteins. Fortunately, studies of the past 15 years have gone a long way toward remedying this problem. In this minireview, we highlight the most important of these advances as they pertain to human meiosis, focusing on temporal aspects of SC assembly, the relationship between the SC and meiotic recombination, and the contribution of SC abnormalities to human infertility.

The case for strategic international alliances to harness nutritional genomics for public and personal health

Nutrigenomics is the study of how constituents of the diet interact with genes, and their products, to alter phenotype and, conversely, how genes and their products metabolise these constituents into nutrients, antinutrients, and bioactive compounds. Results from molecular and genetic epidemiological studies indicate that dietary unbalance can alter gene-nutrient interactions in ways that increase the risk of developing chronic disease. The interplay of human genetic variation and environmental factors will make identifying causative genes and nutrients a formidable, but not intractable, challenge. We provide specific recommendations for how to best meet this challenge and discuss the need for new methodologies and the use of comprehensive analyses of nutrient-genotype interactions involving large and diverse populations. The objective of the present paper is to stimulate discourse and collaboration among nutrigenomic researchers and stakeholders, a process that will lead to an increase in global health and wellness by reducing health disparities in developed and developing countries.

Why do human diversity levels vary at a megabase scale?

Levels of diversity vary across the human genome. This variation is caused by two forces: differences in mutation rates and the differential impact of natural selection. Pertinent to the question of the relative importance of these two forces is the observation that both diversity within species and interspecies divergence increase with recombination rates. This suggests that mutation and recombination are either directly coupled or linked through some third factor. Here, we test these possibilities using the recently generated sequence of the chimpanzee genome and new estimates of human diversity. We find that measures of GC and CpG content, simple-repeat structures, as well as the distance from the centromeres and the telomeres predict diversity as well as divergence. After controlling for these factors, large-scale recombination rates measured from pedigrees are still significant predictors of human diversity and human-chimpanzee divergence. Furthermore, the correlation between human diversity and recombination remains significant even after controlling for human-chimpanzee divergence. Two plausible and non-mutually exclusive explanations are, first, that natural selection has shaped the patterns of diversity seen in humans and, second, that recombination rates across the genome have changed since humans and chimpanzees shared a common ancestor, so that current recombination rates are a better predictor of diversity than of divergence. Because there are indications that recombination rates may have changed rapidly during human evolution, we favor the latter explanation.