TNFRSF13B Genetic variability an anthropological - evolutionary approach to Biomedical Research

In the recent years TNFRSF13B coding variants have been implicated by clinical genetics studies in Common Variable Immunodeficiency (CVID), the most common clinically relevant primary immunodeficiency in individuals of European ancestry, but their functional effects in relation to the development of the disease have not been entirely established. To examine the potential contribution of such variants to CVID, the more comprehensive perspective of an evolutionary approach was applied in this study, underling the belief that evolutionary genetics methods can play a role in dissecting the origin, causes and diffusion of human diseases, representing a powerful tool also in human health research. For this purpose, TNFRSF13B coding region was sequenced in 451 healthy individuals belonging to 26 worldwide populations, in addition to 96 control, 77 CVID and 38 Selective IgA Deficiency (IgAD) individuals from Italy, leading to the first achievement of a global picture of TNFRSF13B nucleotide diversity and haplotype structure and making suggestion of its evolutionary history possible. A slow rate of evolution, within our species and when compared to the chimpanzee, low levels of genetic diversity geographical structure and the absence of recent population specific selective pressures were observed for the examined genomic region, suggesting that geographical distribution of its variability is more plausibly related to its involvement also in innate immunity rather than in adaptive immunity only. This, together with the extremely subtle disease/healthy samples differences observed, suggests that CVID might be more likely related to still unknown environmental and genetic factors, rather than to the nature of TNFRSF13B variants only.

Cytogenetic And Morphologic Approaches Of Hybrids From Experimental Crosses Between Triatoma Lenti Sherlock & Serafim, 1967 And T. Sherlocki Papa Et Al., 2002 (hemiptera: Reduviidae).

The reproductive capacity between Triatoma lenti and Triatoma sherlocki was observed in order to verify the fertility and viability of the offspring. Cytogenetic, morphological and morphometric approaches were used to analyze the differences that were inherited. Experimental crosses were performed in both directions. The fertility rate of the eggs in crosses involving T. sherlocki females was 65% and 90% in F1 and F2 offspring, respectively. In reciprocal crosses, it was 7% and 25% in F1 and F2 offspring, respectively. The cytogenetic analyses of the male meiotic process of the hybrids were performed using lacto-acetic orcein, C-banding and Feulgen techniques. The male F1 offspring presented normal chromosome behavior, a finding that was similar to those reported in parental species. However, cytogenetic analysis of F2 offspring showed errors in chromosome pairing. This post-zygotic isolation, which prevents hybrids in nature, may represent the collapse of the hybrid. This phenomenon is due to a genetic dysregulation that occurs in the chromosomes of F1. The results were similar in the hybrids from both crosses. Morphological features, such as color and size of connexive and the presence of red-orange rings on the femora, were similar to T. sherlocki, while wins size was similar to T. lenti in F1 offspring. The eggshells showed characteristics that were similar to species of origin, whereas the median process of the pygophore resulted in intermediate characteristics in the F1 and a segregating pattern in F2 offspring. Geometric morphometric techniques used on the wings showed that both F1 and F2 offspring were similar to T. lenti. These studies on the reproductive capacity between T. lenti and T. sherlocki confirm that both species are evolutionarily closed; hence, they are included in the brasiliensis subcomplex. The extremely reduced fertility observed in the F2 hybrids confirmed the specific status of the species that were analyzed.

Molecular evolution and phylogeny of elapid snake venom three-finger toxins

Animal venom components are of considerable interest to researchers across a wide variety of disciplines, including molecular biology, biochemistry, medicine, and evolutionary genetics. The three-finger family of snake venom peptides is a particularly interesting and biochemically complex group of venom peptides, because they are encoded by a large multigene family and display a diverse array of functional activities. In addition, understanding how this complex and highly varied multigene family evolved is an interesting question to researchers investigating the biochemical diversity of these peptides and their impact on human health. Therefore, the purpose of our study was to investigate the long-term evolutionary patterns exhibited by these snake venom toxins to understand the mechanisms by which they diversified into a large, biochemically diverse, multigene family. Our results show a much greater diversity of family members than was previously known, including a number of subfamilies that did not fall within any previously identified groups with characterized activities. In addition, we found that the long-term evolutionary processes that gave rise to the diversity of three-finger toxins are consistent with the birth-and-death model of multigene family evolution. It is anticipated that this three-finger toxin toolkit will prove to be useful in providing a clearer picture of the diversity of investigational ligands or potential therapeutics available within this important family.

