Posttranscriptional regulation of FAN1 by miR-124-3p at rs3512 underlies onset-delaying genetic modification in Huntington's disease

Many Mendelian disorders, such as Huntington's disease (HD) and spinocerebellar ataxias, arise from expansions of CAG trinucleotide repeats. Despite the clear genetic causes, additional genetic factors may influence the rate of those monogenic disorders. Notably, genome-wide association studies discovered somewhat expected modifiers, particularly mismatch repair genes involved in the CAG repeat instability, impacting age at onset of HD. Strikingly, FAN1, previously unrelated to repeat instability, produced the strongest HD modification signals. Diverse FAN1 haplotypes independently modify HD, with rare genetic variants diminishing DNA binding or nuclease activity of the FAN1 protein, hastening HD onset. However, the mechanism behind the frequent and the most significant onset-delaying FAN1 haplotype lacking missense variations has remained elusive. Here, we illustrated that a microRNA acting on 3'-UTR (untranslated region) SNP rs3512, rather than transcriptional regulation, is responsible for the significant FAN1 expression quantitative trait loci signal and allelic imbalance in FAN1 messenger ribonucleic acid (mRNA), accounting for the most significant and frequent onset-delaying modifier haplotype in HD. Specifically, miR-124-3p selectively targets the reference allele at rs3512, diminishing the stability of FAN1 mRNA harboring that allele and consequently reducing its levels. Subsequent validation analyses, including the use of antagomir and 3'-UTR reporter vectors with swapped alleles, confirmed the specificity of miR-124-3p at rs3512. Together, these findings indicate that the alternative allele at rs3512 renders the FAN1 mRNA less susceptible to miR-124-3p-mediated posttranscriptional regulation, resulting in increased FAN1 levels and a subsequent delay in HD onset by mitigating CAG repeat instability.

Somatic CAG repeat instability in intermediate alleles of the HTT gene and its potential association with a clinical phenotype

Huntington disease (HD) is a neurodegenerative disorder caused by ≥36 CAGs in the HTT gene. Intermediate alleles (IAs) (27-35 CAGs) are not considered HD-causing, but their potential association with neurocognitive symptoms remains controversial. As HTT somatic CAG expansion influences HD onset, we hypothesised that IAs are somatically unstable, and that somatic CAG expansion may drive phenotypic presentation in some IA carriers. We quantified HTT somatic CAG expansions by MiSeq sequencing in the blood DNA of 164 HD subjects and 191 IA (symptomatic and control) carriers, and in the brain DNA of a symptomatic 33 CAG carrier. We also performed genotype-phenotype analysis. The phenotype of symptomatic IA carriers was characterised by motor (85%), cognitive (27%) and/or behavioural (29%) signs, with a late (58.7 ± 18.6 years), but not CAG-dependent, age at onset. IAs displayed somatic expansion that were CAG and age-dependent in blood DNA, with 0.4% and 0.01% of DNA molecules expanding by CAG and year, respectively. Somatic expansions of +1 and +2 CAGs were detected in the brain of the individual with 33 CAGs, with the highest expansion frequency in the putamen (10.3%) and the lowest in the cerebellum (4.8%). Somatic expansion in blood DNA was not different in symptomatic vs. control IA carriers. In conclusion, we show that HTT IAs are somatically unstable, but we found no association with HD-like phenotypes. It is plausible, however, that some IAs, close to the HD pathological threshold and with a predisposing genetic background, could manifest with neurocognitive symptoms.

Modification of Huntington's disease by short tandem repeats

