Molecular genetics of speciation and human origins

Graft-versus-host disease prophylaxis shapes T cell biology and immune reconstitution after hematopoietic cell transplant

Successful hematopoietic cell transplant requires immunosuppression to prevent graft-versus-host disease (GVHD), a lethal, T-cell-mediated post-transplant complication. The phase 3 BMT CTN 1703 trial demonstrated superior GVHD-free/relapse-free survival for post-transplant cyclophosphamide (PT-Cy)-based GVHD prophylaxis versus tacrolimus/methotrexate (Tac/MTX), but did not improve overall survival. To compare T-cell biology between GVHD prophylaxis regimens, 324 patients were co-enrolled onto BMT CTN 1801 (NCT03959241). We quantified T-cell immune reconstitution using multi-modal analysis, including T-cell receptor (TCR) sequencing of 2,359 longitudinal samples (180,432,350 T-cells). Compared to Tac/MTX, PT-Cy was associated with an early, substantial reduction in TCR diversity that was sustained for 2 years. PT-Cy led to a T-cell reconstitution bottleneck, including reduced thymic output and virus-associated TCRs. Decreased D+14 TCR diversity predicted prevention of chronic GVHD, but also correlated with increased moderate-to-severe infections. This study reveals how distinct immunosuppression strategies have significant effects on the global immune repertoire, underpinning post-transplant clinical outcomes.

Reassessing Hybridisation in Australian Tetragonula Stingless Bees Using Multiple Genetic Markers

We re-examined reports of hybridisation in three cryptic stingless bee species in the genus Tetragonula in South East Queensland, Australia (T. carbonaria, T. davenporti and T. hockingsi). Previous studies on this group using microsatellite markers proposed that hybridisation occasionally takes place. In contrast, we find that using 1745 SNPs we could reliably separate the three species, with no evidence of contemporary (or recent) hybridisation. We found identical amplicon sequences of the nuclear gene EF1alpha across most individuals of the three species, but low and moderate species-specific polymorphisms in the nuclear gene Opsin and the mitochondrial 16S rRNA gene, respectively, with no cases of mito-nuclear discordance at these genes. We confirm that nuclear divergence across these species is low, based on 10-26 kb of non-coding sequence flanking EF1alpha and Opsin (0.7%-1% pairwise difference between species). However, we find mitogenomes to be far more diverged than nuclear genomes (21.6%-23.6% pairwise difference between species). Based on these comprehensive analyses of multiple marker types, we conclude there is no ongoing gene flow among the Tetragonula species of South East Queensland, despite their morphological similarity to one another and the low nuclear divergence among them. The higher resolution provided by multiple SNP markers may lead to lower estimates of contemporary hybridisation more generally.

Characterization of Human-malarial Parasite Species based on DHFR and GST Targets Resulting in Changes in Anti-malarial Drug Binding Conformations

BACKGROUND

In this study, we focused primarily on three anti-malarial drugs, namely chloroquine, mefloquine, and proguanil, and these were tested against two malarial targets DHFR and GST. The species Plasmodium falciparum, Plasmodium malariae, Plasmodium ovale, Plasmodium vivax were used for the study.

OBJECTIVE

The purpose of this study was to determine the sequence and structural similarity of the proteins DHFR and GST among four Plasmodium species as well as to discover the in silico interactions with the aforementioned drug candidates.

METHODS

Bioinformatics databases, such as PDB, UniProt, DrugBank, PubChem, and tools, and software like Phyre 2.0, Clustal O (1.2.4), AutoDock 4, AutoDock Vina, and Discovery Studio Visualizer were used to determine the evolutionary significance of the Plasmodium species.

RESULT

The variations showed a difference in the binding patterns of drugs with our target proteins. Our finding reveals the Plasmodium spp divergence or convergence as well as the structural and sequential similarity or dissimilarity features.

CONCLUSION

Our result suggests that due to the deviation in the sequences and structures, variations in protein-drug binding patterns have emerged.

Capturing Differential Allele-Level Expression and Genotypes of All Classical HLA Loci and Haplotypes by a New Capture RNA-Seq Method

The highly polymorphic human major histocompatibility complex (MHC) also known as the human leukocyte antigen (HLA) encodes class I and II genes that are the cornerstone of the adaptive immune system. Their unique diversity (>25,000 alleles) might affect the outcome of any transplant, infection, and susceptibility to autoimmune diseases. The recent rapid development of new next-generation sequencing (NGS) methods provides the opportunity to study the influence/correlation of this high level of HLA diversity on allele expression levels in health and disease. Here, we describe the NGS capture RNA-Seq method that we developed for genotyping all 12 classical HLA loci (HLA-A, HLA-B, HLA-C, HLA-DPA1, HLA-DPB1, HLA-DQA1, HLA-DQB1, HLA-DRA, HLA-DRB1, HLA-DRB3, HLA-DRB4, and HLA-DRB5) and assessing their allelic imbalance by quantifying their allele RNA levels. This is a target enrichment method where total RNA is converted to a sequencing-ready complementary DNA (cDNA) library and hybridized to a complex pool of RNA-specific HLA biotinylated oligonucleotide capture probes, prior to NGS. This method was applied to 161 peripheral blood mononuclear cells and 48 umbilical cord blood cells of healthy donors. The differential allelic expression of 10 HLA loci (except for HLA-DRA and HLA-DPA1) showed strong significant differences (P < 2.1 × 10-15). The results were corroborated by independent methods. This newly developed NGS method could be applied to a wide range of biological and medical questions including graft rejections and HLA-related diseases.

Computational analysis of functional single nucleotide polymorphisms associated with SLC26A4 gene

