African Genome Research Articles

Evolution of genomic structural variation and genomic architecture in the adaptive radiations of African cichlid fishes.

African cichlid fishes are an ideal system for studying explosive rates of speciation and the origin of diversity in adaptive radiation. Within the last few million years, more than 2000 species have evolved in the Great Lakes of East Africa, the largest adaptive radiation in vertebrates. These young species show spectacular diversity in their coloration, morphology and behavior. However, little is known about the genomic basis of this astonishing diversity. Recently, five African cichlid genomes were sequenced, including that of the Nile Tilapia (Oreochromis niloticus), a basal and only relatively moderately diversified lineage, and the genomes of four representative endemic species of the adaptive radiations, Neolamprologus brichardi, Astatotilapia burtoni, Metriaclima zebra, and Pundamila nyererei. Using the Tilapia genome as a reference genome, we generated a high-resolution genomic variation map, consisting of single nucleotide polymorphisms (SNPs), short insertions and deletions (indels), inversions and deletions. In total, around 18.8, 17.7, 17.0, and 17.0 million SNPs, 2.3, 2.2, 1.4, and 1.9 million indels, 262, 306, 162, and 154 inversions, and 3509, 2705, 2710, and 2634 deletions were inferred to have evolved in N. brichardi, A. burtoni, P. nyererei, and M. zebra, respectively. Many of these variations affected the annotated gene regions in the genome. Different patterns of genetic variation were detected during the adaptive radiation of African cichlid fishes. For SNPs, the highest rate of evolution was detected in the common ancestor of N. brichardi, A. burtoni, P. nyererei, and M. zebra. However, for the evolution of inversions and deletions, we found that the rates at the terminal taxa are substantially higher than the rates at the ancestral lineages. The high-resolution map provides an ideal opportunity to understand the genomic bases of the adaptive radiation of African cichlid fishes.

Increased genetic diversity of ADME genes in African Americans compared with their putative ancestral source populations and implications for pharmacogenomics.

BackgroundAfrican Americans have been treated as a representative population for African ancestry for many purposes, including pharmacogenomic studies. However, the contribution of European ancestry is expected to result in considerable differences in the genetic architecture of African American individuals compared with an African genome. In particular, the genetic admixture influences the genomic diversity of drug metabolism-related genes, and may cause high heterogeneity of drug responses in admixed populations such as African Americans.ResultsThe genomic ancestry information of African-American (ASW) samples was obtained from data of the 1000 Genomes Project, and local ancestral components were also extracted for 32 core genes and 252 extended genes, which are associated with drug absorption, distribution, metabolism, and excretion (ADME) genes. As expected, the global genetic diversity pattern in ASW was determined by the contributions of its putative ancestral source populations, and the whole profiles of ADME genes in ASW are much closer to those in YRI than in CEU. However, we observed much higher diversity in some functionally important ADME genes in ASW than either CEU or YRI, which could be a result of either genetic drift or natural selection, and we identified some signatures of the latter. We analyzed the clinically relevant polymorphic alleles and haplotypes, and found that 28 functional mutations (including 3 missense, 3 splice, and 22 regulator sites) exhibited significantly higher differentiation between the three populations.ConclusionsAnalysis of the genetic diversity of ADME genes showed differentiation between admixed population and its ancestral source populations. In particular, the different genetic diversity between ASW and YRI indicated that the ethnic differences in pharmacogenomic studies are broadly existed despite that African ancestry is dominant in Africans Americans. This study should advance our understanding of the genetic basis of the drug response heterogeneity between populations, especially in the case of population admixture, and have significant implications for evaluating potential inter-population heterogeneity in drug treatment effects.

A targeted next‐generation sequencing toolkit for exon‐based cichlid phylogenomics