Expansions of glutamine-coding CAG trinucleotide repeats cause a number of neurodegenerative diseases, including Huntington's disease and several of spinocerebellar ataxias. In general, age-at-onset of the polyglutamine diseases is inversely correlated with the size of the respective inherited expanded CAG repeat. Expanded CAG repeats are also somatically unstable in certain tissues, and age-at-onset of Huntington's disease corrected for individual HTT CAG repeat length (i.e. residual age-at-onset), is modified by repeat instability-related DNA maintenance/repair genes as demonstrated by recent genome-wide association studies. Modification of one polyglutamine disease (e.g. Huntington's disease) by the repeat length of another (e.g. ATXN3, CAG expansions in which cause spinocerebellar ataxia 3) has also been hypothesized. Consequently, we determined whether age-at-onset in Huntington's disease is modified by the CAG repeats of other polyglutamine disease genes. We found that the CAG measured repeat sizes of other polyglutamine disease genes that were polymorphic in Huntington's disease participants but did not influence Huntington's disease age-at-onset. Additional analysis focusing specifically on ATXN3 in a larger sample set (n = 1388) confirmed the lack of association between Huntington's disease residual age-at-onset and ATXN3 CAG repeat length. Additionally, neither our Huntington's disease onset modifier genome-wide association studies single nucleotide polymorphism data nor imputed short tandem repeat data supported the involvement of other polyglutamine disease genes in modifying Huntington's disease. By contrast, our genome-wide association studies based on imputed short tandem repeats revealed significant modification signals for other genomic regions. Together, our short tandem repeat genome-wide association studies show that modification of Huntington's disease is associated with short tandem repeats that do not involve other polyglutamine disease-causing genes, refining the landscape of Huntington's disease modification and highlighting the importance of rigorous data analysis, especially in genetic studies testing candidate modifiers.

Spanish HTT gene study reveals haplotype and allelic diversity with possible implications for germline expansion dynamics in Huntington disease

We aimed to determine the genetic diversity and molecular characteristics of the Huntington disease (HD) gene (HTT) in Spain. We performed an extended haplotype and exon one deep sequencing analysis of the HTT gene in a nationwide cohort of population-based controls (n = 520) and families with symptomatic individuals referred for HD genetic testing. This group included 331 HD cases and 140 carriers of intermediate alleles. Clinical and family history data were obtained when available. Spanish normal alleles are enriched in C haplotypes (40.1%), whereas A1 (39.8%) and A2 (31.6%) prevail among intermediate and expanded alleles, respectively. Alleles ≥ 50 CAG repeats are primarily associated with haplotypes A2 (38.9%) and C (32%), which are also present in 50% and 21.4%, respectively, of HD families with large intergenerational expansions. Non-canonical variants of exon one sequence are less frequent, but much more diverse, in alleles of ≥27 CAG repeats. The deletion of CAACAG, one of the six rare variants not observed among smaller normal alleles, is associated with haplotype C and appears to correlate with larger intergenerational expansions and early onset of symptoms. Spanish HD haplotypes are characterized by a high genetic diversity, potentially admixed with other non-Caucasian populations, with a higher representation of A2 and C haplotypes than most European populations. Differences in haplotype distributions across the CAG length range support differential germline expansion dynamics, with A2 and C showing the largest intergenerational expansions. This haplotype-dependent germline instability may be driven by specific cis-elements, such as the CAACAG deletion.

Repeat Detector: versatile sizing of expanded tandem repeats and identification of interrupted alleles from targeted DNA sequencing

Targeted DNA sequencing approaches will improve how the size of short tandem repeats is measured for diagnostic tests and preclinical studies. The expansion of these sequences causes dozens of disorders, with longer tracts generally leading to a more severe disease. Interrupted alleles are sometimes present within repeats and can alter disease manifestation. Determining repeat size mosaicism and identifying interruptions in targeted sequencing datasets remains a major challenge. This is in part because standard alignment tools are ill-suited for repetitive and unstable sequences. To address this, we have developed Repeat Detector (RD), a deterministic profile weighting algorithm for counting repeats in targeted sequencing data. We tested RD using blood-derived DNA samples from Huntington's disease and Fuchs endothelial corneal dystrophy patients sequenced using either Illumina MiSeq or Pacific Biosciences single-molecule, real-time sequencing platforms. RD was highly accurate in determining repeat sizes of 609 blood-derived samples from Huntington's disease individuals and did not require prior knowledge of the flanking sequences. Furthermore, RD can be used to identify alleles with interruptions and provide a measure of repeat instability within an individual. RD is therefore highly versatile and may find applications in the diagnosis of expanded repeat disorders and in the development of novel therapies.

A probable cis-acting genetic modifier of Huntington disease frequent in individuals with African ancestry

Huntington disease (HD)is a dominantly inherited neurodegenerative disorder caused by the expansion of a polyglutamine encoding CAG repeat in the huntingtin gene. Recently, it has been established that disease severity in HD is best predicted by the number of pure CAG repeats rather than total glutamines encoded. Along with uncovering DNA repair gene variants as trans-acting modifiers of HD severity, these data reveal somatic expansion of the CAG repeat as a key driver of HD onset. Using high-throughput DNA sequencing, we have determined the precise sequence and somatic expansion profiles of the HTT repeat tract of 68 HD-affected and 158 HD-unaffected African ancestry individuals. A high level of HTT repeat sequence diversity was observed, with three likely African-specific alleles identified. In the most common disease allele (30 out of 68), the typical proline-encoding CCGCCA sequence was absent. This CCGCCA-loss disease allele was associated with an earlier age of diagnosis of approximately 7.1 years and occurred exclusively on haplotype B2. Although somatic expansion was associated with an earlier age of diagnosis in the study overall, the CCGCCA-loss disease allele displayed reduced somatic expansion relative to the typical HTT expansions in blood DNA. We propose that the CCGCCA loss occurring on haplotype B2 is an African cis-acting modifier that appears to alter disease diagnosis of HD through a mechanism that is not driven by somatic expansion. The assessment of a group of individuals from an understudied population has highlighted population-specific differences that emphasize the importance of studying genetically diverse populations in the context of disease.

Genetic modifiers of Huntington disease differentially influence motor and cognitive domains

Genome-wide association studies (GWASs) of Huntington disease (HD) have identified six DNA maintenance gene loci (among others) as modifiers and implicated a two step-mechanism of pathogenesis: somatic instability of the causative HTT CAG repeat with subsequent triggering of neuronal damage. The largest studies have been limited to HD individuals with a rater-estimated age at motor onset. To capitalize on the wealth of phenotypic data in several large HD natural history studies, we have performed algorithmic prediction by using common motor and cognitive measures to predict age at other disease landmarks as additional phenotypes for GWASs. Combined with imputation with the Trans-Omics for Precision Medicine reference panel, predictions using integrated measures provided objective landmark phenotypes with greater power to detect most modifier loci. Importantly, substantial differences in the relative modifier signal across loci, highlighted by comparing common modifiers at MSH3 and FAN1, revealed that individual modifier effects can act preferentially in the motor or cognitive domains. Individual components of the DNA maintenance modifier mechanisms may therefore act differentially on the neuronal circuits underlying the corresponding clinical measures. In addition, we identified additional modifier effects at the PMS1 and PMS2 loci and implicated a potential second locus on chromosome 7. These findings indicate that broadened discovery and characterization of HD genetic modifiers based on additional quantitative or qualitative phenotypes offers not only the promise of in-human validated therapeutic targets but also a route to dissecting the mechanisms and cell types involved in both the somatic instability and toxicity components of HD pathogenesis.

Blood meal analysis of tsetse flies ( Glossina pallidipes: Glossinidae) reveals higher host fidelity on wild compared with domestic hosts

Background: Changes in climate and land use can alter risk of transmission of parasites between domestic hosts and wildlife, particularly when mediated by vectors that can travel between populations. Here we focused on tsetse flies (genus Glossina), the cyclical vectors for both Human African Trypanosomiasis (HAT) and Animal African Trypanosomiasis (AAT). The aims of this study were to investigate three issues related to G. palldipes from Kenya: 1) the diversity of vertebrate hosts that flies fed on; 2) whether host feeding patterns varied in relation to type of hosts, tsetse feeding behaviour, site or tsetse age and sex; and 3) if there was a relationship between trypanosome detection and host feeding behaviours or host types. Methods: Sources of blood meals of Glossina pallidipes were identified by sequencing of the mitochondrial cytochrome b gene and analyzed in relationship with previously determined trypanosome detection in the same flies. Results: In an area dominated by wildlife but with seasonal presence of livestock (Nguruman), 98% of tsetse fed on single wild host species, whereas in an area including a mixture of resident domesticated animals, humans and wildlife (Shimba Hills), 52% of flies fed on more than one host species. Multiple Correspondence Analysis revealed strong correlations between feeding pattern, host type and site but these were resolved along a different dimension than trypanosome status, sex and age of the flies. Conclusions: Our results suggest that individual G. pallidipes in interface areas may show higher feeding success on wild hosts when available but often feed on both wild and domesticated hosts. This illustrates the importance of G. pallidipes as a vector connecting the sylvatic and domestic cycles of African trypanosomes.

