Delimiting species without nuclear monophyly in Madagascar's mouse lemurs.

BACKGROUND: Speciation begins when populations become genetically separated through a substantial reduction in gene flow, and it is at this point that a genetically cohesive set of populations attain the sole property of species: the independent evolution of a population-level lineage. The comprehensive delimitation of species within biodiversity hotspots, regardless of their level of divergence, is important for understanding the factors that drive the diversification of biota and for identifying them as targets for conservation. However, delimiting recently diverged species is challenging due to insufficient time for the differential evolution of characters--including morphological differences, reproductive isolation, and gene tree monophyly--that are typically used as evidence for separately evolving lineages. METHODOLOGY: In this study, we assembled multiple lines of evidence from the analysis of mtDNA and nDNA sequence data for the delimitation of a high diversity of cryptically diverged population-level mouse lemur lineages across the island of Madagascar. Our study uses a multi-faceted approach that applies phylogenetic, population genetic, and genealogical analysis for recognizing lineage diversity and presents the most thoroughly sampled species delimitation of mouse lemur ever performed. CONCLUSIONS: The resolution of a large number of geographically defined clades in the mtDNA gene tree provides strong initial evidence for recognizing a high diversity of population-level lineages in mouse lemurs. We find additional support for lineage recognition in the striking concordance between mtDNA clades and patterns of nuclear population structure. Lineages identified using these two sources of evidence also exhibit patterns of population divergence according to genealogical exclusivity estimates. Mouse lemur lineage diversity is reflected in both a geographically fine-scaled pattern of population divergence within established and geographically widespread taxa, as well as newly resolved patterns of micro-endemism revealed through expanded field sampling into previously poorly and well-sampled regions.

Quantifying Eukaryotic Gene Regulation in Hormone Response and Disease.

Quantifying the function of mammalian enhancers at the genome or population scale has been longstanding challenge in the field of gene regulation. Studies of individual enhancers have provided anecdotal evidence on which many foundational assumptions in the field are based. Genome-scale studies have revealed that the number of sites bound by a given transcription factor far outnumber the genes that the factor regulates. In this dissertation we describe a new method, chromatin immune-enriched reporter assays (ChIP-reporters), and use that approach to comprehensively test the enhancer activity of genomic loci bound by the glucocorticoid receptor (GR). Integrative genomics analyses of our ChIP-reporter data revealed an unexpected mechanism of glucocorticoid (GC)-induced gene regulation. In that mechanism, only the minority of GR bound sites acts as GC-inducible enhancers. Many non-GC-inducible GR binding sites interact with GC-induced sites via chromatin looping. These interactions can increase the activity of GC-induced enhancers. Finally, we describe a method that enables the detection and characterization of the functional effects of non-coding genetic variation on enhancer activity at the population scale. Taken together, these studies yield both mechanistic and genetic evidence that provides context that informs the understanding of the effects of multiple enhancer variants on gene expression.

Sexually coercive male chimpanzees sire more offspring.

In sexually reproducing animals, male and female reproductive strategies often conflict. In some species, males use aggression to overcome female choice, but debate persists over the extent to which this strategy is successful. Previous studies of male aggression toward females among wild chimpanzees have yielded contradictory results about the relationship between aggression and mating behavior. Critically, however, copulation frequency in primates is not always predictive of reproductive success. We analyzed a 17-year sample of behavioral and genetic data from the Kasekela chimpanzee (Pan troglodytes schweinfurthii) community in Gombe National Park, Tanzania, to test the hypothesis that male aggression toward females increases male reproductive success. We examined the effect of male aggression toward females during ovarian cycling, including periods when the females were sexually receptive (swollen) and periods when they were not. We found that, after controlling for confounding factors, male aggression during a female's swollen periods was positively correlated with copulation frequency. However, aggression toward swollen females was not predictive of paternity. Instead, aggression by high-ranking males toward females during their nonswollen periods was positively associated with likelihood of paternity. This indicates that long-term patterns of intimidation allow high-ranking males to increase their reproductive success, supporting the sexual coercion hypothesis. To our knowledge, this is the first study to present genetic evidence of sexual coercion as an adaptive strategy in a social mammal.

Rare hereditary COL4A3/COL4A4 variants may be mistaken for familial focal segmental glomerulosclerosis.