Cichlid fishes (family Cichlidae) are models for evolutionary and ecological research. Massively parallel sequencing approaches have been successfully applied to study relatively recent diversification in groups of African and Neotropical cichlids, but such technologies have yet to be used for addressing larger-scale phylogenetic questions of cichlid evolution. Here, we describe a process for identifying putative single-copy exons from five African cichlid genomes and sequence the targeted exons for a range of divergent (>tens of millions of years) taxa with probes designed from a single reference species (Oreochromis niloticus, Nile tilapia). Targeted sequencing of 923 exons across 10 cichlid species that represent the family's major lineages and geographic distribution resulted in a complete taxon matrix of 564 exons (649,549 bp), representing 559 genes. Maximum likelihood and Bayesian analyses in both species tree and concatenation frameworks yielded the same fully resolved and highly supported topology, which matched the expected backbone phylogeny of the major cichlid lineages. This work adds to the body of evidence that it is possible to use a relatively divergent reference genome for exon target design and successful capture across a broad phylogenetic range of species. Furthermore, our results show that the use of a third-party laboratory coupled with accessible bioinformatics tools makes such phylogenomics projects feasible for research groups that lack direct access to genomic facilities. We expect that these resources will be used in further cichlid evolution studies and hope the protocols and identified targets will also be useful for phylogenetic studies of a wider range of organisms.

Coelacanths

Coelacanths

Reassessment of the Genome Size in Elaeis guineensis and Elaeis oleifera, and Its Interspecific Hybrid.

Aiming at generating a comprehensive genomic database on Elaeis spp., our group is leading several R&D initiatives with Elaeis guineensis (African oil palm) and Elaeis oleifera (American oil palm), including the whole-genome sequencing of the last. Genome size estimates currently available for this genus are controversial, as they indicate that American oil palm genome is about half the size of the African oil palm genome and that the genome of the interspecific hybrid is bigger than both the parental species genomes. We estimated the genome size of three E. guineensis genotypes, five E. oleifera genotypes, and two interspecific hybrids genotypes. On average, the genome size of E. guineensis is 4.32 ± 0.173 pg, while that of E. oleifera is 4.43 ± 0.018 pg. This indicates that both genomes are similar in size, even though E. oleifera is in fact bigger. As expected, the hybrid genome size is around the average of the two genomes, 4.40 ± 0.016 pg. Additionally, we demonstrate that both species present around 38% of GC content. As our results contradict the currently available data on Elaeis spp. genome sizes, we propose that the actual genome size of the Elaeis species is around 4 pg and that American oil palm possesses a larger genome than African oil palm.

The complete sequence of the mitochondrial genome of the African Penguin (Spheniscus demersus)

The complete mitochondrial genome of the African Penguin (Spheniscus demersus) was sequenced. The molecule was sequenced via next generation sequencing and primer walking. The size of the genome is 17,346bp in length. Comparison with the mitochondrial DNA of two other penguin genomes that have so far been reported was conducted namely; Little blue penguin (Eudyptula minor) and the Rockhopper penguin (Eudyptes chrysocome). This analysis made it possible to identify common penguin mitochondrial DNA characteristics. The S. demersus mtDNA genome is very similar, both in composition and length to both the E. chrysocome and E. minor genomes. The gene content of the African penguin mitochondrial genome is typical of vertebrates and all three penguin species have the standard gene order originally identified in the chicken. The control region for S. demersus is located between tRNA-Glu and tRNA-Phe and all three species of penguins contain two sets of similar repeats with varying copy numbers towards the 3′ end of the control region, accounting for the size variance. This is the first report of the complete nucleotide sequence for the mitochondrial genome of the African penguin, S. demersus. These results can be subsequently used to provide information for penguin phylogenetic studies and insights into the evolution of genomes.

Genetic diversity in black South Africans from Soweto

