Haplotype heterogeneity and low linkage disequilibrium reduce reliable prediction of genotypes for the ‑α3.7I form of α-thalassaemia using genome-wide microarray data [version 1; peer review: 1 approved, 1 approved with reservations]

Background: The -α3.7I-thalassaemia deletion is very common throughout Africa because it protects against malaria. When undertaking studies to investigate human genetic adaptations to malaria or other diseases, it is important to account for any confounding effects of α-thalassaemia to rule out spurious associations. Methods: In this study we have used direct α-thalassaemia genotyping to understand why GWAS data from a large malaria association study in Kilifi Kenya did not identify the α-thalassaemia signal. We then explored the potential use of a number of new approaches to using GWAS data for imputing α-thalassaemia as an alternative to direct genotyping by PCR. Results: We found very low linkage-disequilibrium of the directly typed data with the GWAS SNP markers around α-thalassaemia and across the haemoglobin-alpha (HBA) gene region, which along with a complex haplotype structure, could explain the lack of an association signal from the GWAS SNP data. Some indirect typing methods gave results that were in broad agreement with those derived from direct genotyping and could identify an association signal, but none were sufficiently accurate to allow correct interpretation compared with direct typing, leading to confusing or erroneous results. Conclusions: We conclude that going forwards, direct typing methods such as PCR will still be required to account for α-thalassaemia in GWAS studies

A Bayesian Approach to Analyse Genetic Variation within RNA Viral Populations

The development of modern and affordable sequencing technologies has allowed the study of viral populations to an unprecedented depth. This is of particular interest for the study of within-host RNA viral populations, where variation due to error-prone polymerases can lead to immune escape, antiviral resistance and adaptation to new host species. Methods to sequence RNA virus genomes include reverse transcription (RT) and polymerase chain reaction (PCR). RT-PCR is a molecular biology technique widely used to amplify DNA from an RNA template. The method itself relies on the in vitro synthesis of copy DNA from RNA followed by multiple cycles of DNA amplification. However, this method introduces artefactual errors that can act as confounding factors when the sequence data are analysed. Although there are a growing number of published studies exploring the intra- and inter-host evolutionary dynamics of RNA viruses, the complexity of the methods used to generate sequences makes it difficult to produce probabilistic statements about the likely sources of observed sequence variants. This complexity is further compounded as both the depth of sequencing and the length of the genome segment of interest increase. Here we develop a Bayesian method to characterise and differentiate between likely structures for the background viral population. This approach can then be used to identify nucleotide sites that show evidence of change in the within-host viral population structure, either over time or relative to a reference sequence (e.g. an inoculum or another source of infection), or both, without having to build complex evolutionary models. Identification of these sites can help to inform the design of more focussed experiments using molecular biology tools, such as site-directed mutagenesis, to assess the function of specific amino acids. We illustrate the method by applying to datasets from experimental transmission of equine influenza, and a pre-clinical vaccine trial for HIV-1

Recombination Drives Vertebrate Genome Contraction

Selective and/or neutral processes may govern variation in DNA content and, ultimately, genome size. The observation in several organisms of a negative correlation between recombination rate and intron size could be compatible with a neutral model in which recombination is mutagenic for length changes. We used whole-genome data on small insertions and deletions within transposable elements from chicken and zebra finch to demonstrate clear links between recombination rate and a number of attributes of reduced DNA content. Recombination rate was negatively correlated with the length of introns, transposable elements, and intergenic spacer and with the rate of short insertions. Importantly, it was positively correlated with gene density, the rate of short deletions, the deletion bias, and the net change in sequence length. All these observations point at a pattern of more condensed genome structure in regions of high recombination. Based on the observed rates of small insertions and deletions and assuming that these rates are representative for the whole genome, we estimate that the genome of the most recent common ancestor of birds and lizards has lost nearly 20% of its DNA content up until the present. Expansion of transposable elements can counteract the effect of deletions in an equilibrium mutation model; however, since the activity of transposable elements has been low in the avian lineage, the deletion bias is likely to have had a significant effect on genome size evolution in dinosaurs and birds, contributing to the maintenance of a small genome. We also demonstrate that most of the observed correlations between recombination rate and genome contraction parameters are seen in the human genome, including for segregating indel polymorphisms. Our data are compatible with a neutral model in which recombination drives vertebrate genome size evolution and gives no direct support for a role of natural selection in this process

Moving out of the shadows: accomplishing bisexual motherhood

Our qualitative study explored the ways in which bisexual mothers came to identify as such and how they structured their relationships and parenting within hetero-patriarchal society. The experiences of seven self-identified White bisexual women (aged from 28 to 56-years-old) from across England and the Republic of Ireland were investigated through semi-structured interviews. Participants’ children were aged 8 months to 28 years old at the time of their interviews. A thematic narrative analysis highlighted the following issues that participants had encountered in constructing their self-identity: prioritizing children; connecting and disconnecting with others and finessing self-definition; questioning societal relationship expectations. Nevertheless, participants varied considerably in how each of the themes identified were reflected in their lives, in particular depending upon each participant’s interpretation of her local social context. Both motherhood and self-identifying as bisexual gave a sense of meaning and purpose to participants’ life stories, although participants sometimes foregrounded their commitment to their children even at a personal cost to their bisexual identity. Using three different theoretical perspectives from feminist theory, queer theory and life course theory, the narratives analysed revealed ways in which bisexual motherhood not only had been influenced both intentionally and unintentionally by heteronormative expectations but also had directly and indirectly challenged these expectations

The GAA triplet-repeat is unstable in the context of the human FXN locus and displays age-dependent expansions in cerebellum and DRG in a transgenic mouse model