Focal segmental glomerulosclerosis (FSGS) is a histological lesion with many causes, including inherited genetic defects, with significant proteinuria being the predominant clinical finding at presentation. Mutations in COL4A3 and COL4A4 are known to cause Alport syndrome (AS), thin basement membrane nephropathy, and to result in pathognomonic glomerular basement membrane (GBM) findings. Secondary FSGS is known to develop in classic AS at later stages of the disease. Here, we present seven families with rare or novel variants in COL4A3 or COL4A4 (six with single and one with two heterozygous variants) from a cohort of 70 families with a diagnosis of hereditary FSGS. The predominant clinical finding at diagnosis was proteinuria associated with hematuria. In all seven families, there were individuals with nephrotic-range proteinuria with histologic features of FSGS by light microscopy. In one family, electron microscopy showed thin GBM, but four other families had variable findings inconsistent with classical Alport nephritis. There was no recurrence of disease after kidney transplantation. Families with COL4A3 and COL4A4 variants that segregated with disease represent 10% of our cohort. Thus, COL4A3 and COL4A4 variants should be considered in the interpretation of next-generation sequencing data from such patients. Furthermore, this study illustrates the power of molecular genetic diagnostics in the clarification of renal phenotypes.

Olfactory Receptor Subgenomes Linked with Broad Ecological Adaptations in Sauropsida.

Olfactory receptors (ORs) govern a prime sensory function. Extant birds have distinct olfactory abilities, but the molecular mechanisms underlining diversification and specialization remain mostly unknown. We explored OR diversity in 48 phylogenetic and ecologically diverse birds and 2 reptiles (alligator and green sea turtle). OR subgenomes showed species- and lineage-specific variation related with ecological requirements. Overall 1,953 OR genes were identified in reptiles and 16,503 in birds. The two reptiles had larger OR gene repertoires (989 and 964 genes, respectively) than birds (182-688 genes). Overall, birds had more pseudogenes (7,855) than intact genes (1,944). The alligator had significantly more functional genes than sea turtle, likely because of distinct foraging habits. We found rapid species-specific expansion and positive selection in OR14 (detects hydrophobic compounds) in birds and in OR51 and OR52 (detect hydrophilic compounds) in sea turtle, suggestive of terrestrial and aquatic adaptations, respectively. Ecological partitioning among birds of prey, water birds, land birds, and vocal learners showed that diverse ecological factors determined olfactory ability and influenced corresponding olfactory-receptor subgenome. OR5/8/9 was expanded in predatory birds and alligator, suggesting adaptive specialization for carnivory. OR families 2/13, 51, and 52 were correlated with aquatic adaptations (water birds), OR families 6 and 10 were more pronounced in vocal-learning birds, whereas most specialized land birds had an expanded OR family 14. Olfactory bulb ratio (OBR) and OR gene repertoire were correlated. Birds that forage for prey (carnivores/piscivores) had relatively complex OBR and OR gene repertoires compared with modern birds, including passerines, perhaps due to highly developed cognitive capacities facilitating foraging innovations.

Three crocodilian genomes reveal ancestral patterns of evolution among archosaurs.

To provide context for the diversification of archosaurs--the group that includes crocodilians, dinosaurs, and birds--we generated draft genomes of three crocodilians: Alligator mississippiensis (the American alligator), Crocodylus porosus (the saltwater crocodile), and Gavialis gangeticus (the Indian gharial). We observed an exceptionally slow rate of genome evolution within crocodilians at all levels, including nucleotide substitutions, indels, transposable element content and movement, gene family evolution, and chromosomal synteny. When placed within the context of related taxa including birds and turtles, this suggests that the common ancestor of all of these taxa also exhibited slow genome evolution and that the comparatively rapid evolution is derived in birds. The data also provided the opportunity to analyze heterozygosity in crocodilians, which indicates a likely reduction in population size for all three taxa through the Pleistocene. Finally, these data combined with newly published bird genomes allowed us to reconstruct the partial genome of the common ancestor of archosaurs, thereby providing a tool to investigate the genetic starting material of crocodilians, birds, and dinosaurs.

Pathogen-Specific Adaptations to Conserved Signaling Pathways in Cryptococcus neoformans

Cryptococcus neoformans is an opportunistic fungal pathogen that causes significant disease worldwide. Even though this fungus has not evolved specifically to cause human disease, it has a remarkable ability to adapt to many different environments within its infected host. C. neoformans adapts by utilizing conserved eukaryotic and fungal-specific signaling pathways to sense and respond to stresses within the host. Upon infection, two of the most significant environmental changes this organism experiences are elevated temperature and high pH.