BackgroundDue to the unparalleled genetic diversity of its peoples, Africa is attracting growing research attention. Several African populations have been assessed in global initiatives such as the International HapMap and 1000 Genomes Projects. Notably excluded, however, is the southern Africa region, which is inhabited predominantly by southeastern Bantu-speakers, currently suffering under the dual burden of infectious and non-communicable diseases. Limited reference data for these individuals hampers medical research and prevents thorough understanding of the underlying population substructure. Here, we present the most detailed exploration, to date, of genetic diversity in 94 unrelated southeastern Bantu-speaking South Africans, resident in urban Soweto (Johannesburg).ResultsParticipants were typed for ~4.3 million SNPs using the Illumina Omni5 beadchip. PCA and ADMIXTURE plots were used to compare the observed variation with that seen in selected populations worldwide. Results indicated that Sowetans, and other southeastern Bantu-speakers, are a clearly distinct group from other African populations previously investigated, reflecting a unique genetic history with small, but significant contributions from diverse sources. To assess the suitability of our sample as representative of Sowetans, we compared our results to participants in a larger rheumatoid arthritis case–control study. The control group showed good clustering with our sample, but among the cases were individuals who demonstrated notable admixture.ConclusionsSowetan population structure appears unique compared to other black Africans, and may have clinical implications. Our data represent a suitable reference set for southeastern Bantu-speakers, on par with a HapMap type reference population, and constitute a prelude to the Southern African Human Genome Programme.

Complete mitochondrial DNA sequences of the threadfin cichlid (Petrochromis trewavasae) and the blunthead cichlid (Tropheus moorii) and patterns of mitochondrial genome evolution in cichlid fishes.

The cichlid fishes of the East African Great Lakes represent a model especially suited to study adaptive radiation and speciation. With several African cichlid genome projects being in progress, a promising set of closely related genomes is emerging, which is expected to serve as a valuable data base to solve questions on genotype-phenotype relations. The mitochondrial (mt) genomes presented here are the first results of the assembly and annotation process for two closely related but eco-morphologically highly distinct Lake Tanganyika cichlids, Petrochromis trewavasae and Tropheus moorii. The genomic sequences comprise 16,588 bp (P. trewavasae) and 16,590 bp (T. moorii), and exhibit the typical mitochondrial structure, with 13 protein-coding genes, 2 rRNA genes, 22 tRNA genes, and a non-coding control region. Analyses confirmed that the two species are very closely related with an overall sequence similarity of 96%. We analyzed the newly generated sequences in the phylogenetic context of 21 published labroid fish mitochondrial genomes. Consistent with other vertebrates, the D-loop region was found to evolve faster than protein-coding genes, which in turn are followed by the rRNAs; the tRNAs vary greatly in the rate of sequence evolution, but on average evolve the slowest. Within the group of coding genes, ND6 evolves most rapidly. Codon usage is similar among examined cichlid tribes and labroid families; although a slight shift in usage patterns down the gene tree could be observed. Despite having a clearly different nucleotide composition, ND6 showed a similar codon usage. C-terminal ends of Cox1 exhibit variations, where the varying number of amino acids is related to the structure of the obtained phylogenetic tree. This variation may be of functional relevance for Cox1 synthesis.

SP0108 Next Generation Sequencing and Related Technologies, A Way Forward for Arthritis?

The publication of the human genome sequence in 2001 represented one of the major milestones in the history of biology and much effort has since been devoted to describing and...

Exploratory Study of Apolipoprotein E ε4 Genotype and Risk of Alzheimer's Disease in Mexican Hispanics

To the Editor: Hispanics are 1.5 times more likely to develop Alzheimer's disease (AD) than White non-Hispanics1. A community-based study reported that the prevalence of AD is higher in Hispanics than in White non-Hispanics2. Interestingly, a recent study found that the average proportion of Amerindian genetic ancestry in AD cases was a third of that observed in the cognitively intact controls3. These findings suggest that Amerindian ancestry confers a protective effect against disease onset. Hispanics are a mixture of European, African, and Amerindian genomes, with percentages depending on the country of origin4. Individuals of Mexican origin generally are of Amerindian ancestry, while those of Caribbean origin, for example, have more African ancestry5,6. Therefore, the question arises as to whether Mexican Hispanics differ in their susceptibility to AD. The Apolipoprotein E (APOE) gene, with its three common alleles (e2, e3, e4), is involved in cholesterol transport and other biologic pathways potentially relevant to AD7 and has been consistently associated with an increased risk of AD8. Studies indicate that the APOE genetic risk variants may not carry the same risk for Hispanics as for White non-Hispanics. A meta-analysis found an association between the APOE-e4 allele and AD, although the risk conveyed by the APOE-e4 allele frequency was found to be lower in Hispanics, predominately Caribbean origin, than in White non-Hispanics9. Another study of Hispanic Americans of predominantly Caribbean origin showed a smaller relative risk of AD conveyed by the APOE-e4 allele in Hispanics than in White non-Hispanics10. To our knowledge, Mexican Hispanics have not been studied as a single group to determine the associated risk between the APOE-e4 allele and AD onset.

