Contemporary population genetic structure and phylogeography of pink salmon

Thesis (Master's)--University of Washington, 2016-06

Genetic manipulation of Anopheles stephensi immunity to increase plasmodium falciparum salivary gland sporozoite infection levels

Human malaria is responsible for over 700,000 deaths a year. To stay abreast of the threat posed by the parasite, a constant stream of new drugs and vector control methods are required. This study focuses on a vaccine that has the potential to protect against parasite infection, but has been hindered by developmental challenges. In malaria prevention, live, attenuated, aseptic, Plasmodium falciparum sporozoites (PfSPZ) can be administered as a highly protective vaccine. PfSPZ are produced using adult female Anopheles stephensi mosquitoes as bioreactors. Production volume and cost of a PfSPZ vaccine for malaria are expected to be directly correlated with Plasmodium falciparum infection intensity in the salivary glands. The sporogonic development of Plasmodium falciparum in A. stephensi to fully infected salivary gland stage sporozoites is dictated by the activities of several known components of the mosquito’s innate immune system. Here I report on the use of genetic technologies that have been rarely, if ever, used in Anopheles stephensi Sda500 to increase the yield of sporozoites per mosquito and enhance vaccine production. By combining the Gal4/UAS bipartite system with in vivo expression of shRNA gene silencing, activity of the IMD signaling pathway downstream effector LRIM1, an antagonist to Plasmodium development, was reduced in the midgut, fat body, and salivary glands of A. stephensi. In infection studies using P. berghei and P. falciparum these transgenic mosquitoes consistently produced significantly more salivary gland stage sporozoites than wildtype controls, with increases in P. falciparum ranging from 2.5 to 10 fold. Using Plasmodium infection assays and qRT-PCR, two novel findings were identified. First, it was shown that 14 days post Plasmodium infection, transcript abundance of the IMD immune effector genes LRIM1, TEP1 and APL1c are elevated, in the salivary glands of A. stephensi, suggesting the salivary glands may play a role in post midgut defense against the parasite. Second, a non-pathogenic IMD signaling pathway response was observed which could suggest an alternative pathway for IMD activation. The information gained from these studies has significantly increased our knowledge of Plasmodium defense in A. stephensi and moreover could significantly improve vaccine production.

Genetic viability of Nabugabo lakes (LVR satellite lakes) fish species

Natural populations of fish species in Lake Victoria Region (LVR) have under gone dramatic changes including severe reduction in sizes, division of original stocks into disjunct subunits, and segregation into several isolated population units either within a single water body or even worse into separate waters. In addition, these changes have been either preceded or precipitated by introductions of non-indigenous species that out competed the native forms and in case of closely related species genetically swamped them through hybridisation. The latter is especially the case in Nabugabo lakes. Such events lead to fragmentation of populations, which results in reduction in genetic diversity due to genetic drift, inbreeding and reduced or lack of gene flow among independent units. Such phenomena make the continued existence of fisheries stocks in the wild precarious, more so in the face of the competition from exotic species. Species introductions coupled with growing exploitation pressure of the fisheries of these lakes have put the native stocks at risk. Nabugabo lakes harbor cichlid species that are unique to these lakes more so species of the cichlid complex. In this paper the ecological status and genetic viability of key Nabugabo lakes fish species is examined and management options are discussed.

A New Genetic Algorithm for Set Covering Problems

An indirect genetic algorithm for the non-unicost set covering problem is presented. The algorithm is a two-stage meta-heuristic, which in the past was successfully applied to similar multiple-choice optimisation problems. The two stages of the algorithm are an ‘indirect’ genetic algorithm and a decoder routine. First, the solutions to the problem are encoded as permutations of the rows to be covered, which are subsequently ordered by the genetic algorithm. Fitness assignment is handled by the decoder, which transforms the permutations into actual solutions to the set covering problem. This is done by exploiting both problem structure and problem specific information. However, flexibility is retained by a self-adjusting element within the decoder, which allows adjustments to both the data and to stages within the search process. Computational results are presented.

Information and the evolution of codon bias

The informational properties of biological systems are the subject of much debate and research. I present a general argument in favor of the existence and central importance of information in organisms, followed by a case study of the genetic code (specifically, codon bias) and the translation system from the perspective of information. The codon biases of 831 Bacteria and Archeae are analyzed and modeled as points in a 64-dimensional statistical space. The major results are that (1) codon bias evolution does not follow canonical patterns, and (2) the use of coding space in organsims is a subset of the total possible coding space. These findings imply that codon bias is a unique adaptive mechanism that owes its existence to organisms' use of information in representing genes, and that there is a particularly biological character to the resulting biased coding and information use.

Genetic variation in growth and wood-quality traits of Corymbia citriodora subsp. variegata across three sites in south-east Queensland, Australia

Ten growth or wood-quality traits were assessed in three nearby Corymbia citriodora subsp. variegata (CCV) open-pollinated family-within-provenance trials (18 provenances represented by a total of 374 families) to provide information for the development of a breeding program targeting both pulp and solid-wood products. Growth traits (diameter at breast high over bark [DBH], height and conical volume) were assessed at 3 and 7 years of age. Wood-quality traits (density [DEN], Kraft pulp yield [KPY], modulus of elasticity [MoE] and microfibril angle [MfA]) were predicted using near-infrared spectroscopy on wood samples collected from these trials when aged between 10 and 12 years. The high average KPY, DEN and MoE, and low average MfA observed indicates CCV is very suitable for both pulp and timber products. All traits were under moderate to strong genetic control. In across- trials analyses, high (>0.4) heritability estimates were observed for height, DEN, MoE and MfA, while moderate heritability estimates (0.24 to 0.34) were observed for DBH, volume and KPY. Most traits showed very low levels of genotype × site interaction. Estimated age–age genetic correlations for growth traits were strong at both the family (0.97) and provenance (0.99) levels. Relationships among traits (additive genetic correlation estimates) were favourable, with strong and positive estimates between growth traits (0.84 to 0.98), moderate and positive values between growth and wood-quality traits (0.32 to 0.68), moderate and positive between KPY and MoE (0.64), and high and positive between DEN and MoE (0.82). However, negative (but favourable) correlations were detected between MfA and all other evaluated traits (−0.31 to −0.96). The genetic correlation between the same trait expressed on two different sites, at family level, ranged from 0.24 to 0.42 for growth traits, and from 0.29 to 0.53 for wood traits. Therefore simultaneous genetic improvement of growth and wood property traits in CCV for the target environment in south-east Queensland should be possible, given the moderate to high estimates of heritability and favourable correlations amongst all traits studied, unless genotype × site interactions are greater than was evident. © 2016 NISC (Pty) Ltd

Epidemiology and genetic architecture of blood pressure: a family based study of Generation Scotland

Hypertension is a major risk factor for cardiovascular disease and mortality, and a growing global public health concern, with up to one-third of the world’s population affected. Despite the vast amount of evidence for the benefits of blood pressure (BP) lowering accumulated to date, elevated BP is still the leading risk factor for disease and disability worldwide. It is well established that hypertension and BP are common complex traits, where multiple genetic and environmental factors contribute to BP variation. Furthermore, family and twin studies confirmed the genetic component of BP, with a heritability estimate in the range of 30-50%. Contemporary genomic tools enabling the genotyping of millions of genetic variants across the human genome in an efficient, reliable, and cost-effective manner, has transformed hypertension genetics research. This is accompanied by the presence of international consortia that have offered unprecedentedly large sample sizes for genome-wide association studies (GWASs). While GWAS for hypertension and BP have identified more than 60 loci, variants in these loci are associated with modest effects on BP and in aggregate can explain less than 3% of the variance in BP. The aims of this thesis are to study the genetic and environmental factors that influence BP and hypertension traits in the Scottish population, by performing several genetic epidemiological analyses. In the first part of this thesis, it aims to study the burden of hypertension in the Scottish population, along with assessing the familial aggregation and heritialbity of BP and hypertension traits. In the second part, it aims to validate the association of common SNPs reported in the large GWAS and to estimate the variance explained by these variants. In this thesis, comprehensive genetic epidemiology analyses were performed on Generation Scotland: Scottish Family Health Study (GS:SFHS), one of the largest population-based family design studies. The availability of clinical, biological samples, self-reported information, and medical records for study participants has allowed several assessments to be performed to evaluate factors that influence BP variation in the Scottish population. Of the 20,753 subjects genotyped in the study, a total of 18,470 individuals (grouped into 7,025 extended families) passed the stringent quality control (QC) criteria and were available for all subsequent analysis. Based on the BP-lowering treatment exposure sources, subjects were further classified into two groups. First, subjects with both a self-reported medications (SRMs) history and electronic-prescription records (EPRs; n =12,347); second, all the subjects with at least one medication history source (n =18,470). In the first group, the analysis showed a good concordance between SRMs and EPRs (kappa =71%), indicating that SRMs can be used as a surrogate to assess the exposure to BP-lowering medication in GS:SFHS participants. Although both sources suffer from some limitations, SRMs can be considered the best available source to estimate the drug exposure history in those without EPRs. The prevalence of hypertension was 40.8% with higher prevalence in men (46.3%) compared to women (35.8%). The prevalence of awareness, treatment and controlled hypertension as defined by the study definition were 25.3%, 31.2%, and 54.3%, respectively. These findings are lower than similar reported studies in other populations, with the exception of controlled hypertension prevalence, which can be considered better than other populations. Odds of hypertension were higher in men, obese or overweight individuals, people with a parental history of hypertension, and those living in the most deprived area of Scotland. On the other hand, deprivation was associated with higher odds of treatment, awareness and controlled hypertension, suggesting that people living in the most deprived area may have been receiving better quality of care, or have higher comorbidity levels requiring greater engagement with doctors. These findings highlight the need for further work to improve hypertension management in Scotland. The family design of GS:SFHS has allowed family-based analysis to be performed to assess the familial aggregation and heritability of BP and hypertension traits. The familial correlation of BP traits ranged from 0.07 to 0.20, and from 0.18 to 0.34 for parent-offspring pairs and sibling pairs, respectively. A higher correlation of BP traits was observed among first-degree relatives than other types of relative pairs. A variance-component model that was adjusted for sex, body mass index (BMI), age, and age-squared was used to estimate heritability of BP traits, which ranged from 24% to 32% with pulse pressure (PP) having the lowest estimates. The genetic correlation between BP traits showed a high correlation between systolic (SBP), diastolic (DBP) and mean arterial pressure (MAP) (G: 81% to 94%), but lower correlations with PP (G: 22% to 78%). The sibling recurrence risk ratio (λS) for hypertension and treatment were calculated as 1.60 and 2.04 respectively. These findings confirm the genetic components of BP traits in GS:SFHS, and justify further work to investigate genetic determinants of BP. Genetic variants reported in the recent large GWAS of BP traits were selected for genotyping in GS:SFHS using a custom designed TaqMan® OpenArray®. The genotyping plate included 44 single nucleotide polymorphisms (SNPs) that have been previously reported to be associated with BP or hypertension at genome-wide significance level. A linear mixed model that is adjusted for age, age-squared, sex, and BMI was used to test for the association between the genetic variants and BP traits. Of the 43 variants that passed the QC, 11 variants showed statistically significant association with at least one BP trait. The phenotypic variance explained by these variant for the four BP traits were 1.4%, 1.5%, 1.6%, and 0.8% for SBP, DBP, MAP, and PP, respectively. The association of genetic risk score (GRS) that were constructed from selected variants has showed a positive association with BP level and hypertension prevalence, with an average effect of one mmHg increase with each 0.80 unit increases in the GRS across the different BP traits. The impact of BP-lowering medication on the genetic association study for BP traits has been established, with typical practice of adding a fixed value (i.e. 15/10 mmHg) to the measured BP values to adjust for BP treatment. Using the subset of participants with the two treatment exposure sources (i.e. SRMs and EPRs), the influence of using either source to justify the addition of fixed values in SNP association signal was analysed. BP phenotypes derived from EPRs were considered the true phenotypes, and those derived from SRMs were considered less accurate, with some phenotypic noise. Comparing SNPs association signals between the four BP traits in the two model derived from the different adjustments showed that MAP was the least impacted by the phenotypic noise. This was suggested by identifying the same overlapped significant SNPs for the two models in the case of MAP, while other BP traits had some discrepancy between the two sources

Development, characterization and use of genomic SSR markers for assessment of genetic diversity in some Saudi date palm ( Phoenix dactylifera L.) cultivars

Background: The present study was undertaken towards the development of SSR markers and assessing genetic relationships among 32 date palm ( Phoenix dactylifera L.) representing common cultivars grown in different geographical regions in Saudi Arabia. Results: Ninety-three novel simple sequence repeat markers were developed and screened for their ability to detect polymorphism in date palm. Around 71% of genomic SSRs were dinucleotide, 25% tri, 3% tetra and 1% penta nucleotide motives. Twenty-two primers generated a total of 91 alleles with a mean of 4.14 alleles per locus and 100% polymorphism percentage. A 0.595 average polymorphic information content and 0.662 primer discrimination power values were recorded. The expected and observed heterozygosities were 0.676 and 0.763 respectively. Pair-wise similarity values ranged from 0.06 to 0.89 and the overall cultivars averaged 0.41. The UPGMA cluster analysis recovered by principal coordinate analysis illustrated that cultivars tend to group according to their class of maturity, region of cultivation, and fruit color. Analysis of molecular variations (AMOVA) revealed that genetic variation among and within cultivars were 27% and 73%, respectively according to geographical distribution of cultivars. Conclusions: The developed microsatellite markers are additional values to date palm characterization tools that can be used by researchers in population genetics, cultivar identification as well as genetic resource exploration and management. The tested cultivars exhibited a significant amount of genetic diversity and could be suitable for successful breeding program. Genomic sequences generated from this study are available at the National Center for Biotechnology Information (NCBI), Sequence Read Archive (Accession numbers. LIBGSS_039019).

SSR genetic diversity assessment of popular pigeonpea varieties in Malawi reveals unique fingerprints

Background: Pigeonpea ( Cajanus cajan L. Millsp.) is a drought tolerant legume of the Fabaceae family and the only cultivated species in the genus Cajanus. It is mainly cultivated in the semi-arid tropics of Asia and Oceania, Africa and America. In Malawi, it is grown as a source of food and income and for soil improvement in intercropping systems. However, varietal contamination due to natural outcrossing causes significant quality reduction and yield losses. In this study, 48 polymorphic SSR markers were used to assess the diversity among all pigeonpea varieties cultivated in Malawi to determine if a genetic fingerprint could be identified to distinguish the popular varieties. Results: A total of 212 alleles were observed with an average of 5.58 alleles per marker and a maximum of 14 alleles produced by CCttc019 (Marker 40). Polymorphic information content (PIC), ranged from 0.03 to 0.89 with an average of 0.30. A neighbor-joining tree produced 4 clusters. The most commonly cultivated varieties, which include released varieties and cultivated land races, were well-spread across all the clusters observed, indicating that they generally represented the genetic diversity available in Malawi, although substantial variation was evident that can still be exploited through further breeding. Conclusion: Screening of the allelic data associated with the five most popular cultivated varieties, revealed 6 markers – CCB1, CCB7, Ccac035, CCttc003, Ccac026 and CCttc019 – which displayed unique allelic profiles for each of the five varieties. This genetic fingerprint can potentially be applied for seed certification to confirm the genetic purity of seeds that are delivered to Malawi farmers.

Genetic Diversity of Oilseed Rape Fields and Feral Populations in the Context of Coexistence with GM Crops

Despite growing concern about transgenes escaping from fields, few studies have analysed the genetic diversity of crops in an agroecosystem over several years. Accurate information about the dynamics and relationship of the genetic diversity of crops in an agroecosystem is essential for risk assessment and policies concerning the containment of genetically modified crops and their coexistence with crops grown by conventional practices. Here, we analysed the genetic diversity of oilseed rape plants from fields and feral populations over 4 years in an agricultural landscape of 41 km2. We used exact compatibility and maximum likelihood assignment methods to assign these plants to cultivars. Even pure lines and hybrid cultivar seed lots contained several genotypes. The cultivar diversity in fields reflected the conventional view of agroecosystems quite well: that is, there was a succession of cultivars, some grown for longer than others because of their good performance, some used for one year and then abandoned, and others gradually adopted. Three types of field emerged: fields sown with a single cultivar, fields sown with two cultivars, and unassigned fields (too many cultivars or unassigned plants to reliably assign the field). Field plant diversity was higher than expected, indicating the persistence of cultivars that were grown for only one year. The cultivar composition of feral populations was similar to that of field plants, with an increasing number of cultivars each year. By using genetic tools, we found a link between the cultivars of field plants in a particular year and the cultivars of feral population plants in the following year. Feral populations on road verges were more diverse than those on path verges. All of these findings are discussed in terms of their consequences in the context of coexistence with genetically modified crops.

Drinking water disinfection by-products, genetic polymorphisms, and birth outcomes in a european mother-child cohort study

International audience

PREVALENCE, ISOLATION, AND GENETIC CHARACTERIZATION OF TOXOPLASMA GONDII IN CHICKEN

Toxoplasma gondii (T. gondii) is one of the most successful parasites in the world because of its capability of infecting all warm-blooded animals. It has been reported that up to one third of the world population is infected with this parasite. Chickens are recognized as good indicators of the environmental T. gondii oocysts contamination because they obtain food from the ground. Thus, the prevalence of T. gondii in chicken provides more insight related to public health concern from T. gondii. Previous studies have shown a high isolation rate from free-range chickens raised in the United States. The objectives of this study were to evaluate the microbial safety and infection of T. gondii in free-range chickens available at the grocery stores and farms for the consumers to purchase and genotype T. gondii isolates. Chicken hearts were obtained from the local markets and also from the farms raising free- range chickens. Heart juice was obtained from cavities of each heart. Modified agglutination test (MAT) for detection of IgG antibodies was conducted with those heart juice samples with titer of 1:5, 1:25, and 1: 100. Each seropositive heart was pepsin digested and bioassayed into a group of two mice. Six weeks post inoculation (p.i.) mice were bled and euthanized to examine the infection of T. gondii. In addition, multiplex multilocus nested PCR-RFLP was performed to genetically characterize T. gondii isolates with eleven PCR-RFLP markers including SAG1, SAG2, altSAT2, SAG3, BTUB, GRA6, c22-8, c29-a, L358, PK1, and Apico. One hundred fifty from a total of 997 samples (15.0%) were found seropositive for T. gondii. No viable T. gondii was isolated from chicken hearts that were sampled. A total of four genotypes were identified, including one new genotype and three previously identified genotypes. The results suggest that T. gondii oocysts could present in the environment and infect the food animals. T. gondii prevalence in chicken hearts could reflect the environmental contamination of T. gondii and prevalence information can be used to manage T. gondii infection risk.

Spatial patterns of genetic diversity in Mediterranean eagle owl Bubo bubo populations

Little information is available on the patterns of genetic connectivity in owls. We studied the genetic structure of the eagle owl Bubo bubo (Linnaeus, 1758) in southeastern Spain at two different spatial scales. Seven microsatellites previously described for this species were used, although only six loci amplified correctly. The observed low genetic variation could be explained by the short dispersal distance, high mortality rate and high degree of monogamy shown by this large nocturnal predator. As expected, the highest genetic isolation was detected in the geographically most isolated population. Significant genetic differentiation was found among study units separated by less than 50 km. The territorial analysis showed interesting connectivity patterns related with the gene flow and turnover rate of the breeding individuals. The lowest genetic diversity was found in the region with the largest population, which could imply incipient inbreeding.

Prevalence and genetic diversity of Blastocystis in family units living in the United States

The human gut is host to a diversity of microorganisms including the single-celled microbial eukaryote Blastocystis. Although Blastocystis has a global distribution, there is dearth of information relating to its prevalence and diversity in many human populations. The mode of Blastocystis transmission to humans is also insufficiently characterised, however, it is speculated to vary between different populations. Here we investigated the incidence and genetic diversity of Blastocystis in a US population and also the possibility of Blastocystis human-human transmission between healthy individuals using family units (N = 50) living in Boulder, Colorado as our sample-set. Ten of the 139 (~ 7%) individuals in our dataset were positive for Blastocystis, nine of whom were adults and one individual belonging to the children/adolescents group. All positive cases were present in different family units. A number of different Blastocystis subtypes (species) were detected with no evidence of mixed infections. The prevalence of Blastocystis in this subset of the US population is comparatively low relative to other industrialised populations investigated to date; however, subtype diversity was largely consistent with that previously reported in studies of European populations. The distribution of Blastocystis within family units indicates that human-human transmission is unlikely to have occurred within families that participated in this study. It is not unexpected that given the world-wide variation in human living conditions and lifestyles between different populations, both the prevalence of Blastocystis and its mode of transmission to humans may vary considerably.

Favourable areas for expansion and reintroduction of Iberian lynx accounting for distribution trends and genetic diversity of the European rabbit

Although on a local scale Iberian lynx distribution is determined by the availability of prey rabbits, recent modelling analyses have uncovered broad-scale disagreements between these two species’ distribution trends. These analyses showed also that the lynx had become restricted to only a fraction of the rabbit’s genetic diversity, and that this could be jeopardising its survival in the face of environmental hazards and uncertainty. In the present paper, a follow-up was carried out through the building of lynx and rabbit distribution models based on the most recent Spanish mammal atlas. Environmental favourability values for lynx and rabbit were positively correlated within the lynx’s current distribution area, but they were negatively correlated within the total Spanish area where lynx used to occur in the 1980’s. Environmental favourability for rabbits was significantly higher where lynx maintains reproductive populations than where it recently disappeared, indicating that rabbit favourability plays an important role and can be a good predictor of lynx persistence. The lynx and rabbit models were extrapolated to predict favourable areas for both species in Spain as well as in Portugal, on the original scale of the distribution data (10x10 km) and on a 100 times finer spatial resolution (1x1 km). The lynx and rabbit models were also combined through fuzzy logic to forecast the potential for lynx occurrence incorporating information on favourable areas for its main prey. Several areas are proposed as favourable for lynx expansion or re-introduction,