Conserved Rho and Ras family GTPases are central regulators of thermotolerance in C. neoformans. Many GTPases require prenylation to associate with cellular membranes and function properly. Using molecular genetic techniques, microscopy, and infection models, I demonstrated that the prenyltransferase, geranylgeranyl transferase I (GGTase I) is required for thermotolerance and pathogenesis. Using fluorescence microscopy, I found that only a subset of conserved GGTase I substrates requires this enzyme for membrane localization. Therefore, the C. neoformans GGTase I may recognize its substrate in a slightly different manner than other eukaryotic organisms.

The alkaline response transcription factor, Rim101, is a central regulator of stress-response genes important for adapting to the host environment. In particular, Rim101 regulates cell surface alterations involved in immune avoidance. In other fungi, Rim101 is activated by alkaline pH through a conserved signaling pathway, but this pathway had yet been characterized in C. neoformans. Using molecular genetic techniques, I identified and analyzed the conserved members of the Rim pathway. I found that it was only partially conserved in C. neoformans, missing the components that sense pH and initiate pathway activation. Using a genetic screen, I identified a novel Rim pathway component named Rra1. Structural prediction and genetic epistasis experiments suggest that Rra1 may serve as the Rim pathway pH sensor in C. neoformans and other related basidiomycete fungi.

To explore the relevance of Rim pathway signaling in the interaction of C neoformans with its host, I characterized the Rim101-regulated cell wall changes that prevent immune detection. Using HPLC, enzymatic degradation, and cell wall stains, I found that the rim101Δ mutation resulted in increased cell wall chitin exposure. In vitro co-culture assays demonstrated that increased chitin exposure is associated with enhanced activation of macrophages and dendritic cells. To further test this association, I demonstrated that other mutant strains with increased chitin exposure induce macrophage and dendritic cell responses similar to rim101Δ. We used primary macrophages from mutant mouse lines to demonstrate that members of both the Toll-like receptor and C-type lectin receptor families are involved in detecting strains with increased chitin exposure. Finally, in vivo immunological experiments demonstrated that the rim101Δ strain induced a global inflammatory immune response in infected mouse lungs, expanding upon our previous in vivo rim101Δ studies. These results demonstrate that cell wall organization largely determines how fungal cells are detected by the immune system.

The Genetics of Success: How Single-Nucleotide Polymorphisms Associated With Educational Attainment Relate to Life-Course Development.

A previous genome-wide association study (GWAS) of more than 100,000 individuals identified molecular-genetic predictors of educational attainment. We undertook in-depth life-course investigation of the polygenic score derived from this GWAS using the four-decade Dunedin Study (N = 918). There were five main findings. First, polygenic scores predicted adult economic outcomes even after accounting for educational attainments. Second, genes and environments were correlated: Children with higher polygenic scores were born into better-off homes. Third, children's polygenic scores predicted their adult outcomes even when analyses accounted for their social-class origins; social-mobility analysis showed that children with higher polygenic scores were more upwardly mobile than children with lower scores. Fourth, polygenic scores predicted behavior across the life course, from early acquisition of speech and reading skills through geographic mobility and mate choice and on to financial planning for retirement. Fifth, polygenic-score associations were mediated by psychological characteristics, including intelligence, self-control, and interpersonal skill. Effect sizes were small. Factors connecting DNA sequence with life outcomes may provide targets for interventions to promote population-wide positive development.

Evidence for GC-biased gene conversion as a driver of between-lineage differences in avian base composition.

BACKGROUND: While effective population size (Ne) and life history traits such as generation time are known to impact substitution rates, their potential effects on base composition evolution are less well understood. GC content increases with decreasing body mass in mammals, consistent with recombination-associated GC biased gene conversion (gBGC) more strongly impacting these lineages. However, shifts in chromosomal architecture and recombination landscapes between species may complicate the interpretation of these results. In birds, interchromosomal rearrangements are rare and the recombination landscape is conserved, suggesting that this group is well suited to assess the impact of life history on base composition. RESULTS: Employing data from 45 newly and 3 previously sequenced avian genomes covering a broad range of taxa, we found that lineages with large populations and short generations exhibit higher GC content. The effect extends to both coding and non-coding sites, indicating that it is not due to selection on codon usage. Consistent with recombination driving base composition, GC content and heterogeneity were positively correlated with the rate of recombination. Moreover, we observed ongoing increases in GC in the majority of lineages. CONCLUSIONS: Our results provide evidence that gBGC may drive patterns of nucleotide composition in avian genomes and are consistent with more effective gBGC in large populations and a greater number of meioses per unit time; that is, a shorter generation time. Thus, in accord with theoretical predictions, base composition evolution is substantially modulated by species life history.