Friedreich ataxia (FRDA) is caused by homozygosity for FXN alleles containing an expanded GAA triplet-repeat (GAA-TR) sequence. This expanded GAA-TR sequence is unstable in somatic cells of FRDA patients, showing age-dependent expansions in dorsal root ganglia (DRG), the tissue where pathology occurs earliest and is most significant. This is thought to be the basis for the progressive, tissue-specific pathology seen in FRDA, but the mechanism(s) for this somatic instability is unknown. We show that transgenic mice containing the expanded GAA-TR sequence (190 or 82 triplets) in the context of the human FXN locus show tissue-specific and age-dependent somatic instability that mimics the human condition. Small pool PCR analysis, which allows quantitative analysis of instability by assaying individual transgenes in vivo, showed age-dependent expansions specifically in the cerebellum and DRG. The (GAA)190 allele showed some instability by 2 months, progressed at about 0.3 – 0.4 triplets/week, resulting in a significant number of expansions by 12 months. Repeat length determined the age of onset of somatic instability, and the rate and magnitude of expansion. Whereas the GAA-TR was unstable in the context of the human FXN locus, pure GAATR sequences at other genetic loci in the human and murine genomes showed no instability. These data indicate that somatic instability of the GAA-TR sequence in the human FXN gene is determined by a combination of unique cis and trans-acting factors. This mouse model will serve as a useful tool to delineate the mechanism(s) of diseasespecific somatic instability in FRDA

Rare Exonic Minisatellite Alleles in MUC2 Influence Susceptibility to Gastric Carcinoma

BACKGROUND: Mucins are the major components of mucus and their genes share a common, centrally-located region of sequence that encodes tandem repeats. Mucins are well known genes with respect to their specific expression levels; however, their genomic levels are unclear because of complex genomic properties. In this study, we identified eight novel minisatellites from the entire MUC2 region and investigated how allelic variation in these minisatellites may affect susceptibility to gastrointestinal cancer. METHODOLOGY/PRINCIPLE FINDINGS: We analyzed genomic DNA from the blood of normal healthy individuals and multi-generational family groups. Six of the eight minisatellites exhibited polymorphism and were transmitted meiotically in seven families, following Mendelian inheritance. Furthermore, a case-control study was performed that compared genomic DNA from 457 cancer-free controls with DNA from individuals with gastric (455), colon (192) and rectal (271) cancers. A statistically significant association was identified between rare exonic MUC2-MS6 alleles and the occurrence of gastric cancer: odds ratio (OR), 2.56; 95% confidence interval (CI), 1.31-5.04; and p = 0.0047. We focused on an association between rare alleles and gastric cancer. Rare alleles were divided into short (40, 43 and 44) and long (47, 50 and 54), according to their TR (tandem repeats) lengths. Interestingly, short rare alleles were associated with gastric cancer (OR = 5.6, 95% CI: 1.93-16.42; p = 0.00036). Moreover, hypervariable MUC2 minisatellites were analyzed in matched blood and cancer tissue from 28 patients with gastric cancer and in 4 cases of MUC2-MS2, minisatellites were found to have undergone rearrangement. CONCLUSIONS/SIGNIFICANCE: Our observations suggest that the short rare MUC2-MS6 alleles could function as identifiers for risk of gastric cancer. Additionally, we suggest that minisatellite instability might be associated with MUC2 function in cancer cells

Genetic diversity of the African malaria vector Anopheles gambiae

The sustainability of malaria control in Africa is threatened by the rise of insecticide resistance in Anopheles mosquitoes, which transmit the disease1. To gain a deeper understanding of how mosquito populations are evolving, here we sequenced the genomes of 765 specimens of Anopheles gambiae and Anopheles coluzzii sampled from 15 locations across Africa, and identified over 50 million single nucleotide polymorphisms within the accessible genome. These data revealed complex population structure and patterns of gene flow, with evidence of ancient expansions, recent bottlenecks, and local variation in effective population size. Strong signals of recent selection were observed in insecticide-resistance genes, with several sweeps spreading over large geographical distances and between species. The design of new tools for mosquito control using gene-drive systems will need to take account of high levels of genetic diversity in natural mosquito populations

The Impact of Recombination on Nucleotide Substitutions in the Human Genome

Unraveling the evolutionary forces responsible for variations of neutral substitution patterns among taxa or along genomes is a major issue for detecting selection within sequences. Mammalian genomes show large-scale regional variations of GC-content (the isochores), but the substitution processes at the origin of this structure are poorly understood. We analyzed the pattern of neutral substitutions in 1 Gb of primate non-coding regions. We show that the GC-content toward which sequences are evolving is strongly negatively correlated to the distance to telomeres and positively correlated to the rate of crossovers (R2 = 47%). This demonstrates that recombination has a major impact on substitution patterns in human, driving the evolution of GC-content. The evolution of GC-content correlates much more strongly with male than with female crossover rate, which rules out selectionist models for the evolution of isochores. This effect of recombination is most probably a consequence of the neutral process of biased gene conversion (BGC) occurring within recombination hotspots. We show that the predictions of this model fit very well with the observed substitution patterns in the human genome. This model notably explains the positive correlation between substitution rate and recombination rate. Theoretical calculations indicate that variations in population size or density in recombination hotspots can have a very strong impact on the evolution of base composition. Furthermore, recombination hotspots can create strong substitution hotspots. This molecular drive affects both coding and non-coding regions. We therefore conclude that along with mutation, selection and drift, BGC is one of the major factors driving genome evolution. Our results also shed light on variations in the rate of crossover relative to non-crossover events, along chromosomes and according to sex, and also on the conservation of hotspot density between human and chimp