Single Nucleotide Polymorphisms (SNPs) are the most common candidate mutations in human beings that play a vital role in the genetic basis of certain diseases. Previous studies revealed that Solute Carrier Family 26 Member 4 (SLC26A4) being an essential gene of the multi-faceted transporter family SLC26 facilitates reflexive movement of Iodide into follicular lumen through apical membrane of thyrocyte. SLC26A4 gene encodes Pendred protein, a membrane glycoprotein, highly hydrophobic in nature, present at the apical membrane of thyrocyte functioning as transporter of iodide for thyroid cells. A minor genetic variation in SLC26A4 can cause Pendred syndrome, a syndrome associated with thyroid glands and deafness. In this study, we performed in-silico analysis of 674 missense SNPs of SLC26A4 using different computational platforms. The bunch of tools including SNPNEXUS, SNAP-2, PhD-SNP, SNPs&GO, I-Mutant, ConSurf, and ModPred were used to predict 23 highly confident damaging and disease causing nsSNPs (G209V, G197R, L458P, S427P, Q101P, W472R, N392Y, V359E, R409C, Q235R, R409P, G139V, G497S, H723R, D87G, Y127H, F667C, G334A, G95R, S427C, R291W, Q383H and E384G) that could potentially alter the SLC26A4 gene. Moreover, protein structure prediction, protein-ligand docking and Molecular Dynamics simulation were performed to confirm the impact of two evident alterations (Y127H and G334A) on the protein structure and function.

Protective role of HLA-DRB1*11 against juvenile idiopathic arthritis living in North Eastern Iran

Objective(s):

Juvenile idiopathic arthritis (JIA) is one of the most common chronic rheumatic diseases in children. The complex nature of this immune-mediated disease owes itself to several predisposing genes and environmental factors affecting its pathogenesis. Conducted in Iran, this study was originally intended to investigate every possible association between HLA DRB1 alleles and a susceptibility to JIA.

Materials and Methods:

In this case-control study, 45 patients with a definite diagnosis of JIA based on International League against Rheumatism (ILAR) criteria were compared against 46 healthy controls. DNA samples taken from both groups were analyzed using PCR-sequence specific primers (PCR-SSP) method. Data analysis including parametric and nonparametric test and multivariate analysis was undertaken using the SPSS 11.5 software. A P-value< 0.05 was regarded as statistically significant.

Results:

Mean ages in case group and healthy controls were 14.64±6.21 and 13.73±6.39, respectively with no significant difference between the two groups (P=0.515). Sex difference between JIA group and healthy controls was also not significant (P=0.068). The frequency of HLA-DRB1*01 was found the most frequent HLA-RB1 in our patients (33.3%). No significant statistical correlation between various HLA-DRB1 alleles and clinical subtypes of the disease could be established from the data. HLA-DRB1*11 was shown to raise protection to JIA (P=0.035, OR=2.755, 95% CI=0.963-8.055) in northeastern Iran. In addition, we found that HLA-RB1*09 is nominally associated with an increased risk of JIA (P=0.56, OR=2, 05, 95% CI=0.18-23.63).

Conclusion:

HLA-DRB1*11 was shown to raise protection to JIA in northeastern Iran. The disparity of findings in other ethnicities prompts further investigations with larger sample sizes.

Differential allelic expression of SOS1 and hyperexpression of the activating SOS1 c.755C variant in a Noonan syndrome family

Noonan syndrome (NS) is a genetic condition characterized by congenital heart defects, short stature and characteristic facial features. We here present the case of a girl with moderate learning disabilities, delayed language development, craniofacial features and skin anomalies reminiscent of NS. After a mutation screening of the known NS genes PTPN11, SOS1, RAF1, KRAS, GRB2, BRAF and SHOC2 we found the heterozygous c.755T>C variant in SOS1 causing the p.I252T amino-acid substitution, which was considered possibly pathogenetic by bioinformatic predictions. The same variant was present in the proband's mother, displaying some NS features, and maternal grandfather showing no NS traits, but also by a healthy subject in 1000 genomes project database without phenotype informations. The functional analysis revealed that SOS1 c.755C activated the RAS-ERK intracellular pathway, whereas no effects on RAC-JNK cascade have been detected. After a comparison between the sequence of SOS1 cDNA from peripheral blood and SOS1 genomic DNA, we showed for the first time a differential allelic expression of the SOS1 gene in healthy individuals, thus occurring as a physiologic condition. Interestingly, we found that the mutated allele C was 50% more expressed than the wild-type allele T in all familial carriers. The comparable amount of SOS1 mRNA between mutated individuals and the controls indicates that the variant does not affect SOS1 expression. The present study provides a first evidence of allelic imbalance of SOS1 and pinpoints this condition as a possible mechanism underlying a different penetrance of some SOS1-mutated alleles in unrelated carriers.

Major histocompatibility complex class II polymorphisms in forest musk deer (Moschus berezovskii) and their probable association with purulent disease

Genes of the major histocompatibility complex (MHC) family are crucial in immune responses because they present pathogenic peptides to T cells. In this study, we analysed the genetic variation in forest musk deer (Moschus berezovskii) MHC II genes and its potential association with musk deer purulent disease. In total, 53 purulent disease-susceptible and 46 purulent disease-resistant individuals were selected for MHC II exon 2 fragment analysis. Among them, 16 DQ alleles and four additional DR alleles were identified, with DQ exon 2 fragments displaying a low level of polymorphism. The nonsynonymous substitutions exceeded the synonymous substitutions in the peptide-binding sites of DQA2, DQB1 and DQB2. Then, 28 MHC II alleles were used to analyse the distribution patterns of purulent disease between the susceptible and resistant groups. Among them, three alleles (DQA1*01, DQA1*02 and DQA2*04) were found to be resistant, and five alleles (DRB3*07, DQA1*03, DQA1*04, DQA2*05 and DQA2*06) were found to increase susceptibility. Additionally, three haplotypes were found to be putatively associated with musk deer purulent disease. However, these three haplotypes were only found in the resistant or susceptible group, and their frequencies were low. The results from our study support a contributory role of MHC II polymorphisms in the development of purulent disease in forest musk deer.