The African coelacanth genome provides insights into tetrapod evolution

It was a zoological sensation when a living specimen of the coelacanth was first discovered in 1938, as this lineage of lobe-finned fish was thought to have gone extinct 70 million years ago. The modern coelacanth looks remarkably similar to many of its ancient relatives, and its evolutionary proximity to our own fish ancestors provides a glimpse of the fish that first walked on land. Here we report the genome sequence of the African coelacanth, Latimeria chalumnae. Through a phylogenomic analysis, we conclude that the lungfish, and not the coelacanth, is the closest living relative of tetrapods. Coelacanth protein-coding genes are significantly more slowly evolving than those of tetrapods, unlike other genomic features . Analyses of changes in genes and regulatory elements during the vertebrate adaptation to land highlight genes involved in immunity, nitrogen excretion and the development of fins, tail, ear, eye, brain, and olfaction. Functional assays of enhancers involved in the fin-to-limb transition and in the emergence of extra-embryonic tissues demonstrate the importance of the coelacanth genome as a blueprint for understanding tetrapod evolution.

Population Genetic Inference from Personal Genome Data: Impact of Ancestry and Admixture on Human Genomic Variation

Full sequencing of individual human genomes has greatly expanded our understanding of human genetic variation and population history. Here, we present a systematic analysis of 50 human genomes from 11 diverse global populations sequenced at high coverage. Our sample includes 12 individuals who have admixed ancestry and who have varying degrees of recent (within the last 500 years) African, Native American, and European ancestry. We found over 21 million single-nucleotide variants that contribute to a 1.75-fold range in nucleotide heterozygosity across diverse human genomes. This heterozygosity ranged from a high of one heterozygous site per kilobase in west African genomes to a low of 0.57 heterozygous sites per kilobase in segments inferred to have diploid Native American ancestry from the genomes of Mexican and Puerto Rican individuals. We show evidence of all three continental ancestries in the genomes of Mexican, Puerto Rican, and African American populations, and the genome-wide statistics are highly consistent across individuals from a population once ancestry proportions have been accounted for. Using a generalized linear model, we identified subtle variations across populations in the proportion of neutral versus deleterious variation and found that genome-wide statistics vary in admixed populations even once ancestry proportions have been factored in. We further infer that multiple periods of gene flow shaped the diversity of admixed populations in the Americas-70% of the European ancestry in today's African Americans dates back to European gene flow happening only 7-8 generations ago.

Sequential deletion of genes from the African swine fever virus genome using the cre/loxP recombination system

A method has been established to sequentially delete combinations of genes from the ASFV genome to test the effect on virus replication and host responses to infection. Initially the ASFV genes MGF505 2R and MGF505 3R and a truncated MGF360 9L gene were deleted from the genome of the tissue-culture adapted ASFV strain BA71V and replaced with bacteriophage loxP sequences flanking the beta-glucuronidase (GUS) marker gene to create recombinant virus VΔMGF-GUS. Subsequently the GUS gene was removed by site-specific recombination between the two loxP sites involving expression of the bacteriophage Cre recombinase enzyme to create recombinant virus VΔMGFΔGUS. The EP402R and EP153R genes were subsequently deleted from the genome of VΔMGFΔGUS, using the same GUS marker gene, to construct virus VΔMGFΔCD2-Lectin-GUS. These sequential deletions of ASFV genes were shown not to alter virus replication significantly.

Expression sequence tag library derived from peripheral blood mononuclear cells of the chlorocebus sabaeus