Can Genetics Predict Response to Complex Behavioral Interventions? Evidence from a Genetic Analysis of the Fast Track Randomized Control Trial.

Early interventions are a preferred method for addressing behavioral problems in high-risk children, but often have only modest effects. Identifying sources of variation in intervention effects can suggest means to improve efficiency. One potential source of such variation is the genome. We conducted a genetic analysis of the Fast Track randomized control trial, a 10-year-long intervention to prevent high-risk kindergarteners from developing adult externalizing problems including substance abuse and antisocial behavior. We tested whether variants of the glucocorticoid receptor gene NR3C1 were associated with differences in response to the Fast Track intervention. We found that in European-American children, a variant of NR3C1 identified by the single-nucleotide polymorphism rs10482672 was associated with increased risk for externalizing psychopathology in control group children and decreased risk for externalizing psychopathology in intervention group children. Variation in NR3C1 measured in this study was not associated with differential intervention response in African-American children. We discuss implications for efforts to prevent externalizing problems in high-risk children and for public policy in the genomic era.

Functional Evaluation of Causal Mutations Identified in Human Genetic Studies

Human genetics has been experiencing a wave of genetic discoveries thanks to the development of several technologies, such as genome-wide association studies (GWAS), whole-exome sequencing, and whole genome sequencing. Despite the massive genetic discoveries of new variants associated with human diseases, several key challenges emerge following the genetic discovery. GWAS is known to be good at identifying the locus associated with the patient phenotype. However, the actually causal variants responsible for the phenotype are often elusive. Another challenge in human genetics is that even the causal mutations are already known, the underlying biological effect might remain largely ambiguous. Functional evaluation plays a key role to solve these key challenges in human genetics both to identify causal variants responsible for the phenotype, and to further develop the biological insights from the disease-causing mutations.

We adopted various methods to characterize the effects of variants identified in human genetic studies, including patient genetic and phenotypic data, RNA chemistry, molecular biology, virology, and multi-electrode array and primary neuronal culture systems. Chapter 1 is a broader introduction for the motivation and challenges for functional evaluation in human genetic studies, and the background of several genetics discoveries, such as hepatitis C treatment response, in which we performed functional characterization.

Chapter 2 focuses on the characterization of causal variants following the GWAS study for hepatitis C treatment response. We characterized a non-coding SNP (rs4803217) of IL28B (IFNL3) in high linkage disequilibrium (LD) with the discovery SNP identified in the GWAS. In this chapter, we used inter-disciplinary approaches to characterize rs4803217 on RNA structure, disease association, and protein translation.

Chapter 3 describes another avenue of functional characterization following GWAS focusing on the novel transcripts and proteins identified near the IL28B (IFNL3) locus. It has been recently speculated that this novel protein, which was named IFNL4, may affect the HCV treatment response and clearance. In this chapter, we used molecular biology, virology, and patient genetic and phenotypic data to further characterize and understand the biology of IFNL4. The efforts in chapter 2 and 3 provided new insights to the candidate causal variant(s) responsible for the GWAS for HCV treatment response, however, more evidence is still required to make claims for the exact causal roles of these variants for the GWAS association.

Chapter 4 aims to characterize a mutation already known to cause a disease (seizure) in a mouse model. We demonstrate the potential use of multi-electrode array (MEA) system for the functional characterization and drug testing on mutations found in neurological diseases, such as seizure. Functional characterization in neurological diseases is relatively challenging and available systematic tools are relatively limited. This chapter shows an exploratory research and example to establish a system for the broader use for functional characterization and translational opportunities for mutations found in neurological diseases.

Overall, this dissertation spans a range of challenges of functional evaluations in human genetics. It is expected that the functional characterization to understand human mutations will become more central in human genetics, because there are still many biological questions remaining to be answered after the explosion of human genetic discoveries. The recent advance in several technologies, including genome editing and pluripotent stem cells, is also expected to make new tools available for functional studies in human diseases.

Genetic and Molecular Basis of Encapsulation and Capsule Diversity in Kingella kingae