The founding of Mauritian endemic coffee trees by a synchronous long-distance dispersal event

The stochastic process of long-distance dispersal is the exclusive means by which plants colonize oceanic islands. Baker's rule posits that self-incompatible plant lineages are unlikely to successfully colonize oceanic islands because they must achieve a coordinated long-distance dispersal of sufficiently numerous individuals to establish an outcrossing founder population. Here, we show for the first time that Mauritian Coffea species are self-incompatible and thus represent an exception to Baker's rule. The genus Coffea (Rubiaceae) is composed of approximately 124 species with a paleotropical distribution. Phylogenetic evidence strongly supports a single colonization of the oceanic island of Mauritius from either Madagascar or Africa. We employ Bayesian divergence time analyses to show that the colonization of Mauritius was not a recent event. We genotype S-RNase alleles from Mauritian endemic Coffea, and using S-allele gene genealogies, we show that the Mauritian allelic diversity is confined to just seven deeply divergent Coffea S-RNase allelic lineages. Based on these data, we developed an individual-based model and performed a simulation study to estimate the most likely number of founding individuals involved in the colonization of Mauritius. Our simulations show that to explain the observed S-RNase allelic diversity, the founding population was likely composed of fewer than 31 seeds that were likely synchronously dispersed from an ancestral mainland species.

MHC variability supports dog domestication from a large number of wolves: high diversity in Asia

The process of dog domestication is still somewhat unresolved. Earlier studies indicate that domestic dogs from all over the world have a common origin in Asia. So far, major histocompatibility complex (MHC) diversity has not been studied in detail in Asian dogs, although high levels of genetic diversity are expected at the domestication locality. We sequenced the second exon of the canine MHC gene DLA-DRB1 from 128 Asian dogs and compared our data with a previously published large data set of MHC alleles, mostly from European dogs. Our results show that Asian dogs have a higher MHC diversity than European dogs. We also estimated that there is only a small probability that new alleles have arisen by mutation since domestication. Based on the assumption that all of the currently known 102 DLA-DRB1 alleles come from the founding wolf population, we simulated the number of founding wolf individuals. Our simulations indicate an effective population size of at least 500 founding wolves, suggesting that the founding wolf population was large or that backcrossing has taken place.

Trans-species polymorphism and selection in the MHC class II DRA genes of domestic sheep

Highly polymorphic genes with central roles in lymphocyte mediated immune surveillance are grouped together in the major histocompatibility complex (MHC) in higher vertebrates. Generally, across vertebrate species the class II MHC DRA gene is highly conserved with only limited allelic variation. Here however, we provide evidence of trans-species polymorphism at the DRA locus in domestic sheep (Ovis aries). We describe variation at the Ovar-DRA locus that is far in excess of anything described in other vertebrate species. The divergent DRA allele (Ovar-DRA*0201) differs from the sheep reference sequences by 20 nucleotides, 12 of which appear non-synonymous. Furthermore, DRA*0201 is paired with an equally divergent DRB1 allele (Ovar-DRB1*0901), which is consistent with an independent evolutionary history for the DR sub-region within this MHC haplotype. No recombination was observed between the divergent DRA and B genes in a range of breeds and typical levels of MHC class II DR protein expression were detected at the surface of leukocyte populations obtained from animals homozygous for the DRA*0201, DRB1*0901 haplotype. Bayesian phylogenetic analysis groups Ovar-DRA*0201 with DRA sequences derived from species within the Oryx and Alcelaphus genera rather than clustering with other ovine and caprine DRA alleles. Tests for Darwinian selection identified 10 positively selected sites on the branch leading to Ovar-DRA*0201, three of which are predicted to be associated with the binding of peptide antigen. As the Ovis, Oryx and Alcelaphus genera have not shared a common ancestor for over 30 million years, the DRA*0201 and DRB1*0901 allelic pair is likely to be of ancient origin and present in the founding population from which all contemporary domestic sheep breeds are derived. The conservation of the integrity of this unusual DR allelic pair suggests some selective advantage which is likely to be associated with the presentation of pathogen antigen to T-cells and the induction of protective immunity.

Phylogenetic reconstruction and the identification of ancient polymorphism in the Bovini tribe (Bovidae, Bovinae)

BACKGROUND

The Bovinae subfamily incorporates an array of antelope, buffalo and cattle species. All of the members of this subfamily have diverged recently. Not surprisingly, a number of phylogenetic studies from molecular and morphological data have resulted in ambiguous trees and relationships amongst species, especially for Yak and Bison species. A partial phylogenetic reconstruction of 13 extant members of the Bovini tribe (Bovidae, Bovinae) from 15 complete or partially sequenced autosomal genes is presented.

RESULTS

We identified 3 distinct lineages after the Bovini split from the Boselaphini and Tragelaphini tribes, which has lead to the (1) Buffalo clade (Bubalus and Syncerus species) and a more recent divergence leading to the (2) Banteng, Gaur and Mithan and (3) Domestic cattle clades. A fourth lineage may also exist that leads to Bison and Yak. However, there was some ambiguity as to whether this was a divergence from the Banteng/Gaur/Mithan or the Domestic cattle clade. From an analysis of approximately 30,000 sites that were amplified in all species 133 sites were identified with ambiguous inheritance, in that all trees implied more than one mutation at the same site. Closer examination of these sites has identified that they are the result of ancient polymorphisms that have subsequently undergone lineage sorting in the Bovini tribe, of which 53 have remained polymorphic since Bos and Bison species last shared a common ancestor with Bubalus between 5-8 million years ago (MYA).