Association Analysis of Chromosome X to Identify Genetic Modifiers of Huntington's Disease

BACKGROUND

Huntington's disease (HD) is caused by an expanded (>35) CAG trinucleotide repeat in huntingtin (HTT). Age-at-onset of motor symptoms is inversely correlated with the size of the inherited CAG repeat, which expands further in brain regions due to somatic repeat instability. Our recent genetic investigation focusing on autosomal SNPs revealed that age-at-onset is also influenced by genetic variation at many loci, the majority of which encode genes involved in DNA maintenance/repair processes and repeat instability.

OBJECTIVE

We performed a complementary association analysis to determine whether variants in the X chromosome modify HD.

METHODS

We imputed SNPs on chromosome X for ∼9,000 HD subjects of European ancestry and performed an X chromosome-wide association study (XWAS) to test for association with age-at-onset corrected for inherited CAG repeat length.

RESULTS

In a mixed effects model XWAS analysis of all subjects (males and females), assuming random X-inactivation in females, no genome-wide significant onset modification signal was found. However, suggestive significant association signals were detected at Xq12 (top SNP, rs59098970; p-value, 1.4E-6), near moesin (MSN), in a region devoid of DNA maintenance genes. Additional suggestive signals not involving DNA repair genes were observed in male- and female-only analyses at other locations.