Poor value of pulsed-field gel electrophoresis to investigate long-term scale epidemiology of methicillin-resistant Staphylococcus aureus.

Pulsed-field gel electrophoresis (PFGE) is widely used for epidemic investigations of methicillin-resistant Staphylococcus aureus (MRSA). In the present study, we evaluated its use in a long-term epidemiological setting (years to few decades, country to continent level). The clustering obtained from PFGE patterns after SmaI digestion of the DNA of 20 strains was compared to that obtained using a phylogenetic typing method (multiprimer RAPD). The results showed that the analysis of small PFGE bands (10-85kb) correlates better with multiprimer RAPD than the analysis of large PFGE bands (>85-700kb), suggesting that the analysis of small bands would be more suitable for the investigation of long-term epidemiological setting. However, given the technical difficulties to obtain a good resolution of these bands and the putative presence of plasmids among them, PFGE does not appear to be a method of choice for the long-term epidemiology analysis of MRSA.

Integrated Genetic Epidemiology of Infectious Diseases: The Chagas Model

Genetic typing of pathogenic agents and of vectors has known impressive developments in the last 10 years, thanks to the progresses of molecular biology, and to the contribution of the concepts of evolutionary genetics. Moreover, we know more and more on the genetic susceptibility of man to infectious diseases. I propose here to settle a new, synthetic field of research, which I call `integrated genetic epidemiology of infectious diseases' (IGEID). I aim at evaluating, by an evolutionary genetic approach, the respective impact, on the transmission and pathogenicity of infectious diseases, of the host's, the pathogen's and the vector's genetic diversity, and their possible interactions (co-evolution phenomena). Chagas' disease constitutes a fine model to develop the IGEID methodology, by both field and experimental studies.

Genetic isolation of insular populations of the Maghrebian bat, Myotis punicus, in the Mediterranean Basin

We investigate the population genetic structure of the Maghrebian bat, Myotis punicus, between the mainland and islands to assess the island colonization pattern and current gene flow between nearby islands and within the mainland. Location North Africa and the Mediterranean islands of Corsica and Sardinia. Methods We sequenced part of the control region (HVII) of 79 bats across 11 colonies. The phylogeographical pattern was assessed by analysing molecular diversity indices, examining differentiation among populations and estimating divergence time. In addition, we genotyped 182 bats across 10 colonies at seven microsatellite loci. We used analysis of molecular variance and a Bayesian approach to infer nuclear population structure. Finally, we estimated sex-specific dispersal between Corsica and Sardinia. Results Mitochondrial analyses indicated that colonies between Corsica, Sardinia and North Africa are highly differentiated. Within islands there was no difference between colonies, while at the continental level Moroccan and Tunisian populations were highly differentiated. Analyses with seven microsatellite loci showed a similar pattern. The sole difference was the lack of nuclear differentiation between populations in North Africa, suggesting a male-biased dispersal over the continental area. The divergence time of Sardinian and Corsican populations was estimated to date back to the early and mid-Pleistocene. Main conclusions Island colonization by the Maghrebian bats seems to have occurred in a stepping-stone manner and certainly pre-dated human colonization. Currently, open water seems to prevent exchange of bats between the two islands, despite their ability to fly and the narrowness of the strait of Bonifacio. Corsican and Sardinian populations are thus currently isolated from any continental gene pool and must therefore be considered as different evolutionarily significant units (ESU).

Family Size Evolution in Drosophila Chemosensory Gene Families: A Comparative Analysis with a Critical Appraisal of Methods

Gene turnover rates and the evolution of gene family sizes are important aspects of genome evolution. Here, we use curated sequence data of the major chemosensory gene families from Drosophila-the gustatory receptor, odorant receptor, ionotropic receptor, and odorant-binding protein families-to conduct a comparative analysis among families, exploring different methods to estimate gene birth and death rates, including an ad hoc simulation study. Remarkably, we found that the state-of-the-art methods may produce very different rate estimates, which may lead to disparate conclusions regarding the evolution of chemosensory gene family sizes in Drosophila. Among biological factors, we found that a peculiarity of D. sechellia's gene turnover rates was a major source of bias in global estimates, whereas gene conversion had negligible effects for the families analyzed herein. Turnover rates vary considerably among families, subfamilies, and ortholog groups although all analyzed families were quite dynamic in terms of gene turnover. Computer simulations showed that the methods that use ortholog group information appear to be the most accurate for the Drosophila chemosensory families. Most importantly, these results reveal the potential of rate heterogeneity among lineages to severely bias some turnover rate estimation methods and the need of further evaluating the performance of these methods in a more diverse sampling of gene families and phylogenetic contexts. Using branch-specific codon substitution models, we find further evidence of positive selection in recently duplicated genes, which attests to a nonneutral aspect of the gene birth-and-death process.

Patterns of nucleotide diversity at the regions encompassing the Drosophila insulin-like peptide (dilp) genes: demography vs positive selection in Drosophila melanogaster.

In Drosophila, the insulin-signaling pathway controls some life history traits, such as fertility and lifespan, and it is considered to be the main metabolic pathway involved in establishing adult body size. Several observations concerning variation in body size in the Drosophila genus are suggestive of its adaptive character. Genes encoding proteins in this pathway are, therefore, good candidates to have experienced adaptive changes and to reveal the footprint of positive selection. The Drosophila insulin-like peptides (DILPs) are the ligands that trigger the insulin-signaling cascade. In Drosophila melanogaster, there are several peptides that are structurally similar to the single mammalian insulin peptide. The footprint of recent adaptive changes on nucleotide variation can be unveiled through the analysis of polymorphism and divergence. With this aim, we have surveyed nucleotide sequence variation at the dilp1-7 genes in a natural population of D. melanogaster. The comparison of polymorphism in D. melanogaster and divergence from D. simulans at different functional classes of the dilp genes provided no evidence of adaptive protein evolution after the split of the D. melanogaster and D. simulans lineages. However, our survey of polymorphism at the dilp gene regions of D. melanogaster has provided some evidence for the action of positive selection at or near these genes. The regions encompassing the dilp1-4 genes and the dilp6 gene stand out as likely affected by recent adaptive events.

Comparative genomics of Hedgehog loci in chordates and the origin of Shh regulatory novelties

The origin and evolution of the complex regulatory landscapes of some vertebrate developmental genes, often spanning hundreds of Kbp and including neighboring genes, remain poorly understood. The Sonic Hedgehog (Shh) genomic regulatory block (GRB) is one of the best functionally characterized examples, with several discrete enhancers reported within its introns, vast upstream gene-free region and neighboring genes (Lmbr1 and Rnf32). To investigate the origin and evolution of this GRB, we sequenced and characterized the Hedgehog (Hh) loci from three invertebrate chordate amphioxus species, which share several early expression domains with Shh. Using phylogenetic footprinting within and between chordate lineages, and reporter assays in zebrafish probing >30 Kbp of amphioxus Hh, we report large sequence and functional divergence between both groups. In addition, we show that the linkage of Shh to Lmbr1 and Rnf32, necessary for the unique gnatostomate-specific Shh limb expression, is a vertebrate novelty occurred between the two whole-genome duplications.

Broken colinearity of the amphioxus Hox cluster.

Background In most eumetazoans studied so far, Hox genes determine the identity of structures along the main body axis. They are usually linked in genomic clusters and, in the case of the vertebrate embryo, are expressed with spatial and temporal colinearity. Outside vertebrates, temporal colinearity has been reported in the cephalochordate amphioxus (the least derived living relative of the chordate ancestor) but only for anterior and central genes, namely Hox1 to Hox4 and Hox6. However, most of the Hox gene expression patterns in amphioxus have not been reported. To gain global insights into the evolution of Hox clusters in chordates, we investigated a more extended expression profile of amphioxus Hox genes. Results Here we report an extended expression profile of the European amphioxus Branchiostoma lanceolatum Hox genes and describe that all Hox genes, except Hox13, are expressed during development. Interestingly, we report the breaking of both spatial and temporal colinearity for at least Hox6 and Hox14, which thus have escaped from the classical Hox code concept. We show a previously unidentified Hox6 expression pattern and a faint expression for posterior Hox genes in structures such as the posterior mesoderm, notochord, and hindgut. Unexpectedly, we found that amphioxus Hox14 had the most divergent expression pattern. This gene is expressed in the anterior cerebral vesicle and pharyngeal endoderm. Amphioxus Hox14 expression represents the first report of Hox gene expression in the most anterior part of the central nervous system. Nevertheless, despite these divergent expression patterns, amphioxus Hox6 and Hox14 seem to be still regulated by retinoic acid. Conclusions Escape from colinearity by Hox genes is not unusual in either vertebrates or amphioxus and we suggest that those genes escaping from it are probably associated with the patterning of lineage-specific morphological traits, requiring the loss of those developmental constraints that kept them colinear.

Prototypical Arthropod Gene Content and Genome Organisation in the Centipede Strigamia maritima

Myriapods (e.g., centipedes and millipedes) display a simple homonomous body plan relative to other arthropods. All members of the class are terrestrial, but they attained terrestriality independently of insects. Myriapoda is the only arthropod class not represented by a sequenced genome. We present an analysis of the genome of the centipede Strigamia maritima. It retains a compact genome that has undergone less gene loss and shuffling than previously sequenced arthropods, and many orthologues of genes conserved from the bilaterian ancestor that have been lost in insects. Our analysis locates many genes in conserved macro-synteny contexts, and many small-scale examples of gene clustering. We describe several examples where S. maritima shows different solutions from insects to similar problems. The insect olfactory receptor gene family is absent from S. maritima, and olfaction in air is likely effected by expansion of other receptor gene families. For some genes S. maritima has evolved paralogues to generate coding sequence diversity, where insects use alternate splicing. This is most striking for the Dscam gene, which in Drosophila generates more than 100,000 alternate splice forms, but in S. maritima is encoded by over 100 paralogues. We see an intriguing linkage between the absence of any known photosensory proteins in a blind organism and the additional absence of canonical circadian clock genes. The phylogenetic position of myriapods allows us to identify where in arthropod phylogeny several particular molecular mechanisms and traits emerged. For example, we conclude that juvenile hormone signalling evolved with the emergence of the exoskeleton in the arthropods and that RR-1 containing cuticle proteins evolved in the lineage leading to Mandibulata. We also identify when various gene expansions and losses occurred. The genome of S. maritima offers us a unique glimpse into the ancestral arthropod genome, while also displaying many adaptations to its specific life history.

Analysis of a Plant Complex Resistance Gene Locus Underlying Immune-Related Hybrid Incompatibility and Its Occurrence in Nature

Mechanisms underlying speciation in plants include detrimental (incompatible) genetic interactions between parental alleles that incur a fitness cost in hybrids. We reported on recessive hybrid incompatibility between an Arabidopsis thaliana strain from Poland, Landsberg erecta (Ler), and many Central Asian A. thaliana strains. The incompatible interaction is determined by a polymorphic cluster of Toll/interleukin-1 receptor-nucleotide binding-leucine rich repeat (TNL) RPP1 (Recognition of Peronospora parasitica1)-like genes in Ler and alleles of the receptor-like kinase Strubbelig Receptor Family 3 (SRF3) in Central Asian strains Kas-2 or Kond, causing temperature-dependent autoimmunity and loss of growth and reproductive fitness. Here, we genetically dissected the RPP1-like Ler locus to determine contributions of individual RPP1-like Ler (R1R8) genes to the incompatibility. In a neutral background, expression of most RPP1-like Ler genes, except R3, has no effect on growth or pathogen resistance. Incompatibility involves increased R3 expression and engineered R3 overexpression in a neutral background induces dwarfism and sterility. However, no individual RPP1-like Ler gene is sufficient for incompatibility between Ler and Kas-2 or Kond, suggesting that co-action of at least two RPP1-like members underlies this epistatic interaction. We find that the RPP1-like Ler haplotype is frequent and occurs with other Ler RPP1-like alleles in a local population in Gorzów Wielkopolski (Poland). Only Gorzów individuals carrying the RPP1-like Ler haplotype are incompatible with Kas-2 and Kond, whereas other RPP1-like alleles in the population are compatible. Therefore, the RPP1-like Ler haplotype has been maintained in genetically different individuals at a single site, allowing exploration of forces shaping the evolution of RPP1-like genes at local and regional population scales.

High homologous gene conservation despite extreme autopolyploid redundancy in sugarcane

P>Modern sugarcane (Saccharum spp.) is the leading sugar crop and a primary energy crop. It has the highest level of `vertical` redundancy (2n = 12x = 120) of all polyploid plants studied to date. It was produced about a century ago through hybridization between two autopolyploid species, namely S. officinarum and S. spontaneum. In order to investigate the genome dynamics in this highly polyploid context, we sequenced and compared seven hom(oe)ologous haplotypes (bacterial artificial chromosome clones). Our analysis revealed a high level of gene retention and colinearity, as well as high gene structure and sequence conservation, with an average sequence divergence of 4% for exons. Remarkably, all of the hom(oe)ologous genes were predicted as being functional (except for one gene fragment) and showed signs of evolving under purifying selection, with the exception of genes within segmental duplications. By contrast, transposable elements displayed a general absence of colinearity among hom(oe)ologous haplotypes and appeared to have undergone dynamic expansion in Saccharum, compared with sorghum, its close relative in the Andropogonea tribe. These results reinforce the general trend emerging from recent studies indicating the diverse and nuanced effect of polyploidy on genome dynamics.

Intraspecific evolution of Rhizoctonia solani AG-1 IA associated with soybean and rice in Brazil based on polymorphisms at the ITS-5.8S rDNA operon

Rhizoctonia solani AG-1 IA causes leaf blight on soybean and rice. Despite the fact that R. solani AG-1 IA is a major pathogen affecting soybean and rice in Brazil and elsewhere in the world, little information is available on its genetic diversity and evolution. This study was an attempt to reveal the origin, and the patterns of movement and amplification of epidemiologically significant genotypes of R. solani AG-1 IA from soybean and rice in Brazil. For inferring intraspecific evolution of R. solani AG-1 IA sampled from soybean and rice, networks of ITS-5.8S rDNA sequencing haplotypes were built using the statistical parsimony algorithm from Clement et al. (2000) Molecular Ecology 9: 1657-1660. Higher haplotype diversity (Nei M 1987, Molecular Evolutionary Genetics Columbia University Press, New york: 512p.) was observed for the Brazilian soybean sample of R. solani AG-1 IA (0.827) in comparison with the rest of the world sample (0.431). Within the south-central American clade (3-2), four haplotypes of R. solani AG-1 IA from Mato Grosso, one from Tocantins, one from Maranhao, and one from Cuba occupied the tips of the network, indicating recent origin. The putative ancestral haplotypes had probably originated either from Mato Grosso or Maranhao States. While 16 distinct haplotypes were found in a sample of 32 soybean isolates of the pathogen, the entire rice sample (n=20) was represented by a single haplotype (haplotype 5), with a worldwide distribution. The results from nested-cladistic analysis indicated restricted gene flow with isolation by distance (or restricted dispersal by distance in nonsexual species) for the south-central American clade (3-2), mainly composed by soybean haplotypes.

Male-specific contributions to the Brazilian population of Espirito Santo

Y chromosome markers have been widely studied due to their various applications in the fields of forensic and evolutionary genetics. In this study, 35 Y-SNPs and 17 Y-STRs were genotyped in 253 males from the State of Espirito Santo, Brazil. A total of 18 haplogroups and 243 haplotypes were detected; the haplogroup and haplotype diversities were 0.7794 and 0.9997, respectively. Genetic distance analysis using the Y-STR data showed no statistically significant differences between Espirito Santo and other admixed populations from Brazil. The classification of paternal lineages based on haplogroups showed a predominant European contribution (85.88 %), followed by African (11.37 %) and Amerindian (2.75 %) contributions.