CONCLUSION

Uncertainty arises in our phylogenetic reconstructions because many species in the Bovini diverged over a short period of time. It appears that a number of sites with ambiguous inheritance have been maintained in subsequent populations by chance (lineage sorting) and that they have contributed to an association between Yak and Domestic cattle and an unreliable phylogenetic reconstruction for the Bison/Yak clade. Interestingly, a number of these aberrant sites are in coding sections of the genome and their identification may have important implications for studying the neutral rate of mutation at nonsynonymous sites. The presence of these sites could help account for the apparent contradiction between levels of polymorphism and effective population size in domesticated cattle.

Full-length sequence analysis of the HLA-DRB1 locus suggests a recent origin of alleles

The HLA region harbors some of the most polymorphic loci in the human genome. Among them is the class II locus HLA-DRB1, with more than 400 known alleles. The age of the polymorphism and the rate at which new alleles are generated at HLA loci has caused much controversy over the years. Previous studies have mostly been restricted to the 270 base pairs that constitute the second exon and represent the most variable part of the gene. Here, we investigate the evolutionary history of the HLA-DRB1 locus on the basis of an analysis of 15 genomic full-length alleles (10-15 kb). In addition, the variation in 49 complete coding sequences and 322 exon 2 sequences were analyzed. When excluding exon 2 from the analysis, the diversity at the synonymous sites was found to be similar to the intron diversity. The overall diversity in noncoding region was also similar to the genome average. The DRB1*03 lineage has been found in human, chimpanzee, bonobo, gorilla, and orangutan. An ancestral "proto HLA-DRB1*03 lineage" appeared to have diverged in the last 5 million years into the human-specific lineages *08, *11, *13, and *14. With exception to exon 2, both the coding- and the noncoding diversity suggests a recent origin (<1 million years ago) for most of the alleles at the HLA-DRB1 locus. Sites encoding for amino acids involved in antigen binding [antigen recognizing sites (ARS)] appear to have a more ancient origin. Taken together, the recent origin of most alleles, the high diversity between allelic lineages, and the ancient origin of sequence motifs in exon 2, is consistent with a relatively rapid generation of novel alleles by gene conversion like events.

Extreme context specificity in differential allelic expression

Variability in cis-regulation of gene expression has been implicated in the phenotypic manifestation of complex traits including common, multifactorial diseases. The differential expression of alleles due to polymorphism in cis-regulatory elements is common in the human genome, but there is a paucity of information about the context specificity of these control elements. In this study, we examined the differential allelic expression (DAE) of BMP5 in human mesenchymal tissues obtained from 16 donors undergoing joint replacement for treatment of osteoarthritis. We observed significant differences in BMP5 allelic output, with allelic ratios greater than 4:1 (P < 10(-20)) in the tissues of some donors. We also discovered a significant variability in allelic expression within the different tissues of donors. For 12 of our donors, we examined the allelic expression of BMP5 in two different regions of cartilage: cartilage adjacent to the site of the osteoarthritic lesion and cartilage distal from the lesion. Five of these 12 donors demonstrated highly significant differences (P < or = 10(-8)) in allelic expression between the different regions of their cartilage. Using DAE as a phenotype, we attempted to map tissue-specific cis-regulatory polymorphisms, and we identified a single nucleotide polymorphism located downstream of BMP5, which was significantly associated with DAE in some but not all of the examined tissues. These findings suggest that allelic expression can be highly context specific and that when interrogating the cis-regulatory control of a particular gene, one cannot necessarily assume that allelic expression is conserved across different tissues or even across different regions of the same tissue.

Conserved extended haplotypes of the major histocompatibility complex: further characterization

Since the complete sequencing of a human major histocompatibility complex (MHC) haplotype, interest in non-human leucocyte antigen (HLA) genes encoded in the MHC has been growing. Non-HLA genes, which outnumber the HLA genes, may contribute to or account for HLA and disease associations. Most information on non-HLA genes has been obtained in separate studies of individual loci. To comprehensively address polymorphisms of relevant non-HLA genes in 'conserved extended haplotypes' (CEH), we investigated 101 International Histocompatibility Workshop reference cell lines and nine additional anonymous samples representing all 37 unambiguously characterized CEHs at MICA, NFKBIL1, LTA, NCR3, AIF1, HSPA1A, HSPA1B, BF, NOTCH4 and a single nucleotide polymorphism (SNP) at HLA-DQA1 as well as MICA, NOTCH4, HSPA1B and all five tumour necrosis factor short tandem repeat (STR) polymorphisms. This work (1) provides an extensive catalogue of MHC polymorphisms in all CEHs, (2) unravels interrelationships between HLA and non-HLA haplotypical lineages, (3) resolves reported typing ambiguities and (4) describes haplospecific markers for a number of CEHs. Analysis also identified a DQA1 SNP and segments containing MHC class III polymorphisms that corresponded with class II (DRB3 and DRB4) lineages. These results portray the MHC where lineages containing non-HLA and HLA variants in linkage disequilibrium may operate in concert and can guide more thorough design and interpretation of HLA-disease relationships.

Rapid evolution of major histocompatibility complex class I genes in primates generates new disease alleles in humans via hitchhiking diversity