CONCLUSION

Although not genome-wide significant, potentially due to small effect size compared to the power of the current study, our data leave open the possibility of modification of HD by a non-DNA repair process. Our XWAS results are publicly available at the updated GEM EURO 9K website hosted at https://www.hdinhd.org/ for browsing, pathway analysis, and data download.

Approaches to Sequence the HTT CAG Repeat Expansion and Quantify Repeat Length Variation

Background:

Huntington’s disease (HD) is an autosomal dominant neurodegenerative disorder caused by the expansion of the HTT CAG repeat. Affected individuals inherit ≥36 repeats and longer alleles cause earlier onset, greater disease severity and faster disease progression. The HTT CAG repeat is genetically unstable in the soma in a process that preferentially generates somatic expansions, the proportion of which is associated with disease onset, severity and progression. Somatic mosaicism of the HTT CAG repeat has traditionally been assessed by semi-quantitative PCR-electrophoresis approaches that have limitations (e.g., no information about sequence variants). Genotyping-by-sequencing could allow for some of these limitations to be overcome.

Objective:

To investigate the utility of PCR sequencing to genotype large (>50 CAGs) HD alleles and to quantify the associated somatic mosaicism.

Methods:

We have applied MiSeq and PacBio sequencing to PCR products of the HTT CAG repeat in transgenic R6/2 mice carrying ∼55, ∼110, ∼255 and ∼470 CAGs. For each of these alleles, we compared the repeat length distributions generated for different tissues at two ages.

Results:

We were able to sequence the CAG repeat full length in all samples. However, the repeat length distributions for samples with ∼470 CAGs were biased towards shorter repeat lengths.

Conclusion:

PCR sequencing can be used to sequence all the HD alleles considered, but this approach cannot be used to estimate modal allele size or quantify somatic expansions for alleles ⪢250 CAGs. We review the limitations of PCR sequencing and alternative approaches that may allow the quantification of somatic contractions and very large somatic expansions.

Genetic Risk Underlying Psychiatric and Cognitive Symptoms in Huntington's Disease

BACKGROUND

Huntington's disease (HD) is an inherited neurodegenerative disorder caused by an expanded CAG repeat in the HTT gene. It is diagnosed following a standardized examination of motor control and often presents with cognitive decline and psychiatric symptoms. Recent studies have detected genetic loci modifying the age at onset of motor symptoms in HD, but genetic factors influencing cognitive and psychiatric presentations are unknown.

METHODS

We tested the hypothesis that psychiatric and cognitive symptoms in HD are influenced by the same common genetic variation as in the general population by 1) constructing polygenic risk scores from large genome-wide association studies of psychiatric and neurodegenerative disorders and of intelligence and 2) testing for correlation with the presence of psychiatric and cognitive symptoms in a large sample (n = 5160) of patients with HD.

RESULTS

Polygenic risk score for major depression was associated specifically with increased risk of depression in HD, as was schizophrenia risk score with psychosis and irritability. Cognitive impairment and apathy were associated with reduced polygenic risk score for intelligence.

CONCLUSIONS

Polygenic risk scores for psychiatric disorders, particularly depression and schizophrenia, are associated with increased risk of the corresponding psychiatric symptoms in HD, suggesting a common genetic liability. However, the genetic liability to cognitive impairment and apathy appears to be distinct from other psychiatric symptoms in HD. No associations were observed between HD symptoms and risk scores for other neurodegenerative disorders. These data provide a rationale for treatments effective in depression and schizophrenia to be used to treat depression and psychotic symptoms in HD.