BackgroundAfrican Green Monkeys (AGM) are amongst the most frequently used nonhuman primate models in clinical and biomedical research, nevertheless only few genomic resources exist for this species. Such information would be essential for the development of dedicated new generation technologies in fundamental and pre-clinical research using this model, and would deliver new insights into primate evolution.ResultsWe have exhaustively sequenced an Expression Sequence Tag (EST) library made from a pool of Peripheral Blood Mononuclear Cells from sixteen Chlorocebus sabaeus monkeys. Twelve of them were infected with the Simian Immunodeficiency Virus. The mononuclear cells were or not stimulated in vitro with Concanavalin A, with lipopolysacharrides, or through mixed lymphocyte reaction in order to generate a representative and broad library of expressed sequences in immune cells. We report here 37,787 sequences, which were assembled into 14,410 contigs representing an estimated 12% of the C. sabaeus transcriptome. Using data from primate genome databases, 9,029 assembled sequences from C. sabaeus could be annotated. Sequences have been systematically aligned with ten cDNA references of primate species including Homo sapiens, Pan troglodytes, and Macaca mulatta to identify ortholog transcripts. For 506 transcripts, sequences were quasi-complete. In addition, 6,576 transcript fragments are potentially specific to the C. sabaeus or corresponding to not yet described primate genes.ConclusionsThe EST library we provide here will prove useful in gene annotation efforts for future sequencing of the African Green Monkey genomes. Furthermore, this library, which particularly well represents immunological and hematological gene expression, will be an important resource for the comparative analysis of gene expression in clinically relevant nonhuman primate and human research.

Genomic sovereignty and the African promise: mining the African genome for the benefit of Africa

Scientific interest in genomics in Africa is on the rise with a number of funding initiatives aimed specifically at supporting research in this area. Genomics research on material of African...

Genomics in Africa: Avoiding Past Pitfalls

A landmark genomics project is taking shape in Africa that shifts the power and prominence to local scientists. If successful, the program will offer valuable insights into the inheritance of common diseases and reshape the paradigm of foreign-funded research.

Biomedical Cloud Computing With Amazon Web Services

In this overview to biomedical computing in the cloud, we discussed two primary ways to use the cloud (a single instance or cluster), provided a detailed example using NGS mapping, and highlighted the associated costs. While many users new to the cloud may assume that entry is as straightforward as uploading an application and selecting an instance type and storage options, we illustrated that there is substantial up-front effort required before an application can make full use of the cloud's vast resources. Our intention was to provide a set of best practices and to illustrate how those apply to a typical application pipeline for biomedical informatics, but also general enough for extrapolation to other types of computational problems. Our mapping example was intended to illustrate how to develop a scalable project and not to compare and contrast alignment algorithms for read mapping and genome assembly. Indeed, with a newer aligner such as Bowtie, it is possible to map the entire African genome using one m2.2xlarge instance in 48 hours for a total cost of approximately $48 in computation time. In our example, we were not concerned with data transfer rates, which are heavily influenced by the amount of available bandwidth, connection latency, and network availability. When transferring large amounts of data to the cloud, bandwidth limitations can be a major bottleneck, and in some cases it is more efficient to simply mail a storage device containing the data to AWS (http://aws.amazon.com/importexport/). More information about cloud computing, detailed cost analysis, and security can be found in references.

Structural variation in two human genomes mapped at single-nucleotide resolution by whole genome de novo assembly.

Here we use whole-genome de novo assembly of second-generation sequencing reads to map structural variation (SV) in an Asian genome and an African genome. Our approach identifies small- and intermediate-size homozygous variants (1-50 kb) including insertions, deletions, inversions and their precise breakpoints, and in contrast to other methods, can resolve complex rearrangements. In total, we identified 277,243 SVs ranging in length from 1-23 kb. Validation using computational and experimental methods suggests that we achieve overall <6% false-positive rate and <10% false-negative rate in genomic regions that can be assembled, which outperforms other methods. Analysis of the SVs in the genomes of 106 individuals sequenced as part of the 1000 Genomes Project suggests that SVs account for a greater fraction of the diversity between individuals than do single-nucleotide polymorphisms (SNPs). These findings demonstrate that whole-genome de novo assembly is a feasible approach to deriving more comprehensive maps of genetic variation.