A plausible explanation for many MHC-linked diseases is lacking. Sequencing of the MHC class I region (coding units or full contigs) in several human and nonhuman primate haplotypes allowed an analysis of single nucleotide variations (SNV) across this entire segment. This diversity was not evenly distributed. It was rather concentrated within two gene-rich clusters. These were each centered, but importantly not limited to, the antigen-presenting HLA-A and HLA-B/-C loci. Rapid evolution of MHC-I alleles, as evidenced by an unusually high number of haplotype-specific (hs) and hypervariable (hv) (which could not be traced to a single species or haplotype) SNVs within the classical MHC-I, seems to have not only hitchhiked alleles within nearby genes, but also hitchhiked deleterious mutations in these same unrelated loci. The overrepresentation of a fraction of these hvSNV (hv1SNV) along with hsSNV, as compared to those that appear to have been maintained throughout primate evolution (trans-species diversity; tsSNV; included within hv2SNV) tends to establish that the majority of the MHC polymorphism is de novo (species specific). This is most likely reminiscent of the fact that these hsSNV and hv1SNV have been selected in adaptation to the constantly evolving microbial antigenic repertoire.

Population bottlenecks in quasispecies dynamics

The characteristics of natural populations result from different stochastic and deterministic processes that include reproduction with error, selection, and genetic drift. In particular, population fluctuations constitute a stochastic process that may play a very relevant role in shaping the structure of populations. For example, it is expected that small asexual populations will accumulate mutations at a higher rate than larger ones. As a consequence, in any population the fixation of mutations is accelerated when environmental conditions cause population bottlenecks. Bottlenecks have been relatively frequent in the history of life and it is generally accepted that they are highly relevant for speciation. Although population bottlenecks can occur in any species, their effects are more noticeable in organisms that form large and heterogeneous populations, such as RNA viral quasispecies. Bottlenecks can also positively select and isolate particles that still keep the ability to infect cells from a disorganized population created by crossing the error threshold.

The AGT gene in Africa: a distinctive minor allele haplotype, a polymorphism (V326I), and a novel PH1 mutation (A112D) in Black Africans

We describe a novel missense mutation (A112D) and polymorphism (V326I) in the human AGT gene in two black African patients with primary hyperoxaluria type 1, an autosomal recessive disease resulting from a deficiency of the liver peroxisomal enzyme alanine:glyoxylate aminotransferase (AGT; EC 2.6.1.44). V326I was found in DNA from normal control Blacks with an allele frequency of 3%. Expression studies confirmed that A112D reduced AGT enzyme activity by 95% while V326I had no effect. Both A112D and V326I were homozygous in both patients and lie on a variant of the minor allele of the AGT gene. This variant haplotype, Mi(A), includes an intron 1 duplication and intron 4 VNTR (38 repeat) but lacks the P11L and I340M normally associated with the minor allele in Caucasians. Among the South African Blacks tested, the Mi(A) haplotype had an allele frequency of 12% compared to 3 % for the Caucasian-type minor allele haplotype.

Is a Multivariate Consensus Representation of Genetic Relationships Among Populations Always Meaningful?

To determine the relationships among closely related populations or species, two methods are commonly used in the literature: phylogenetic reconstruction or multivariate analysis. The aim of this article is to assess the reliability of multivariate analysis. We describe a method that is based on principal component analysis and Mantel correlations, using a two-step process: The first step consists of a single-marker analysis and the second step tests if each marker reveals the same typology concerning population differentiation. We conclude that if single markers are not congruent, the compromise structure is not meaningful. Our model is not based on any particular mutation process and it can be applied to most of the commonly used genetic markers. This method is also useful to determine the contribution of each marker to the typology of populations. We test whether our method is efficient with two real data sets based on microsatellite markers. Our analysis suggests that for closely related populations, it is not always possible to accept the hypothesis that an increase in the number of markers will increase the reliability of the typology analysis.

Increased heterozygosity for MHC class II lineages in newborn males

In plants, fungi and marine invertebrates, there are genetic compatibility systems to ensure diversity in the offspring. The importance of genetic compatibility in gametic union and selective abortion in vertebrate animals has also been appreciated recently. There have been suggestions that the major histocompatibility complex (HLA in humans) may be a compatibility system in vertebrates. HLA class II haplotypes often contain a second expressed DRB locus which can be either DRB3, DRB4 or DRB5. These encode the supertypical specificities and mark the ancestral lineages. The members of each lineage have related DNA sequences at the main class II locus HLA-DRB1. We analysed 415 newborns at all expressed DRB loci by PCR analysis to seek evidence for sex-specific prenatal selection events. While there was no significant change in heterozygosity rates between males and females at DRB1, the proportion of males carrying two DRB1 specificities from different ancestral lineages was significantly increased (53.7% in males vs 39.3% in females, P = 0.003). The genotypes consisting of phylogenetically most distinct ones, namely the DRB3 and DRB4 haplotypes, showed the most striking difference between sexes (P = 0.007). These results suggested a more favourable outcome for male concepti heterozygous for supertypical haplotypes. Heterozygosity for most divergent haplotypical families ensures the highest degree of functional heterozygosity at the main HLA class II locus DRB1 while increasing the likelihood of heterozygosity also at other MHC loci. Our observations agree with the previously reported heterozygote excess in male newborn rats and mice. Correlations between MHC class II heterozygosity and advertised male quality in deer and pheasant as well as increased reproductive success in MHC class II heterozygous male macaques are examples of postnatal benefits of heterozygosity in males that may be behind the development of prenatal selection mechanisms. The MHC-mediated prenatal selection of males may also be one of the selective events suggested by the very high primary (male-to-female) sex ratio at fertilization reaching close to unity at birth in humans. These results provide an appealing working hypothesis for further studies in humans and other vertebrates.

The genetical history of humans and the great apes

When and where did modern humans evolve? How did our ancestors spread over the world? Traditionally, answers to questions such as these have been sought in historical, archaeological, and fossil records. However, increasingly genetic data provide information about the evolution of our species. In this review, we focus on the comparison of the variation in the human gene pool to that of our closest evolutionary relatives, the great apes, because this provides a relevant perspective on human genetical evolution. For instance, comparisons to the great apes show that humans are unique in having little genetic variation as well as little genetic structure in their gene pool. Furthermore, genetic data indicate that humans, but not the great apes, have experienced a period of dramatic growth in their early history.