CAG Repeat Not Polyglutamine Length Determines Timing of Huntington's Disease Onset

Variable, glutamine-encoding, CAA interruptions indicate that a property of the uninterrupted HTT CAG repeat sequence, distinct from the length of huntingtin's polyglutamine segment, dictates the rate at which Huntington's disease (HD) develops. The timing of onset shows no significant association with HTT cis-eQTLs but is influenced, sometimes in a sex-specific manner, by polymorphic variation at multiple DNA maintenance genes, suggesting that the special onset-determining property of the uninterrupted CAG repeat is a propensity for length instability that leads to its somatic expansion. Additional naturally occurring genetic modifier loci, defined by GWAS, may influence HD pathogenesis through other mechanisms. These findings have profound implications for the pathogenesis of HD and other repeat diseases and question the fundamental premise that polyglutamine length determines the rate of pathogenesis in the "polyglutamine disorders."

MSH3 modifies somatic instability and disease severity in Huntington's and myotonic dystrophy type 1

The mismatch repair gene MSH3 has been implicated as a genetic modifier of the CAG·CTG repeat expansion disorders Huntington's disease and myotonic dystrophy type 1. A recent Huntington's disease genome-wide association study found rs557874766, an imputed single nucleotide polymorphism located within a polymorphic 9 bp tandem repeat in MSH3/DHFR, as the variant most significantly associated with progression in Huntington's disease. Using Illumina sequencing in Huntington's disease and myotonic dystrophy type 1 subjects, we show that rs557874766 is an alignment artefact, the minor allele for which corresponds to a three-repeat allele in MSH3 exon 1 that is associated with a reduced rate of somatic CAG·CTG expansion (P = 0.004) and delayed disease onset (P = 0.003) in both Huntington's disease and myotonic dystrophy type 1, and slower progression (P = 3.86 × 10-7) in Huntington's disease. RNA-Seq of whole blood in the Huntington's disease subjects found that repeat variants are associated with MSH3 and DHFR expression. A transcriptome-wide association study in the Huntington's disease cohort found increased MSH3 and DHFR expression are associated with disease progression. These results suggest that variation in the MSH3 exon 1 repeat region influences somatic expansion and disease phenotype in Huntington's disease and myotonic dystrophy type 1, and suggests a common DNA repair mechanism operates in both repeat expansion diseases.

Sodalis glossinidius presence in wild tsetse is only associated with presence of trypanosomes in complex interactions with other tsetse-specific factors

Background

Susceptibility of tsetse flies (Glossina spp.) to trypanosomes of both humans and animals has been associated with the presence of the endosymbiont Sodalis glossinidius. However, intrinsic biological characteristics of the flies and environmental factors can influence the presence of both S. glossinidius and the parasites. It thus remains unclear whether it is the S. glossinidius or other attributes of the flies that explains the apparent association. The objective of this study was to test whether the presence of Trypanosoma vivax, T. congolense and T. brucei are related to the presence of S. glossinidius in tsetse flies when other factors are accounted for: geographic location, species of Glossina, sex or age of the host flies.

Results

Flies (n = 1090) were trapped from four sites in the Shimba Hills and Nguruman regions in Kenya. Sex and species of tsetse (G. austeni, G. brevipalpis, G. longipennis and G. pallidipes) were determined based on external morphological characters and age was estimated by a wing fray score method. The presence of trypanosomes and S. glossinidius was detected using PCR targeting the internal transcribed spacer region 1 and the haemolysin gene, respectively. Sequencing was used to confirm species identification. Generalised Linear Models (GLMs) and Multiple Correspondence Analysis (MCA) were applied to investigate multivariable associations. The overall prevalence of trypanosomes was 42.1%, but GLMs revealed complex patterns of associations: the presence of S. glossinidius was associated with trypanosome presence but only in interactions with other factors and only in some species of trypanosomes. The strongest association was found for T. congolense, and no association was found for T. vivax. The MCA also suggested only a weak association between the presence of trypanosomes and S. glossinidius. Trypanosome-positive status showed strong associations with sex and age while S. glossinidius-positive status showed a strong association with geographic location and species of fly.

Conclusions