Kingella kingae is a bacterial pathogen that is increasingly recognized as an etiology of septic arthritis, osteomyelitis, bacteremia, and endocarditis in young children. The pathogenesis of K. kingae disease starts with bacterial adherence to the respiratory epithelium of the posterior pharynx. Previous work has identified type IV pili and a trimeric autotransporter protein called Knh (Kingella NhhA homolog) as critical factors for adherence to human epithelial cells. Additional studies established that the presence of a polysaccharide capsule interferes with Knh-mediated adherence. Given the inhibitory role of capsule during adherence we sought to uncover the genes involved in capsule expression to understand how capsule is elaborated on the cell surface. Additionally, this work aimed to further characterize capsule diversity among K. kingae clinical isolates and to investigate the relationship between capsule type and site of isolation.

We first set out to identify the carbohydrates present in the K. kingae capsule present in the prototype strain 269-492. Glycosyl composition and NMR analysis of surface extractable polysaccharides demonstrated two distinct polysaccharides, one consisting of GalNAc and Kdo with the structure →3)-β-GalpNAc-(1→5)-β-Kdop-(2→ and the other containing galactose alone with the structure →5)-β-Galf-(1→.

To discern the two polysaccharides we disrupted the ctrA gene required for surface localization of the K. kingae polysaccharide capsule and observed a loss of GalNAc and Kdo but no effect on the presence of Gal in bacterial surface extracts. In contrast, deletion of the pamABCDE locus involved in production of a reported galactan exopolysaccharide eliminated Gal but had no effect on the presence of GalNAc and Kdo in surface extracts. These results established that K. kingae strain KK01 produces a polysaccharide capsule with the structure →3)-β-GalpNAc-(1→5)-β-Kdop-(2→ and a separate exopolysaccharide with the structure →5)-β-Galf-(1→.

Having established that K. kingae produces a capsule comprised of GalNAc and Kdo, we next set out to identify the genetic determinants of capsule through a transposon mutagenesis screen. In addition to the previously identified ctrABCD operon, lipA, lipB, and a putative glycosyltransferase termed csaA (capsule synthesis region A gene A) were found to be essential for the production of surface-localized capsule. The ctr operon, lipA, lipB, and csaA were found to be present at unlinked locations throughout the genome, which is atypical for gram-negative organisms that elaborate a capsule dependent on an ABC-type transporter for surface localization. Through examining capsule localization in the ctrA, lipA, lipB, and csaA mutant strains, we determined that the ctrABCD, lipA/lipB, and csaA gene products respectively function in capsule export, assembly, and synthesis, respectively. The GalNAc transferase and Kdo transferase domains found in CsaA further support its role in catalyzing the synthesis of the GalNAc-Kdo capsule in the K. kingae prototype strain.

To investigate the capsule diversity that exists in K. kingae we screened a panel of strains isolated from patients with invasive disease or healthy carriers for the csaA capsule synthesis locus. We discovered that Kingella kingae expresses one of 4 capsule synthesis loci (csa, csb, csc, or csd) associated with a capsule consisting of Kdo and GalNAc (type a), Kdo and GlcNAc (type b), Kdo and ribose (type c), and GlcNAc and galactose (type d), respectively. Cloning of the csa, csb, csc, or csd locus into the empty flanking gene region in a non-encapsulated mutant (creation of an isogenic capsule swap) was sufficient to produce either the type a, type b, or type c capsule, respectively, further supporting the role of these loci in expression of a specific polysaccharide linkage. Capsule type a and capsule type b accounted for 96% of invasive strains. Conversely, capsule type c and capsule type d were found disproportionately among carrier isolates, suggesting that capsule type is important in promoting invasion and dissemination.

In conclusion, we discovered that Kingella kingae expresses a polysaccharide capsule and an exopolysaccharide on its surface that require distinct genetic loci for surface localization. Further investigation into genetic determinants of encapsulation revealed the loci ctrABCD, lipA/lipB, and a putative glycosyltransferase are required for capsule expression, with the gene products having roles in capsule export, assembly, and synthesis, respectively. The putative glycosyltransferase CsaA was determined to be a bifunctional enzyme with both GalNAc-transferase and Kdo-transferase activity. Furthermore, we discovered a total of 4 capsule types expressed in clinical isolates of K. kingae, each with a distinct capsule synthesis locus. The variation in the proportion of capsule types found between invasive strains and carriage strains suggest that capsule type is important in promoting invasion and dissemination. Taken together, this work expands our knowledge of the capsule types expressed among K. kingae carrier and invasive isolates and provides insights into the common genetic determinants of capsule expression. These contributions may lead to selecting clinically relevant capsule types to develop into a capsule based vaccine to prevent K. kingae colonization.