The genetics of hyphal fusion and vegetative incompatibility in filamentous ascomycete fungi

Filamentous fungi grow as a multicellular, multinuclear network of filament-shaped cells called hyphae. A fungal individual can be viewed as a fluid, dynamic system that is characterized by hyphal tip growth, branching, and hyphal fusion (anastomosis). Hyphal anastomosis is especially important in such nonlinear systems for the purposes of communication and homeostasis. Filamentous fungi can also undergo hyphal fusion with different individuals to form heterokaryons. However, the viability of such heterokaryons is dependent upon genetic constitution at heterokaryon incompatibility (het) loci. If hyphal fusion occurs between strains that differ in allelic specificity at het loci, vegetative incompatibility, which is characterized by hyphal compartmentation and cell lysis, is induced. This review covers microscopic and genetic analysis of hyphal fusion and the molecular and genetic analysis of the consequence of hyphal fusion between individuals that differ in specificity at het loci in filamentous ascomycetes.

Mining for SNPs: putting the common variants--common disease hypothesis to the test

Classical molecular genetic strategies have succeeded in identifying mutations responsible for numerous rare diseases with Mendelian patterns of inheritance, but have been largely unsuccessful in unravelling the (genetic basis of complex medical conditions like cardiovascular disease' diabetes and mental illness. These common disorders are shaped by multiple genes that exert weak allelic effects in the setting of confounding environmental variables. Association study designs provide statistical povwer to reveal the modest contributions of weak alleles, and evidence is mounting that common genetic polymorphisms play a role in complex diseases. Cataloguing genetic variation in human populations is a prerequisite for further validation of the 'common variants-common disease' hypothesis, and polymorphism discovery has begun in earnest in the academic and private sector. We will review several strategies for high-throughput polymorphism discovery and discuss the implications of early results from polymorphism screens for future genetic studies.

Cellular and molecular mechanisms of sexual incompatibility in plants and fungi

Plants and fungi show an astonishing diversity of mechanisms to promote outbreeding, the most widespread of which is sexual incompatibility. Sexual incompatibility involves molecular recognition between mating partners. In fungi and algae, highly polymorphic mating-type loci mediate mating through complementary interactions between molecules encoded or regulated by different mating-type haplotypes, whereas in flowering plants polymorphic self-incompatibility loci regulate mate recognition through oppositional interactions between molecules encoded by the same self-incompatibility haplotypes. This subtle mechanistic difference is a consequence of the different life cycles of fungi, algae, and flowering plants. Recent molecular and biochemical studies have provided fascinating insights into the mechanisms of mate recognition and are beginning to shed light on evolution and population genetics of these extraordinarily polymorphic genetic systems of incompatibility.

Characterization of single-nucleotide polymorphisms in coding regions of human genes

A major goal in human genetics is to understand the role of common genetic variants in susceptibility to common diseases. This will require characterizing the nature of gene variation in human populations, assembling an extensive catalogue of single-nucleotide polymorphisms (SNPs) in candidate genes and performing association studies for particular diseases. At present, our knowledge of human gene variation remains rudimentary. Here we describe a systematic survey of SNPs in the coding regions of human genes. We identified SNPs in 106 genes relevant to cardiovascular disease, endocrinology and neuropsychiatry by screening an average of 114 independent alleles using 2 independent screening methods. To ensure high accuracy, all reported SNPs were confirmed by DNA sequencing. We identified 560 SNPs, including 392 coding-region SNPs (cSNPs) divided roughly equally between those causing synonymous and non-synonymous changes. We observed different rates of polymorphism among classes of sites within genes (non-coding, degenerate and non-degenerate) as well as between genes. The cSNPs most likely to influence disease, those that alter the amino acid sequence of the encoded protein, are found at a lower rate and with lower allele frequencies than silent substitutions. This likely reflects selection acting against deleterious alleles during human evolution. The lower allele frequency of missense cSNPs has implications for the compilation of a comprehensive catalogue, as well as for the subsequent application to disease association.

Distribution of haplotypes from a chromosome 21 region distinguishes multiple prehistoric human migrations

Despite mounting genetic evidence implicating a recent origin of modern humans, the elucidation of early migratory gene-flow episodes remains incomplete. Geographic distribution of haplotypes may show traces of ancestral migrations. However, such evolutionary signatures can be erased easily by recombination and mutational perturbations. A 565-bp chromosome 21 region near the MX1 gene, which contains nine sites frequently polymorphic in human populations, has been found. It is unaffected by recombination and recurrent mutation and thus reflects only migratory history, genetic drift, and possibly selection. Geographic distribution of contemporary haplotypes implies distinctive prehistoric human migrations: one to Oceania, one to Asia and subsequently to America, and a third one predominantly to Europe. The findings with chromosome 21 are confirmed by independent evidence from a Y chromosome phylogeny. Loci of this type will help to decipher the evolutionary history of modern humans.

Evidence for balancing selection operating at the het-c heterokaryon incompatibility locus in a group of filamentous fungi

In filamentous fungi, het loci (for heterokaryon incompatibility) are believed to regulate self/nonself-recognition during vegetative growth. As filamentous fungi grow, hyphal fusion occurs within an individual colony to form a network. Hyphal fusion can occur also between different individuals to form a heterokaryon, in which genetically distinct nuclei occupy a common cytoplasm. However, heterokaryotic cells are viable only if the individuals involved have identical alleles at all het loci. One het locus, het-c, has been characterized at the molecular level in Neurospora crassa and encodes a glycine-rich protein. In an effort to understand the role of this locus in filamentous fungi, we chose to study its evolution by analyzing het-c sequence variability in species within Neurospora and related genera. We determined that the het-c locus was polymorphic in a field population of N. crassa with close to equal frequency of each of the three allelic types. Different species and even genera within the Sordariaceae shared het-c polymorphisms, indicating that these polymorphisms originated in an ancestral species. Finally, an analysis of the het-c specificity region shows a high occurrence of nonsynonymous substitution. The persistence of allelic lineages, the nearly equal allelic distribution within populations, and the high frequency of nonsynonymous substitutions in the het-c specificity region suggest that balancing selection has operated to maintain allelic diversity at het-c. Het-c shares this particular evolutionary characteristic of departing from neutrality with other self/nonself-recognition systems such as major histocompatibility complex loci in mammals and the S (self-incompatibility) locus in angiosperms.

Genetic polymorphism and natural selection in the malaria parasite Plasmodium falciparum

We have studied the genetic polymorphism at 10 Plasmodium falciparum loci that are considered potential targets for specific antimalarial vaccines. The polymorphism is unevenly distributed among the loci; loci encoding proteins expressed on the surface of the sporozoite or the merozoite (AMA-1, CSP, LSA-1, MSP-1, MSP-2, and MSP-3) are more polymorphic than those expressed during the sexual stages or inside the parasite (EBA-175, Pfs25, PF48/45, and RAP-1). Comparison of synonymous and nonsynonymous substitutions indicates that natural selection may account for the polymorphism observed at seven of the 10 loci studied. This inference depends on the assumption that synonymous substitutions are neutral, which we test by analyzing codon bias and G+C content in a set of 92 gene loci. We find evidence for an overall trend towards increasing A+T richness, but no evidence for mutation bias. Although the neutrality of synonymous substitutions is not definitely established, this trend towards an A+T rich genome cannot explain the accumulation of substitutions at least in the case of four genes (AMA-1, CSP, LSA-1, and PF48/45) because the Gleft and right arrow C transversions are more frequent than expected. Moreover, the Tajima test manifests positive natural selection for the MSP-1 and, less strongly, MSP-3 polymorphisms; the McDonald-Kreitman test manifests natural selection at LSA-1 and PF48/45. We conclude that there is definite evidence for positive natural selection in the genes encoding AMA-1, CSP, LSA-1, MSP-1, and Pfs48/45. For four other loci, EBA-175, MSP-2, MSP-3, and RAP-1, the evidence is limited. No evidence for natural selection is found for Pfs25.

Recent origin of HLA-DRB1 alleles and implications for human evolution

The HLA class I and class II loci are the most highly polymorphic coding regions in the human genome. Based on the similarity of the coding sequences of alleles between species, it has been claimed that the HLA polymorphism is ancient and predates the separation of human (Homo) and chimpanzee (Pan), 4-7.4 Myr ago. Analysis of intron sequences, however, provides support for a more recent origin and for rapid generation of alleles at the HLA class II DRB1 locus. The human DRB1 alleles can be divided into groups (allelic lineages); most of these lineages have diverged from each other before the separation of Homo and Pan. Alleles within such a lineage, however, appear to be, on average, 250,000 years old, implying that the vast majority (greater than 90%) of the more than 135 contemporary human DRB1 alleles have been generated after the separation of Homo and Pan. The coalescence time of alleles within allelic lineages indicates that the effective population size (Ne) for early hominids (over the last 1 Myr) was approximately 10(4) individuals, similar to estimates based on other nuclear loci and mitochondrial DNA. With a single exception, the genetic mechanisms (gene conversion and point mutation) that have diversified the exon-2 sequences do not appear to extend into the adjacent intron sequences. The part of exon 2 encoding the beta-sheet evolves in concert with the surrounding introns, while the alpha-helix appears to have been subjected to gene conversion-like events, suggesting that such exchange events are highly localised and occur over extremely short sequence tracts.

Analysis of intron sequences at the class II HLA-DRB1 locus: implications for the age of allelic diversity

Analyses of the coding sequences of HLA class II alleles have revealed high similarity between species, indicating that much of the polymorphism predates the separation of human (Homo) and chimpanzee (Pan), 4-7.4 million years ago. Recent studies of the intron sequences of alleles provide support for a much more recent origin and rapid generation of HLA alleles. At the DRB1 locus, intron analysis indicates that most of the allelic lineages have diverged from each other before the separation of Homo and Pan, consistent with the exon analysis. However, the intron sequences of alleles within lineages are almost identical, indicating that many of the alleles have been generated after the divergence of the Homo and Pan lineages.

Departure from neutrality at the mitochondrial NADH dehydrogenase subunit 2 gene in humans, but not in chimpanzees

To test whether patterns of mitochondrial DNA (mtDNA) variation are consistent with a neutral model of molecular evolution, nucleotide sequences were determined for the 1041 bp of the NADH dehydrogenase subunit 2 (ND2) gene in 20 geographically diverse humans and 20 common chimpanzees. Contingency tests of neutrality were performed using four mutational categories for the ND2 molecule: synonymous and nonsynonymous mutations in the transmembrane regions, and synonymous and nonsynonymous mutations in the surface regions. The following three topological mutational categories were also used: intraspecific tips, intraspecific interiors, and interspecific fixed differences. The analyses reveal a significantly greater number of nonsynonymous polymorphisms within human transmembrane regions than expected based on interspecific comparisons, and they are inconsistent with a neutral equilibrium model. This pattern of excess nonsynonymous polymorphism is not seen within chimpanzees. Statistical tests of neutrality, such as TAJIMA's D test, and the D and F tests proposed by FU and LI, indicate an excess of low frequency polymorphisms in the human data, but not in the chimpanzee data. This is consistent with recent directional selection, a population bottleneck or background selection of slightly deleterious mutations in human mtDNA samples. The analyses further support the idea that mitochondrial genome evolution is governed by selective forces that have the potential to affect its use as a "neutral" marker in evolutionary and population genetic studies.

Allelic specificity at the het-c heterokaryon incompatibility locus of Neurospora crassa is determined by a highly variable domain

In filamentous fungi, the ability to form a productive heterokaryon with a genetically dissimilar individual is controlled by specific loci termed het loci. Only strains homozygous for all het loci can establish a heterokaryon. In Neurospora crassa, 11 loci, including the mating-type locus, regulate the capacity to form heterokaryons. An allele of the het-c locus (het-cOR) of N. crassa has been previously characterized and encodes a nonessential 966 amino acid glycine-rich protein. Herein, we describe the genetic and molecular characterization of two hei-c alleles, het-cPA and het-cOR, that have a different specificity from that of het-cOR, showing that vegetative incompatibility is mediated by multiple alleles at het-c. By constructing chimeric alleles, we show that het-c specificity is determined by a highly variable domain of 34-48 amino acids in length. In this regard, het-c is similar to loci that regulate recognition in other species, such as the (S) self-incompatibility locus in plants, the sexual compatibility locus in basidiomycetes and the major histocompatibility complex (MHC) genes in vertebrates.

Allelic disequilibrium and allele frequency distribution as a function of social and demographic history

Allelic disequilibrium between closely linked genes is a common observation in human populations and often gives rise to speculation concerning the role of selective forces. In a previous treatment, we have developed a population model of the expected distribution of rare variants (including private polymorphisms) in Amerindians and have argued that, because of the great expansion of Amerindian numbers with the advent of agriculture, most of these rare variants are of relatively recent origin. Many other populations have similar histories of striking recent expansions. In this treatment, we demonstrate that, in consequence of this fact, a high degree of linkage disequilibrium between two nonhomologous alleles <0.5 cM apart is the "normal" expectation, even in the absence of selection. This expectation is enhanced by the previous subdivision of human populations into relatively isolated tribes characterized by a high level of endogamy and inbreeding. We also demonstrate that the alleles associated with a recessive disease phenotype are expected to exist in a population in very variable frequencies: there is no need to postulate positive selection with respect to the more common disease-associated alleles for such entities as phenylketonuria or cystic fibrosis.

Population biology of antigen presentation by MHC class I molecules

In principle, the function of major histocompatibility complex (MHC) molecules is simple: to bind a peptide and engage a T cell. In practice, placing this function within the context of the immune response begs questions of population biology; How does the immune response emerge from the interactions among populations of peptides, T cells and MHC molecules? Within a population of vertebrates, how does MHC polymorphism stamp individuality on the response? Does polymorphism confer differential advantages in responding to parasites? How are the pressures on the MHC reflected in turnover of alleles? The role of mutation, recombination, selection, and drift in the generation and maintenance of MHC class 1 polymorphism are considered.

Global patterns of linkage disequilibrium at the CD4 locus and modern human origins

Haplotypes consisting of alleles at a short tandem repeat polymorphism (STRP) and an Alu deletion polymorphism at the CD4 locus on chromosome 12 were analyzed in more than 1600 individuals sampled from 42 geographically dispersed populations (13 African, 2 Middle Eastern, 7 European, 9 Asian, 3 Pacific, and 8 Amerindian). Sub-Saharan African populations had more haplotypes and exhibited more variability in frequencies of haplotypes than the Northeast African or non-African populations. The Alu deletion was nearly always associated with a single STRP allele in non-African and Northeast African populations but was associated with a wide range of STRP alleles in the sub-Saharan African populations. This global pattern of haplotype variation and linkage disequilibrium suggests a common and recent African origin for all non-African human populations.

The myth of Eve: molecular biology and human origins

It has been proposed that modern humans descended from a single woman, the "mitochondrial Eve" who lived in Africa 100,000 to 200,000 years ago. The human immune system DRB1 genes are extremely polymorphic, with gene lineages that coalesce into an ancestor who lived around 60 million years ago, a time before the divergence of the apes from the Old World monkeys. The theory of gene coalescence suggests that, throughout the last 60 million years, human ancestral populations had an effective size of 100,000 individuals or greater. Molecular evolution data favor the African origin of modern humans, but the weight of the evidence is against a population bottleneck before their emergence. The mitochondrial Eve hypothesis emanates from a confusion between gene genealogies and individual genealogies.

New HLA-DPB1 alleles generated by interallelic gene conversion detected by analysis of sperm

The rate at which allelic diversity at the HLA loci evolves has been the subject of considerable controversy. The patchwork pattern of sequence polymorphism within the second exon of the HLA class II loci, particularly in the DPB1 locus, may have been generated by segmental exchange (gene conversion). We have analysed the frequency of variant DPB1 sequences that have been created by interallelic gene conversion in the germline by screening pools of sperm using PCR amplification and oligonucleotide probe typing. Our results indicate that about 1/10,000 sperm represents a new DPB1 sequence generated by short tracts of segmental exchange (gene conversion) within the second exon, suggesting that interallelic gene conversion may have an important role in generating the extensive allelic diversity at the HLA loci.

Absence of polymorphism at the ZFY locus on the human Y chromosome

DNA polymorphism in the Y chromosome, examined at a 729-base pair intron located immediately upstream of the ZFY zinc-finger exon, revealed no sequence variation in a worldwide sample of 38 human males. This finding cannot be explained by global constraint on the intron sequence, because interspecific comparisons with other nonhuman primates revealed phylogenetically informative sequence changes. The invariance likely results from either a recent selective sweep, a recent origin for modern Homo sapiens, recurrent male population bottlenecks, or historically small effective male population sizes. A coalescence model predicts an expected time to a most recent common ancestral male lineage of 270,000 years (95 percent confidence limits: 0 to 800,000 years).