We suggest that previous conclusions about the presence of endosymbionts increasing probability of trypanosome presence in tsetse flies may have been confounded by other factors, such as community composition of the tsetse flies and the specific trypanosomes found in different regions.

Electronic supplementary material

The online version of this article (10.1186/s12866-018-1285-6) contains supplementary material, which is available to authorized users.

Morphological re-description and molecular identification of Tabanidae (Diptera) in East Africa

Biting flies of the family Tabanidae are important vectors of human and animal diseases across continents. However, records of Africa tabanids are fragmentary and mostly cursory. To improve identification, documentation and description of Tabanidae in East Africa, a baseline survey for the identification and description of Tabanidae in three eastern African countries was conducted. Tabanids from various locations in Uganda (Wakiso District), Tanzania (Tarangire National Park) and Kenya (Shimba Hills National Reserve, Muhaka, Nguruman) were collected. In Uganda, octenol baited F-traps were used to target tabanids, while NG2G traps baited with cow urine and acetone were employed in Kenya and Tanzania. The tabanids were identified using morphological and molecular methods. Morphologically, five genera (Ancala, Tabanus, Atylotus, Chrysops and Haematopota) and fourteen species of the Tabanidae were identified. Among the 14 species identified, six belonged to the genus Tabanus of which two (T. donaldsoni and T. guineensis) had not been described before in East Africa. The greatest diversity of tabanid species were collected from the Shimba Hills National Reserve, while collections from Uganda (around the shores of Lake Victoria) had the fewest number of species. However, the Ancala genus was found in Uganda, but not in Kenya or Tanzania. Maximum likelihood phylogenies of mitochondrial cytochrome c oxidase 1 (COI) genes sequenced in this study show definite concordance with morphological species identifications, except for Atylotus. This survey will be critical to building a complete checklist of Tabanidae prevalent in the region, expanding knowledge of these important vectors of human and animal diseases.

Molecular identification of different trypanosome species and subspecies in tsetse flies of northern Nigeria

Background

Animal African Trypanosomiasis (AAT) is caused by several species of trypanosomes including Trypanosoma congolense, T. vivax, T. godfreyi, T. simiae and T. brucei. Two of the subspecies of T. brucei also cause Human African Trypanosomiasis. Although some of them can be mechanically transmitted by biting flies; these trypanosomes are all transmitted by tsetse flies which are the cyclical vectors of Trypanosoma congolense, T. godfreyi, T. simiae and T. brucei. We present here the first report assessing the prevalence of trypanosomes in tsetse flies in Nigeria using molecular tools.

Methods

488 tsetse flies of three species, Glossina palpalis palpalis, G. tachinoides and G. morsitans submorsitans were collected from Wuya, Niger State and Yankari National Park, Bauchi State in 2012. Trypanosomes were detected and identified using an ITS1 PCR assay on DNA purified from the ‘head plus proboscis’ (H + P) and abdomen (ABD) parts of each fly.

Results

T. vivax and T. congolense Savannah were the major parasites detected. Trypanosomes prevalence was 7.1 % in G. p. palpalis, 11.9 % in G. tachinoides and 13.5 % in G. m. submorsitans. Prevalences of T. congolense Savannah ranged from 2.5 to 6.7 % and of T. vivax were approximately 4.5 %. Trypanosoma congolense Forest, T. godfreyi and T. simiae were also detected in the site of Yankari. The main biological and ecological determinants of trypanosome prevalence were the fly sex, with more trypanosomes found in females than males, and the site, with T. congolense subspp. being more abundant in Yankari than in Wuya. As expected, the trypanosome species diversity was higher in Yankari National Park than in the more agricultural site of Wuya where vertebrate host species diversity is lower.

Conclusions

Our results show that T. congolense Savannah and T. vivax are the main species of parasite potentially causing AAT in the two study sites and that Yankari National Park is a potential reservoir of trypanosomes both in terms of parasite abundance and species diversity.

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-016-1585-3) contains supplementary material, which is available to authorized users.

Laboratory colonisation and genetic bottlenecks in the tsetse fly Glossina pallidipes

Background

The IAEA colony is the only one available for mass rearing of Glossina pallidipes, a vector of human and animal African trypanosomiasis in eastern Africa. This colony is the source for Sterile Insect Technique (SIT) programs in East Africa. The source population of this colony is unclear and its genetic diversity has not previously been evaluated and compared to field populations.