Frequency of HLA B∗5701 allele carriers in abacavir treated-HIV infected patients and controls from northeastern Brazil

In clinical practice, using pharmacogenetics has risen in importance since the technologies for genomic variation searches for different responses to drugs became widespread, easier, more rapid, and affordable for a specialized laboratory. The main pharmacogenetics studies have mostly involved drug resistance monitoring in oncology patients;1 however, other applications are being developed to detect genetic molecular markers associated with resistance or hypersensitivity/adverse reactions to antiretroviral drug treatment in Human Immunodeficiency Virus (HIV)-infected patients.2 HIV treatment is known to be limited by adverse drug reactions and the development of resistance. One significant example is the strong association of the Abacavir hypersensitivity reaction with HLA-B*5701 in HIV-positive patients.3 Abacavir, a common drug for treating HIV-infected patients, is an efficient nucleoside reverse-transcriptase inhibitor with a beneficial long-term toxicity profile, often used with other antiretroviral agents. However, Abacavir can yield adverse effects such as immunologically-driven hypersensitivity in 5 to 8% of HIV-positive subjects during the first six weeks of use.2 Hypersensitivity symptoms immediately reverse after the interruption of Abacavir.4 The hypersensitivity reaction to Abacavir has been reportedly associated with the presence of the major histocompatibility complex class I allele HLA-B*5701 and this association has been confirmed in several replication studies in different ethnic groups of HIV-positive patients.3 In our study, we searched for the presence of HLA B*5701 in 96 HIV-positive patients treated with Abacavir and in 243 healthy subjects from Northeastern Brazil (Recife, Pernambuco) to verify the percentage of HLA B*5701 allele carriers in HIV patients and in the general population from Northeast Brazil. This area is known to harbor a tri-hybrid population resulting from contributing African (44%), Caucasian (34%), and native American (22%) genomes.5

Origin and taxonomic status of the Palearctic population of the stem borer Sesamia nonagrioides (Lefèbvre) (Lepidoptera: Noctuidae)

The major pest of maize in Mediterranean Europe, the stem borer Sesamia nonagrioides (Lefèbvre) (Lepidoptera: Noctuidae), has a fragmented distribution, north and south of the Sahara. The present study aimed: (1) to clarify the uncertain taxonomic status of the Palearctic and sub-Saharan populations which were first considered as different species and later on as subspecies (Sesamia nonagrioides nonagrioides and Sesamia nonagrioides botanephaga) and (2) to investigate the origin of the Palearctic population which extends from Spain to Iran, outside what is considered typical for this mainly tropical genus. We reconstructed the evolutionary history of both populations using one nuclear and two mitochondrial genes. The sub-Saharan taxon was fragmented in two isolated populations (West and East) whose mitochondrial genes were distant by 2.3%. The Palearctic population was included in the East African clade and its genes were close or identical to those of a population from Central Ethiopia, where the species was discovered for the first time. Similarly, in Africa, the alleles of the nuclear gene were distributed mainly in two West and East clades, whereas some Palearctic alleles belonged to the West clade. The Palearctic population originated therefore from East and West Africa and is the progeny of the cross between these two African populations. The main species concepts were in agreement, leading to the conclusion that the three populations are still conspecific. In the surveyed regions, the species therefore does not include two subspecies but three isolated populations. The Palearctic population suffered from severe bottlenecks that resulted in the fixation of one East African mitochondrial genome and the large reduction in its genetic diversity compared to the African populations. The data suggest that natural colonization of the Palearctic region was more plausible than human introduction. The allelic distribution of the Palearctic population was similar to that of species that survived the last glaciation. It is concluded that the African populations expanded during the last interglacial, crossed the Sahara and mixed in North Africa where fixation of the East mitochondrial genome occurred. The species then colonized Europe westward through only one eastern entrance. The coalescent-based estimate of the time to the ancestor of the Palearctic population was 108 000 years, which is consistent with this scenario.