Methodology/Principal Findings

We examined the genetic variation within and between the IAEA colony and its potential source populations in north Zimbabwe and the Kenya/Uganda border at 9 microsatellites loci to retrace the demographic history of the IAEA colony. We performed classical population genetics analyses and also combined historical and genetic data in a quantitative analysis using Approximate Bayesian Computation (ABC). There is no evidence of introgression from the north Zimbabwean population into the IAEA colony. Moreover, the ABC analyses revealed that the foundation and establishment of the colony was associated with a genetic bottleneck that has resulted in a loss of 35.7% of alleles and 54% of expected heterozygosity compared to its source population. Also, we show that tsetse control carried out in the 1990's is likely reduced the effective population size of the Kenya/Uganda border population.

Conclusions/Significance

All the analyses indicate that the area of origin of the IAEA colony is the Kenya/Uganda border and that a genetic bottleneck was associated with the foundation and establishment of the colony. Genetic diversity associated with traits that are important for SIT may potentially have been lost during this genetic bottleneck which could lead to a suboptimal competitiveness of the colony males in the field. The genetic diversity of the colony is lower than that of field populations and so, studies using colony flies should be interpreted with caution when drawing general conclusions about G. pallidipes biology.

There is only one mass reared laboratory colony of Glossina pallidipes, a vector of human African trypanosomiasis and arguably the main vector of animal African trypanosomiasis in eastern Africa. This colony is the main one used for basic research on this species and is intended to be used for Sterile Insect Technique (SIT) programs for control of field populations. The origins of this colony are not clear and the extent to which it is genetically representative of the species is unknown. Using population genetics analyses to compare the current colony with two potential source populations we have shown that the colony is from the Kenya/Uganda border and that its foundation and establishment in the laboratory were associated with a genetic bottleneck, i.e. reduction of genetic variation due to increased genetic drift in a population of reduced size. As a consequence, the genetic diversity of the colony is lower than that of G. pallidipes field populations.

Secondary contact and admixture between independently invading populations of the western corn rootworm, Diabrotica virgifera virgifera in Europe

The western corn rootworm, Diabrotica virgifera virgifera (Coleoptera: Chrysomelidae), is one of the most destructive pests of corn in North America and is currently invading Europe. The two major invasive outbreaks of rootworm in Europe have occurred, in North-West Italy and in Central and South-Eastern Europe. These two outbreaks originated from independent introductions from North America. Secondary contact probably occurred in North Italy between these two outbreaks, in 2008. We used 13 microsatellite markers to conduct a population genetics study, to demonstrate that this geographic contact resulted in a zone of admixture in the Italian region of Veneto. We show that i) genetic variation is greater in the contact zone than in the parental outbreaks; ii) several signs of admixture were detected in some Venetian samples, in a bayesian analysis of the population structure and in an approximate bayesian computation analysis of historical scenarios and, finally, iii) allelic frequency clines were observed at microsatellite loci. The contact between the invasive outbreaks in North-West Italy and Central and South-Eastern Europe resulted in a zone of admixture, with particular characteristics. The evolutionary implications of the existence of a zone of admixture in Northern Italy and their possible impact on the invasion success of the western corn rootworm are discussed.

Stratified dispersal and increasing genetic variation during the invasion of Central Europe by the western corn rootworm, Diabrotica virgifera virgifera

Invasive species provide opportunities for investigating evolutionary aspects of colonization processes, including initial foundations of populations and geographic expansion. Using microsatellite markers and historical information, we characterized the genetic patterns of the invasion of the western corn rootworm (WCR), a pest of corn crops, in its largest area of expansion in Europe: Central and South-Eastern (CSE) Europe. We found that the invaded area probably corresponds to a single expanding population resulting from a single introduction of WCR and that gene flow is geographically limited within the population. In contrast to what is expected in classical colonization processes, an increase in genetic variation was observed from the center to the edge of the outbreak. Control measures against WCR at the center of the outbreak may have decreased effective population size in this area which could explain this observed pattern of genetic variation. We also found that small remote outbreaks in southern Germany and north-eastern Italy most likely originated from long-distance dispersal events from CSE Europe. We conclude that the large European outbreak is expanding by stratified dispersal, involving both continuous diffusion and discontinuous long-distance dispersal. This latter mode of dispersal may accelerate the expansion of WCR in Europe in the future.